Changeset 237 for trunk/SOURCES/flottab2-0.7.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/flottab2-0.7.f90 (modified) (12 diffs)

trunk/SOURCES/flottab2-0.7.f90

@@ -122 +122 @@
        where ( mk_init(:,:).eq.1)                 ! pose
           flot(:,:) = .False.
-          H(:,:)=max(H(:,:),(10.+sealevel-Bsoc(:,:))*row/ro)  ! pour avoir archim=10
+          H(:,:)=max(H(:,:),(10.+sealevel_2d(:,:)-Bsoc(:,:))*row/ro)  ! pour avoir archim=10
           S(:,:) = Bsoc(:,:) + H(:,:) 
        elsewhere
@@ -175 +175 @@
           uys1(i,j)=uybar(i,j)
 
-          archim = Bsoc(i,j)+H(i,j)*ro/row -sealevel
+          archim = Bsoc(i,j)+H(i,j)*ro/row -sealevel_2d(i,j)
           !      if ((i.eq.132).and.(j.eq.183)) print*,'archim=',archim
 
@@ -189 +189 @@
                 if (igrdline.eq.1) then   ! en cas de grounding line prescrite
                    flot(i,j)=.false.
-                   H(i,j)=(10.+sealevel-Bsoc(i,j))*row/ro  ! pour avoir archim=10
+                   H(i,j)=(10.+sealevel_2d(i,j)-Bsoc(i,j))*row/ro  ! pour avoir archim=10
                    new_flot_point(i,j)=.false.
                 endif
@@ -202 +202 @@
                    if (igrdline.eq.1) then   ! en cas de grounding line prescrite
                       flot(i,j)=.false.
-                      H(i,j)=(10.+sealevel-Bsoc(i,j))*row/ro  ! pour avoir archim=10
+                      H(i,j)=(10.+sealevel_2d(i,j)-Bsoc(i,j))*row/ro  ! pour avoir archim=10
                       new_flot_point(i,j)=.false.
                    endif
@@ -228 +228 @@
              flot(i,j)=.true.  !cdc points ocean sont flot meme sans glace
              H(i,j)=0.
-             S(i,j)=sealevel
+             S(i,j)=sealevel_2d(i,j) !afq -- WARNING: est-ce qu'on veut vraiment mettre S a la valeur locale du niveau marin?
              B(i,j)=S(i,j)-H(i,j)
           endif arch
@@ -239 +239 @@
 
              surnet=H(i,j)*(1.-ro/row)
-             S(i,j)=surnet+sealevel
+             S(i,j)=surnet+sealevel_2d(i,j) !afq -- WARNING: est-ce qu'on veut vraiment mettre S a la valeur locale du niveau marin?
              B(i,j)=S(i,j)-H(i,j)
           end if
@@ -246 +246 @@
     end do
     !$OMP END DO
-    !~ print*,'flottab 2',S(132,183),H(132,183),BSOC(132,183),B(132,183),sealevel
-    !~ print*,'flottab 2',flot(132,183),ice(132,183)
-    !~ if (S(132,183).LT.sealevel) print*,'BUGGGGGGGGGGGGG !!!!!!!!!!!!!!!!'
        
 !!$ do i=1,nx
@@ -262 +259 @@
 
 
-
     !-----------------------------------------------------------------------
     !$OMP DO
@@ -298 +294 @@
           !                 sur le demi noeud condition archim < 100 m
 
-          archim=(Bsoc(i,j)+Bsoc(i-1,j))*0.5-sealevel+ro/row*Hmx(i,j)
+          archim=(Bsoc(i,j)+Bsoc(i-1,j))*0.5-(sealevel_2d(i,j)+sealevel_2d(i-1,j))*0.5+ro/row*Hmx(i,j)
           gzmx(i,j)=gzmx(i,j).and.(archim.le.100.) 
           cotemx(i,j)=gzmx(i,j)
@@ -337 +333 @@
           !                 sur le demi noeud condition archim > 100 m
 
-          archim=(Bsoc(i,j)+Bsoc(i,j-1))*0.5-sealevel+ro/row*Hmy(i,j)
+          archim=(Bsoc(i,j)+Bsoc(i,j-1))*0.5-(sealevel_2d(i,j)+sealevel_2d(i,j-1))*0.5+ro/row*Hmy(i,j)
           gzmy(i,j)=gzmy(i,j).and.(archim.le.100.) 
           cotemy(i,j)=gzmy(i,j)
@@ -605 +601 @@
           ice(:,:)=0
           H(:,:)=0.
-          S(:,:)=H(:,:)*(1.-ro/row) + sealevel
+          S(:,:)=H(:,:)*(1.-ro/row) + sealevel_2d(:,:)
           B(:,:)=S(:,:) - H(:,:)
        end where
@@ -790 +786 @@
             flgzmy(smax_i,smax_j)=.false. ; flgzmx(smax_i,smax_j+1)=.false.
      
-            if (Smax_.le.sealevel) then
+            if (Smax_.le.sealevel) then  ! afq -- WARNING: en fait avec le niveau marin local cest plus complique que ca!
+            !if (Smax_.le.sealevel_2d(i,j)) then ! afq --WARNING: il faudrait regarder point par point sur la tache...
               write(num_tracebug,*)'Attention, une ile avec la surface sous l eau'
               write(num_tracebug,*)'time=',time,'   coord:',smax_i,smax_j