Changeset 237 for trunk/SOURCES/Greeneem_files/lect-clim-act-greeneem_mois_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/Greeneem_files/lect-clim-act-greeneem_mois_mod.f90 (modified) (3 diffs)

trunk/SOURCES/Greeneem_files/lect-clim-act-greeneem_mois_mod.f90

@@ -74 +74 @@
       do i=1,nx
          read (20,*) S0CLIM(i,j)
-         S0CLIM(i,j)=max(S0CLIM(i,j),slv(i,j))   ! pour etre au niveau des mers
+         S0CLIM(i,j)=max(S0CLIM(i,j),sealevel_2d(i,j))   ! pour etre au niveau des mers
       end do
    end do
@@ -120 +120 @@
       do i=1,nx
 ! Zs est l'altitude de la surface qu'elle soit mer, glace ou lac
-         ZS(i,j)=max(slv(i,j),S0CLIM(I,J)) !slv=sealev pour nolake
+         ZS(i,j)=max(sealevel_2d(i,j),S0CLIM(I,J)) !sealevel_2d=sealev pour nolake
 
          do mo=1,mois
@@ -261 +261 @@
 do j=1,ny
    do i=1,nx
-      ZS(i,j)=max(slv(i,j),S(i,j))
+      ZS(i,j)=max(sealevel_2d(i,j),S(i,j))
       Tann(i,j)=d_ann+gamma_ann*S(i,j)+c_ann*ylat(i,j)+kappa_ann*(360-xlong(i,j))
       Tjuly(i,j)=d_july+gamma_july*S(i,j)+c_july*ylat(i,j)+kappa_july*(360-xlong(i,j))