43 integer(i4b),
intent(in) :: i, j
44 real(dp),
intent(in) :: z_sl
46 real(dp) :: rhosw_rho_ratio, h_ice, h_sea
48 rhosw_rho_ratio = rho_sw/rho
52 if ( (maske(j,i) == 1_i2b).or.(maske(j,i) == 2_i2b) )
then
54 if (zl(j,i) > z_sl)
then
66 if (zl(j,i) > z_sl)
then
73 h_ice = zs(j,i)-zb(j,i)
76 if ( h_ice < (rhosw_rho_ratio*h_sea) )
then
78 #if ( MARGIN==1 || ( MARGIN==2 && MARINE_ICE_FORMATION==1 ) )
81 #elif ( MARGIN==2 && MARINE_ICE_FORMATION==2 )
integer(i2b) function mask_update(z_sl, i, j)
Update of the topography mask due to changes of the sea level.
Declarations of kind types for SICOPOLIS.
Declarations of global variables for SICOPOLIS (for the ANT domain).
Declarations of global variables for SICOPOLIS.