Re: Regridding Daymet Joined Tiles with ESMF_Regridding on NCL 6.1.0

The wget had 12100. That was a mistake. It should be 12296
as indicated in the plots associated with each example.

Th WWW site has been updated. Make sure you reload all links.

===
Regarding tile 12100 ... The tile has no useable data for year 2010

On 10/31/12 2:22 PM, Ping Yang wrote:
> Hi NCL,
>
> I followed the http://www.ncl.ucar.edu/__Applications/
> <http://www.ncl.ucar.edu/Applications/> (http://www.ncl.ucar.edu/Applications/Scripts/daymet_2.ncl)
> for Daymet data, and I found problems, actually the file name should be
> prcp.2010_12100.nc <http://prcp.2010_12100.nc>.
>
> *I didn't see the **dstFile generated *and the plotting is wrong:
> Inline image 1
>
> The
>
> There were warning information:
>
> warning:ContourPlotInitialize: scalar field is constant; ContourPlot not
> possible:[errno=1102]
> warning:ContourPlotSetValues: Data values out of range of levels set by
> MANUALLEVELS mode
> warning:ContourPlotInitialize: scalar field is constant; ContourPlot not
> possible:[errno=1102]
> warning:ContourPlotSetValues: Data values out of range of levels set by
> MANUALLEVELS mode
>
> I also have questions on the code creating the destination grid:
>
> ;---Create the destination lat/lon grid
> NLAT = nlat/4 ; 4 is arbitrary
> MLON = mlon/3 ; 3 is arbitrary
> lat = fspan( min(lat2d),max(lat2d),NLAT)
> lon = fspan( min(lon2d),max(lon2d),MLON)
>
> Opt@DstGridType = "rectilinear"
> Opt@DstGridLat = lat
> Opt@DstGridLon = lon
>
>
> I saw the number 4 and 3 are arbitrary, why? If I plan to regrid this 1
> KM to 4 KM, are they should be
> *nlat/2 *and *mlon/2*?
>
> Looking forward to hearing from you.
>
> Regards,
>
> Ping
>
>
> On Wed, Oct 31, 2012 at 2:09 PM, Dennis Shea <shea@ucar.edu
> <mailto:shea@ucar.edu>> wrote:
>
> Dave [, Ping,]
>
> http://www.ncl.ucar.edu/__Applications/
> <http://www.ncl.ucar.edu/Applications/>
>
> Click Daymet under Datasets
>
> Comments?
>
> D
>
>
>
>
> On 10/30/12 3:51 PM, Ping Yang wrote:
>
> 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
> <http://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
> <http://curv_to_3minute.nc>
> <http://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
> <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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 31 14:35:49 2012