Changeset 237 for trunk/SOURCES/Hudson_files/eaubasale-0.5_hudson_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/Hudson_files/eaubasale-0.5_hudson_mod.f90 (modified) (4 diffs)

trunk/SOURCES/Hudson_files/eaubasale-0.5_hudson_mod.f90

@@ -139 +139 @@
         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
@@ -182 +182 @@
            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
@@ -308 +308 @@
 
               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.
@@ -320 +320 @@
 
               ! Attention le bloc suivant est pour le run 20 
-              !           Zflot=row/ro*(sealevel-Bsoc(i,j))-10.
+              !           Zflot=row/ro*(sealevel_2d(i,j)-Bsoc(i,j))-10.
               Zflot=H(i,j)*rog/rofreshg