Changeset 1724 for LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Feb 7, 2013, 3:56:38 PM (11 years ago)

Author:

musat

Message:

Ajout coefficient beta (=Evap/Evap? potentielle) fonction de la latitude
et de la saison (voir these AbdulKhadre? Traore).

Frederic Hourdin
IM

File:

• LMDZ5/trunk/libf/phylmd/conf_phys_m.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -1 @@
-
 !
 ! $Id: conf_phys.F90 1668 2012-10-12 10:47:37Z idelkadi $
@@ -80 +79 @@
 
   character (len = 6),SAVE  :: type_ocean_omp, version_ocean_omp, ocean_omp
+  character (len = 10),SAVE  :: type_veget_omp
   CHARACTER(len = 8),SAVE   :: aer_type_omp
-  logical,SAVE              :: ok_veget_omp, ok_newmicro_omp
+  logical,SAVE              :: ok_newmicro_omp
   logical,SAVE        :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp        
   logical,SAVE        :: ok_LES_omp   
@@ -217 +217 @@
 !Config Help = Type de modele de vegetation utilise
 !
-  ok_veget_omp = .false.
-  call getin('VEGET', ok_veget_omp)
+  type_veget_omp ='orchidee'
+  call getin('VEGET', type_veget_omp)
 !
 !Config Key  = OK_journe
@@ -289 +289 @@
   aerosol_couple_omp = .false.
   CALL getin('aerosol_couple',aerosol_couple_omp)
-
 !
 !Config Key  = flag_aerosol
@@ -1675 +1674 @@
     type_ocean = type_ocean_omp
     version_ocean = version_ocean_omp
-    ok_veget = ok_veget_omp
+
+    ok_veget=.true.
+    type_veget=type_veget_omp
+    if (type_veget=='n' .or. type_veget=='bucket' .or. type_veget=='betaclim' ) &
+         & then
+         ok_veget=.false.
+    endif
+
     ok_newmicro = ok_newmicro_omp
     ok_journe = ok_journe_omp
@@ -1796 +1802 @@
   write(lunout,*)' Type ocean = ', type_ocean
   write(lunout,*)' Version ocean = ', version_ocean
-  write(lunout,*)' Config veget = ', ok_veget
+  write(lunout,*)' Config veget = ', ok_veget,type_veget
   write(lunout,*)' Sortie journaliere = ', ok_journe
   write(lunout,*)' Sortie haute frequence = ', ok_hf