[SICOPOLIS:113] Regional time-series (scalar) output

Ralf Greve greve at lowtem.hokudai.ac.jp
Mon Mar 22 10:23:44 JST 2021

Dear all,

during the last days, I implemented a new function that allows writing 
regional (rather than only ice-sheet-wide) scalar output to NetCDF 
files. To do so, a region mask must be defined that marks the regions by 
integers. For instance, for Antarctica, this may be 1 for East 
Antarctica, 2 for West Antarctica and 3 for the Peninsula. The region 
mask file can be provided either in ASCII or NetCDF format, and its name 
is defined in the run-specs header by the parameter MASK_REGION_FILE. 
For the case of three regions, scalar output for these regions is 
written to files runname_ser_region{01,02,03}.nc. Up to 99 regions are 

For Antarctica, the latest version of the input-file repository contains 
the region mask as an additional variable 'mask_region' in the NetCDF 
versions of the IMBIE basin files 
ant{08,10,16,20,32,40,64}_imbie2016_basins_extrapolated.nc (run the 
script get_input_files.sh to download the update). So, if you define one 
of these files (depending on resolution) by MASK_REGION_FILE, the 
regions will be defined and the corresponding scalar output produced.

Actually, for Antarctica, even without specifying MASK_REGION_FILE, 
three regions will be defined automatically, namely East Antarctica, 
West Antarctica and the Peninsula. However, they are by bounded by 
latitudes and longitudes, which is a bit rough. Therefore, it is better 
to define the regions by the IMBIE basin files as explained above.

I hope you'll find this useful. Please let me know if you have any 

Best regards,


Dr. Ralf Greve
Professor, Glacier and Ice Sheet Research
Institute of Low Temperature Science (ILTS)
Hokkaido University
Kita-19, Nishi-8, Kita-ku
Sapporo 060-0819, Japan
greve at lowtem.hokudai.ac.jp

More information about the SICOPOLIS mailing list