controlled by the new preprocessor switch NETCDF4_ENABLED (0/1).
Default is NetCDF-3 Classic as before.

relative to xmin, xmax, ymin, ymax of the computational domain.
If a defined ice core is out of the domain, execution does not stop
any more, only a warning is written.

now only performed for physical-parameter files in ASCII format.

netCDF commands to facilitate error detection.

sico_specs_v5_grl20_ss25ka.h replaced by sico_specs_v5_grl16_bm4_ss25ka.h
(SIA, 16 km resolution, using the new BedMachine V4 topography).
Counterpart for Antarctica renamed
sico_specs_v5_ant40_ss25ka.h -> sico_specs_v5_ant40_b2_ss25ka.h.

changed from ASCII (not available any more in the repository) to NetCDF.

(max./min. number of nonlinear iterations for the SSA/SStA solver)
changed from (10,1) to (25,2).

Settings of TOL_ITER_SSA and N_ITER_SSA changed.

minimum number of iterations N_ITER_SSA_MIN (default 1).

in the run-specs headers.

Parameter CHECK_RES_IMAX_JMAX mentioned under "Domain XYZ".

If undefined, the default value 1 a = 31556925.445 s
(IUPAC-IUGS year for epoch 2000.0) is used.
Value stored in global variable year2sec (previously a parameter).
New global variable sec2year = 1/year2sec introduced.

added for extracting the climatology and reference thickness
from time-slice output files.

If 1 or undefined, compatibility check between horizontal resolution
and number of grid points carried out as before.
If 0, compatibility check not carried out.

Instead, the name of the simulation is extracted directly from the
file name of the run-specs header.
New global variable run_name introduced.
Length of run_name and anfdatname changed from 100 to 256.

in addition to the ASCII file for the isostatically relaxed
bedrock topography.

changed for better comprehensibility: Load due to ice and sea water
(load_ice_water) and steady-state displacement of the lithosphere
(wss) calculated explicitly for all cases.

New horizontal resolutions 6.25 km and 12.5 km implemented.

Script get_input_files.sh updated accordingly.

It now downloads the topography files for the new, combined
EISMINT domain.

and EMTSHELF (modified EISMINT, with rectangular ocean basin,
under XYZ) combined into single domain EISMINT.
Topographies are read from files rather than computed within the
module sico_init_m.

If zl == z_sl (bedrock and sea level exactly equal), ice-free land
or ocean points always remain ice-free land or ocean points.

It now downloads also the new topography files for the
EISMINT-Shelf domain.

Run-specs header for simulation v5_emtshelf25_expH updated.

such that the topography is now read from files rather than computed
within the module. Allows more flexibility!

Coments for options TSURFACE==6, ACCSURFACE==6 and ABLSURFACE==6
updated.

* Possibility added to prescribe spatio-temporally constant values for
the surface-temperature vertical gradient (topographic lapse rate)
and the SMB vertical gradient.
* If the SMB anomaly and the SMB vertical gradient are read from
netCDF files, two different units are now allowed:
kg/(m2*s) or m/a ice equiv., [kg/(m2*s)]/m or [m/a ice equiv.]/m, resp.

changed to make it more robust.

that variable 'year2sec' is not present in SICOPOLIS output files.

Variable dlimdt (total ice mass change) added to output.