Re: Regridding Daymet Joined Tiles with ESMF_Regridding on NCL 6.1.0

Dear NCL,

I installed a new version of gcc and located the library directory to
DYLD_LIBRARY_PATH,
and now it the libgomp.1.dylib issue went away. Thank you so much.

Can you have a look at this script(Regrid a curvilinear grid from 1km to
rectilinear grid at 4km, I was trying to find an example on NCL website,
but there is one regriding from WRF grid to the same resolution and another
template regridding to 1 degree, so I intended to compile them together):

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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
begin

*;---Data file containing source grid*
* src_file = "~/Data/Daymet_NE_20100418_prcp.nc" *
* sfile = addfile(src_file,"r")*
* dst_file = "~/Data/Daymet_NE_20100418_Regridded_prcp.nc"*
* dfile = addfile(dst_file,"c")*
*;---Get variable to regrid*
* varname = "prcp" *
* var = sfile->$varname$*
* src_lat = sfile->lat ;;---Change (maybe)*
* src_lon = sfile->lon ;;---Change (maybe)*
*;---Set up regridding options*
* Opt = True*
* Opt@SrcFileName = src_file ; Name of source and*
* Opt@DstFileName = dst_file ; destination files*
*;---"bilinear" is the default. "patch" and "conserve" are other options.*
* Opt@InterpMethod = "bilinear"*
* Opt@WgtFileName = "~/Data/curv_to_3minute.nc"*
*
*
* Opt@SrcGridLat = src_lat ; source grid*
* Opt@SrcGridLon = src_lon*
* Opt@SrcRegional = True ;;--Change (maybe)*
* Opt@DstRegional = True*
*
*
* newlon = fspan(-84.025,-64.875,384)*
* newlat = fspan(34.875,50.125,306)*
*
*
*;---Create the destination lat/lon grid*
*; lat = fspan( 53.2420, 85.4022,nlat)*
*; lon = fspan(135.7750,258.1880,nlon)*
*
*
* Opt@DstGridType = "rectilinear"*
* Opt@DstGridLat = newlat*
* Opt@DstGridLon = newlon*
*
*
* Opt@SrcRegional = True*
* Opt@DstRegional = True*
*
*
* Opt@ForceOverwrite = True*
* Opt@PrintTimings = True*
* Opt@Debug = True*
* *
* var_regrid = ESMF_regrid(var,Opt) ; Do the regridding*
* printVarSummary(var_regrid)*
*
*
*;----------------------------------------------------------------------*
*; Plotting section*
*;*
*; This section creates filled contour plots of both the original*
*; data and the regridded data, and panels them.*
*;----------------------------------------------------------------------*
* var@lat2d = src_lat ; Needed for plotting.*
* var@lon2d = src_lon*
*
*
* wks = gsn_open_wks("pdf","curv_to_3min")*
* xwks = gsn_open_wks("x11","curv_to_3min")*
*
*
* res = True*
*
*
* res@gsnMaximize = True*
*
*
* res@gsnDraw = False*
* res@gsnFrame = False*
*
*
* res@cnFillOn = True*
* res@cnLinesOn = False*
* res@cnLineLabelsOn = False*
* res@cnFillMode = "RasterFill"*
*
*
* res@lbLabelBarOn = False ; Turn on later in panel*
*
*
* res@mpMinLatF = min(src_lat)*
* res@mpMaxLatF = max(src_lat)*
* res@mpMinLonF = min(src_lon)*
* res@mpMaxLonF = max(src_lon)*
*
*
*;;--Change (likely)*
*; res@cnLevelSelectionMode = "ManualLevels"*
*; res@cnMinLevelValF = 0*
*; res@cnMaxLevelValF = 200*
*; res@cnLevelSpacingF = 10*
*
*
*;---Resources for plotting regridded data*
* res@gsnAddCyclic = False ;;---Change (maybe)*
* res@tiMainString = "3 min grid (" + Opt@InterpMethod + ")"*
*
*
* plot_regrid_p = gsn_csm_contour_map(wks,var_regrid(0,:,:),res)*
* plot_regrid_x = gsn_csm_contour_map(xwks,var_regrid(0,:,:),res)*
*
*
*;---Resources for plotting original data*
* dims = tostring(dimsizes(var))*
* res@gsnAddCyclic = False ;;---Change (maybe)*
* res@tiMainString = "Original curvilinear grid (" + \*
* str_join(dims," x ") + ")"*
*
*
* plot_orig_p = gsn_csm_contour_map(wks,var,res) *
* plot_orig_x = gsn_csm_contour_map(xwks,var,res)*
*;---Compare the plots in a panel*
* pres = True*
* pres@gsnMaximize = True*
* pres@gsnPanelLabelBar = True*
*
*
* gsn_panel(wks,(/plot_orig_p,plot_regrid_p/),(/2,1/),pres)*
* gsn_panel(xwks,(/plot_orig_x,plot_regrid_x/),(/2,1/),pres)*
*end*

The information for the variable to be regrided:
*Variable: x*
*Type: float*
*Total Size: 11613888 bytes*
* 2903472 values*
*Number of Dimensions: 3*
*Dimensions and sizes: [time | 1] x [y | 1872] x [x | 1551]*
*Coordinates: *
* time: [11055.5..11055.5]*
*Number Of Attributes: 6*
* long_name : daily total precipitation*
* units : mm/day*
* coordinates : lon lat*
* _FillValue : -9999*
* valid_range : ( 0, 200 )*
* cell_methods : area: sum time: sum*

It is running forever without seeing the end. I am not quite understanding
the example on the defining of the resource grid and destination grid, Is
the resource grid of the source file? and also about the destination
grid(is it will generate automatically?) this seems much complicated
comparing to the precious version.

Looking forward to hearing from you.

Regards,

Ping

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 30 15:52:09 2012