 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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the ANT domain) |
   sub_ice_shelf_melting_param.F90 | Sub-ice-shelf melting parameterisation for Antarctica |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the ASF domain) |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the EMTP2SGE domain) |
   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 |
   calc_dxyz.F90 | Computation of all components of the strain-rate tensor, the full effective strain rate and the shear fraction |
   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_gia.F90 | Computation of the glacial isostatic adjustment of the lithosphere surface |
   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_const.F90 | Isothermal mode: Setting of the temperature and age to constant values |
   calc_temp_enth.F90 | Computation of temperature, water content and age with the enthalpy method |
   calc_temp_enth_1.F90 | Computation of temperature and age for a cold ice column with the enthalpy method |
   calc_temp_enth_2.F90 | Computation of temperature and age for an ice column with a temperate base with the enthalpy method |
   calc_temp_enth_ssa.F90 | Computation of temperature and age for ice shelves (floating ice) with the enthalpy method |
   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_1.F90 | Computation of the ice topography (including gradients) for SIA-only dynamics of ice sheets without ice shelves |
   calc_top_2.F90 | Computation of the ice topography (including gradients) for hybrid SIA/SStA dynamics of ice sheets without ice shelves |
   calc_top_3.F90 | Computation of the ice topography (including gradients) for coupled SIA/SSA or SIA/SStA/SSA dynamics of ice sheets with ice shelves |
   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_vxy_static.F90 | Computation of the horizontal velocity vx, vy, the horizontal volume flux qx, qy etc. for static ice |
   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_vz_static.F90 | Computation of the vertical velocity vz for static ice |
   calc_water_bas.F90 | Computation of the thickness of the water column under the ice base |
   creep.F90 | Creep response function for ice |
   enth_temp_omega.F90 | Conversion from temperature (temp) and water content (omega) to enthalpy (enth) and vice versa |
   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 |
   ice_material_quantities.F90 | Material quantities of ice: Rate factor, heat conductivity, specific heat (heat capacity) |
   init_temp_water_age.F90 | Initial temperature, water content and age |
   make_zl0.F90 | Generation of an isostatically relaxed lithosphere surface topography for either the rigid lithosphere, the local lithosphere or the elastic lithosphere model (depending on the setting of the parameter REBOUND). This routine is not to be used regularly, and it is only executed if the parameter EXEC_MAKE_ZL0 is defined in the header file |
   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 |
   nc_check.F90 | NetCDF error capturing |
   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) |
   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_sle_solvers.F90 | Solvers for systems of linear equations used by SICOPOLIS |
   sico_types.F90 | Declarations of kind types for SICOPOLIS |
   sico_variables.F90 | Declarations of global variables for SICOPOLIS |
   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 |
   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) or both (for the enthalpy method) |
  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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the GRL domain) |
   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) |
  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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the NHEM domain) |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the NMARS domain) |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the SCAND domain) |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the SMARS domain) |
   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 |
   sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the TIBET domain) |
   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) |
  xyz | |
   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 |
    sico_vars.F90 | Declarations of global variables for SICOPOLIS (for the HEINO domain) |
    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) |
 sicopolis.F90 | Main program file of SICOPOLIS |