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

[Ralf Greve]

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 
allowed.

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 
questions.

Best regards,

Ralf

-- 
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
(+81)-(0)11-706-6891
greve at lowtem.hokudai.ac.jp
http://wwwice.lowtem.hokudai.ac.jp/~greve/
https://twitter.com/RalfGreve