56 integer(i4b),
intent(in) :: i, j
57 real(dp),
intent(in) :: z_sl
59 integer(i2b) :: mask_update_sea_level
60 real(dp) :: rhosw_rho_ratio, H_ice, H_sea
66 if ( (
maske(j,i) == 1_i2b).or.(
maske(j,i) == 2_i2b) )
then 68 if (
zl(j,i) > z_sl)
then 69 mask_update_sea_level = 1_i2b
72 mask_update_sea_level = 2_i2b
80 if (
zl(j,i) > z_sl)
then 82 mask_update_sea_level = 0_i2b
87 h_ice =
zs(j,i)-
zb(j,i)
90 if ( h_ice < (rhosw_rho_ratio*h_sea) )
then 92 #if (MARGIN==1 || (MARGIN==2 && MARINE_ICE_FORMATION==1)) 93 mask_update_sea_level = 2_i2b
95 #elif (MARGIN==2 && MARINE_ICE_FORMATION==2) 96 mask_update_sea_level = 0_i2b
98 mask_update_sea_level = 3_i2b
104 mask_update_sea_level = 0_i2b
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 ...
integer(i2b) function, public mask_update_sea_level(z_sl, i, j)
Main function of mask_update_m: Update of the ice-land-ocean mask due to changes of the sea level...
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).
Declarations of kind types for SICOPOLIS.
real(dp), dimension(0:jmax, 0:imax) zs
zs(j,i): Coordinate z of the surface topography
real(dp) rho_sw
RHO_SW: Density of sea water.
real(dp) rho
RHO: Density of ice.
Update of the ice-land-ocean mask due to changes of the sea level.
real(dp), dimension(0:jmax, 0:imax) zb
zb(j,i): Coordinate z of the ice base
Declarations of global variables for SICOPOLIS.