Re: ESMF_regrid

Hi Laura,

The ESMF_regrid function got confused because you had set both DstGridType
and DstGridLat/DstGridLon. I will fix the code so that it exits more
gracefully in this case. In your case, I think you can fix the problem by
not setting DstGridType.

Let me know if that doesn't help.

--Mary


On Sun, Jun 8, 2014 at 10:54 AM, Laura Fowler <laura@ucar.edu> wrote:

> Hi:
> I am trying to regrid data from my MPAS mesh to some NNRP monthly data
> using ESMF_regrid, and to create a "weight" file that I can reuse
> later. My end goal is to make difference plots between my MPAS
> regridded outputs and NNRP data.
>
> I am getting the following error message and do not know how to solve
> this:
> (0) x1.40962.output.2005-05-01_00.00.00.nc
> (0) nTime = 124
> (0) nCells = 40962
> (0) get_src_grid_info: source lat dims = (40962)
> (0) get_src_grid_info: source lon dims = (40962)
> (0) get_src_grid_info: source grid type is 'unstructured'
> (0) unstructured_to_ESMF: triangulating the data ... this can be slow
> (0) min/max ElementVertices = 0/40961
> (0) unstructured_to_ESMF: total number of elements created: 81920
> (0) unstructured_to_ESMF: Element Area: min:0.0001103995352902398
> max:0.0001587103318516636
> (0) get_dst_grid_info: destination lat dims = (73)
> (0) get_dst_grid_info: destination lon dims = (144)
> fatal:Variable (latlon_dims) is undefined
> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 3191
> in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> The corresponding source code is:
> opt = True
> opt@Debug = True
> opt@SrcFileName = "source.x1_40962.nc"
> opt@DstFileName = "destination.nnrp_73x144rectangular.nc"
> opt@WgtFileName =
> "weights_onCells.40962_to_nnrp.73x144rectangular.nc"
> opt@ForceOverwrite = True
> opt@SrcGridLat = latCell
> opt@SrcGridLon = lonCell
> opt@DstGridLat = lat_nnrp
> opt@DstGridLon = lon_nnrp
> opt@DstGridType = "2.5deg" ; destination grid
> opt@DstTitle = "NNRP Rectangular Grid 2.5 degree resolution"
> ;opt@DstLLCorner = (/90.0d, 0.0d /)
> ;opt@DstURCorner = (/-90.0d,357.5d/)
>
> latCell and lonCell contains the latitude and longitude of every MPAS
> cells. lat_nnrp and lon_nnrp are the latitudes and longitudes of the
> NNRP data:
>
> lat_2 = 90, 87.5, 85, 82.5, 80, 77.5, 75, 72.5, 70, 67.5, 65, 62.5,
> 60,
> 57.5, 55, 52.5, 50, 47.5, 45, 42.5, 40, 37.5, 35, 32.5, 30, 27.5,
> 25,
> 22.5, 20, 17.5, 15, 12.5, 10, 7.5, 5, 2.5, 0, -2.5, -5, -7.5,
> -10, -12.5,
> -15, -17.5, -20, -22.5, -25, -27.5, -30, -32.5, -35, -37.5, -40,
> -42.5,
> -45, -47.5, -50, -52.5, -55, -57.5, -60, -62.5, -65, -67.5, -70,
> -72.5,
> -75, -77.5, -80, -82.5, -85, -87.5, -90 ;
>
> lon_2 = 0, 2.5, 5, 7.5, 10, 12.5, 15, 17.5, 20, 22.5, 25, 27.5, 30,
> 32.5,
> 35, 37.5, 40, 42.5, 45, 47.5, 50, 52.5, 55, 57.5, 60, 62.5, 65,
> 67.5, 70,
> 72.5, 75, 77.5, 80, 82.5, 85, 87.5, 90, 92.5, 95, 97.5, 100,
> 102.5, 105,
> 107.5, 110, 112.5, 115, 117.5, 120, 122.5, 125, 127.5, 130,
> 132.5, 135,
> 137.5, 140, 142.5, 145, 147.5, 150, 152.5, 155, 157.5, 160,
> 162.5, 165,
> 167.5, 170, 172.5, 175, 177.5, 180, 182.5, 185, 187.5, 190,
> 192.5, 195,
> 197.5, 200, 202.5, 205, 207.5, 210, 212.5, 215, 217.5, 220,
> 222.5, 225,
> 227.5, 230, 232.5, 235, 237.5, 240, 242.5, 245, 247.5, 250,
> 252.5, 255,
> 257.5, 260, 262.5, 265, 267.5, 270, 272.5, 275, 277.5, 280,
> 282.5, 285,
> 287.5, 290, 292.5, 295, 297.5, 300, 302.5, 305, 307.5, 310,
> 312.5, 315,
> 317.5, 320, 322.5, 325, 327.5, 330, 332.5, 335, 337.5, 340,
> 342.5, 345,
> 347.5, 350, 352.5, 355, 357.5 ;
>
> I really appreciate your help with this. The entire source code can
> also be found on
>
> /glade/scratch/laura/MPAS.projects/MPAS.nrcm.2014-04-20/NNRP.accum_diagnostics.gwdo/PLOTS
>
> Laura
>
> !----------------------------------------------------
> Laura D. Fowler
> Mesoscale and Microscale Meteorology Division (MMM)
> National Center for Atmospheric Research
> P.O. Box 3000, Boulder CO 80307-3000
>
> e-mail: laura@ucar.edu
> phone : 303-497-1628
>
> !----------------------------------------------------
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

Received on Thu Jun 12 06:16:33 2014