56 real(dp),
intent(in) :: z_sl
58 real(dp) :: year_sec_inv
59 real(dp) :: rhosw_rho_ratio
60 real(dp) :: calv_uw_coeff, r1_calv_uw, r2_calv_uw
61 real(dp),
dimension(0:JMAX,0:IMAX) :: H, H_sea
65 year_sec_inv = 1.0_dp/year_sec
71 #if (defined(CALV_UW_COEFF)) 72 calv_uw_coeff = calv_uw_coeff * year_sec_inv
74 stop
' >>> calving_underwater_ice: CALV_UW_COEFF undefined!' 77 #if (defined(R1_CALV_UW)) 78 r1_calv_uw = r1_calv_uw
80 stop
' >>> calving_underwater_ice: R1_CALV_UW undefined!' 83 #if (defined(R2_CALV_UW)) 84 r2_calv_uw = r2_calv_uw
86 stop
' >>> calving_underwater_ice: R2_CALV_UW undefined!' 92 h_sea = max(z_sl -
zl, 0.0_dp)
96 where ( (
maske == 0_i2b).and.(h < rhosw_rho_ratio*h_sea) )
97 calv_uw_ice = calv_uw_coeff * h**r1_calv_uw * h_sea**r2_calv_uw
integer(i2b), dimension(0:jmax, 0:imax) maske
maske(j,i): Ice-land-ocean mask. 0: grounded ice, 1: ice-free land, 2: ocean, 3: floating ice ...
real(dp), dimension(0:jmax, 0:imax) zl
zl(j,i): Coordinate z of the lithosphere surface
Declarations of global variables for SICOPOLIS (for the ANT domain).
real(dp), dimension(0:jmax, 0:imax) calv_uw_ice
Calving of "underwater ice".
Declarations of kind types for SICOPOLIS.
real(dp), dimension(0:jmax, 0:imax) h_c
H_c(j,i): Thickness of ice in the upper (kc) domain (thickness of the cold-ice layer for POLY...
subroutine calving_underwater_ice(z_sl)
Main routine: Calving of "underwater ice".
real(dp) rho_sw
RHO_SW: Density of sea water.
real(dp) rho
RHO: Density of ice.
real(dp), dimension(0:jmax, 0:imax) h_t
H_t(j,i): Thickness of ice in the lower (kt) domain (thickness of the temperate layer for POLY...
Declarations of global variables for SICOPOLIS.