(no subject)

Hi:
I am trying to write a template to regrid an unstructured MPAS grid to
an other higher or lower resolution unstructured MPAS grid using ESMF.
I am using the function ESMF_regrid_gen_weights so that I can reuse
the file containing the weights. See the ncl script below.

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

begin

r2d = 180.0d/(atan(1)*4.0d) ; conversion from radians to degrees.

src_file = "x1.40962.grid.nc"
sfile = addfile(src_file,"r")

src_file2 = "x1.10242.grid.nc"
sfile2 = addfile(src_file2,"r")

opt = True
opt@Debug = True
opt@InterpMethod = "bilinear"
opt@ForceOverwrite = True
opt@SrcESMF = True
opt@DstESMF = True

opt@SrcGridLon = sfile2->lonCell*r2d ; source grid lon input.
opt@SrcGridLat = sfile2->latCell*r2d ; source grid lat input.
opt@DstGridLon = sfile->lonCell*r2d ; destination grid lon grid.
opt@DstGridLat = sfile->latCell*r2d ; destination grid lat grid.

srcGridName = src_file2
dstGridName = "x1.10242_to_40962.nc"
wgtGridName = "weights.10242_to_40962.nc"
ESMF_regrid_gen_weights(srcGridName,dstGridName,wgtGridName,opt)

end

I get the following error message:
(0) x1.10242.grid.nc
(0) ESMF_regrid_gen_weights: number of processors used: 1
(0) --------------------------------------------------
(0) ESMF_regrid_gen_weights: the following command is about to be
executed on the system:
(0) 'ESMF_RegridWeightGen --source x1.10242.grid.nc --destination
x1.10242_to_40962.nc --weight weights.10242_to_40962.nc --method
bilinear --src_type ESMF --dst_type ESMF -i --64bit_offset'
(0) --------------------------------------------------
  ESMF_LogWrite(): Deprecated: Use logmsgFlag instead of logmsgList.
(0) ESMF_regrid_gen_weights: output from 'ESMF_RegridWeightGen':
(0) Starting weight generation with these inputs:
(1) Source File: x1.10242.grid.nc
(2) Destination File: x1.10242_to_40962.nc
(3) Weight File: weights.10242_to_40962.nc
(4) Source File is in ESMF format
(5) Source Grid is a global grid
(6) Source Grid is an unstructured grid
(7) Destination File is in ESMF format
(8) Destination Grid is a global grid
(9) Destination Grid is an unstructured grid
(10) Regrid Method: bilinear
(11) Pole option: ALL
(12) Ignore unmapped destination points
(13) Output weight file in 64bit offset NetCDF file format
(14)
(15) NetCDF Error: Dimension nodeCount in x1.10242.grid.nc :
NetCDF: Invalid dimension ID or name
(16) ERROR: Problem on processor 0 . Please see the
PET*.RegridWeightGen.Log files for a traceback.
(0) --------------------------------------------------
(0) ESMF_regrid_gen_weights: 'ESMF_RegridWeightGen' was not
successful.

and my PET file looks like

20131009 143041.025 INFO PET0 Running with ESMF Version 5.2.0rp2
20131009 143041.120 ERROR PET0 ESMF_IOScrip.F90:1963
ESMF_GetMeshFromFile netCDF Status Return Error
20131009 143041.120 ERROR PET0 ESMF_Mesh.F90:1192
ESMF_MeshCreateFromUnstruct() Failure - Internal subroutine call
returned Error
20131009 143041.120 ERROR PET0 ESMF_Mesh.F90:1075
ESMF_MeshCreateFromFile() Failure - Internal subroutine call returned
Error
20131009 144050.758 INFO PET0 Running with ESMF Version 5.2.0rp2
20131009 144050.763 ERROR PET0 ESMF_IOScrip.F90:1963
ESMF_GetMeshFromFile netCDF Status Return Error
20131009 144050.763 ERROR PET0 ESMF_Mesh.F90:1192
ESMF_MeshCreateFromUnstruct() Failure - Internal subroutine call
returned Error
20131009 144050.763 ERROR PET0 ESMF_Mesh.F90:1075
ESMF_MeshCreateFromFile() Failure - Internal subroutine call returned
Error

I did not find any information in the ESMF documentation regarding
that error message. I would really like a hint on what to do next.

Thanks,
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 Wed Oct 9 14:42:52 2013