Re: ESMF Regridding

Before proceeding, could you please do a

%> ncdump -h air.2m.mon.mean.nc | less

What numeric type is the variable "air" ?

If it is type 'short', then you *must* unpack prior to use.
Use 'short2flt' located in 'contributed.ncl'
http://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml

     datAir = short2flt( sfile->air )
     printVarSummary(datAir)
     print("datAir: min="+min(datAir)+" max="+max(datAir))

Perhaps you were getting the "large errors" cuz you did not unpack.

On 2/19/13 5:20 PM, Travis, Katherine wrote:
> Hi,
> Could anyone who uses the ESMF_regridding.ncl routine tell me how to properly set the following options?
> I am trying to regrid sparse emissions data on a regional (US) 12km by 12km curvilinear grid to a 1x1degree grid, and the bilinear and path methods give me large errors.
> I am trying to use the conserve method, but I get the error that my grid corners are bad. Here is how I am currently trying to set them up:
>
> srcFile = "air.2m.mon.mean.nc"
> sfile = addfile(srcFile,"r")
> datAir = sfile->air
> lat = sfile->lat
> lon = sfile->lon
>
> Opt@SrcGridCornerLat = (/min(lat(0,:)),min(lat(:,0)), max(lat(0,:)),max(lat(:,0))/)
> Opt@SrcGridCornerLon = (/min(lon(0,:)),min(lon(:,0)), max(lon(0,:)),max(lon(:,0))/)
> Opt@DstGridCornerLat = (/1,1,11,47/)
> Opt@DstGridCornerLon = (/-144,-179,-68, 179/)
>
> Thank you,
> --
> Katherine R. Travis
> PhD Student
> Harvard University
> SEAS
> ktravis@fas.harvard.edu
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 19 22:57:21 2013