Changeset 187 for branches/iLoveclim/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

Timestamp:

03/15/18 14:18:24 (6 years ago)

Author:

aquiquet

Message:

Grisli-iloveclim branch merged to trunk at revision 185

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90 (modified) (12 diffs)

branches/iLoveclim/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

@@ -62 +62 @@
 real,dimension(:),allocatable :: tpert          !> temperature for climate forcing
 real,dimension(:),allocatable :: spert          !> sea surface perturbation
+real,dimension(:),allocatable :: bpert          !> basal melt index perturbation
 real :: coefT                    !> pour modifier l'amplitude de la perturb. T
-character(len=120) :: filforc    !> nom du fichier forcage
+character(len=150) :: filforc    !> nom du fichier forcage
 integer :: pertsmb               !> boolean: do we modify the smb?
 real    :: rapsmb                !> if we modify the smb this is the equivalent of rappact
+integer :: pertbmb               !> boolean: do we modify the bmelt under ice shelves?
+real    :: coefbmb               !> if we modify the bmelt, we read a perturbation (based on dD, insolation, etc.) and apply this coefficient
 
 
@@ -99 +102 @@
   !aurel for Tafor:
   character(len=8) :: control      !label to check clim. forc. file (filin) is usable
-  character(len=80):: filin
+  character(len=100):: filin
 
   namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file
@@ -129 +132 @@
   !   bm(:,:) = bm(:,:) - 5.                     ! pour faire un masque a l'exterieur du Groenland actuel
   !end where
-  where ((H(:,:).lt.1.).and.(Bsoc(:,:).gt.0.).and.(bm(:,:).gt.0.))
-     bm(:,:) = -0.1
-  end where
+  !where ((H(:,:).lt.1.).and.(Bsoc(:,:).gt.0.).and.(bm(:,:).gt.0.))
+  !   bm(:,:) = -0.1
+  !end where
 !cdc test debug Hemin15 et Greeneem15
 !  where (bm(:,:).lt.-1000) bm(:,:)=0.
@@ -278 +281 @@
   end if
   
+  if (.not.allocated(bpert)) then
+     allocate(bpert(nft),stat=err)
+     if (err/=0) then
+        print *,"erreur a l'allocation du tableau Bpert",err
+        stop 4
+     end if
+  end if
+  
   do i=1,nft
-     read(num_forc,*) tdate(i),spert(i),tpert(i)
+     if (pertbmb==1) then
+        read(num_forc,*) tdate(i),spert(i),tpert(i),bpert(i)
+     else
+        read(num_forc,*) tdate(i),spert(i),tpert(i)
+        bpert(i)=0.
+     endif
   end do
   close(num_forc)
   
   tpert(:)=tpert(:)*coefT
-  
+  bpert(:)=max(1.+bpert(:)*coefbmb,0.01) !freezing is not allowed
 
   
@@ -325 +341 @@
   implicit none
   namelist/lapse_rates/T_lapse_rate
-  namelist/clim_pert_massb/coefT,filforc,pertsmb,rapsmb
+  namelist/clim_pert_massb/coefT,filforc,pertsmb,rapsmb,pertbmb,coefbmb
   
   rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
@@ -353 +369 @@
   write(num_rep_42,*) 'pertsmb        = ', pertsmb
   write(num_rep_42,*) 'rapsmb         = ', rapsmb
+  write(num_rep_42,*) 'pertbmb        = ', pertbmb
+  write(num_rep_42,*) 'coefbmb         = ', coefbmb
   write(num_rep_42,*)'/'
   write(num_rep_42,*)
@@ -388 +406 @@
            tafor=tpert(1)
            sealevel=spert(1)
+           coefbmshelf=bpert(1)
            ift=1
 
@@ -393 +412 @@
            tafor=tpert(nft)
            sealevel=spert(nft)
+           coefbmshelf=bpert(nft)
            ift=nft
            
@@ -402 +422 @@
                  sealevel=spert(i)+(spert(i+1)-spert(i))*    &
                       (time-tdate(i))/(tdate(i+1)-tdate(i))
+                 coefbmshelf=bpert(i)+(bpert(i+1)-bpert(i))*    &
+                      (time-tdate(i))/(tdate(i+1)-tdate(i))
                  ift=i
                  goto 100
@@ -409 +431 @@
 100     continue
 
-        
         Tann (:,:) = Ta0 (:,:) + T_lapse_rate * (S(:,:)-S0(:,:)) +Tafor
         Ts(:,:)    = Tann(:,:)
@@ -420 +441 @@
            end if
         end if
+
+        !afq -- as in Pollard and Deconto, melt is more efficient at high temperature... ->
+        if (coefbmshelf.gt.1.) then
+           coefbmshelf = 1. + ( coefbmshelf - 1. ) * 1.5 !afq.... tuning for the degla?...
+        end if
         
-        ! coefmshelf est un coefficient qui fait varier bmgrz et bmshelf en fonction de tafor
-        coefbmshelf=(1.+tafor/10.)       ! coefbmshelf=0 pour tafor=-10deg
-        coefbmshelf=max(coefbmshelf,0.)  
-        coefbmshelf=min(coefbmshelf,2.)
+!afq --old:        ! coefmshelf est un coefficient qui fait varier bmgrz et bmshelf en fonction de tafor
+!afq --old:        coefbmshelf=(1.+tafor/10.)       ! coefbmshelf=0 pour tafor=-10deg
+!afq --old:        coefbmshelf=max(coefbmshelf,0.)  
+!afq --old:        coefbmshelf=min(coefbmshelf,2.)
         
 !        coefbmshelf=1.  !afq for calibration Greenland
@@ -441 +467 @@
   implicit none
   integer             :: k_snap               ! pour calculer les indices de temps
-  real                :: time_bis             ! pour repliquer les dernieres annees
   integer             :: k