ESMCI_DistGrid.C:5101 ESMCI::DistGrid::getSequenceInde Invalid argument - SeqIndex type mismatch detected

[axelschweiger]

Requirements

• Reviewed ESMF Reference Manual
• Searched GitHub Discussions

Affiliation(s)

University of Washington

ESMF Version

v8.4.0

Issue

I’m trying to use ESMF_RegridWeightGen to regrid the GEBCO 2024 Bathymetry data set to a regional grid using a curvilinear coordinate system specified in lat,lon coordinates. The Pole is in inside the North American Continent (Alaska).

Here is the command I am using:
ESMF_RegridWeightGen -s 2024/GEBCO_2024.nc --dst_regional --ignore-unmapped -d capmas_static_scalar.temp.nc -w Gebco2024x864x640.script.nc -m bilinear

The generated Log file is getting very large (I’ve tried –no_log) with messages:
20241125 170059.675 ERROR PET0 ESMCI_DistGrid.C:5101 ESMCI::DistGrid::getSequenceInde Invalid argument - SeqIndex type mismatch detected.

Running this for testing in single processor mode for now.

This is the output from the script.
Starting weight generation with these inputs:
Source File: 2024/GEBCO_2024.nc
Destination File: capmas_static_scalar.temp.nc
Weight File: Gebco2024x864x640.script.nc
Source File is in CF Grid format
Source Grid is a global grid
Source Grid is a logically rectangular grid
Use the center coordinates of the source grid to do the regrid
Destination File is in CF Grid format
Destination Grid is a regional grid
Destination Grid is a logically rectangular grid
Use the center coordinates of the destination grid to do the regrid
Regrid Method: bilinear
Pole option: ALL
Line Type: cartesian
Norm Type: dstarea
Extrap. Method: none

ncdump -h on output grid show this:

[](dimensions:
              x = 864 ;
              y = 640 ;
variables:
              float lat(y, x) ;
                             lat:_FillValue = 1.e+20f ;
                             lat:units = "degrees_north" ;
              float lon(y, x) ;
                             lon:_FillValue = 1.e+20f ;
                             lon:units = "degrees_east" ;

// global attributes:
                             :history = "Thu Nov 21 13:09:26 2024: ncatted -a units,lat,o,c,degrees_north -a units,lon,o,c,degrees_east capmas_static_scalar.20241121.nc -o test.nc" ;
                             :NCO = "netCDF Operators version 5.1.4 (Homepage = http://nco.sf.net/, Code = http://github.com/nco/nco)" ;
}
)

Input grid: 
dimensions:
              lon = 86400 ;
              lat = 43200 ;
variables:
              double lon(lon) ;
                             lon:standard_name = "longitude" ;
                             lon:long_name = "longitude" ;
                             lon:units = "degrees_east" ;
                             lon:axis = "X" ;
                             lon:sdn_parameter_urn = "SDN:P01::ALONZZ01" ;
                             lon:sdn_parameter_name = "Longitude east" ;
                             lon:sdn_uom_urn = "SDN:P06::DEGE" ;
                             lon:sdn_uom_name = "Degrees east" ;
              double lat(lat) ;
                             lat:standard_name = "latitude" ;
                             lat:long_name = "latitude" ;
                             lat:units = "degrees_north" ;
                             lat:axis = "Y" ;
                             lat:sdn_parameter_urn = "SDN:P01::ALATZZ01" ;
                             lat:sdn_parameter_name = "Latitude north" ;
                             lat:sdn_uom_urn = "SDN:P06::DEGN" ;
                             lat:sdn_uom_name = "Degrees north" ;
              char crs ;
                             crs:grid_mapping_name = "latitude_longitude" ;
                             crs:epsg_code = "EPSG:4326" ;
                             crs:inverse_flattening = 298.257223563 ;
                             crs:semi_major_axis = 6378137. ;
              short elevation(lat, lon) ;
                             elevation:standard_name = "height_above_mean_sea_level" ;
                             elevation:long_name = "Elevation relative to sea level" ;
                             elevation:units = "m" ;
                             elevation:grid_mapping = "crs" ;
                             elevation:sdn_parameter_urn = "SDN:P01::BATHHGHT" ;
                             elevation:sdn_parameter_name = "Sea floor height (above mean sea level) {bathymetric height}" ;
                             elevation:sdn_uom_urn = "SDN:P06::ULAA" ;
                             elevation:sdn_uom_name = "Metres" ;

// global attributes:
                             :title = "The GEBCO_2024 Grid - a continuous terrain model for oceans and land at 15 arc-second intervals"
;
                             :summary = "The GEBCO_2024 Grid is a continuous, global terrain model for ocean and land with a spatial re
(

Autotag

@oehmke