using ESMF

Hi,
I'm trying out the ESMF with some data:

      Opt = True
      Opt@Mask2D = where(.not.ismissing(sstc),1,0)
      Opt@GridCornerLat = fhc->lat_vertices
      Opt@GridCornerLon = fhc->lon_vertices
      Opt@PrintTimings = True
      Opt@ForceOverwrite = True
      Opt@Debug = True
      curvilinear_to_SCRIP("A1p.nc",lat_A1pg,lon_A1pg,Opt)
      delete(Opt)

      Opt = True
      Opt@ForceOverwrite = True
      Opt@LLCorner = (/-89.75d, 0.00d /)
      Opt@URCorner = (/ 89.75d, 359.75d /)
      Opt@Title = "World Grid 1-degree Resolution"
      Opt@PrintTimings = True
      Opt@Debug = True
      latlon_to_SCRIP("Had.nc","1x1",Opt)
      delete(Opt)

      Opt = True
      Opt@ForceOverwrite = True
      Opt@SrcESMF = False
      Opt@DstESMF = False
      Opt@InterpMethod = "bilinear"
      Opt@PrintTimings = True
      ESMF_regrid_gen_weights("A1p.nc", "Had.nc", "weight.nc", Opt)
      delete(Opt)

but when it gets to the ESMF_regrid_gen_weights I get the error message:

(0) curvilinear_to_SCRIP: using grid corners provided by user...
(0) =====> CPU Elapsed Time: curvilinear_to_SCRIP: 0.029996 <=====
(0) curvilinear_to_SCRIP: calculating grid corners...
(0) calc_SCRIP_corners
(0) min/max original lat: -89.75/89.75
(0) min/max original lon: 0/359.75
(0) calc_SCRIP_corners
(0) min/max grid_corner_lat2d: -89.75/89.75
(0) min/max grid_corner_lon2d: 0/359.75
(0) calc_SCRIP_corners
(0) min/max grid_corner_lat: -89.75/89.75
(0) min/max grid_corner_lon: 0/359.75
(0) =====> CPU Elapsed Time: latlon_to_SCRIP: 8.17876 <=====
ESMF_RegridWeightGen: error while loading shared libraries: libgfortran.so.3: cannot open shared object file: No such file or directory
fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
fatal:["Execute.c":7556]:Execute: Error occurred at or near line 1972 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl

fatal:["Execute.c":7556]:Execute: Error occurred at or near line 110 in file plot_sstbias.ncl

I have normally not needed to make any additional environment variables to get NCL to work properly. Is this a bug or something wrong with my setup or code?

[mac599@dcc Steve]$ ncl -V
6.1.0-beta

[mac599@dcc Steve]$ uname -a
Linux dcc 2.6.18-274.7.1.el5 #1 SMP Thu Oct 20 16:21:01 EDT 2011 x86_64 x86_64 x86_64 GNU/Linux

The new interpolation interface looks very promising.

I could not see any examples going from a rectlinear to curvlinear grid, this is the case that I am working on. It is quite typical to generate a model bias on the model grid, and therefore observations need to be interpolated to a model grid.

Many thanks,
Mark.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 25 00:49:41 2012