interpolation

Dear all,

I did a test interpolate WRF output to regular lon-lat grid using
rcm2rgrid and interpolate back to the original WRF
grid using rgrid2rcm. The figure of output seems not the same as the
original one. Many grids with no data. what¢s
wrong? The following is my script. Any suggetion will be appreciated.

.....

Lon = fspan (230,290,60)
Lat = fspan (25,55,30)
U=f->U(:,:)
Lat2d =f->lat(:,:)
Lon2d=f->lon(:,:)
Xgrd_U = rcm2rgrid(lat2d,lon2d,U,lat,lon,0 )
U_back = rgrid2rcm(lat,lon,U,lat2d,lon2d,0 )

filout = "output.nc"
fout = addfile (filout, "c")
fout->prate = U
fout->xgrd_UQ =xgrd_U
fout->temp =U_back
 

Yanhong

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed May 19 18:08:08 2010