7 !! Update of the topography mask due to changes of the sea level.
11 !! Copyright 2009-2013 Ralf Greve
15 !! This file is part of SICOPOLIS.
17 !! SICOPOLIS is free software: you can redistribute it and/or modify
18 !! it under the terms of the GNU General Public License as published by
19 !! the Free Software Foundation, either version 3 of the License, or
20 !! (at your option) any later version.
22 !! SICOPOLIS is distributed in the hope that it will be useful,
23 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
24 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
25 !! GNU General Public License for more details.
27 !! You should have received a copy of the GNU General Public License
28 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
30 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
33 !> Update of the topography mask due to changes of the sea level.
34 !<------------------------------------------------------------------------------
42 integer(i4b),
intent(in) :: i, j
43 real(dp),
intent(in) :: z_sl
45 real(dp) :: rhosw_rho_ratio, h_ice, h_sea
47 rhosw_rho_ratio = rho_sw/rho
51 if ( (maske(j,i) == 1_i2b).or.(maske(j,i) == 2_i2b) )
then
53 if (zl(j,i) > z_sl)
then
65 if (zl(j,i) > z_sl)
then
72 h_ice = zs(j,i)-zb(j,i)
75 if ( h_ice < (rhosw_rho_ratio*h_sea) )
then
77 #if ( MARGIN==1 || ( MARGIN==2 && MARINE_ICE_FORMATION==1 ) )
80 #elif ( MARGIN==2 && MARINE_ICE_FORMATION==2 )