Changeset 237 for trunk/SOURCES/Greeneem_files/output_greeneem_mod-0.4.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/Greeneem_files/output_greeneem_mod-0.4.f90 (modified) (1 diff)

trunk/SOURCES/Greeneem_files/output_greeneem_mod-0.4.f90

@@ -141 +141 @@
      do j=1,ny
         do i=1,nx
-           if (mask_cal(i,j,k).and..not.flot(i,j).and.H(i,j).GT.1..and.((sealevel-B(i,j)).LE.0.)) then
+           if (mask_cal(i,j,k).and..not.flot(i,j).and.H(i,j).GT.1..and.((sealevel_2d(i,j)-B(i,j)).LE.0.)) then
               isvolf(kk) = isvolf(kk) + H(i,j)
            else
-              isvolf(kk) = isvolf(kk) + H(i,j)-row/ro*(sealevel-B(i,j))
+              isvolf(kk) = isvolf(kk) + H(i,j)-row/ro*(sealevel_2d(i,j)-B(i,j))
            endif
         enddo