Changeset 173 for trunk/SOURCES

Timestamp:

01/08/18 17:15:03 (7 years ago)

Author:

aquiquet

Message:

Basal melt rates depend on distance from limit of the continental shelf (tof&afq)

File:

• trunk/SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90 (modified) (9 diffs)

trunk/SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90

@@ -23 +23 @@
   real,dimension(nx,ny) ::  bmshelf     !< fusion basale sous shelf
 
-  REAL,dimension(nx,ny) :: dist_talu    !< distance du point au talu continental
   REAL,dimension(nx,ny) :: typeshelf    !< Type de shelf Ronne->1 Ross ->2 ....
 
@@ -38 +37 @@
   real,dimension(nx,ny) ::  bmgrz_0               !< fusion basale a la grounding zone
   real,dimension(nx,ny) ::  bmshelf_0             !< fusion basale sous shelf
+  integer :: flag_dist
+  real,dimension(nx,ny) :: coef_bmelt_dist
   
 CONTAINS
@@ -48 +49 @@
 
     real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer
+    real :: bmelt_dist0
     character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves
     character(len=100) :: file_bmelt_anom
+    real,dimension(nx,ny) :: dist_talus  ! distance au talus
+    character(len=100) :: file_dist_talus ! fichier distance au talus
     
     namelist/bmelt_ant_reg_initmip/bmelt_regions,bmgrz_regions,  &  
-         bmelt_talus,bmgrz_talus,bmelt_coef,file_number_shelves  
+         bmelt_talus,bmgrz_talus,bmelt_coef,file_number_shelves, &
+         flag_dist,file_dist_talus,bmelt_dist0  
     namelist/bmelt_anom_initMIP/file_bmelt_anom, bmelt_time  
 
@@ -71 +76 @@
     write(num_rep_42,*) 'bmelt_coef      = ', bmelt_coef
     write(num_rep_42,*) 'file_numer_shelves = ', file_number_shelves
+    write(num_rep_42,*) 'flag_dist = ', flag_dist
+    write(num_rep_42,*) 'file_dist_talus = ', file_dist_talus
+    write(num_rep_42,*) 'bmelt_dist0 = ', bmelt_dist0
     write(num_rep_42,*)'/'                            
     write(num_rep_42,428) '! bmelt_regions    :  fonte basale sous shelves pour les 18 regions initMIP'
@@ -77 +85 @@
     write(num_rep_42,428) '! bmelt_coef                 :  coef fonte (1 pour conserver val)'
     write(num_rep_42,428) '! file_numer_ice-shelves     : fichier zones ice shelves'
+    write(num_rep_42,428) '! flag_dist : flag pour bmelt fnct distance talus'
+    write(num_rep_42,428) '! file_dist_talus : fichier de distance talus'
+    write(num_rep_42,428) '! bmelt_dist0 : coef de bmelt au talus'
     write(num_rep_42,*)
 
@@ -87 +98 @@
     typeshelf(:,:)  = tab(:,:)
 
+
+    if (flag_dist.eq.1) then
+       file_dist_talus=TRIM(DIRNAMEINP)//trim(file_dist_talus)
+       call Read_Ncdf_var('z',file_dist_talus,tab)
+       dist_talus(:,:)  = min(max(tab(:,:),0.),400.)
+       coef_bmelt_dist(:,:)= 1. / ((1/bmelt_dist0)*(1+(bmelt_dist0-1)*(dist_talus(:,:)/400.)**2))
+    else
+       coef_bmelt_dist(:,:)= 1.
+    endif
 
     bms_init:    do j=1,ny
@@ -172 +192 @@
           shelf:    if (flot(i,j)) then    ! partie flottante
 
-             bmelt(i,j)=coef_talus*bmshelf(i,j)
+             bmelt(i,j)=coef_talus*bmshelf(i,j)*coef_bmelt_dist(i,j)
 
              if (fbm(i,j)) then 
-                bmelt(i,j)=coef_talus*bmgrz(i,j)
+                bmelt(i,j)=coef_talus*bmgrz(i,j)*coef_bmelt_dist(i,j)
              endif
 
@@ -197 +217 @@
              !   en fonction du nombre de points flottants
 
-             bmelt(i,j)= ngr/4.*bmgrz(i,j)*coef_talus+(1.-ngr/4.)*bmelt(i,j)
+             bmelt(i,j)= (ngr/4.*bmgrz(i,j)*coef_talus)*coef_bmelt_dist(i,j)+(1.-ngr/4.)*bmelt(i,j)
 
 
@@ -204 +224 @@
        end do
     end do
-