6 #define DATE '2013-05-19'
10 !! @section Description
12 !! SICOPOLIS (SImulation COde for POLythermal Ice Sheets) is a 3-d
13 !! dynamic/thermodynamic model which simulates the evolution of large ice
14 !! sheets. It was originally created as a part of the doctoral thesis by
15 !! Greve (1995) in a version for the Greenland ice sheet. Since then,
16 !! SICOPOLIS has been developed continuously and applied to problems of
17 !! past, present and future glaciation of Greenland, Antarctica, the
18 !! entire northern hemisphere and also the polar ice caps of the planet
21 !! The model is based on the shallow ice approximation for grounded ice
22 !! and the shallow shelf approximation for floating ice (e.g., Greve and
23 !! Blatter 2009). It is coded in Fortran 90 and uses finite-difference
24 !! discretization on a staggered grid, the velocity components being
25 !! taken between grid points. Its particularity is the detailed treatment
26 !! of basal temperate layers (that is, regions with a temperature at the
27 !! pressure melting point), which are positioned by fulfilling a
28 !! Stefan-type jump condition at the interface to the cold ice
29 !! regions. Within the temperate layers, the water content is computed,
30 !! and its influence on the ice viscosity is taken into account.
32 !! Required model forcing:
33 !! @li Surface mass balance (precipitation, evaporation, runoff).
34 !! @li Mean annual air temperature above the ice.
35 !! @li Eustatic sea level.
36 !! @li Geothermal heat flux.
38 !! Output (as functions of position and time):
39 !! @li Extent and thickness of the ice sheet.
40 !! @li Velocity field.
41 !! @li Temperature field.
42 !! @li Water-content field (temperate regions).
43 !! @li Age of the ice.
44 !! @li Isostatic displacement and temperature of the lithosphere.\n
47 !! @li Greve, R. 1995.\n
48 !! Thermomechanisches Verhalten polythermer Eisschilde - Theorie, Analytik,
50 !! Doctoral thesis, Department of Mechanics, Darmstadt University of
51 !! Technology, Germany, 226 pp.\n
52 !! ISBN: 3-8265-0999-4.
53 !! @li Greve, R. and H. Blatter. 2009.\n
54 !! Dynamics of Ice Sheets and Glaciers.\n
55 !! Springer, Berlin, Germany etc., 287 pp. ISBN: 978-3-642-03414-5.
56 !! @li SICOPOLIS website: http://sicopolis.greveweb.net/
60 !! Copyright 2009-2013 Ralf Greve\n
61 !! (with contributions by Reinhard Calov, Thorben Dunse, Thomas Goelles,
62 !! Philipp Hancke, Nina Kirchner, Sascha Knell, Alex Robinson, Tatsuru Sato,
67 !! SICOPOLIS is free software: you can redistribute it and/or modify
68 !! it under the terms of the GNU General Public License as published by
69 !! the Free Software Foundation, either version 3 of the License, or
70 !! (at your option) any later version.
72 !! SICOPOLIS is distributed in the hope that it will be useful,
73 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
74 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
75 !! GNU General Public License for more details.
77 !! You should have received a copy of the GNU General Public License
78 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
82 !! Main program file of SICOPOLIS.
86 !! Copyright 2009-2013 Ralf Greve
90 !! This file is part of SICOPOLIS.
92 !! SICOPOLIS is free software: you can redistribute it and/or modify
93 !! it under the terms of the GNU General Public License as published by
94 !! the Free Software Foundation, either version 3 of the License, or
95 !! (at your option) any later version.
97 !! SICOPOLIS is distributed in the hope that it will be useful,
98 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
99 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
100 !! GNU General Public License for more details.
102 !! You should have received a copy of the GNU General Public License
103 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
106 ! *****************************************************************
108 ! *****************************************************************
109 ! V1.0: Modular version of SICOPOLIS, applicable to any
111 ! -----------------------------------------------------------
112 ! Based on sicopolis_nhem.F95 V4.0
113 ! (northern-hemisphere version of SICOPOLIS).
114 ! -----------------------------------------------------------------
115 ! V1.1: Optional new boundary conditions for the ice
116 ! surface at the margin of the computational domain
117 ! (parameter BC_MARG).
118 ! -----------------------------------------------------------------
119 ! V1.2: Optional first-order-upstream discretization
120 ! of vertical advection terms without artificial
121 ! diffusion in the age equation
122 ! (parameter ADV_VERT).
123 ! -----------------------------------------------------------
124 ! Optional exponential coupling between accumulation
125 ! and delta_ts (parameter ACCSURFACE==3).
126 ! -----------------------------------------------------------------
127 ! V1.3: New module for the southern hemisphere included
128 ! (parameter DOMAIN==4).
129 ! -----------------------------------------------------------
130 ! Improved reading of ice-core-temperature and
131 ! sea-level data, no longer restricted to the
132 ! interval from 250 kyr BP to today.
133 ! -----------------------------------------------------------
134 ! New model domain for Scandinavia included
135 ! (parameter DOMAIN==2).
136 ! -----------------------------------------------------------------
137 ! V1.4: New module for Antarctica included.
138 ! -----------------------------------------------------------
139 ! Bugfix in subroutines sico_init
140 ! (for grl, shem, nmars):
141 ! Parameter list in calls of subroutine boundary
142 ! for parameter ANF_DAT==3 corrected.
143 ! -----------------------------------------------------------
144 ! Mars module: Degree-day model replaced by simpler
145 ! accumulation-ablation parameterization.
146 ! -----------------------------------------------------------------
147 ! V2.0: All file extensions changed from '.F95' to '.F90'.
148 ! -----------------------------------------------------------
149 ! COMMON blocks in input file 'sico.common' replaced
150 ! by module 'sico_variables'.
151 ! -----------------------------------------------------------
152 ! Character variable 'runname' (name of current run)
153 ! can have an arbitrary length of less than
155 ! -----------------------------------------------------------
156 ! Variables d_help_c, d_help_t, h_diff,
157 ! maske, maske_help, n_cts, n_cts_neu, kc_cts,
159 ! now declared globally in module 'sico_variables'.
160 ! -----------------------------------------------------------
161 ! The header file 'sicoheader.F95' does no longer
162 ! exist; however, 'sicopolis.F90' is compiled
164 ! -----------------------------------------------------------
165 ! Header file with simulation specifications must
166 ! be called 'sico_specs.h'.
167 ! -----------------------------------------------------------
168 ! Explicit definition of paths INPATH (for input
169 ! files), OUTPATH (for output files), ANFDATPATH
170 ! (for initial-value file) and Z_SLE_PATH (for
171 ! previously computed time-series file with
172 ! sea-level equivalents) in specifications file.
173 ! -----------------------------------------------------------
174 ! Parameter DOMAIN removed. Instead, the
175 ! computational domain is now defined by a macro
176 ! ANT, GRL, SCAND, NHEM, SHEM, EMTP2SGE or NMARS.
177 ! -----------------------------------------------------------
178 ! Parameter BC_MARG (see version 1.1) removed, back
179 ! to zero-thickness boundary condition at the margin
180 ! of the numerical domain as the only possibility.
181 ! -----------------------------------------------------------
182 ! Geothermal heat flux (Q_GEO), flow-enhancement factor
183 ! (ENH) and minimum ice thickness interpreted as glaciation
184 ! (H_MIN) now defined in header file sico_specs.h.
185 ! -----------------------------------------------------------
186 ! Physical parameters now read from a file by the
187 ! new subroutine 'phys_para'.
188 ! -----------------------------------------------------------
189 ! Indexing of arrays RF(T), KAPPA(T) and C(T)
190 ! changed such that the index T now corresponds to
191 ! the temperature in deg C.
192 ! -----------------------------------------------------------
193 ! Computation of the horizontal volume flux qx, qy
194 ! moved from subroutine 'calc_top' to subroutine
196 ! -----------------------------------------------------------
197 ! Forcing quantities delta_ts and z_sl added to
198 ! the output for the time-slice files produced by
199 ! subroutine 'output1'.
200 ! -----------------------------------------------------------
201 ! Bugfix in subroutine 'calc_top':
202 ! In the diffusive smoothing scheme for zs_neu, the
203 ! denominator of the correction factor 'korrfakt'
204 ! can no longer become zero due to the insertion
205 ! of a volume offset eps_vol.
206 ! -----------------------------------------------------------------
207 ! V2.0.1: Maximum surface velocity now written to the time-series
209 ! -----------------------------------------------------------------
210 ! V2.0.2: Module for the southern hemisphere deleted.
211 ! ---------------------------------------------------------
212 ! New Greenland topography data based on Letreguilly and
213 ! ETOPO5 with slightly smoothed relaxed bedrock used.
214 ! -----------------------------------------------------------------
215 ! V2.0.3: Limitation of the computed horizontal velocities to an
216 ! interval [-VH_MAX, VH_MAX] to be prescribed in the
218 ! -----------------------------------------------------------------
219 ! V2.0.4: Implementation of alternative flow laws
220 ! (Goldsby-Kohlstedt, Durham) in addition to the usual
221 ! Glen's flow law (parameter FLOW_LAW).
222 ! ---------------------------------------------------------
223 ! For the north polar cap of Mars, the average dust content
224 ! of the ice can be specified (parameter FRAC_DUST).
225 ! -----------------------------------------------------------------
226 ! V2.1: Spatially variable geothermal heat flux implemented, can
227 ! be chosen by new parameter Q_GEO_MOD.
228 ! -----------------------------------------------------------
229 ! Output of time-series file for ice-core locations
230 ! prescribed in subroutine sico_init implemented
231 ! (parameter OUTSER==3).
232 ! -----------------------------------------------------------
233 ! New module for the south-polar cap of Mars included
235 ! -----------------------------------------------------------
236 ! For the Martian ice caps (NMARS, SMARS), the mass balance
237 ! is now converted from water equiv. to (ice+dust) equiv.
238 ! (previously the dust was ignored).
239 ! -----------------------------------------------------------------
240 ! V2.2: Sub-melt basal sliding implemented (parameter GAMMA_SLIDE
241 ! in phys_para file; 0: no sub-melt sliding).
242 ! -----------------------------------------------------------
244 ! - New Bamber topography implemented, horizontal resolution
246 ! - Surface temperatures (mean annual, mean July)
247 ! parameterised by Ritz et al. (1997).
248 ! - Precipitation instead of solid accumulation used as
250 ! - Reeh's degree day model with statistical fluctuations,
251 ! supplemented by explicit consideration of rainfall
252 ! and different ice/snow degree-day factors for warm
253 ! and cold conditions (parameters BETA1_LT_0, BETA1_HT_0,
254 ! BETA2_LT_0, BETA2_HT_0 in phys_para file).
255 ! - Net surface mass balance computed as precipitation minus
256 ! runoff (previously solid accumulation minus melting).
257 ! - Optional climate forcing by glacial index and prescribed
258 ! LGM anomalies for surface temperature and precipitation
259 ! (parameters TSURFACE==5 and ACCSURFACE==5).
260 ! -----------------------------------------------------------
261 ! Martian ice caps (NMARS, SMARS):
262 ! - Accumulation input either in water equiv.
263 ! or in (ice+dust) equiv. (parameter ACC_UNIT).
264 ! - Lower limit of -125/-128 deg C (sublimation temperature
265 ! of CO2 over the NPC/SPC) introduced for the ice surface
267 ! -----------------------------------------------------------------
268 ! V2.3: Implementation of Bjoern Grieger's surface-temperature
269 ! parameterisation for the Martian ice caps (parameter
271 ! -----------------------------------------------------------
272 ! Greenland: Optional linear or exponential interpolation
273 ! function for glacial-index-based precipitation forcing
274 ! (parameter PRECIP_ANOM_INTERPOL).
275 ! -----------------------------------------------------------------
276 ! V2.4: Implementation of an elastic lithosphere (ELRA model),
277 ! based on Sascha Knell's diploma thesis
278 ! (parameter REBOUND==2)
279 ! -----------------------------------------------------------
280 ! Implementation of Bernd Muegge's conservative first-order
281 ! upstream scheme for the age equation
282 ! (parameters ADV_HOR==3 and ADV_VERT==3).
283 ! -----------------------------------------------------------
284 ! Antarctica, Scandinavia (Greenland see V2.2):
285 ! - Precipitation instead of solid accumulation used as
287 ! - Reeh's degree day model with statistical fluctuations,
288 ! supplemented by explicit consideration of rainfall.
289 ! - Net surface mass balance computed as precipitation minus
290 ! runoff (previously solid accumulation minus melting).
291 ! - Optional climate forcing by glacial index and prescribed
292 ! LGM anomalies for surface temperature and precipitation
293 ! (parameters TSURFACE==5 and ACCSURFACE==5).
294 ! -----------------------------------------------------------
295 ! Northern-hemisphere module no longer up-to-date, therefore
296 ! deleted for the time being.
297 ! -----------------------------------------------------------
298 ! Compatibility check between version number of SICOPOLIS
299 ! and version number in header file.
300 ! -----------------------------------------------------------------
301 ! V2.4.1: Antarctica, Scandinavia, Greenland:
302 ! - Optional computation of surface melting (ablation)
303 ! by the positive-degree-day (PDD) method
304 ! or the linear-temperature-index (LTI) method
305 ! (parameter ABLSURFACE).
306 ! - Computation of positive degree days (subroutine pdd)
307 ! changed from annual sinus cycle of the surface
308 ! temperature to direct use of monthly-mean temperatures.
309 ! ---------------------------------------------------------
311 ! - Climate input (surface temperature, precipitation)
312 ! changed from seasonal to monthly-mean values.
313 ! - Optional horizontal resolutions 40, 20 or 10 km.
314 ! -----------------------------------------------------------------
315 ! V2.5: Basal sliding now computed in new subroutine calc_vxy_b.
316 ! Parameters of sliding law are set in the header file
317 ! (previously in phys_para file). New options: Weertman
318 ! exponents can be chosen freely, alternative use of
319 ! full or reduced pressure in sliding law.
320 ! -----------------------------------------------------------
321 ! Antarctica, Greenland, Scandinavia:
322 ! Two variants of PDD, either with instantaneous runoff of
323 ! rainfall (ABLSURFACE==1), or with contribution of rainfall
324 ! to formation of superimposed ice (ABLSURFACE==2).
325 ! Linear-temperature-index (LTI) method is now chosen by
327 ! -----------------------------------------------------------
328 ! North- and south-polar cap of Mars:
329 ! In addition to the density and the heat conductivity,
330 ! now the specific heat is also computed as a
331 ! volume-fraction-weighed mean of the ice and dust values,
332 ! respectively. Function c_val changed accordingly.
333 ! -----------------------------------------------------------
334 ! EISMINT Phase 2 Simplified Geometry Experiment:
335 ! Input file surf_para_emtp2sge.dat and subroutine surf_para
336 ! deleted. Parameters now defined in header file.
337 ! -----------------------------------------------------------
338 ! Minimum bedrock elevation for extent of marine ice
339 ! defined by new variable z_mar. Can be chosen as a constant
340 ! (MARGIN==2, old option), as proportional to the sea-level
341 ! stand (MARGIN==3, new option) or as a piecewise linear
342 ! function of sea level (MARGIN==4, new option; as proposed
343 ! by Zweck and Huybrechts, 2005).
344 ! -----------------------------------------------------------------
345 ! V2.6: New module for ISMIP HEINO included (defined by macro
347 ! -----------------------------------------------------------
348 ! Additional output of the basal frictional heating in the
349 ! time-series file for the ice-core locations.
350 ! -----------------------------------------------------------
351 ! All constants throughout the program defined consistently
352 ! as double-precision constants.
353 ! -----------------------------------------------------------
354 ! New global parameters pi, pi_inv, pi_180, pi_180_inv, eps.
355 ! -----------------------------------------------------------------
356 ! V2.7: Special treatment for the chasms of the Martian NPC (Chasma
357 ! Borealis) and SPC (Chasma Australis) implemented (optional,
358 ! parameter CHASM==2). Chasm regions defined in additional
359 ! mask file -> larger geothermal heat flux and wind erosion
360 ! rate active within prescribed time interval (specified in
362 ! -----------------------------------------------------------
363 ! ISMIP HEINO: Time-series output for the sediment region
364 ! ("Hudson Bay, Hudson Strait") in new file runname.sed.
365 ! -----------------------------------------------------------
366 ! Improved restart: Time derivatives dzs_dtau, dzm_dtau,
367 ! dzb_dtau, dH_c_dtau and dH_t_dtau written to time-slice
368 ! files and read by restart routine topography3. Further,
369 ! Q_b_tot is computed as the sum of Q_bm and Q_tld in
370 ! initialization routine sico_init.
371 ! -----------------------------------------------------------
372 ! In addition to Q_b_tot (formerly called Q_mean), now also
373 ! limitations of Q_bm and Q_tld with the same lower and upper
374 ! limits, QB_MIN and QB_MAX, respectively.
375 ! -----------------------------------------------------------
376 ! In subroutine output2 Reinhard Calov's method for computing
377 ! the ice equivalent of the water volume below the sea level
378 ! (V_ges_redu) implemented.
379 ! -----------------------------------------------------------------
380 ! V2.7.1: Martian ice caps (NMARS, SMARS):
381 ! New topography inputs based on the MOLA MEGDR data with
382 ! 64 pix/deg resolution. Optional horizontal resolutions
384 ! -----------------------------------------------------------------
385 ! V2.7.2: x and y coordinates of the origin point (i,j) = (0,0)
386 ! now defined in header file as X0 and Y0.
387 ! ---------------------------------------------------------
388 ! Optional computation of the evolution of the free
389 ! surface as usual (parameter ZS_EVOL==1),
390 ! or keeping the free surface fixed on the initial
391 ! topography (parameter ZS_EVOL==0).
392 ! ---------------------------------------------------------
394 ! Update of the ice-core positions according to Dahl-Jensen
395 ! (pers. comm. 2006) and inclusion of the planned NEEM
397 ! ---------------------------------------------------------
398 ! Martian ice caps (NMARS, SMARS):
399 ! New version of module instemp for the LIT scheme
401 ! -----------------------------------------------------------------
403 ! Optional special sliding law for the North-East Greenland
404 ! Ice Stream NEGIS (parameter ICE_STREAM==2).
405 ! Surface meltwater effect on basal sliding implemented.
406 ! -----------------------------------------------------------
408 ! Optional 25-km resolution implemented (in addition to the
409 ! standard resolution of 50 km).
410 ! -----------------------------------------------------------
412 ! Re-formatting and re-naming of the input data (now same
413 ! style and units like for Greenland). Topography input newly
414 ! created from RAMPDEM V2 and BEDMAP data-sets, with an
415 ! additional 20-km resolution implemented.
416 ! -----------------------------------------------------------------
417 ! V2.8.1: Northern-hemisphere module reintroduced, topography
418 ! files newly created for horizontal resolutions of
420 ! ---------------------------------------------------------
421 ! Explicit introduction of the evaporation (variable evap).
422 ! In the modules for the terrestrial ice sheets, evap is
423 ! set to zero, whereas in the modules for the Martian
424 ! ice caps, evap accounts for the entire ablation.
425 ! ---------------------------------------------------------
427 ! Introduction of the option ANF_DAT==1, which defines the
428 ! initial topography as a thin ice layer (with a thickness
429 ! of 2*H_MIN) everywhere on the land area.
430 ! -----------------------------------------------------------------
431 ! V2.9: New module for Tibet included (defined by macro TIBET).
432 ! -----------------------------------------------------------
434 ! Parameterization for the surface meltwater effect on basal
435 ! sliding changed; additional dependence on ice thickness
436 ! and possible nonlinear dependence on runoff.
437 ! -----------------------------------------------------------
438 ! Topography files (for zs, zb, zb0, mask) can now be
439 ! specified in the header file.
440 ! -----------------------------------------------------------------
441 ! V3.0: Handling of ice shelf and marine ice dynamics.
442 ! -----------------------------------------------------------
443 ! SeaRISE set-up implemented.
444 ! -----------------------------------------------------------
445 ! Changing water load now considered for the LLRA model
446 ! (REBOUND==1) in the same way as for the ELRA model.
447 ! -----------------------------------------------------------
448 ! Optional output of time-slice files in NetCDF format
449 ! (controlled by parameter NETCDF).
450 ! -----------------------------------------------------------
451 ! Name of the physical-parameter file is now specified in
453 ! -----------------------------------------------------------
454 ! New module 'sico_types' for declarations of kind types
456 ! -----------------------------------------------------------
457 ! Subroutine 'calc_top':
458 ! Diffusive smoothing scheme for zs_neu deleted.
459 ! -----------------------------------------------------------
460 ! Introduction of module Austfonna (ASF) available at 1, 2
461 ! and 4 km horizontal resolution - Thorben Dunse.
462 ! Output5.F90 writes time-series data at specific surface
464 ! -----------------------------------------------------------
465 ! Stereographic projection changed from a spherical planet
466 ! to an ellipsoidal planet.
467 ! *****************************************************************
469 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
473 #include "sico_specs.h"
477 #include "subroutines/general/sico_types.F90"
478 #include "subroutines/general/sico_variables.F90"
484 #include "subroutines/nmars/instemp.f90"
486 #include "subroutines/smars/instemp.f90"
490 !> Main program of SICOPOLIS.
491 !<------------------------------------------------------------------------------
561 integer(i4b) :: ndat2d, ndat3d
562 integer(i4b) :: n_output
563 integer(i4b),
dimension((IMAX+1)*(JMAX+1)) :: ii, jj
564 integer(i4b),
dimension(0:JMAX,0:IMAX) :: nn
565 real(dp) :: delta_ts, glac_index
566 real(dp) :: sum_accum, mean_accum, mean_accum_inv
567 real(dp) :: dtime, dtime_temp, dtime_wss, dtime_out, dtime_ser
568 real(dp) :: time, time_init, time_end
569 real(dp),
dimension(100) :: time_output
570 real(dp) :: dxi, deta, dzeta_c, dzeta_t, dzeta_r
571 real(dp) :: z_sl, dzsl_dtau, z_mar
572 character (len=100) :: runname
577 sum_accum, mean_accum, mean_accum_inv, &
578 dtime, dtime_temp, dtime_wss, dtime_out, dtime_ser, &
579 time, time_init, time_end, time_output, &
580 dxi, deta, dzeta_c, dzeta_t, dzeta_r, &
581 z_sl, dzsl_dtau, z_mar, &
583 ndat2d, ndat3d, n_output, &
589 sum_accum, mean_accum, mean_accum_inv, &
590 dtime, dtime_temp, dtime_wss, dtime_out, dtime_ser, &
591 time, time_init, time_end, time_output, &
592 dxi, deta, dzeta_c, dzeta_t, dzeta_r, &
593 z_sl, dzsl_dtau, z_mar, &
595 ndat2d, ndat3d, n_output, &
607 #include "subroutines/nhem/sico_init.F90"
609 #include "subroutines/scand/sico_init.F90"
611 #include "subroutines/grl/sico_init.F90"
613 #include "subroutines/nmars/sico_init.F90"
615 #include "subroutines/smars/sico_init.F90"
616 #elif defined(EMTP2SGE)
617 #include "subroutines/emtp2sge/sico_init.F90"
619 #include "subroutines/heino/sico_init.F90"
621 #include "subroutines/ant/sico_init.F90"
623 #include "subroutines/tibet/sico_init.F90"
625 #include "subroutines/asf/sico_init.F90"
628 #include "subroutines/general/sico_main_loop.F90"
629 #include "subroutines/general/sico_end.F90"
632 #include "subroutines/nhem/phys_para.F90"
634 #include "subroutines/scand/phys_para.F90"
636 #include "subroutines/grl/phys_para.F90"
638 #include "subroutines/nmars/phys_para.F90"
640 #include "subroutines/smars/phys_para.F90"
641 #elif defined(EMTP2SGE)
642 #include "subroutines/emtp2sge/phys_para.F90"
644 #include "subroutines/heino/phys_para.F90"
646 #include "subroutines/ant/phys_para.F90"
648 #include "subroutines/tibet/phys_para.F90"
650 #include "subroutines/asf/phys_para.F90"
653 #include "subroutines/general/read_phys_para_value.F90"
654 #include "subroutines/general/read_kei.F90"
657 #include "subroutines/nhem/topography2.F90"
658 #include "subroutines/nhem/topography3.F90"
660 #include "subroutines/nhem/topography3_nc.F90"
663 #include "subroutines/scand/topography2.F90"
664 #include "subroutines/scand/topography3.F90"
666 #include "subroutines/scand/topography3_nc.F90"
669 #include "subroutines/grl/topography1.F90"
670 #include "subroutines/grl/topography2.F90"
671 #include "subroutines/grl/topography3.F90"
673 #include "subroutines/grl/topography3_nc.F90"
675 #include "subroutines/grl/read_target_topo_nc.F90"
679 #include "subroutines/nmars/topography1.F90"
680 #include "subroutines/nmars/topography2.F90"
681 #include "subroutines/nmars/topography3.F90"
683 #include "subroutines/nmars/topography3_nc.F90"
686 #include "subroutines/smars/topography1.F90"
687 #include "subroutines/smars/topography2.F90"
688 #include "subroutines/smars/topography3.F90"
690 #include "subroutines/smars/topography3_nc.F90"
692 #elif defined(EMTP2SGE)
693 #include "subroutines/emtp2sge/topography2.F90"
694 #include "subroutines/emtp2sge/topography3.F90"
696 #include "subroutines/emtp2sge/topography3_nc.F90"
699 #include "subroutines/heino/topography1.F90"
700 #include "subroutines/heino/topography2.F90"
701 #include "subroutines/heino/topography3.F90"
703 #include "subroutines/heino/topography3_nc.F90"
706 #include "subroutines/ant/topography1.F90"
707 #include "subroutines/ant/topography2.F90"
708 #include "subroutines/ant/topography3.F90"
710 #include "subroutines/ant/topography3_nc.F90"
712 #include "subroutines/ant/read_target_topo_nc.F90"
716 #include "subroutines/tibet/topography2.F90"
717 #include "subroutines/tibet/topography3.F90"
719 #include "subroutines/tibet/topography3_nc.F90"
722 #include "subroutines/asf/topography1.F90"
723 #include "subroutines/asf/topography2.F90"
724 #include "subroutines/asf/topography3.F90"
726 #include "subroutines/asf/topography3_nc.F90"
731 #include "subroutines/nhem/boundary.F90"
733 #include "subroutines/scand/boundary.F90"
735 #include "subroutines/grl/boundary.F90"
737 #include "subroutines/nmars/boundary.F90"
739 #include "subroutines/smars/boundary.F90"
740 #elif defined(EMTP2SGE)
741 #include "subroutines/emtp2sge/boundary.F90"
743 #include "subroutines/heino/boundary.F90"
745 #include "subroutines/ant/boundary.F90"
747 #include "subroutines/tibet/boundary.F90"
749 #include "subroutines/asf/boundary.F90"
752 #include "subroutines/general/mask_update.F90"
753 #include "subroutines/general/pdd.F90"
754 #include "subroutines/general/erfcc.F90"
755 #include "subroutines/general/calc_temp_poly.F90"
756 #include "subroutines/general/calc_temp_cold.F90"
757 #include "subroutines/general/calc_temp1.F90"
758 #include "subroutines/general/calc_temp2.F90"
759 #include "subroutines/general/calc_temp3.F90"
760 #include "subroutines/general/calc_temp_r.F90"
761 #include "subroutines/general/shift_cts_upward.F90"
762 #include "subroutines/general/shift_cts_downward.F90"
763 #include "subroutines/general/calc_enhance.F90"
764 #include "subroutines/general/calc_vxy_b_sia.F90"
765 #include "subroutines/general/calc_vxy_sia.F90"
766 #include "subroutines/general/calc_vz_sia.F90"
768 #include "subroutines/general/calc_temp_ssa.F90"
769 #include "subroutines/general/calc_vxy_ssa.F90"
770 #include "subroutines/general/calc_vxy_ssa_matrix.F90"
771 #include "subroutines/general/calc_vis_ssa.F90"
772 #include "subroutines/general/calc_vz_ssa.F90"
773 #include "subroutines/general/viscosity.F90"
775 #include "subroutines/general/calc_qbm.F90"
776 #include "subroutines/general/calc_water_bas.F90"
777 #include "subroutines/general/calc_top.F90"
778 #include "subroutines/general/calc_elra.F90"
779 #include "subroutines/general/sor_sprs.F90"
780 #include "subroutines/general/topograd_1.F90"
781 #include "subroutines/general/topograd_2.F90"
782 #include "subroutines/general/calc_temp_melt.F90"
783 #include "subroutines/general/calc_temp_bas.F90"
784 #include "subroutines/general/creep.F90"
785 #include "subroutines/general/ratefac.F90"
786 #include "subroutines/general/ratefac_t.F90"
787 #include "subroutines/general/kappa_val.F90"
788 #include "subroutines/general/c_val.F90"
789 #include "subroutines/general/metric.F90"
790 #include "subroutines/general/geo_coord.F90"
791 #include "subroutines/general/num_coord.F90"
792 #include "subroutines/general/tri_sle.F90"
793 #include "subroutines/general/output1.F90"
794 #include "subroutines/general/output2.F90"
795 #include "subroutines/general/output3.F90"
796 #include "subroutines/general/output4.F90"
799 #include "subroutines/asf/output5.F90"
803 #include "subroutines/general/output_nc.F90"
804 #include "subroutines/general/check.F90"
807 #include "subroutines/general/borehole.F90"
809 #if (ADV_HOR==3 && ADV_VERT==3)
810 #include "subroutines/general/calc_age_c.F90" /* Algorithm by B. Muegge */
811 #include "subroutines/general/calc_age_t.F90" /* Algorithm by B. Muegge */