Re: ESMF regrid

Hi Noel,

Can you provide me with your script and data? It looks like there's some problem creating the source description file.
The error makes me think that the file is too large, but the SrcLargeFile should have handled this.

--Mary

On Mar 5, 2013, at 4:34 PM, Noel Aloysius wrote:

> Hi Mary,
>
> I am using NCL 6.1.0.
>
> I modified the script as per you suggestion, but still get the following error message.
>
> It works when I select a smaller region though.
>
> ncendef: ncid 131072: NetCDF: One or more variable sizes violate format constraints
> fatal:NetCDF: Operation not allowed in define mode: error attempting to write variable (grid_dims) to file (/home/fas/saiers/nra23/scratch/GPP_NPP_MODIS/GPP_src_SCRIP.nc)
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 1168 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 1277 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 2976 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 3089 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> fatal:["Execute.c":7741]:Execute: Error occurred at or near line 90
>
> Thanks,
> Noel
>
>
>
> On Tue, Mar 5, 2013 at 5:25 PM, Mary Haley <haley@ucar.edu> wrote:
> Hi Noel,
>
> What version of NCL are you using? It looks like a large file issue, so I think this should fix the problem:
>
> Opt@SrcLargeFile = True
> Opt@DstLargeFile = True
>
> These are the defaults for Version 6.1.2.
>
> --Mary
>
> On Mar 5, 2013, at 1:01 PM, Noel Aloysius wrote:
>
> > Hi NCL-talk,
> >
> > I want to regrid a high-resolution (~1km) MODIS GPP to a coarse resolution (1x1) grid using ESMF_regrid. The script is appended below.
> > I get the following error messages,
> >
> > ncendef: ncid 131072: NetCDF: One or more variable sizes violate format constraints
> > fatal:NetCDF: Operation not allowed in define mode: error attempting to write variable (grid_dims) to file (/home/fas/scratch/GPP_NPP_MODIS/src_SCRIP.nc)
> > fatal:["Execute.c":7743]:Execute: Error occurred at or near line 1168 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> > fatal:["Execute.c":7743]:Execute: Error occurred at or near line 1277 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> > fatal:["Execute.c":7743]:Execute: Error occurred at or near line 2976 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> > fatal:["Execute.c":7743]:Execute: Error occurred at or near line 3089 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
> > fatal:["Execute.c":7741]:Execute: Error occurred at or near line 27
> >
> > Thanks in advance for your help to resolve this issue,
> > Noel
> >
> >
> > ;**************************************************************************
> > srcDir = "~/"
> > srcFile = "MOD17A3_Science_GPP_2000.nc"
> > srcFileName = str_concat((/srcDir,srcFile/))
> >
> > sfile = addfile(srcFileName,"r") ; source grid
> > gpp = int2flt (sfile->MOD17A3_Science_GPP_2000)
> > gpp = gpp(::-1,:) ; South -> North
> >
> > Opt = True
> >
> > Opt@SrcFileName = srcDir + "src_SCRIP.nc"
> > Opt@DstFileName = srcDir + "dst_SCRIP.nc"
> > Opt@WgtFileName = srcDir + "GPP1km_2_halfddeg.nc"
> > Opt@ForceOverwrite = True
> >
> > Opt@DstGridType = "1x1"
> > Opt@DstLLCorner = (/ -59.5, -179.5/)
> > Opt@DstURCorner = (/ 79.5, 179.5/)
> >
> > gpp_regrid = ESMF_regrid(gpp,Opt)
> > ;**************************************************************************
> >
> > printVarSummary(gpp),
> >
> > Variable: gpp
> > Type: float
> > Total Size: 2903040000 bytes
> > 725760000 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [lat | 16800] x [lon | 43200]
> > Coordinates:
> > lat: [-59.99583333334029..79.99583333333334]
> > lon: [-179.9958333333333..179.9958333333346]
> > Number Of Attributes: 6
> > _FillValue : 1e+20
> > long_name : MOD17A3_Science_GPP_2000
> > esri_pe_string : GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
> > coordinates : lon lat
> > units : Degree
> > missing_value : 0
> >
> > ncl 12> print(min(gpp))
> > (0) 1
> > ncl 13> print(max(gpp))
> > (0) 65535
> >
> >
> > _______________________________________________
> > 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Mar 7 07:56:30 2013