Changeset 108 for trunk/libf/dyn3d/conf_gcm.F

Timestamp:

Apr 12, 2011, 11:16:02 AM (14 years ago)

Author:

slebonnois

Message:

Sebastien Lebonnois: sponge layer et dissip horizontale.

File:

• trunk/libf/dyn3d/conf_gcm.F (modified) (1 diff)

trunk/libf/dyn3d/conf_gcm.F

@@ -303 +303 @@
 ! Parametres controlant la variation sur la verticale des constantes de
 ! dissipation.
-! Pour le moment actifs uniquement dans la version a 39 niveaux
-! avec ok_strato=y
-
-       dissip_factz=4.
-       dissip_deltaz=10.
-       dissip_zref=30.
-       CALL getin('dissip_factz',dissip_factz )
+! Actifs uniquement avec ok_strato=y
+
+       dissip_fac_mid=2.
+       dissip_fac_up=10.
+       dissip_deltaz=10.! Intervalle (km) pour le changement mid / up
+       dissip_hdelta=5. ! scale height (km) dans la zone de la transition(m)
+       dissip_pupstart=1.e3  ! pression (Pa) au bas la transition mid / up
+       CALL getin('dissip_fac_mid',dissip_fac_mid )
+       CALL getin('dissip_fac_up',dissip_fac_up )
        CALL getin('dissip_deltaz',dissip_deltaz )
-       CALL getin('dissip_zref',dissip_zref )
-
+       CALL getin('dissip_hdelta',dissip_hdelta )
+       CALL getin('dissip_pupstart',dissip_pupstart )
+
+! Parametres controlant la couche eponge
+! Actifs uniquement avec ok_strato=y
+!   mode = 0 : pas de sponge
+!   mode = 1 : u et v -> 0
+!   mode = 2 : u et v -> moyenne zonale
+!   mode = 3 : u, v et h -> moyenne zonale
        iflag_top_bound=1
+       mode_top_bound=3
        tau_top_bound=1.e-5
        CALL getin('iflag_top_bound',iflag_top_bound)
+       CALL getin('mode_top_bound',mode_top_bound)
        CALL getin('tau_top_bound',tau_top_bound)