subroutines | |
ant | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
read_target_topo_nc.F90 | Reading of the target-topography file (in NetCDF format) |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
asf | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
output5.F90 | Writing of time-series data for all defined surface points on file in ASCII format. Modification of Tolly's output7 by Thorben Dunse |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
emtp2sge | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
general | |
borehole.F90 | Computation of an arbitrary field quantity for a given borehole position x_pos, y_pos by weighed averaging of the corresponding gridded 2-d field |
c_val.F90 | Specific heat of ice: Linear interpolation of tabulated values in C(.) |
calc_age_c.F90 | Computation of age for a cold ice column |
calc_age_t.F90 | Computation of age for a cold / temperate ice column |
calc_elra.F90 | Computation of the isostatic steady-state displacement of the lithosphere (wss) for the ELRA model |
calc_enhance.F90 | Computation of the flow enhancement factor |
calc_qbm.F90 | Computation of the basal melting rate Q_bm. Summation of Q_bm and Q_tld |
calc_temp1.F90 | Computation of temperature and age for a cold ice column |
calc_temp2.F90 | Computation of temperature and age for an ice column with a temperate base overlain by cold ice |
calc_temp3.F90 | Computation of temperature, water content and age for an ice column with a temperate base overlain by a temperate-ice layer |
calc_temp_bas.F90 | Computation of the basal temperatures |
calc_temp_cold.F90 | Computation of temperature and age in cold-ice mode |
calc_temp_melt.F90 | Computation of the melting temperatures |
calc_temp_poly.F90 | Computation of temperature, water content and age in polythermal mode |
calc_temp_r.F90 | Computation of temperature for an ice-free column |
calc_temp_ssa.F90 | Computation of temperature and age for ice shelves (floating ice) |
calc_top.F90 | Computation of the topography (including gradients) |
calc_vis_ssa.F90 | Computation of the depth-integrated viscosity vis_int_g in the shallow shelf approximation |
calc_vxy_b_sia.F90 | Computation of the basal horizontal velocity vx_b, vy_b in the shallow ice approximation |
calc_vxy_sia.F90 | Computation of the shear stresses txz, tyz, the effective shear stress sigma, the depth-averaged fluidity flui_ave_sia, the horizontal velocity vx, vy and the horizontal volume flux qx, qy in the shallow ice approximation |
calc_vxy_ssa.F90 | Computation of the horizontal velocity vx, vy, the horizontal volume flux qx, qy and the flux across the grounding line q_gl_g in the shallow shelf approximation |
calc_vxy_ssa_matrix.F90 | Solution of the system of linear equations for the horizontal velocities vx_m, vy_m in the shallow shelf approximation |
calc_vz_sia.F90 | Computation of the vertical velocity vz in the shallow ice approximation |
calc_vz_ssa.F90 | Computation of the vertical velocity vz in the shallow shelf approximation |
calc_water_bas.F90 | Computation of the effective thickness of subglacial water |
check.F90 | NetCDF error capturing |
creep.F90 | Creep response function for ice |
erfcc.F90 | Computation of the complementary error function erfc(x) = 1-erf(x) with a fractional error everywhere less than 1.2 x 10^(-7) (formula by Press et al., 'Numerical Recipes in Fortran 77') |
geo_coord.F90 | Computation of longitude lambda and latitude phi for position (x,y) in the numerical domain |
kappa_val.F90 | Heat conductivity of ice: Linear interpolation of tabulated values in KAPPA(.) |
mask_update.F90 | Update of the topography mask due to changes of the sea level |
metric.F90 | Definition of the components g11 and g22 of the metric tensor of the applied coordinates |
num_coord.F90 | Computation of position (x,y) in the numerical domain for longitude lambda and latitude phi (in rad) |
output1.F90 | Writing of time-slice files in binary format |
output2.F90 | Writing of time-series data on file in ASCII format |
output3.F90 | Writing of time-series data of the ice thickness field on file in ASCII format |
output4.F90 | Writing of time-series data of the deep ice cores on file in ASCII format |
output_nc.F90 | Writing of time-slice files in NetCDF format |
pdd.F90 | Computation of the positive degree days (PDD) with statistical temperature fluctuations; based on semi-analytical solution by Calov and Greve (2005). Note that this subroutine uses years as time unit (as opposed to the seconds which are usually used by SICOPOLIS) |
ratefac.F90 | Rate factor for cold ice: Linear interpolation of tabulated values in RF(.) |
ratefac_t.F90 | Rate factor for temperate ice |
read_kei.F90 | Reading of the tabulated kei function |
read_phys_para_value.F90 | Reading of a value of a physical parameter from the phys_para file |
shift_cts_downward.F90 | Downward shifting of the CTS |
shift_cts_upward.F90 | Upward shifting of the CTS |
sico_end.F90 | Ending of SICOPOLIS |
sico_main_loop.F90 | Main loop of SICOPOLIS |
sico_types.F90 | Declarations of kind types for SICOPOLIS |
sico_variables.F90 | Declarations of global variables for SICOPOLIS |
sor_sprs.F90 | SOR solver for a system of linear equations lgs_a*lgs_x=lgs_b [matrix in compressed row storage CRS, represented by arrays lgs_a_value(values), lgs_a_index (indices) and lgs_a_ptr (pointers)] |
topograd_1.F90 | Calculation of topography gradients on the staggered grid and on the grid points (the latter by second-order discretization). Length rescaling with the corresponding components of the metric tensor is already included |
topograd_2.F90 | Calculation of topography gradients on the staggered grid and on the grid points (the latter by fourth-order discretization). Length rescaling with the corresponding components of the metric tensor is already included |
tri_sle.F90 | Solution of a system of linear equations Ax=b with tridiagonal matrix A |
viscosity.F90 | Ice viscosity as a function of the effective strain rate and the temperature (in cold ice) or the water content (in temperate ice) |
grl | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
read_target_topo_nc.F90 | Reading of the target-topography file (in NetCDF format) |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
heino | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For an initial topography with a thin ice layer (thickness 2*H_MIN) everywhere on the land area |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
nhem | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
nmars | |
boundary.F90 | Mars Atmosphere-Ice Coupler MAIC-1.5: Computation of the surface temperature (must be less than 0 deg C) and of the accumulation-ablation rate for the north polar cap of Mars. Computation of the geothermal heat flux |
instemp.f90 | Computation of the daily mean surface temperature of Mars based on obliquity, eccentricity and the anomaly of vernal equinox (local insolation temperature scheme = LIS scheme) |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
scand | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
smars | |
boundary.F90 | Mars Atmosphere-Ice Coupler MAIC-1.5: Computation of the surface temperature (must be less than 0 deg C) and of the accumulation-ablation rate for the south polar cap of Mars. Computation of the geothermal heat flux |
instemp.f90 | Computation of the daily mean surface temperature of Mars based on obliquity, eccentricity and the anomaly of vernal equinox (local insolation temperature scheme = LIS scheme) |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography1.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For present-day initial topography |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
tibet | |
boundary.F90 | Computation of the surface temperature (must be less than 0 deg C!) and of the accumulation-ablation function |
phys_para.F90 | Reading of physical parameters |
sico_init.F90 | Initialisations for SICOPOLIS |
topography2.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For ice-free initial topography with relaxed lithosphere |
topography3.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in binary format) |
topography3_nc.F90 | Definition of the initial surface and bedrock topography (including gradients) and of the horizontal grid spacings dxi, deta. For initial topography from previous simulation (read in NetCDF format) |
sicopolis.F90 | Main program file of SICOPOLIS |