Adding new surface data error: "gridmap_check ERROR from mksectorwater: mapping areas not conserved"

[sabin_taranu]

Hi everybody,

While trying to add new surface data I run into the following issue:
gridmap_check ERROR from mksectorwater: mapping areas not conserved
global sum output field = 0.0000000000E+00
global sum input field = -0.5100487376E+13

Does anyone know what can be the problem and possible solutions?

To reproduce the error:
git clone GitHub - TaranuDev/CTSM: Community Terrestrial Systems Model (includes the Community Land Model of CESM) CTSM_input_test_ST
cd CTSM_input_test_ST
git checkout read_input_data2
./manage_externals/checkout_externals
cd tools/mksurfdata_map/src
gmake
cd ..
./mksurfdata.pl -r 0.9x1.25 -y 2000

In addition to see all files I modified (assuming you are in CTSM_test_ST/):
grep -crl "sectorwater"

 

[oleson]

Since the global sum input field is a large negative number, I suspect that your input mask is all missing values.
I think you need a LANDMASK variable on your input file.

 

[sabin_taranu]

Hi Keith, you were right!
I added the LANDMASK variable and it worked like magic ;)

Thanks for your help on previous question too!

Have a nice day!
Sabin