From greve at lowtem.hokudai.ac.jp Fri Mar 12 15:57:22 2021 From: greve at lowtem.hokudai.ac.jp (Ralf Greve) Date: Fri, 12 Mar 2021 15:57:22 +0900 Subject: [SICOPOLIS:112] Reading 2D input (topography etc.) from NetCDF files Message-ID: <511d84f1-2939-5594-247e-75a7f885a0fd@lowtem.hokudai.ac.jp> Dear all, during the last days, I made a number of commits that allow providing 2D (plan-view) input, such as topography, surface temperature or precipitation, in NetCDF files. Previously, only ASCII was supported. The file types are recognized automatically by the respective extension ('nc' for NetCDF, otherwise ASCII is assumed), so that no further settings are needed. I hope you'll find this new feature useful. 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 From greve at lowtem.hokudai.ac.jp Mon Mar 22 10:23:44 2021 From: greve at lowtem.hokudai.ac.jp (Ralf Greve) Date: Mon, 22 Mar 2021 10:23:44 +0900 Subject: [SICOPOLIS:113] Regional time-series (scalar) output Message-ID: 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 From greve at lowtem.hokudai.ac.jp Mon Mar 22 16:22:48 2021 From: greve at lowtem.hokudai.ac.jp (Ralf Greve) Date: Mon, 22 Mar 2021 16:22:48 +0900 Subject: [SICOPOLIS:114] Region mask files for Greenland Message-ID: Dear all, I also made available region mask files in NetCDF format for Greenland. Similar to Antarctica, the required variable 'mask_region' is contained in the latest versions of the IMBIE basin files grl{05,10,20,40}_imbie2016_basins_extrapolated.nc. To download them, run the script get_input_files.sh that is part of the SICOPOLIS package. The regions defined are the six IMBIE basins NO, NE, SE, SW, CW and NW (, scroll down to "Rignot Basins"). So, to create regional scalar output for these six regions, get the latest revision of SICOPOLIS and define #define MASK_REGION_FILE 'grlxx_imbie2016_basins_extrapolated.nc' (xx to be replaced by the chosen resolution) in the run-specs header. The rest will happen automatically. 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