[ncl-talk] Using NCL to Regrid NetCDF data for ArcGIS

Matthew Shultz matthew.shultz at yale.edu
Tue Jun 23 15:31:20 MDT 2015


Here’s the output I get when I run the command.  Line 40 in my file is just the call to  ESMF_regrid
Thanks for you assistance


(0) get_src_grid_info: source lat dims = (277,349)
(0) get_src_grid_info: source lon dims = (277,349)
(0) get_src_grid_info: source grid type is 'curvilinear'
fatal:Eq: Dimension size, for dimension number 0, of operands does not match, can't continue

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 1614 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 3853 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 3983 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":8575]:Execute: Error occurred at or near line 40 in file regrid3.ncl




On Jun 23, 2015, at 1:22 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:

> What line number does the error message refer to?
> 
> 
> 
> On Mon, Jun 22, 2015 at 7:42 PM, Matthew Shultz <matthew.shultz at yale.edu> wrote:
> Greetings,
> I am very new to netcdf  and need some assistance.  I am trying to calculate zonal statistics (sums and averages within watershed shapes) using some tools in ArcGIS for netcdf.  The tools require netcdf files to be in a regular grid.  My source data is NARR data, which isn't being recognized as a regular grid (I believe it is curvilinear).  So I am trying to regrid using templates provided on the ucar website.  Attached is the code I am currently using.  It gives me the following error.  Despite search archives of this list, other responses have not resolved my problem. Help is very much appreciated, clearly there is something I am not understanding.
> 
> 
> 
> fatal:Eq: Dimension size, for dimension number 0, of operands does not match, can't continue
> 
> 
> 
> 
> 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/esmf/ESMF_regridding.ncl"
> 
> ;---Data file containing source grid
>     src_file = "air.sfc.2009.nc"                    ;;---Change (likely)
>     sfile    = addfile(src_file,"r")
> 
> ;---Get variable to regrid
>     varname = "air"                    ;;---Change (likely)
>     var     = sfile->air
>     src_lat = sfile->lat               ;;---Change (maybe)
>     src_lon = sfile->lon               ;;---Change (maybe)
> 
> ;---Set up regridding options
>     Opt                   = True
> 
> ;---"bilinear" is the default. "patch" and "conserve" are other options.
>     Opt at InterpMethod      = "bilinear"        ;;---Change (maybe)
> 
>     Opt at WgtFileName       = "curv_to_1deg.nc"
> 
>     Opt at SrcGridLat        = src_lat           ; source grid
>     Opt at SrcGridLon        = src_lon
>     Opt at SrcRegional       = True             ;;--Change (maybe)
>     Opt at SrcInputFileName  = src_file          ; optional, but good idea
>     Opt at SrcMask2D         = where(.not.ismissing(var),1,0) ; Necessary if has
>                                                            ; missing values.
> 
>     Opt at DstGridType       = "0.5x0.5"            ; Destination grid
>     ;Opt at DstLLCorner       = (/-89.75d,   0.00d /)   ;;--Change (likely)
>     ;Opt at DstURCorner       = (/ 89.75d, 359.75d /)   ;;--Change (likely)
>     Opt at DstRegional       = True                 ;;--Change (maybe)
> 
>     Opt at ForceOverwrite    = True
>     Opt at PrintTimings      = True
>     Opt at Debug             = True
> 
>     var_regrid = ESMF_regrid(var,Opt)     ; Do the regridding
> 
>     printVarSummary(var_regrid)
> 
> 
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150623/02a3b153/attachment.html 


More information about the ncl-talk mailing list