Changeset 237 for trunk/SOURCES/isostasie_mod-0.3.f90

Timestamp:

01/09/19 17:09:26 (5 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/isostasie_mod-0.3.f90 (modified) (2 diffs)

trunk/SOURCES/isostasie_mod-0.3.f90

@@ -103 +103 @@
     do i=1,nx
        do j=1,ny
-          if (ro*h0(i,j).ge.-row*(BSOC0(i,j)-sealevel)) then
+          if (ro*h0(i,j).ge.-row*(BSOC0(i,j)-sealevel_2d(i,j))) then
 
              !        glace ou terre
@@ -110 +110 @@
           else
              !        ocean
-             charge(i,j)=-rowg*(Bsoc0(i,j)-sealevel)
+             charge(i,j)=-rowg*(Bsoc0(i,j)-sealevel_2d(i,j))
           endif
        end do