Changeset 476

Timestamp:

02/07/24 16:31:29 (5 months ago)

Author:

aquiquet

Message:

Cleaning branch: ablation_bord with explicit arguments

Location:

branches/GRISLIv3/SOURCES

Files:

• ablation_bord.f90 (modified) (2 diffs)
• main3D-0.4-40km.f90 (modified) (2 diffs)
• steps_time_loop.f90 (modified) (2 diffs)

branches/GRISLIv3/SOURCES/ablation_bord.f90

@@ -1 +1 @@
 !> \file ablation_bord.f90
-!! Clacul de l'ablation sur le bord
+!! Evaluation of the ablation at the edge (bord) of the ice sheet 
 !<
       
 
 !> SUBROUTINE: Ablation_bord
-!! Calcule l'ablation pour les points du bord 
+!! Evaluation of the ablation at the edge (bord) of the ice sheet
 !! \author ...
 !! \date ....
-!! @note cas des bords de calotte : un point du bord peut perdre en plus la glace 
-!! amenée par l'écoulement.
-!! @note Calcule l'ablation pour les points du bord avec H=0 mais apport par l'amont
-!! diagnostic pour la fermeture du bilan d'eau
-!! par convention ablation negative dans la zone d'ablation 
-!! @note Used modules:
-!! @note  - use module3D_phy
+!! @note ice edge case: it can loose ice brought by ice transport
+!! @note Evaluate the ablation at the edge with H=0 but ice brought by mass transport
+!! diagnosed for water conservation
+!! ablation is defined as negative when ice is removed
 !>
 
-subroutine ablation_bord
+subroutine ablation_bord(ice,bm,bmelt,ablbord,ablbord_dtt,dt)
 
+  implicit none
 
-  use module3D_phy, only: ice,ablbord,bmelt,bm,dt,ablbord_dtt
-  use geography, only: nx,ny
-  implicit none
-  integer :: ip1,im1,jp1,jm1   ! pour chercher les voisins
+  real                  , intent(in)    :: dt
+  integer,dimension(:,:), intent(in)    :: ice
+  real,   dimension(:,:), intent(in)    :: bm,bmelt
+  real,   dimension(:,:), intent(inout) :: ablbord,ablbord_dtt
+  
+  integer :: nxloc,nyloc
+  integer :: ip1,im1,jp1,jm1   ! to look for neighbours
   integer :: somm
   integer :: i,j
 
-  ! Dans les zones ou le calcul icethick a donne un H negatif : ablbord=H
+  nxloc=ubound(ice,1)
+  nyloc=ubound(ice,2)
 
-  ! pour être au bord, il suffit que la somme des ice voisins soit > 0
-  do j=1,ny
-     do i=1,nx
+  ! When icethick has given negative values for H: ablbord=H
+
+  ! to be bord (edge), the sum of neighbouring ice points should be > 0
+  do j=1,nyloc
+     do i=1,nxloc
         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)
+           ip1=min(i+1,nxloc)
+           jp1=min(j+1,nyloc)
 
            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
+           if (somm.gt.0.and.(ablbord(i,j).LT.0.)) then        ! glaciated neighbours
+
               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
@@ -57 +58 @@
   ablbord_dtt(:,:) = ablbord_dtt(:,:) + ablbord(:,:)
 
-  !cdc deplace dans bilan_flux_output_mod ablbord(:,:)=0.
-
 end subroutine ablation_bord

branches/GRISLIv3/SOURCES/main3D-0.4-40km.f90

@@ -132 +132 @@
 
   USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen, &
+                          ice,bm,bmelt,ablbord,ablbord_dtt,dt,    &
                           s,h,b,bsoc,flot,mk,mk0,uxbar,uybar,hwater,time,timemax,ndebug,ndebug_max
   use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,runname,itracebug
@@ -238 +239 @@
            call flottab()
            call calving
-           call ablation_bord
+           call ablation_bord(ice,bm,bmelt,ablbord,ablbord_dtt,dt)
            call flottab
            call Neffect()

branches/GRISLIv3/SOURCES/steps_time_loop.f90

@@ -15 +15 @@
 
   use module3d_phy, only: ispinup,isynchro,timemax,time,flot,S,B,Bsoc,H,sealevel_2d,hmx,hmy,&
+       ice,bm,bmelt,ablbord,ablbord_dtt,dt, &
        uxbar,uybar,hdot,dtmax,dtmin,iout
   use runparam, only: itracebug,nt,num_tracebug
@@ -60 +61 @@
      call determin_tache
      call calving
-     call ablation_bord
+     call ablation_bord(ice,bm,bmelt,ablbord,ablbord_dtt,dt)
      call bilan_eau
      call bilan_flux_output