Changeset 237 for trunk/SOURCES/eaubasale-0.5_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/eaubasale-0.5_mod.f90 (modified) (3 diffs)

trunk/SOURCES/eaubasale-0.5_mod.f90

@@ -158 +158 @@
         if (flot(i,j))then  ! points flottants
            klimit(i,j)=1
-           limit_hw(i,j)=(sealevel-Bsoc(i,j))*rowg/rofreshg
+           limit_hw(i,j)=(sealevel_2d(i,j)-Bsoc(i,j))*rowg/rofreshg
 
         else if (IBASE(I,J).eq.1) then ! base froide
@@ -200 +200 @@
            pot_w(I,J)=pot_w(I,J)+rog*H(I,J)
 
-           pot_f(I,J)=rowg*(sealevel-S(i,j)+H(I,J))  ! pression a la base de l'ice shelf
+           pot_f(I,J)=rowg*(sealevel_2d(i,j)-S(i,j)+H(I,J))  ! pression a la base de l'ice shelf
         enddo
      enddo
@@ -314 +314 @@
 
               if (flot(i,j)) then  ! if flot > hwater=hwater in ocean 
-                 hwater(i,j)= sealevel - bsoc(i,j)
+                 hwater(i,j)= sealevel_2d(i,j) - bsoc(i,j)
                  !             hwater(i,j)= max(0.,hwater(i,j))
                  if (hwater(i,j).lt.0.) hwater(i,j)=0.