;===============================================================================
; open a restart file
; change units on salt
; write back out
;===============================================================================

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

;-------------------------------------------------------------------------------
begin

fileName = "BEN1.c10b05nrcm.B1850CN.f09_g16.00.pop.r.0001-01-03-0000_ic.nc"


fileHandle = addfile(fileName,"w")

field    = "salt"
data = fileHandle->$field$
printVarSummary(data)
print("pre  max = "+max(data(:,:,:,:)))
print("pre  min = "+min(data(:,:,:,:)))
data(:,:,:,:) = data(:,:,:,:)*1000.0000000000000000000000
;do i = 0,ni-1
;do j = 0,nj-1
   ;factor = 1000.0
   ;data(k,j,i) = data(k,j,i)*factor
;end do
;end do
print("post max = " + max(data(:,:,:,:) ))
print("post min = " + min(data(:,:,:,:) ))
fileHandle->$field$ = data
delete(data)

end

;===============================================================================