question for lhtfl regrid

Dear NCL:
       Recently ,I perform areal conservative remapping from gaussian grid(192¡Á94) to one fixed grid (144¡Á72 grid) using area_conserve_remap_Wrap .However,when I use matlab to read the result,I found a little program.So ,I need your help .Can you tell me the principle of area_conserve_remap_Wrap ?I read the follow USL http://www.ncl.ucar.edu/Document/Functions/Contributed/area_conserve_remap_Wrap.shtml,
 but only the input and output are introduced. Could you help me check my code?
 "lhtfl.sfc.gauss.2012.nc" is source data,and "lhtfl_2007regrid.nc" is the data after the interpolation.
 I look forward to hearing from you.
  
  
 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
 p= addfile("/home/administrator/lhtfl.sfc.gauss.2012.nc","r")
f=p->lhtfl
f=f(:,::-1,:)
 opt = False
   NLAT = 72
   MLON = 144
   LAT = latGlobeFo(NLAT, "LAT", "latitude" , "degrees_north")
   LON = lonGlobeFo(MLON, "LON", "longitude", "degrees_east" )
   X = area_conserve_remap_Wrap(f&lon, f&lat, f,LON, LAT, opt)
   X=X(:,::-1,:)
 out=addfile("lhtfl_2007newregrid.nc","c")
out->X=X
 end


´ÓQQÓÊÏä·¢À´µÄ³¬´ó¸½¼þ

data.tar.gz (20.14M, 2013Äê12ÔÂ25ÈÕ 10:28 µ½ÆÚ)½øÈëÏÂÔØÒ³Ã棺http://mail.qq.com/cgi-bin/ftnExs_download?k=7830633901878acbb81827071e39051a0152500d5d0e010c150952000d1403035c544e000d0f5618005307080a5b04060b00530e38323751594402174c58451b5f4a6304&t=exs_ftn_download&code=80c98975

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Nov 24 19:29:43 2013