Changeset 237 for trunk/SOURCES/3D-physique-gen_mod.f90

Timestamp:

01/09/19 17:09:26 (6 years ago)

Author:

aquiquet

Message:

Sealevel is now treated as a 2D variable (sealevel_2d while sealevel remains the eustatic sea level), results should remain identical as sealevel_2d is equal to sealevel in this revision.

File:

• trunk/SOURCES/3D-physique-gen_mod.f90 (modified) (6 diffs)

trunk/SOURCES/3D-physique-gen_mod.f90

@@ -141 +141 @@
   real ::  SECYEAR                      !< for relation an/seconds
   !real ::  S22                          ! ct for PDD calculation
-  real ::  sealevel                     !< niveau des mers 
+  real ::  sealevel                     !< niveau des mers, afq: now used as eustatic sea level
   ! real ::  SIGMA                        ! variabilite Tday
   real ::  SURF                         !< 
@@ -235 +235 @@
   integer,dimension(nx,ny) :: Mk_init      !< initial mask (with islands, outcrops, ...
   integer,dimension(nx,ny) :: MK           !< masks (ice sheet, max, above water, below water, 1)
-  integer,dimension(nx,ny) :: M_sealev     !< masks for sea level calculation (present ice sheet)
   integer,dimension(nx,ny) :: MNEG         !< masks (ice sheet, max, above water, below water, 1)
   integer,dimension(nx,ny) :: IBASE     !< type de base (froide, temperee)
@@ -245 +244 @@
   real,dimension(nx,ny) :: ACQUA        !< Surface des surfaces en eau
   real,dimension(nx,ny) :: B            !< Altitude de la base de la glace  'o'
-  real,dimension(nx,ny) :: B_sealev     !< Altitude de la base de la glace pour le
-  !< calcul du nivx de la mer (/present)
   real,dimension(nx,ny) :: BDOT         !< derivee de B / t
   real,dimension(nx,ny) :: B1           !< 
@@ -310 +307 @@
 
 
-  real,dimension(nx,ny) :: H_sealev     !< ice thickness for sea level calculation        ! (present ice sheet)
   real,dimension(nx,ny) :: HDOT         !< ice thickness derivee / t
   real,dimension(nx,ny) :: HDOTWATER
@@ -334 +330 @@
   real,dimension(nx,ny) :: SW           !< for bedrock isostasy
   real,dimension(nx,ny) :: S            !< altitude of ice sheet surface
-  real,dimension(nx,ny) :: S_sealev     !< altitude of ice sheet surface for sea
-  !< level calculation (present ice sheet)
+  real,dimension(nx,ny) :: sealevel_2d  !< local sea surface elevation
   real,dimension(nx,ny) :: SLOPE        !< slope 'o'
   real,dimension(nx,ny) :: SDX          !< slope derivee / x '>'
@@ -344 +339 @@
   real,dimension(nx,ny) :: SLOPE2my     !< = Sdy**2 + Sdxmy**2 '^'
   real,dimension(nx,ny) :: S0           !< altitude actuelle de la surface
-  real,dimension(nx,ny) :: slv          !< niveau de flottaison (sealevel et lakes)
+!afq -- replaced by sealevel_2D:  real,dimension(nx,ny) :: slv          !< niveau de flottaison (sealevel et lakes)
   real,dimension(nx,ny) :: TJULY        !< Ground air temperature July
   real,dimension(nx,ny) :: TANN         !< Ground air temperature annual