Changeset 2126 for LMDZ5/trunk/libf/phylmd/lmdz1d.F90

Timestamp:

Oct 15, 2014, 2:03:57 AM (10 years ago)

Author:

fhourdin

Message:

Introduction du cas Dice couplé atmosphere/surface
+ nouveau paramètre de contrôle f_ri_cd_min, seuil minimum
sur la fonction f(Ri) en facteur du coefficient de traîné neutre.
Par défaut : f_ri_cd_min=0.1 (comme avant)

Introduction of the coupled atmosphere/surface dice case.
+ new parameter f_ri_cd_min, minimum threshold on the f(Ri) fonction
that multiply the neutral drag coefficient.

File:

• LMDZ5/trunk/libf/phylmd/lmdz1d.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/lmdz1d.F90

@@ -112 +112 @@
         logical :: forcing_twpice  = .false.
         logical :: forcing_amma    = .false.
+        logical :: forcing_dice    = .false.
         logical :: forcing_GCM2SCM = .false.
         logical :: forcing_GCSSold = .false.
@@ -178 +179 @@
 !---------------------------------------------------------------------
       logical :: ok_writedem =.true.
+      real :: sollw_in = 0.
+      real :: solsw_in = 0.
       
 !---------------------------------------------------------------------
@@ -262 +265 @@
 !             initial profiles from AMMA nc file 
 !             LS convergence, omega and surface fluxes imposed from AMMA file  
+!forcing_type = 7 ==> forcing_dice = .true.
+!             initial profiles and large scale forcings in dice_driver.nc
+!             Different stages: soil model alone, atm. model alone
+!             then both models coupled
 !forcing_type = 40 ==> forcing_GCSSold = .true.
 !             initial profile from GCSS file
@@ -296 +303 @@
       elseif (forcing_type .eq.6) THEN
        forcing_amma = .true.
+      elseif (forcing_type .eq.7) THEN
+       forcing_dice = .true.
       elseif (forcing_type .eq.40) THEN
        forcing_GCSSold = .true.
@@ -319 +328 @@
 
         type_ts_forcing = 0
-        if (forcing_toga.or.forcing_sandu.or.forcing_astex)                 &
+        if (forcing_toga.or.forcing_sandu.or.forcing_astex .or. forcing_dice)                 &
      &    type_ts_forcing = 1
 
@@ -366 +375 @@
 ! Special case for amma which lasts less than one day : 64800s !! (MPL 20120216)
       IF(forcing_type .EQ. 6) fnday=64800./86400.
+!     IF(forcing_type .EQ. 6) fnday=50400./86400.
       annee_ref = anneeref
       mois = 1
@@ -391 +401 @@
      & (year_ini_amma,mth_ini_amma,day_ini_amma,heure_ini_amma              &
      & ,day_ju_ini_amma)
-
+      ELSEIF (forcing_type .eq.7) THEN
+! Convert the initial date of DICE to Julian day
+      call ymds2ju                                                         &
+     & (year_ini_dice,mth_ini_dice,day_ini_dice,heure_ini_dice             & 
+     & ,day_ju_ini_dice)
       ELSEIF (forcing_type .eq.59) THEN
 ! Convert the initial date of Sandu case to Julian day
@@ -405 +419 @@
 
       ELSEIF (forcing_type .eq.61) THEN
-
 ! Convert the initial date of Arm_cu case to Julian day
       call ymds2ju                                                          &