Changeset 237 for trunk/SOURCES/MISMIP3D_files/lect-mismip3d_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/MISMIP3D_files/lect-mismip3d_mod.f90 (modified) (1 diff)

trunk/SOURCES/MISMIP3D_files/lect-mismip3d_mod.f90

@@ -111 +111 @@
   do j=1,ny
      do i=1,nx
-        if ((bsoc(i,j)+h(i,j)*ro/row -sealevel).lt.0.) then
+        if ((bsoc(i,j)+h(i,j)*ro/row -sealevel_2d(i,j)).lt.0.) then
            flot(i,j)=.true.
         else