ESMF_regrid

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 Sun Jun 08 04:54:55 2014