[ncl-talk] error ESMF_regrid

Mary Haley haley at ucar.edu
Tue May 16 10:16:57 MDT 2017 

There was no map projection provided on the files that you gave me, so I
just (lazily) used a C.E. projection.

I should have plotted the data over the north pole.  See the attached. I
also doubled the number of contour levels.

--Mary




On Tue, May 16, 2017 at 6:54 AM, Vanúcia Schumacher <
vanucia-schumacher at hotmail.com> wrote:

> Thanks for your help.
>
>
> I wounder if there is a way to fix the projection of my source data for
> don't look funny.
>
> When I make a ncview in the source data, it does not have this projection.
>
>
>
>
> Thanks for collaboration,
> ------------------------------
> *De:* Mary Haley <haley at ucar.edu>
> *Enviado:* terça-feira, 16 de maio de 2017 02:57:46
>
> *Para:* Vanúcia Schumacher
> *Cc:* ncl-talk at ucar.edu
> *Assunto:* Re: [ncl-talk] error ESMF_regrid
>
> Vanúcia,
>
> I was able to download your files and run your script. Thanks for
> providing these.
>
> There was main issue with your script: You are subsetting the source lat /
> lon grid with:
>
>  latSrc   = g->XLAT(::6,::6)
>  lonSrc   = g->XLONG(::6,::6)
>
> but you are not subsetting the source variable in the same way.
>
> In order to correctly regrid a data array, you must provide the lat/lon
> point for every element of your data array. To fix this, you either need to
> NOT subset the latSrc/lonSrc grids, or you need to subset varSrc0 with
> (::6,::6) as well. I chose the first option.
>
> Once I made this change, the regridding seemed to work fine.
>
> I created a panel plot so you could compare the plots from the first time.
> I don't know the correct projection of your source data, so I just plotted
> the data in a basic lat/lon plot.  The shape looks funny, but the
> contouring is correct.
>
> See the attached modified version of your script and the graphic. I also
> included a version of the script that creates an animation.
>
> --Mary
>
>
>
>
> On Mon, May 15, 2017 at 9:42 PM, Mary Haley <haley at ucar.edu> wrote:
>
>> ​Hi Vanúcia,
>>
>> ​
>> Just an FYI about what's going on here:
>>
>> With the destination grid, if you provide one-dimensional arrays for lat
>> and lon that are both the same length (let's called the length DLEN), the
>> regridding code has no idea if your intention is to regrid the data to a
>> rectilinear grid that is DLEN x DLEN, or to an unstructured set of random
>> points of length DLEN.  This is why you
>> ​ are​
>> get the warning:
>>
>>  get_dst_grid_info: Error: the lat/lon arrays are both 1D and the same
>> length.
>>    This could be a rectilinear or unstructured grid. Please set
>> 'grid_type'
>>
>>
>> Since you
>> ​said you wanted
>>  to regrid to a 1 degree grid, you did the right thing by setting:
>>
>> Opt at DstGridType = "rectilinear"
>>
>>
>> ​I should point out, however, that if you want to regrid to a 1 degree
>> grid, you could have also just set:
>>
>> Opt at DstGridType = "1deg"
>>
>>
>> and then you don't need to set:
>>
>> Opt at DstGridLat     = latDst
>> Opt at DstGridLon     = lonDst
>>
>>
>> ​How​ever, setting DstGridLat and DstGridLon yourself should be okay.
>>
>> I'm not sure what you mean by "the out_file It's a mess, with the
>> curvilinear grid."
>> ​  I am downloading the files you provided and will take a look.
>>
>> --Mary
>>
>>
>> On Mon, May 15, 2017 at 7:35 PM, Vanúcia Schumacher <
>> vanucia-schumacher at hotmail.com> wrote:
>>
>>> Hi
>>>
>>>
>>> I sent via ftp site the complete script, and
>>>
>>> the 2_file and 3_file used. Sorry, but I can't to reduze the size of the
>>> 1_file (shf_asr-2000-classic.nc) to send via ftp, but I attached in
>>> this e-mail (https://1drv.ms/f/s!AhH-sSGcLfkLj0WsaPkz5a_nbuo4)
>>> [image: imageData708542]
>>> <https://1drv.ms/f/s!AhH-sSGcLfkLj0WsaPkz5a_nbuo4>
>>> help <https://1drv.ms/f/s!AhH-sSGcLfkLj0WsaPkz5a_nbuo4>
>>> Compartilhado via OneDrive
>>>
>>>
>>> I include the line Opt at DstGridType = "rectilinear" and  the script
>>> works, but the out_file It's a mess, with the curvilinear grid.
>>>
>>> I would like to regress the file_1 (shf_asr-2000-classic.nc) to
>>> resolution of 1ºx1º with the grid rectilinear. I can not do this correctly.
>>>
>>>
>>> The names of the files:
>>>
>>>
>>> file_1 used: shf_asr-2000-classic.nc
>>>
>>>
>>> file_2: asr15km.static.2d.2007070100.nc
>>>
>>> file_3: shf_era-1x1.nc
>>>
>>>
>>> Thanks
>>> ------------------------------
>>> *De:* Kevin Hallock <hallock at ucar.edu>
>>> *Enviado:* segunda-feira, 15 de maio de 2017 21:53:51
>>> *Para:* Vanúcia Schumacher
>>> *Cc:* ncl-talk at ucar.edu
>>> *Assunto:* Re: [ncl-talk] error ESMF_regrid
>>>
>>> Hi,
>>>
>>> Could you please provide the complete script you are trying to run?
>>>
>>> Thanks,
>>> Kevin
>>>
>>> On May 15, 2017, at 10:57 AM, Vanúcia Schumacher <
>>> vanucia-schumacher at hotmail.com> wrote:
>>>
>>> Hi NCL users,
>>>
>>> I am trying to interpolate a data to a resolution of 1x1, but it is
>>> returning the following error:
>>>
>>> Thank you any collaboration
>>>
>>> (0) get_src_grid_info: source lat dims = (120,120)
>>> (0) get_src_grid_info: source lon dims = (120,120)
>>> (0) get_src_grid_info: source grid type is 'curvilinear'
>>> (0) curvilinear_to_SCRIP: calculating grid corners...
>>> (0) curvilinear_to_SCRIP: no lat values are at the poles, so
>>> (0)       calculating grid corners using
>>> (0)       calc_SCRIP_corners_noboundaries...
>>> (0) calc_SCRIP_corners_noboundaries
>>> (0)     min/max original lat: 24.6335/89.8972
>>> (0)     min/max original lon: -179.998/179.975
>>> (0) calc_SCRIP_corners_noboundaries
>>> (0)     min/max Extlat2d: 23.76/89.8972
>>> (0)     min/max Extlon2d: -179.998/179.975
>>> (0) calc_SCRIP_corners_noboundaries
>>> (0)     min/max ExtGridCenter_lat: 24.1982/89.2938
>>> (0)     min/max ExtGridCenter_lon: -179.399/179.417
>>> (0) get_dst_grid_info: destination lat dims = (360)
>>> (0) get_dst_grid_info: destination lon dims = (360)
>>> (0) get_dst_grid_info: Error: the lat/lon arrays are both 1D and the
>>> same length.
>>> (0)   This could be a rectilinear or unstructured grid. Please set
>>> 'grid_type'
>>> (0)   to 'rectilinear' or 'unstructured'.
>>>
>>>
>>> Variable: varSrc0
>>> Type: float
>>> Total Size: 323481600 bytes
>>>             80870400 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes:  [Time | 156] x [XLAT | 720] x [XLONG | 720]
>>> Coordinates:
>>>             Time: [868188..981394.5]
>>> Number Of Attributes: 6
>>>   FieldType : 104
>>>   MemoryOrder :  XY
>>>   description :  LATENT HEAT FLUX AT THE SURFACE
>>>   units : W m-2
>>>   stagger :
>>>   coordinates :  XLONG XLAT
>>>
>>> Opt                = True
>>> Opt at SrcGridLat     = latSrc
>>> Opt at SrcGridLon     = lonSrc
>>> Opt at DstGridLat     = latDst
>>> Opt at DstGridLon     = lonDst
>>> Opt at SrcRegional    = True ; because the source grid (ASR)  is a
>>> regional grid
>>> Opt at DstRegional    = False ; if regriddin to global; set to True if
>>> smoothing ASR
>>> Opt at InterpMethod   = "bilinear" ; this is the default method
>>> Opt at NoPETLog       = True
>>> Opt at RemoveSrcFile  = True
>>> Opt at RemoveDstFile  = True
>>> Opt at RemoveWgtFile  = True
>>> Opt at Debug          = True ; set this to False to make it less verbose
>>> Opt at CopyVarCoords  = True ; if True, the lat-lon data will be attached
>>> as
>>> Opt at ForceOverwrite = True
>>>
>>>
>>> ; Regridding function
>>>
>>> varDst0 = ESMF_regrid(varSrc0,Opt)
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>>
>>> _______________________________________________
>>> 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/20170516/d032de11/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: OutlookEmoji-imageData708542.png
Type: image/png
Size: 1785 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170516/d032de11/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: regrid_esmf_mod.ncl
Type: application/octet-stream
Size: 5212 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170516/d032de11/attachment-0001.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: asr_era_polar.png
Type: image/png
Size: 284694 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170516/d032de11/attachment-0003.png

More information about the ncl-talk mailing list