Changeset 844 for trunk/LMDZ.GENERIC_2.0/deftank

Timestamp:

Nov 18, 2012, 2:03:20 PM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. Added a new folder for the interface with LMDZ.COMMON and LMDZ5

Location:

trunk/LMDZ.GENERIC_2.0

Files:

• . (added)
• deftank (added)
• deftank/gray_lmdz.generic_lmdz5 (copied) (copied from trunk/LMDZ.GENERIC/deftank/gray_lmdz.generic_lmdz5)
• deftank/gray_lmdz.generic_lmdz5/run.def (modified) (10 diffs)
• deftank/gray_lmdz.generic_lmdz5/run.def.128.96 (added)
• deftank/linkgen.sh (copied) (copied from trunk/LMDZ.GENERIC/libf/phystd/for_lmdz5/linkgen.sh)
• deftank/z2sig.earth (added)

trunk/LMDZ.GENERIC_2.0/deftank/gray_lmdz.generic_lmdz5/run.def

@@ -1 +1 @@
 
+#########
+######### THIS IS INTENDED TO WORK BOTH WITH LMDZ.GENERIC + old LMDZ core
+#########                                AND LMDZ.GENERIC + new LMDZ core
+#########
+
+### A INCLURE ?
+## calend=venus  ## calendrier specifique
+## resetvarc=n   ## reinit des variables de controle
+## iflag_trac=0  ## avec ou sans traceurs
+
+################################################################
 #### LMDZ5
 #### IN GCM.DEF
+#### -- those parameters are ignored when using old LMDZ core
+################################################################
 
 ## period appel a la dissipation (multiple iperiod, 0 si calcul automatique)
+## (remplace ce qui etait idissip dans l'ancienne dynamique)
 dissip_period = 0
 dissip_period = 2 
@@ -21 +35 @@
 read_start=y 
 
-##  Avec ou sans strato
+##  Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horizontale
 ok_strato=n
+
+### Dissipation horizontale
+#dissip_fac_mid=2.
+#dissip_fac_up=2.
+## deltaz et hdelta en km
+#dissip_deltaz=10.
+#dissip_hdelta=5.
+## pupstart en Pa
+#dissip_pupstart=1.e3
 
 #  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
 iflag_top_bound=0
+
+########################################
+## LMDZ.COMMON
+##    1: dans les 4 derniers niveaux
+##    2: dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+#iflag_top_bound=1
+## Mode Couche eponge
+##   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
+#mode_top_bound=1
+#######################################
+
 #  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=1.e-5
 tau_top_bound=5.e-5
 
+################################################################
 #### LMDZ5
 #### IN RUN.DEF
+#### -- those parameters are ignored when using old LMDZ core
+################################################################
 
 ## Type de calendrier utilise
@@ -40 +81 @@
 ## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
 dayref=10120
-dayref=40
+#dayref=40
 
 ##  Annee de l'etat  initial (   avec  4  chiffres   )
@@ -72 +113 @@
 iecri=120
 iecri=480
-iecri=4800
+#iecri=4800
+#iecri=9600
 
 ## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.)
@@ -83 +125 @@
 #### AND NOW THE OLD run.def ####
 #################################
-
+## -----> REMOVED
+## was present and conflicting with LMDZ5
+## (although useless in old OLD LMDz DYN CORE
+## - iecri [periode d'ecriture du fichier histoire (en jour)]
+## - periodav [periode de stockage fichier histmoy (en jour) dynzon?]
+#################################
+#################################
 
 #----------------------------------------------------------------------- 
@@ -91 +139 @@
 # Nombre de jours d'integration                                         
      nday = 720
-     nday = 20
-     nday = 100
-     nday = 1800
+#     nday = 20
+#     nday = 100
+#     nday = 1800
 
 # nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )      
  day_step = 240
- day_step = 480
+# day_step = 480
+# day_step = 960
 
 # periode pour le pas Matsuno (en pas)                                  
@@ -105 +154 @@
   iconser=120
 
-############# VOIR CI-DESSUS ne sert pas dans generic
-############# ... mais sert dans LMDZ5 !!!!!!
-############# ... donc on commente ...
-## periode d'ecriture du fichier histoire (en jour)                      
-#    iecri=1
-#
-## periode de stockage fichier histmoy (en jour) dynzon?                        
-# periodav=10.
-############# VOIR CI-DESSUS
-############# VOIR CI-DESSUS
-
-############# VOIR CI-DESSUS dissip_period
-## periode de la dissipation (en pas)                                    
-#  idissip=5
-############# VOIR CI-DESSUS
+# periode de la dissipation (en pas)                                    
+# -- this is ignored by LMDZ5 which rely on dissip_period instead
+  idissip=5
 
 # choix de l'operateur de dissipation (star ou  non star )              
- lstardis=.true. #ou y ???
+ lstardis=.true. 
 
 # avec ou sans coordonnee hybrides
@@ -137 +174 @@
 
 # temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
- tetagdiv=1000. 2500. 4000.
+ tetagdiv=2500. 4000.
 
 # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
- tetagrot=1000. 2500. 4000.
+ tetagrot=2500. 4000.
 
 # temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
- tetatemp=1000. 2500. 4000.
+ tetatemp=2500. 4000.
 
 # coefficient pour gamdissip                                            
@@ -149 +186 @@
 
 # choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
-  purmats=.false. # ou n ???
+  purmats=.false. 
 
 # avec ou sans physique                                                 
@@ -197 +234 @@
   ysinus= .false.
 
-###################### MARCHE SEULEMENT AVEC L ANCIENNE DYNAMIQUE
-## Avec sponge layer
-#  callsponge  = .true.
-# 
-## Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
-#  mode_sponge= 2
-#
-## Sponge:  hauteur de sponge (km)
-#  hsponge= 50
-#
-## Sponge:  tetasponge (secondes)
-#  tetasponge = 50000
-###################### MARCHE SEULEMENT AVEC L ANCIENNE DYNAMIQUE
+###################### SPONGE LAYER IN OLD DYNAMICS
+###################### -- this is ignored by LMDZ5
+# Avec sponge layer
+  callsponge  = .true.
+# Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
+  mode_sponge= 2
+# Sponge:  hauteur de sponge (km)
+  hsponge= 50
+# Sponge:  tetasponge (secondes)
+  tetasponge = 50000
+###################### 
 
 # some definitions for the physics, in file 'callphys.def'