continents contours and corresponding variable values don't overlap

Hi,

I have a binary data file with global temperature values. I want to mask out the continents and have the temperature over the ocean only. I use mask function and all works but when I plot the temperature the continents contours and the temperature values are not corresponding - they are somehow translated ( i.e. the contour of Australia is over Africa temperature values). Could anyone please tell me how to fix this problem ?

The code I used is the following :

; **********************
; read the files
  path = ("/data/oldfarm2/colfi/cam3.1.p2/CCSM/IEt42_test1/temp/ts_interp.b30C.dat")
  nlat = 64
  nlon = 128
  nrec=1560
  ts = new((/nrec,nlat,nlon/),float)
  
  do ir=1,1559
  rec=ir
  ts(ir,:,:)=fbindirread(path,rec,(/nlat,nlon/),"float")
  end do

  ins=addfile("/data/oldfarm2/colfi/cam3.1.p2/CCSM/IEt42_test1/scripts/IEt42_test1.cam.1870-01.nc","r")
  oro = ins->LANDFRAC(0,:,:)

;************************************************
 mask out land
  
   ocean_only = ts
   ocean_only = mask(ts,oro,0) ; Return the values of ts where oro=0.8

Thanks,

Ioana
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Feb 11 11:58:48 2011