Changeset 55 for trunk/SOURCES/GrIce2sea_files

Timestamp:

04/05/16 14:52:40 (8 years ago)

Author:

aquiquet

Message:

SMB forced: T perturbation from ice cores which may affect SMB (rappact)

Location:

trunk/SOURCES/GrIce2sea_files

Files:

• climat_GrIce2sea_years_mod.f90 (modified) (2 diffs)
• climat_GrIce2sea_years_perturb_mod.f90 (modified) (7 diffs)

trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90

@@ -116 +116 @@
 !cdc test debug Hemin15 et Greeneem15
 !  where (bm(:,:).lt.-1000) bm(:,:)=0.
-
+  where (bm(:,:).eq.0) bm(:,:)=-10. !afq
 
   acc(:,:) = 0.
@@ -138 +138 @@
 !cdc test debug Hemin15
 !  where (Tann(:,:).lt.-100.) Tann(:,:)=10.
-
+  where (Tann(:,:).lt.-100.) Tann(:,:)=10. !afq
 
   ta0(:,:)   = Tann(:,:)

trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

@@ -58 +58 @@
 
 ! aurel, pour climat type perturb:
-integer nft             !< nombre de lignes a lire dans le fichier forcage climatique
-real,dimension(:),allocatable :: tdate          !< time for climate forcing
-real,dimension(:),allocatable :: tpert          !< temperature for climate forcing
-real,dimension(:),allocatable :: spert          !< sea surface perturbation
-real :: coefT                   !< pour modifier l'amplitude de la perturb. T
-character(len=120) :: filforc    !< nom du fichier forcage
+integer nft             !> nombre de lignes a lire dans le fichier forcage climatique
+real,dimension(:),allocatable :: tdate          !> time for climate forcing
+real,dimension(:),allocatable :: tpert          !> temperature for climate forcing
+real,dimension(:),allocatable :: spert          !> sea surface perturbation
+real :: coefT                    !> pour modifier l'amplitude de la perturb. T
+character(len=120) :: 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
 
 
@@ -130 +132 @@
 !cdc test debug Hemin15 et Greeneem15
 !  where (bm(:,:).lt.-1000) bm(:,:)=0.
-
+  where (bm(:,:).eq.0) bm(:,:)=-10. !afq
+  
   acc(:,:) = 0.
   abl(:,:) = 0.
@@ -151 +154 @@
 !cdc test debug Hemin15
 !  where (Tann(:,:).lt.-100.) Tann(:,:)=10.
-
+  where (Tann(:,:).lt.-100.) Tann(:,:)=10. !afq
 
   ta0(:,:)   = Tann(:,:)
@@ -228 +231 @@
 
 ! copie la valeur de reference dans bm_0
-     bm_0(:,:) = bm(:,:)
+!     bm_0(:,:) = bm(:,:)
   endif
-
+! copie la valeur de reference dans bm_0
+  bm_0(:,:) = bm(:,:) !afq marion dufresne: we do it even if we don't read any snapshot
+  
   filin=trim(dirforcage)//trim(filforc)
 
@@ -318 +323 @@
   implicit none
   namelist/lapse_rates/T_lapse_rate
-  namelist/clim_pert_massb/coefT,filforc
+  namelist/clim_pert_massb/coefT,filforc,pertsmb,rapsmb
   
   rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
@@ -344 +349 @@
   write(num_rep_42,*) 'coefT        = ', coefT
   write(num_rep_42,'(A,A)') ' filforc      = ', filforc
+  write(num_rep_42,*) 'coefT        = ', pertsmb
+  write(num_rep_42,*) 'coefT        = ', rapsmb
   write(num_rep_42,*)'/'
   write(num_rep_42,*)
@@ -403 +410 @@
         Tann (:,:) = Ta0 (:,:) + T_lapse_rate * (S(:,:)-S0(:,:)) +Tafor
         Ts(:,:)    = Tann(:,:)
+
+        ! aurel marion dufresne: we might want to decrease the SMB during glacials..?
+        if (pertsmb.eq.1) then
+           !bm(:,:) = bm_0(:,:) * exp( rapsmb *(Tann(:,:)-Ta0(:,:)))
+           ! afq, to match present day smb, I don't want to take S-elev feedb. into account
+           bm(:,:) = bm_0(:,:) * exp( rapsmb * Tafor )
+        end if
+        
      case(1)
         call massb_Ice2sea_RCM