Changeset 844


Ignore:
Timestamp:
Nov 18, 2012, 2:03:20 PM (13 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:
158 added
1 edited
8 copied

Legend:

Unmodified
Added
Removed

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

    r838 r844  
    11
     2#########
     3######### THIS IS INTENDED TO WORK BOTH WITH LMDZ.GENERIC + old LMDZ core
     4#########                                AND LMDZ.GENERIC + new LMDZ core
     5#########
     6
     7### A INCLURE ?
     8## calend=venus  ## calendrier specifique
     9## resetvarc=n   ## reinit des variables de controle
     10## iflag_trac=0  ## avec ou sans traceurs
     11
     12################################################################
    213#### LMDZ5
    314#### IN GCM.DEF
     15#### -- those parameters are ignored when using old LMDZ core
     16################################################################
    417
    518## period appel a la dissipation (multiple iperiod, 0 si calcul automatique)
     19## (remplace ce qui etait idissip dans l'ancienne dynamique)
    620dissip_period = 0
    721dissip_period = 2
    … …  
    2135read_start=y
    2236
    23 ##  Avec ou sans strato
     37##  Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horizontale
    2438ok_strato=n
     39
     40### Dissipation horizontale
     41#dissip_fac_mid=2.
     42#dissip_fac_up=2.
     43## deltaz et hdelta en km
     44#dissip_deltaz=10.
     45#dissip_hdelta=5.
     46## pupstart en Pa
     47#dissip_pupstart=1.e3
    2548
    2649#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
    2750iflag_top_bound=0
     51
     52########################################
     53## LMDZ.COMMON
     54##    1: dans les 4 derniers niveaux
     55##    2: dans les couches de pression plus faible que 100 fois la pression de la derniere couche
     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=1
     63#######################################
     64
    2865#  Coefficient pour la couche eponge (valeur derniere couche)
     66tau_top_bound=1.e-5
    2967tau_top_bound=5.e-5
    3068
     69################################################################
    3170#### LMDZ5
    3271#### IN RUN.DEF
     72#### -- those parameters are ignored when using old LMDZ core
     73################################################################
    3374
    3475## Type de calendrier utilise
    … …  
    4081## Jour de l'etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
    4182dayref=10120
    42 dayref=40
     83#dayref=40
    4384
    4485##  Annee de l'etat  initial (   avec  4  chiffres   )
    … …  
    72113iecri=120
    73114iecri=480
    74 iecri=4800
     115#iecri=4800
     116#iecri=9600
    75117
    76118## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.)
    … …  
    83125#### AND NOW THE OLD run.def ####
    84126#################################
    85 
     127## -----> REMOVED
     128## was present and conflicting with LMDZ5
     129## (although useless in old OLD LMDz DYN CORE
     130## - iecri [periode d'ecriture du fichier histoire (en jour)]
     131## - periodav [periode de stockage fichier histmoy (en jour) dynzon?]
     132#################################
     133#################################
    86134
    87135#-----------------------------------------------------------------------
    … …  
    91139# Nombre de jours d'integration                                         
    92140     nday = 720
    93      nday = 20
    94      nday = 100
    95      nday = 1800
     141#     nday = 20
     142#     nday = 100
     143#     nday = 1800
    96144
    97145# nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )     
    98146 day_step = 240
    99  day_step = 480
     147# day_step = 480
     148# day_step = 960
    100149
    101150# periode pour le pas Matsuno (en pas)                                 
    … …  
    105154  iconser=120
    106155
    107 ############# VOIR CI-DESSUS ne sert pas dans generic
    108 ############# ... mais sert dans LMDZ5 !!!!!!
    109 ############# ... donc on commente ...
    110 ## periode d'ecriture du fichier histoire (en jour)                     
    111 #    iecri=1
    112 #
    113 ## periode de stockage fichier histmoy (en jour) dynzon?                       
    114 # periodav=10.
    115 ############# VOIR CI-DESSUS
    116 ############# VOIR CI-DESSUS
    117 
    118 ############# VOIR CI-DESSUS dissip_period
    119 ## periode de la dissipation (en pas)                                   
    120 #  idissip=5
    121 ############# VOIR CI-DESSUS
     156# periode de la dissipation (en pas)                                   
     157# -- this is ignored by LMDZ5 which rely on dissip_period instead
     158  idissip=5
    122159
    123160# choix de l'operateur de dissipation (star ou  non star )             
    124  lstardis=.true. #ou y ???
     161 lstardis=.true.
    125162
    126163# avec ou sans coordonnee hybrides
    … …  
    137174
    138175# temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    139  tetagdiv=1000. 2500. 4000.
     176 tetagdiv=2500. 4000.
    140177
    141178# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    142  tetagrot=1000. 2500. 4000.
     179 tetagrot=2500. 4000.
    143180
    144181# temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
    145  tetatemp=1000. 2500. 4000.
     182 tetatemp=2500. 4000.
    146183
    147184# coefficient pour gamdissip                                           
    … …  
    149186
    150187# choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
    151   purmats=.false. # ou n ???
     188  purmats=.false.
    152189
    153190# avec ou sans physique                                                 
    … …  
    197234  ysinus= .false.
    198235
    199 ###################### MARCHE SEULEMENT AVEC L ANCIENNE DYNAMIQUE
    200 ## Avec sponge layer
    201 #  callsponge  = .true.
    202 #
    203 ## Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
    204 #  mode_sponge= 2
    205 #
    206 ## Sponge:  hauteur de sponge (km)
    207 #  hsponge= 50
    208 #
    209 ## Sponge:  tetasponge (secondes)
    210 #  tetasponge = 50000
    211 ###################### MARCHE SEULEMENT AVEC L ANCIENNE DYNAMIQUE
     236###################### SPONGE LAYER IN OLD DYNAMICS
     237###################### -- this is ignored by LMDZ5
     238# Avec sponge layer
     239  callsponge  = .true.
     240# Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
     241  mode_sponge= 2
     242# Sponge:  hauteur de sponge (km)
     243  hsponge= 50
     244# Sponge:  tetasponge (secondes)
     245  tetasponge = 50000
     246######################
    212247
    213248# some definitions for the physics, in file 'callphys.def'

  • trunk/LMDZ.GENERIC_2.0/libf/phygeneric/iniphysiq.F

    r843 r844  
    88     $           pdayref,ptimestep,
    99     $           plat,plon,parea,pcu,pcv,
    10      $           prad,pg,pr,pcpp)
     10     $           prad,pg,pr,pcpp,iflag_phys)
    1111      USE dimphy
    1212      USE mod_grid_phy_lmdz
    … …  
    7777      INTEGER pdayref
    7878      INTEGER :: ibegin,iend,offset
     79      INTEGER,INTENT(IN) :: iflag_phys ! type of physics to be called
    7980
    8081      INTEGER :: ngrid_sub
    … …  
    187188      !!print*, plat,plon,parea,prad,pg,pr,pcpp
    188189
    189       RETURN
    190 9999  CONTINUE
    191       abort_message ='Cette version demande les fichier rnatur.dat
    192      & et surf.def'
    193       CALL abort_gcm (modname,abort_message,1)
     190!      RETURN
     191!9999  CONTINUE
     192!      abort_message ='Cette version demande les fichier rnatur.dat
     193!     & et surf.def'
     194!      CALL abort_gcm (modname,abort_message,1)
    194195
    195196      END
Note: See TracChangeset for help on using the changeset viewer.