Changeset 237 for trunk/SOURCES/Hemin40_files/bmelt-hemin40-depth_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/Hemin40_files/bmelt-hemin40-depth_mod.f90 (modified) (2 diffs)

trunk/SOURCES/Hemin40_files/bmelt-hemin40-depth_mod.f90

@@ -85 +85 @@
          do i=1,nx
 
-            if ((sealevel-bsoc(i,j)).lt.600) then
+            if ((sealevel_2d(i,j)-bsoc(i,j)).lt.600) then
                     bmshelf(i,j)=0.2 ! 0.45!   0.65  
                     bmgrz(i,j)=0.2 ! 0.45 !  0.65
@@ -129 +129 @@
 shelf:    if (flot(i,j)) then    ! partie flottante
 
-            if ((sealevel-bsoc(i,j)).lt.600) then
+            if ((sealevel_2d(i,j)-bsoc(i,j)).lt.600) then
                     if ((j.lt.75).and.(i.lt.145) ) then  !Atlantique 
                     bmelt(i,j)=(.6+.4*coefbmshelf)*bmshelf(i,j)