Changeset 490

Timestamp:

02/21/24 16:45:10 (4 months ago)

Author:

dumas

Message:

depth_shelf_min : minimal depth (-200m) used to compute basal melting in bmelt-beckmann-gcm_mod

Location:

branches/GRISLIv3

Files:

• SOURCES/bmelt-beckmann-gcm_mod.f90 (modified) (2 diffs)
• launching_area/Master-files/Param/anteis1_param_list_lhs-nolin.dat (modified) (1 diff)

branches/GRISLIv3/SOURCES/bmelt-beckmann-gcm_mod.f90

@@ -81 +81 @@
     integer :: nmax ! to avoid inifite loop 
     integer :: err
-
-    namelist/bmelt_beckmann_gcm_mod/K_t,TO_file,SO_file,Bassin_file,nbassins,deltaT_bool,deltaT_file,bmelt_empty,nonlocal_bool,grlmelt_bool
+    real :: depth_shelf_min ! minimal depth used to compute bmelt
+    integer :: ksup
+
+    namelist/bmelt_beckmann_gcm_mod/K_t,TO_file,SO_file,Bassin_file,nbassins,deltaT_bool,deltaT_file,bmelt_empty,nonlocal_bool,grlmelt_bool,depth_shelf_min
     
     rewind(num_param)        ! loop back at the beginning of the param_list.dat file
@@ -326 +328 @@
    enddo
    !$OMP END DO
+
+   ! limitation de la temperature a une profondeur depth_shelf_min pour eviter tx fonte tres faible au front des ice-shelves
+   ! recherche de l'indice de profondeur 
+   do k=1,nzoc
+      if ( zoc(k).GT.depth_shelf_min) ksup = k + 1
+   enddo
+   if (ksup.LT.2) print*,'bmelt-beckmann-gcm_mod, check depth_shelf_min : ksup = ',ksup
+
+   do k=1,ksup - 1
+      !$OMP WORKSHARE
+      temp_ocean(:,:,k) = temp_ocean(:,:,ksup)
+      salinity_ocean(:,:,k) = salinity_ocean(:,:,ksup)
+      !$OMP END WORKSHARE
+   enddo
+            
+   
    
    !$OMP WORKSHARE

branches/GRISLIv3/launching_area/Master-files/Param/anteis1_param_list_lhs-nolin.dat

@@ -231 +231 @@
 nonlocal_bool = 1
 grlmelt_bool = 0
+depth_shelf_min = -200.
 /
 !---------------------------------------------------------