[ncl-talk] ESMF_RegridWeightGen: issue with conservative mapping

[Michele Petrini - CITG]

Dear helpdesk,

I am trying to generate an ensemble of mapping files between different model grids within the CESM2.0 model framework.

Whereas everything works fine between global grids (e.g., for CAM, POP, CLM and RTM), I am having problems between regional-global grids;

- CISM regional grid of Northern Hemisphere (lat between ~7 and ~90 deg N, lon between ~0 and ~360)
- CLM global grid

When I perform a bilinear interpolation, no problems occur - whereas when the interpolation method is ‘aave’, that is conservative mapping, I get the following:

20181030 085614.748 INFO             PET0 Running with ESMF Version 7.0.0
20181030 085614.764 INFO             PET0 Running with ESMF Version 7.0.0
20181030 085616.153 ERROR            PET0 ESMF_Grid.F90:5120 ESMF_GridCreate Wrong argument specified  - - Bad corner array in SCRIP file
20181030 085616.153 ERROR            PET0 ESMF_Grid.F90:6307 ESMF_GridCreateFrmScrip Wrong argument specified  - Internal subroutine call returned Error
20181030 085616.154 ERROR            PET0 ESMF_Grid.F90:5976 ESMF_GridCreateFrmNCFile Wrong argument specified  - Internal subroutine call returned Error
20181030 085616.154 ERROR            PET0 ESMF_RegridWeightGen.F90:1027 ESMF_RegridWeightGenFile Wrong argument specified  - Internal subroutine call returned Error

This happens in both directions. I checked for missing values in the CISM SCRIP file, but everything looks okay - I am wondering if having local grid in lat direction and global in lon direction might be a problem, or if the issue is with the Pole option - do you haver any clue why this happens?

Here you can find the message I get during runtime

Creating map_lgm21ka_cism4km_TO_lgm21ka_clm1deg_aave.181030.nc...
regional source
cmd is  ESMF_RegridWeightGen --ignore_unmapped -m conserve -w map_lgm21ka_cism4km_TO_lgm21ka_clm1deg_aave.181030.nc -s SCRIPgrid_lgm21ka_cism.nc -d SCRIPgrid_lgm21ka_clm.nc  --src_regional
   ESMF_VERSION_STRING:       7.0.0

 Starting weight generation with these inputs:
   Source File: SCRIPgrid_lgm21ka_cism.nc
   Destination File: SCRIPgrid_lgm21ka_clm.nc
   Weight File: map_lgm21ka_cism4km_TO_lgm21ka_clm1deg_aave.181030.nc
   Source File is in SCRIP format
   Source Grid is a regional grid
   Source Grid is a logically rectangular grid
   Use the center coordinates of the source grid to do the regrid
   Destination File is in SCRIP format
   Destination Grid is a global grid
   Destination Grid is a logically rectangular grid
   Use the center coordinates of the destination grid to do the regrid
   Regrid Method: conserve
   Pole option: NONE
   Ignore unmapped destination points
   Line Type: greatcircle
   Norm Type: dstarea



 ERROR: Problem on processor            0
 . Please see the PET*.RegridWeightGen.Log files for a traceback.
Error status returned from create_ESMF_map script

thanks!

Michele

******************************************************************
Michele Petrini, PhD

Department of Geoscience and Remote Sensing, TUDelft
Delft, The Netherlands

M.Petrini at tudelft.nl<mailto:M.Petrini at tudelft.nl>
+39 3398367372
michele.petrins (Skype alias)



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181030/25485557/attachment.html>