  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 |
1.8.1
