Changeset 146 for branches/iLoveclim/SOURCES/bilan_eau_mod.f90

Timestamp:

10/20/17 09:31:39 (7 years ago)

Author:

aquiquet

Message:

Grisli-iLoveclim branch: merged to trunk at revision 145

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/bilan_eau_mod.f90 (modified) (4 diffs)

branches/iLoveclim/SOURCES/bilan_eau_mod.f90

@@ -31 +31 @@
 real,dimension(nx,ny) :: Bmelt_dtt    !< basal melting on ice points accumulated during dtt
 real,dimension(nx,ny) :: calv_dtt     !< calving sur dtt (pour calcul bilan d'eau)
+real,dimension(nx,ny) :: archimtab    ! point pose si > 0
+real,dimension(nx,ny) :: grline_dtt   ! grounding line flux during dtt
 
 real :: sum_H_old
 real :: diff_H
+real :: alpha_flot                      ! ro/row 
 
 real,dimension(nx,ny) :: bm_dt,bmelt_dt
@@ -42 +45 @@
 subroutine init_bilan_eau
  ! initialisation des variables
-        diff_H=0.
-        sum_H_old = sum(H(2:nx-1,2:ny-1),mask=ice(2:nx-1,2:ny-1)==1)
-        tot_water(:,:)=0.
-        bm_dt(:,:)=0.
-        bmelt_dt(:,:)=0.
+  diff_H=0.
+  sum_H_old = sum(H(2:nx-1,2:ny-1),mask=ice(2:nx-1,2:ny-1)==1)
+  tot_water(:,:)=0.
+  bm_dt(:,:)=0.
+  bmelt_dt(:,:)=0.
+  alpha_flot=ro/row
 ! iLOVECLIM initialisation of water conservation related variables
-        trendWAC=0.
-        smbWAC(:,:)=0.
-        bmeltWAC(:,:)=0.
-        calvingWAC(:,:)=0.
+  trendWAC=0.
+  smbWAC(:,:)=0.
+  bmeltWAC(:,:)=0.
+  calvingWAC(:,:)=0.
 end subroutine init_bilan_eau
 
-        
+
 
 
@@ -88 +92 @@
 
 where (ice(2:nx-1,2:ny-1).eq.1)
-        Bm_dtt(2:nx-1,2:ny-1) = Bm_dtt(2:nx-1,2:ny-1) + Bm_dt(2:nx-1,2:ny-1) !* dt ! somme Bm sur dt
-        bmelt_dtt(2:nx-1,2:ny-1) = bmelt_dtt(2:nx-1,2:ny-1) + bmelt_dt(2:nx-1,2:ny-1) ! * dt ! somme bmelt sur dt
+   Bm_dtt(2:nx-1,2:ny-1) = Bm_dtt(2:nx-1,2:ny-1) + Bm_dt(2:nx-1,2:ny-1) !* dt ! somme Bm sur dt
+   bmelt_dtt(2:nx-1,2:ny-1) = bmelt_dtt(2:nx-1,2:ny-1) + bmelt_dt(2:nx-1,2:ny-1) ! * dt ! somme bmelt sur dt
 endwhere
+
+archimtab(:,:) = Bsoc(:,:)+H(:,:)*alpha_flot - sealevel
+gr_line(:,:)=0
+do j=1,ny
+  do i=1,nx
+     !afq    if ((H(i,j).gt.0.).and.(archimtab(i,j).GE.0.).and.(Bsoc(i,j).LE.sealevel)) then ! grounded with ice
+     if ((H(i,j).gt.0.).and.(archimtab(i,j).GE.0.)) then ! grounded with ice
+      if (archimtab(i-1,j).LT.0..and.Uxbar(i,j).LT.0..and..not.flot_marais(i-1,j)) gr_line(i,j)=1
+      if (archimtab(i+1,j).LT.0..and.Uxbar(i+1,j).GT.0..and..not.flot_marais(i+1,j)) gr_line(i,j)=1
+      if (archimtab(i,j-1).LT.0..and.Uybar(i,j).LT.0..and..not.flot_marais(i,j-1)) gr_line(i,j)=1
+      if (archimtab(i,j+1).LT.0..and.Uybar(i,j+1).GT.0..and..not.flot_marais(i,j+1)) gr_line(i,j)=1
+    endif
+  enddo
+enddo
+
+where (gr_line(:,:).eq.1)
+!~   grline_dtt(:,:)= (((uxbar(:,:)+eoshift(uxbar(:,:),shift=1,boundary=0.,dim=1))**2+ &
+!~                     (uybar(:,:)+eoshift(uybar(:,:),shift=1,boundary=0.,dim=2))**2)**0.5)*0.5 &
+!~                     *H(:,:) + grline_dtt(:,:) 
+   grline_dtt(:,:)= - sqrt(                                                                        &
+                           ( (uxbar(:,:)+eoshift(uxbar(:,:),shift=1,boundary=0.,dim=1))*dy/2. )**2 +  &
+                           ( (uybar(:,:)+eoshift(uybar(:,:),shift=1,boundary=0.,dim=2))*dx/2. )**2 )  &
+                           * H(:,:) * dt + grline_dtt(:,:)
+endwhere
+
 
 if (isynchro.eq.1) then
@@ -103 +132 @@
    
 ! bilan d'eau sur la grille :
-        water_bilan=sum(tot_water(:,:))
-        diff_H = diff_H/dtt
+  water_bilan=sum(tot_water(:,:))
+  diff_H = diff_H/dtt
 
 !999 format(f0.2,1x,e15.8,1x,i10,8(1x,e15.8))
 !       write(6,999),time,sum_H,count(ice(:,:)==1),diff_H,water_bilan,sum(calv_dtt(:,:))/dtt,sum(ablbord_dtt(:,:))/dtt,sum(bmelt_dtt(:,:),mask=ice(:,:)==1)/dtt,sum(bm(:,:),mask=ice(:,:)==1),sum(Bm_dtt(:,:))/dtt,sum(bmelt_dtt(:,:))/dtt
-        diff_H_water_bilan(2:nx-1,2:ny-1)=tot_water(2:nx-1,2:ny-1)-diff_H_2D(2:nx-1,2:ny-1)
+  diff_H_water_bilan(2:nx-1,2:ny-1)=tot_water(2:nx-1,2:ny-1)-diff_H_2D(2:nx-1,2:ny-1)
 
 endif