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

Timestamp:

Jul 14, 2017, 10:23:40 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Parametrization of drag by copses
Need version 4465 of ORCHIDEE at least

6. Cheruy

Location:

LMDZ5/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/conf_phys_m.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -189 +189 @@
     INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
     INTEGER,SAVE :: iflag_pbl_split_omp
+!FC
+    INTEGER,SAVE :: ifl_pbltree_omp
+    REAL,SAVE :: Cd_frein_omp
+!FC
     INTEGER,SAVE :: iflag_order2_sollw_omp
     INTEGER, SAVE :: lev_histins_omp, lev_histLES_omp 
@@ -1404 +1408 @@
     iflag_pbl_omp = 1
     CALL getin('iflag_pbl',iflag_pbl_omp)
+
+!FC
+    !Config Key  = ifl_pbltree
+    !Config Desc = drag from trees 0 no activated
+    !Config Def  = 0
+    !Config Help =
+    !
+    ifl_pbltree_omp = 0
+    CALL getin('ifl_pbltree',ifl_pbltree_omp)
+!FC
+    !Config Key  = Cd_frein
+    !Config Desc = drag from trees 
+    !Config Def  = 7.5E-02 (valeur Masson mais fait planter avec des LAI eleves)
+    !Config Help =
+    !
+    Cd_frein_omp = 7.5E-02
+    CALL getin('Cd_frein',Cd_frein_omp)
+
     !
     !Config Key  = iflag_pbl_split
@@ -2199 +2221 @@
     iflag_pbl = iflag_pbl_omp
     iflag_pbl_split = iflag_pbl_split_omp
+!FC
+    ifl_pbltree = ifl_pbltree_omp
+    Cd_frein    =Cd_frein_omp
     iflag_order2_sollw = iflag_order2_sollw_omp
     lev_histhf = lev_histhf_omp
@@ -2440 +2465 @@
        CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1)
     ENDIF
+
+    ! ORCHIDEE must be activated for ifl_pbltree=1
+    IF (.NOT. ok_veget .AND. ifl_pbltree==1) THEN
+       WRITE(lunout,*)' ORCHIDEE must be activated for ifl_pbltree=1'
+       CALL abort_physic('conf_phys','ok_veget and ifl_pbltree not coherent',1)
+    END IF
 
     !$OMP MASTER
@@ -2566 +2597 @@
     write(lunout,*)' freq_calNMC = ',freq_calNMC
     write(lunout,*)' iflag_pbl = ', iflag_pbl
+!FC
+    write(lunout,*)' ifl_pbltree = ', ifl_pbltree
+    write(lunout,*)' Cd_frein = ', Cd_frein
     write(lunout,*)' iflag_pbl_split = ', iflag_pbl_split
     write(lunout,*)' iflag_order2_sollw = ', iflag_order2_sollw