Changeset 67

Timestamp:

05/27/16 17:08:36 (8 years ago)

Author:

aquiquet

Message:

Small bug correction in shelf melting for clim-pert-like experiments

Location:

trunk/SOURCES

Files:

• GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90 (modified) (4 diffs)
• bmelt-seuil-profondeur_mod.f90 (modified) (1 diff)

trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

@@ -16 +16 @@
 
 
-use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel
+use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel,coefbmshelf
 !use lect_climref_Ice2sea
 use netcdf
@@ -130 +130 @@
   !end where
   where ((H(:,:).lt.1.).and.(Bsoc(:,:).gt.0.).and.(bm(:,:).gt.0.))
-     bm(:,:) = -0.01
+     bm(:,:) = -0.1
   end where
 !cdc test debug Hemin15 et Greeneem15
@@ -380 +380 @@
 
   integer i,ift
-  
+
   select case (massb_time)
      case(0)
@@ -421 +421 @@
         end if
         
+        coefbmshelf=(1.+tafor/10.)
+        coefbmshelf=max(coefbmshelf,0.)  
+        coefbmshelf=min(coefbmshelf,2.)
+  
      case(1)
         call massb_Ice2sea_RCM

trunk/SOURCES/bmelt-seuil-profondeur_mod.f90

@@ -118 +118 @@
 
 talus_nochange : if (Bsoc0(i,j).lt.depth_talus) then 
+           coef_talus = 1         ! pas de changement au dessus de l'ocecan profond
+        else
            coef_talus = coefbmshelf
-        else
-           coef_talus = 1         ! pas de changement au dessus de l'ocecan profond
         endif talus_nochange