Changeset 119

Timestamp:

06/22/17 17:03:39 (7 years ago)

Author:

aquiquet

Message:

S<sealevel bug correction

Location:

trunk/SOURCES

Files:

• flottab2-0.7.f90 (modified) (3 diffs)
• steps_time_loop.f90 (modified) (1 diff)

trunk/SOURCES/flottab2-0.7.f90

@@ -215 +215 @@
                 FLOT(I,J)=.false.
                 BOOST=.false.
-
              endif
              !cdc correction topo pour suivre  variations sealevel
@@ -227 +226 @@
 !cdc  1m           H(i,j)=min(1.,max(0.,(sealevel - Bsoc(i,j))*row/ro-0.01))
              H(i,j)=0.
-             surnet=H(i,j)*(1.-ro/row)
-             S(i,j)=surnet+sealevel
+             S(i,j)=sealevel
              B(i,j)=S(i,j)-H(i,j)
           endif arch
@@ -249 +247 @@
     !~ print*,'flottab 2',flot(132,183),ice(132,183)
     !~ if (S(132,183).LT.sealevel) print*,'BUGGGGGGGGGGGGG !!!!!!!!!!!!!!!!'
-
+       
 !!$ do i=1,nx
 !!$    do j=1,ny

trunk/SOURCES/steps_time_loop.f90

@@ -81 +81 @@
      where(.not.flot(:,:))
         S(:,:)=Bsoc(:,:)+H(:,:)
+        S(:,:)=max(S(:,:),sealevel)
         B(:,:)=Bsoc(:,:)
      end where