Changeset 1171 for trunk


Ignore:
Timestamp:
Jan 30, 2014, 12:16:55 PM (11 years ago)
Author:
aslmd
Message:

LMDZ.UNIVERSAL. clean and tidied reference run.def

Location:
trunk/LMDZ.UNIVERSAL
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.UNIVERSAL/README.aymeric

    r1105 r1171  
     1
     2
     3makelmdz_fcm -s 1 -d 64x48x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm
     4makelmdz_fcm -s 1 -d 128x96x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm
     5#####################################
     6## REF SANDRINE
     7makelmdz_fcm -cpp NODYN -s 1 -d 8x128x64 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm
     8#####################################
     9
     10
     11
     12
     13
     14makelmdz_fcm -s 1 -d 64x48x64 -b 17x23 -arch GNOMEp -parallel mpi -p generic gcm
     15
     16test rings
     17makelmdz_fcm -cpp NODYN -s 1 -d 8x48x32 -b 17x23 -arch GNOMEp -parallel mpi -p generic gcm
     18makelmdz_fcm -cpp NODYN -s 1 -d 8x48x32 -b 20x30 -arch GNOMEp -parallel mpi -p generic gcm
     19
     20-----------------------------------
    121
    222
  • trunk/LMDZ.UNIVERSAL/deftank/PROTO_Saturn_Article/run.def

    r1165 r1171  
     1###########################################################################
     2### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ### RUN.DEF ###
     3###########################################################################
     4### A MERGED VERSION TO BE USED WITH OLD OR RECENT LMDZ CORE
     5### i.e. compiles gcm.def+run.def in LMDZ5 with run.def in old planeto LMDZ
     6### (*) indicates same parameters in old and recent LMDZ core
     7###########################################################################
    18
    2 #########
    3 ######### THIS IS INTENDED TO WORK BOTH WITH LMDZ.GENERIC + old LMDZ core
    4 #########                                AND LMDZ.GENERIC + new LMDZ core
    5 #########
     9###########################################################################
     10### NUMERICAL INTEGRATIONS
     11###########################################################################
     12### nday --> number of days of integration (*)
     13### iperiod --> period for the Matsuno step (in dynamical steps) (*)
     14### day_step --> number of steps per day (multiple of iperiod) (*)
     15### purmats --> temporal integration: Matsuno or Matsuno-leapfrog? (*)
     16### grireg --> regular grid or not (*)
     17### adjust --> flag to activate load optimization
     18### use_filtre_fft --> flag to activate FFT filtering
     19###########################################################################
     20nday = 65000
     21iperiod = 5
     22day_step = 70
     23purmats = .false.
     24grireg = .true.
     25adjust = n
     26use_filtre_fft = n
    627
    7 ### A INCLURE ?
    8 ## calend=venus  ## calendrier specifique
    9 ## resetvarc=n   ## reinit des variables de controle
    10 ## iflag_trac=0  ## avec ou sans traceurs
     28###########################################################################
     29### WHAT RELATES TO PHYSICS
     30###########################################################################
     31### iphysiq --> period for calling physics (in dynamical steps) (*)
     32### ecritphy --> frequency for outputs in diagfi (in dynamical steps) (*)
     33### planet_type --> the planet
     34### disvert_type --> type of vertical coordinates (2 for planeto)
     35### hybrid --> use hybrid coordinates or not? (*)
     36### iflag_phys=0 --> no physics (e.g. shallow water)
     37### iflag_phys=1 --> physics (e.g. phylmd)
     38### iflag_phys=2 --> simple newtonian relaxation
     39### physic --> flag with or without physics (*)
     40###########################################################################
     41iphysiq = 14
     42ecritphy = 7000
     43planet_type = generic
     44disvert_type = 2
     45hybrid = .false.
     46iflag_phys = 1
     47physic = .true.
    1148
    12 ################################################################
    13 #### LMDZ5
    14 #### IN GCM.DEF
    15 #### -- those parameters are ignored when using old LMDZ core
    16 ################################################################
     49###########################################################################
     50### INITIALIZATION and NUDGING
     51###########################################################################
     52### read_start --> with or without start files ?
     53###                (without: initialization by iniacademic)
     54### raz_date=0 --> use the date in start file
     55### raz_date=1 --> erase the date and start from scratch
     56### ok_guide --> flag for nudging towards winds (needs relax.nc)
     57### facwind --> time constant for relaxation
     58###             5.e8 medium 1.e7-5.e6 strong 1.e6 very fast
     59###########################################################################
     60read_start = y
     61raz_date = 0
     62ok_guide = n
     63facwind = 1.e6
    1764
    18 ## period appel a la dissipation (multiple iperiod, 0 si calcul automatique)
    19 ## (remplace ce qui etait idissip dans l'ancienne dynamique)
    20 #dissip_period = 0
    21 #dissip_period = 2
     65###########################################################################
     66### HORIZONTAL DISSIPATION
     67###########################################################################
     68### dissip_period --> period call to dissipation (in dynamical steps)
     69### - must be a multiple of iperiod
     70### - 0 means automatic computation
     71### - use idissip instead if old LMDZ core (idissip ignored in LMDZ5)
     72### dissip_fac_mid --> enhancement factor in middle atmosphere
     73### - set to 1 to prevent dissipation enhancement
     74### lstardis --> star or not star dissipation operator (*)
     75### nitergdiv --> number of iterations for gradiv operator (*)
     76### nitergrot --> number of iterations for nxgradrot operator (*)
     77### niterh --> number of iterations for divgrad operator (*)
     78### tetagdiv --> dissipation timescale of smallest wvl: u,v (gradiv) (*)
     79### tetagrot --> dissipation timescale of smallest wvl: u,v (nxgradrot) (*)
     80### tetatemp --> dissipation timescale of smallest wvl: h (divgrad) (*)
     81### coefdis --> coefficient for gamdissip (*)
     82###########################################################################
    2283dissip_period = 5
     84dissip_fac_mid = 1.
     85lstardis = .true.
     86nitergdiv = 1
     87nitergrot = 2
     88niterh = 2
     89tetagdiv = 5400. 8000. 160000. 1000. 2500. 4000.
     90tetagrot = 7200. 8000. 160000. 1000. 2500. 4000.
     91tetatemp = 7200. 8000. 160000. 1000. 2500. 4000.
    2392
    24 ## avec ou sans physique
    25 ## 0: pas de physique (e.g. en mode Shallow Water)
    26 ## 1: avec physique (e.g. physique phylmd)
    27 ## 2: avec rappel newtonien dans la dynamique                                         
    28 iflag_phys=2
    29 iflag_phys=1
     93###########################################################################
     94### ENHANCEMENT OF HORIZONTAL DISSIPATION IN UPPERMOST LEVELS
     95###########################################################################
     96### - not activated if ok_strato=n
     97### - dissip_fac_up --> enhancement factor
     98### - dissip_deltaz --> altitude range of transition (km)
     99### - dissip_hdelta --> scale height at transition (km)
     100### - dissip_pupstart --> pressure at bottom of transition region (Pa)
     101###########################################################################
     102ok_strato = n
     103dissip_fac_up = 5
     104dissip_deltaz = 150.
     105dissip_hdelta = 50.
     106dissip_pupstart = 50.
    30107
    31 ## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ?
    32 ## (sans fichiers de demarrage, initialisation des champs par iniacademic
    33 ##  dans la dynamique)
    34 #read_start=n
    35 read_start=y
     108###########################################################################
     109### DAMPING ("SPONGE") LAYER AT MODEL TOP
     110###########################################################################
     111### iflag_top_bound=0 for no sponge
     112### iflag_top_bound=1 for sponge over 4 topmost layers
     113### iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
     114### mode_top_bound=0 for no sponge
     115### mode_top_bound=1 for u,v --> 0
     116### mode_top_bound=2 for u,v --> zonal mean
     117### mode_top_bound=3 for u,v,h --> zonal mean
     118### tau_top_bound --> damping frequency at last layer
     119### e-5 medium / e-4 strong yet reasonable / e-3 very strong
     120###########################################################################
     121iflag_top_bound = 1
     122mode_top_bound = 1
     123tau_top_bound = 1.e-4
    36124
    37 ##  Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horizontale
    38 ok_strato=n
     125###########################################################################
     126### OUTPUT FROM DYNAMICS
     127###########################################################################
     128### prt_level --> level of control prints (1-10)
     129### ok_dyn_ins --> flag instantaneous outputs (dyn_hist.nc)
     130### iecri --> frequency of instantaneous outputs (in dynamical steps)
     131### ok_dyn_ave --> flag averaged outputs (dyn_hist_ave.nc)
     132### periodav --> period of output of averages in dynamics (in days)
     133### ok_dynzon --> flag outputs dynzon
     134### NB: iecri and periodav used to be in old LMDZ core although useless
     135### iconser --> frequency of outputs of control variables (in dyn steps)(*)
     136###########################################################################
     137prt_level = 1
     138ok_dyn_ins = n
     139iecri = 480
     140ok_dyn_ave = n
     141periodav = 1
     142ok_dynzon = 0
     143iconser = 120
    39144
    40 ### Dissipation horizontale
    41 dissip_fac_mid=1. # for no var dissip (see inidissip)
    42 dissip_fac_up=1. # for no var dissip (see inidissip)
    43 ## deltaz et hdelta en km
    44 #dissip_deltaz=10.
    45 #dissip_hdelta=5.
    46 ## pupstart en Pa
    47 #dissip_pupstart=1.e3
     145###########################################################################
     146### ZOOMED GRID (same parameters in old and recent LMDZ core)
     147###########################################################################
     148### clon --> longitude of center of zoom (deg) (*)
     149### clat --> latitude of center of zoom (deg) (*)
     150### grossismx --> zoom factor in longitude (*)
     151### grossismy --> zoom factor in latitude (*)
     152### dzoomx --> longitude extent of zoom zone (fraction of total area) (*)
     153### dzoomy --> latitude extent of zoom zone (fraction of total area) (*)
     154### taux --> "stifness" of zoom in longitude (*)
     155### tauy --> "stifness" of zoom in latitude (*)
     156### fxyhypb --> f(y) hyperbolic (T) or sinusoidal (F) (*)
     157### ysinus --> f(y) = sin(lat) (T) or lat (F) (*)
     158###########################################################################
     159clon = 0.
     160clat = 0.
     161grossismx = 1.
     162grossismy = 1.
     163dzoomx = 0.
     164dzoomy = 0.
     165taux = 2.
     166tauy = 2.
     167fxyhypb = .false.
     168ysinus = .false.
    48169
    49 ########################################
    50 ## LMDZ.COMMON
    51 ## iflag_top_bound=0 for no sponge
    52 ## iflag_top_bound=1 for sponge over 4 topmost layers
    53 ## iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
    54 iflag_top_bound=0
    55 iflag_top_bound=2
    56 iflag_top_bound=1
    57 ## Mode Couche eponge
    58 ##   mode = 0 : pas de sponge
    59 ##   mode = 1 : u et v -> 0
    60 ##   mode = 2 : u et v -> moyenne zonale
    61 ##   mode = 3 : u, v et h -> moyenne zonale
    62 mode_top_bound=0
    63 mode_top_bound=2
    64 mode_top_bound=1
    65 #######################################
    66 
    67 
    68 #  Coefficient pour la couche eponge (valeur derniere couche)
    69 tau_top_bound=1.e-5
    70 #tau_top_bound=5.e-5
    71 #tau_top_bound=1.e-4
    72 
    73 ## guidage par vent dans relax.nc
    74 ok_guide = y
    75 ok_guide = n
    76 facwind = 2.e7
    77 facwind = 5.e8
    78 
    79 ################################################################
    80 #### LMDZ5
    81 #### IN RUN.DEF
    82 #### -- those parameters are ignored when using old LMDZ core
    83 ################################################################
    84 
    85 ## Type de calendrier utilise
    86 ## valeur possible: earth_360d (defaut), earth_365d, earth_366d
    87 ip_ebil_phy=1
    88 ip_ebil_dyn=1
    89 calend=earth_360d
    90 
    91 ### Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
    92 #dayref=10120
    93 ###  Annee de l'etat  initial (   avec  4  chiffres   )
    94 #anneeref=2000
    95 
    96 ## Remise a zero de la date initiale
    97 raz_date=0
    98 ### for restart runs
    99 # use 1 for 'from scratch' runs
    100 
    101 ## flag de sortie dynzon
    102 ok_dynzon=n
    103 ## activation du calcul d equilibrage de charge
    104 adjust=n
    105 ## activation du filtre fft
    106 use_filtre_fft=n
    107 ## niveau d'impression de controle
    108 prt_level=10
    109 prt_level=1
    110 
    111 ## coordonnees type planeto
    112 disvert_type = 2
    113 
    114 ## the generic physics
    115 planet_type = generic
    116 
    117 ## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.)
    118 ok_dyn_ins=n
    119 #ok_dyn_ins=y
    120 ## periode d'ecriture des sorties instantanees dans la dynamique
    121 ## (en pas dynamiques)
    122 #iecri=30
    123 #iecri=120
    124 iecri=480
    125 #iecri=960
    126 #iecri=4800
    127 #iecri=9600
    128 
    129 ## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.)
    130 ok_dyn_ave=n
    131 ## periode de stockage des moyennes dans la dynamique (en jours)
    132 periodav=1
    133 
    134 
    135 #################################
    136 #### AND NOW THE OLD run.def ####
    137 #################################
    138 ## -----> REMOVED
    139 ## was present and conflicting with LMDZ5
    140 ## (although useless in old OLD LMDz DYN CORE
    141 ## - iecri [periode d'ecriture du fichier histoire (en jour)]
    142 ## - periodav [periode de stockage fichier histmoy (en jour) dynzon?]
    143 #################################
    144 #################################
    145 
    146 #-----------------------------------------------------------------------
    147 #Parametres de controle du run:                                         
    148 #------------------------------                                         
    149 
    150 # Nombre de jours d'integration                                         
    151 #     nday = 720
    152 #     nday = 20
    153 #     nday = 100
    154 #     nday = 1800
    155 #     nday = 250000
    156      nday = 65000
    157 
    158 # nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )     
    159 # day_step = 25
    160 # day_step = 50
    161 day_step = 70
    162 # day_step = 240
    163 # day_step = 480
    164 # day_step = 960
    165 
    166 # periode pour le pas Matsuno (en pas)                                 
    167   iperiod=5
    168 
    169 # periode de sortie des variables de controle (en pas)                 
    170   iconser=120
    171 
    172 # periode de la dissipation (en pas)                                   
    173 # -- this is ignored by LMDZ5 which rely on dissip_period instead
    174   idissip=5
    175 
    176 # choix de l'operateur de dissipation (star ou  non star )             
    177  lstardis=.true.
    178 
    179 # avec ou sans coordonnee hybrides
    180  hybrid=.false. # calculs 1D sandrine #.true.
    181 
    182 # nombre d'iterations de l'operateur de dissipation   gradiv           
    183 nitergdiv=1
    184 
    185 # nombre d'iterations de l'operateur de dissipation  nxgradrot         
    186 nitergrot=2
    187 
    188 # nombre d'iterations de l'operateur de dissipation  divgrad           
    189    niterh=2
    190 
    191 # temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    192  tetagdiv= 5400. 8000. 160000. 1000. 2500. 4000.
    193 
    194 # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    195  tetagrot= 7200. 8000. 160000. 1000. 2500. 4000.
    196 
    197 # temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
    198  tetatemp= 7200. 8000. 160000. 1000. 2500. 4000.
    199 
    200 # coefficient pour gamdissip                                           
    201   coefdis=0.
    202 
    203 # choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
    204   purmats=.false.
    205 
    206 # avec ou sans physique                                                 
    207    physic=.true.
    208 
    209 # periode de la physique (en pas) (en pas dynamiques, n'a de sens que si iflag_phys=1)                                     
    210   iphysiq=5
    211   iphysiq=14
    212 #  iphysiq=1
    213 
    214 # choix d'une grille reguliere                                         
    215   grireg=.true.
    216 
    217 # frequence (en pas) de l'ecriture du fichier diagfi               
    218   ecritphy=500
    219   ecritphy=7000
    220 
    221 # longitude en degres du centre du zoom                                 
    222    clon=0.
    223 
    224 # latitude en degres du centre du zoom                                 
    225    clat=0.
    226 
    227 # facteur de grossissement du zoom,selon longitude                     
    228   grossismx=1.
    229 
    230 # facteur de grossissement du zoom ,selon latitude                     
    231  grossismy=1.
    232 
    233 #  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
    234   fxyhypb=.false.
    235 
    236 # extension en longitude  de la zone du zoom  ( fraction de la zone totale)
    237    dzoomx= 0.
    238 
    239 # extension en latitude de la zone  du zoom  ( fraction de la zone totale)
    240    dzoomy=0.
    241 
    242 #  raideur du zoom en  X
    243     taux=2.
    244 
    245 #  raideur du zoom en  Y
    246     tauy=2.
    247 
    248 #  Fonction  f(y) avec y = Sin(latit.) si = .TRUE. ,  Sinon  y = latit.
    249   ysinus= .false.
    250 
    251 ###################### SPONGE LAYER IN OLD DYNAMICS
    252 ###################### -- this is ignored by LMDZ5 !!! NO LONGER !!!
     170###########################################################################
     171### LEFTOVERS
     172### calend=earth_360d ## specific calendar
     173### resetvarc=n   ## reinit des variables de controle
     174### iflag_trac=0  ## avec ou sans traceurs
     175### dayref=10120  ## initial date
     176### anneeref=2000 ## year of initial state
     177### --- old version of spong layer
    253178## Avec sponge layer
    254179#  callsponge  = .true.
     
    259184## Sponge:  tetasponge (secondes)
    260185#  tetasponge = 50000
    261 #######################
    262186
    263 # some definitions for the physics, in file 'callphys.def'
     187###########################################################################
     188### INCLUDE OTHER DEF FILES (physics, specific settings, etc...)
     189###########################################################################
    264190INCLUDEDEF=callphys.def
    265 #INCLUDEDEF=modif.def
     191INCLUDEDEF=modif.def
Note: See TracChangeset for help on using the changeset viewer.