Changeset 377

Timestamp:

03/20/23 18:42:16 (13 months ago)

Author:

dumas

Message:

use only in subroutine ablation_bord

File:

• branches/GRISLIv3/SOURCES/ablation_bord.f90 (modified) (1 diff)

branches/GRISLIv3/SOURCES/ablation_bord.f90

@@ -20 +20 @@
 
 
-use module3D_phy
-implicit none
-integer :: ip1,im1,jp1,jm1   ! pour chercher les voisins
-integer :: somm
+  use module3D_phy, only:nx,ny,ice,ablbord,bmelt,bm,dt,ablbord_dtt
+  implicit none
+  integer :: ip1,im1,jp1,jm1   ! pour chercher les voisins
+  integer :: somm
+  integer :: i,j
 
-! Dans les zones ou le calcul icethick a donne un H negatif : ablbord=H
+  ! Dans les zones ou le calcul icethick a donne un H negatif : ablbord=H
 
-! pour être au bord, il suffit que la somme des ice voisins soit > 0
-do j=1,ny
-   do i=1,nx
-      if (ice(i,j).eq.0) then   ! seulment pour les noeuds sans glace
-         im1=max(i-1,1)
-         jm1=max(j-1,1)
-         ip1=min(i+1,nx)
-         jp1=min(j+1,ny)
+  ! pour être au bord, il suffit que la somme des ice voisins soit > 0
+  do j=1,ny
+     do i=1,nx
+        if (ice(i,j).eq.0) then   ! seulment pour les noeuds sans glace
+           im1=max(i-1,1)
+           jm1=max(j-1,1)
+           ip1=min(i+1,nx)
+           jp1=min(j+1,ny)
 
-         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
-!~           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
+           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
+              !~             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
-end do
+           ablbord(i,j)=0.
+        end if
+     end do
+  end do
 
-ablbord_dtt(:,:) = ablbord_dtt(:,:) + ablbord(:,:)
+  ablbord_dtt(:,:) = ablbord_dtt(:,:) + ablbord(:,:)
 
-!cdc deplace dans bilan_flux_output_mod ablbord(:,:)=0.
+  !cdc deplace dans bilan_flux_output_mod ablbord(:,:)=0.
 
 end subroutine ablation_bord