Changeset 250 for trunk

Timestamp:

01/28/19 17:12:22 (5 years ago)

Author:

aquiquet

Message:

Fast atmospheric variability taken into account in climat_forcage using a glacial index

File:

• trunk/SOURCES/climat_forcage_mod.f90 (modified) (5 diffs)

trunk/SOURCES/climat_forcage_mod.f90

@@ -28 +28 @@
   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
+  real    :: r_atmvar              !< a ratio to account fast atmos variability between the 2 snapshots (in %: 0=>no fast var,1=>only fast var)
 
 
@@ -190 +191 @@
 
 ! Warning: when using climat_forcage we should in principle use a tpert which is a glacial index (ranging from about 0 to 1)
-!          here we can still use a coefT to scale the signal if needed
+!          here we should use the coefT to scale the signal if needed
 !          the coefbmb now is copy and paste from what we have done for Antarctica in Quiquet et al. GMD 2018
 
@@ -205 +206 @@
 
 
-  namelist/clim_forcage/clim_ref_file,ttr,forcage_file1,forcage_file2,typerun,lapserate,rappact,filforc,pertbmb,coeft,coefbmb
+  namelist/clim_forcage/clim_ref_file,ttr,forcage_file1,forcage_file2,typerun,lapserate,rappact,filforc,pertbmb,coeft,coefbmb,r_atmvar
 
   rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
@@ -224 +225 @@
   write(num_rep_42,*) 'coefT   = ', coefT
   write(num_rep_42,*) 'coefbmb = ', coefbmb
+  write(num_rep_42,*) 'r_atmvar = ', r_atmvar
   write(num_rep_42,*)'/'                            
   write(num_rep_42,*)
@@ -344 +346 @@
   endif   ! fin du test avant,apres,milieu
 
+  if (abs(Tafor).gt.2.) then
+     write(*,*) "A glacial index greater than 2? Really??"
+     STOP
+  endif
+  coeftime = (1-r_atmvar)*coeftime + r_atmvar*(1-Tafor)
 
 !if (geoplace(1:5).eq.'hemin') then