Changeset 102 for trunk/SOURCES/ablation_bord.f90

Timestamp:

05/05/17 16:44:01 (7 years ago)

Author:

dumas

Message:

Update closed water cycle | H = 1m on sea suppressed

File:

• trunk/SOURCES/ablation_bord.f90 (modified) (3 diffs)

trunk/SOURCES/ablation_bord.f90

@@ -32 +32 @@
 ! mise a jour du masque "ice"
 where (flot(:,:))        ! points flottants, sera éventuellement réévalué dans flottab
-   where(H(:,:).gt.max(Hmin,Hmin+BM(:,:)-Bmelt(:,:)))
+!cdc 1m  where(H(:,:).gt.max(Hmin,Hmin+BM(:,:)-Bmelt(:,:)))
+   where(H(:,:).gt.0.)
       ice(:,:)=1
    elsewhere
       ice(:,:)=0
-      H(:,:)=max(1.,min(0.,(sealevel - Bsoc(:,:))*row/ro-0.01))
+!cdc 1m      H(:,:)=max(1.,min(0.,(sealevel - Bsoc(:,:))*row/ro-0.01))
+                        H(:,:)=0.
    endwhere
 elsewhere                ! points posés
@@ -47 +49 @@
 
 
+
+
+
 where (ice(:,:).eq.1)    !  quand glace : abl=bm-acc 
    ablbord(:,:)=0.
 endwhere
+!~ where (calv(:,:).LT.0.) ! il ne peut pas y avoir calving et ablbord sur un meme point
+!~       ablbord(:,:)=0.
+!~ endwhere
 
 ! pour être au bord, il suffit que la somme des ice voisins soit > 0
@@ -62 +70 @@
          somm=ice(ip1,j)+ice(im1,j)+ice(i,jp1)+ice(i,jm1)
          if (somm.gt.0.and.(ablbord(i,j).LT.0.)) then                               ! voisins englaces
-!         if (somm.gt.0.and.((ablbord(i,j).LT.0.).OR.(flot(ip1,j).OR.flot(im1,j).OR.flot(i,jp1).OR.flot(i,jm1)))) then   ! voisins englaces
-!            ablbord(i,j)= ablbord(i,j) - (Acc(i,j) * dt)     !  Ablvrai=bmpdd-Acc-Hcalc/dt
-            ablbord(i,j)= ablbord(i,j) + ((Bmelt(i,j) - Bm(i,j)) * dt)  ! version Tof : ablbord = Hcons_mass - Bm - Bmelt
+!~ !         if (somm.gt.0.and.((ablbord(i,j).LT.0.).OR.(flot(ip1,j).OR.flot(im1,j).OR.flot(i,jp1).OR.flot(i,jm1)))) then   ! voisins englaces
+!~ !            ablbord(i,j)= ablbord(i,j) - (Acc(i,j) * dt)     !  Ablvrai=bmpdd-Acc-Hcalc/dt
+!~           if (ablbord(i,j).LT.0..and.front(i,j).LE.4.and.front(i,j).GE.0) then
+            ablbord(i,j)= max(0.,ablbord(i,j) + ((Bmelt(i,j) - Bm(i,j)) * dt))  ! version Tof : ablbord = Hcons_mass - Bm - Bmelt
 
-         else
-! version cat            ablbord(i,j)=hdot(i,j)-acc(i,j)   ! on ne fond que ce qui est tombe + existant
-                                                ablbord(i,j)=0.
-         endif
-
+          else
+!~ ! version cat            ablbord(i,j)=hdot(i,j)-acc(i,j)   ! on ne fond que ce qui est tombe + existant
+                                                ablbord(i,j)=0.
+          endif
+        else
+                      ablbord(i,j)=0.
       end if
    end do