Ignore:
Timestamp:
Oct 17, 2023, 9:40:14 AM (15 months ago)
Author:
slebonnois
Message:

BdeBatz? : Cleans microphysics and makes few corrections for physics

Location:
trunk/LMDZ.TITAN/deftank
Files:
1 deleted
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.TITAN/deftank/callphys.def

    r3085 r3090  
    9090# Chemistry is computed every "ichim" phys. timestep
    9191ichim      = 20
    92 # Infinite tank of CH4 ?
     92# Infinite tank of CH4 ? (.false. == lakes at the poles)
    9393resCH4_inf = .false.
    9494
    … …  
    107107opt4clouds = .true.
    108108# Cloudy fraction in the Clear / Dark column method ?
    109 Fcloudy    = 0.15
     109Fcloudy    = 0.2
    110110# Adjustment of ssca albedo for cloudy dropplet ?
    111111Fssa       = 0.90

  • trunk/LMDZ.TITAN/deftank/diagfi.def

    r1681 r3090  
    1 ISR
    2 ASR
    3 pphi
    4 Ls
    5 tsurf
    6 ps
    7 temp
    8 u
    9 v
    10 w
    11 p
    12 q2
    13 sensibFlux
    14 zdtsw
    15 zdtlw
    16 zdtdyn
    17 zdudif
    18 zdudyn

  • trunk/LMDZ.TITAN/deftank/gases.def

    r1681 r3090  
    1 # gases
     1# Gases
    2210
    330.985

  • trunk/LMDZ.TITAN/deftank/profile.def

    r2241 r3090  
    1 #profile (Reference CH4 vertical profile (mol/mol) : surface->top)
     1# Profile (Reference CH4 vertical profile (mol/mol) : Surface ---> Top)
    220.0565
    330.0567

  • trunk/LMDZ.TITAN/deftank/run.def.11

    r1681 r3090  
    1 #
    2 #-----------------------------------------------------------------------
    3 #Parametres de controle du run:                                         
    4 #------------------------------                                         
     1# -----------------------------------------------------------------------
     2# Parametres de controle du run :
     3# -----------------------------------------------------------------------
    54
    65# Planete
    … …  
    87
    98# Nombre de jours d'integration                                         
    10      nday=11
     9nday = 11
    1110
    12 # nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )     
    13  day_step = 8000
     11# Nombre de pas par jour (multiple de iperiod) (ici pour dt = 1 min)
     12day_step = 8000
    1413
    15 # periode pour le pas Matsuno (en pas)                                 
    16   iperiod=5
     14# Periode pour le pas Matsuno (en pas)                       
     15iperiod = 5
    1716
    18 # periode de sortie des variables de controle (en pas)                 
    19   iconser=500
     17# Periode de sortie des variables de controle (en pas)
     18iconser = 1000
    2019
    21 # periode d'ecriture du fichier histoire (en jour)                     
    22     iecri=100
     20# Periode d'ecriture du fichier histoire (en jour)
     21iecri = 100
    2322
    24 # periode de stockage fichier histmoy (en jour)                         
    25  periodav=60.
     23# Periode de stockage fichier histmoy (en jour)
     24periodav = 60.
    2625
    27 # periode de la dissipation (en pas)                                   
    28   dissip_period=5
     26# Periode de la dissipation (en pas)
     27dissip_period = 5
    2928
    30 # choix de l'operateur de dissipation (star ou  non star )             
    31  lstardis=.true.
     29# Choix de l'operateur de dissipation (star ou non star)
     30lstardis = .true.
    3231
    33 # avec ou sans coordonnee hybrides
    34  hybrid=.true.
     32# Avec ou sans coordonnee hybrides
     33hybrid = .true.
    3534
    36 # nombre d'iterations de l'operateur de dissipation   gradiv           
    37 nitergdiv=1
     35# Nombre d'iterations de l'operateur de dissipation gradiv
     36nitergdiv = 1
    3837
    39 # nombre d'iterations de l'operateur de dissipation  nxgradrot         
    40 nitergrot=2
     38# Nombre d'iterations de l'operateur de dissipation nxgradrot
     39nitergrot = 2
    4140
    42 # nombre d'iterations de l'operateur de dissipation  divgrad           
    43    niterh=2
     41# Nombre d'iterations de l'operateur de dissipation divgrad
     42niterh = 2
    4443
    45 # temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    46  tetagdiv=50000.
     44# Temps de dissipation des plus petites long. d'ondes pour u,v (gradiv)
     45tetagdiv = 25000.
    4746
    48 # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    49  tetagrot=50000.
     47# Temps de dissipation des plus petites long. d'ondes pour u,v (nxgradrot)
     48tetagrot = 25000.
    5049
    51 # temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
    52  tetatemp=50000.
     50# Temps de dissipation des plus petites long. d'ondes pour h (divgrad)
     51tetatemp = 25000.
    5352
    5453# Facteur d'amplification de la dissip en haut du modele
    55   dissip_fac_up = 30.
     54dissip_fac_up = 30.
    5655
    5756# Intervalle (km) pour le changement mid / up de la dissip
    58 #   dissip_deltaz = 30.
     57# dissip_deltaz = 30.
    5958
    6059# Pression (Pa) au bas la transition mid / up
    61    dissip_pupstart = 1.e1
     60dissip_pupstart = 1.e1
    6261
    63 # coefficient pour gamdissip                                           
    64   coefdis=0.
     62# Coefficient pour gamdissip
     63coefdis = 0.
    6564
    66 # choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
    67   purmats=.false.
     65# Choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
     66purmats = .false.
    6867
    69 # avec ou sans physique                                                 
    70    physic=.true.
     68# Avec ou sans physique
     69physic = .true.
    7170
    72 # periode de la physique (en pas)                                       
    73   iphysiq=40
     71# Periode de la physique (en pas)
     72iphysiq = 40
    7473
    75 # Titan: marees gravitationnelles
    76   tidal=.true.
     74# Titan : marees gravitationnelles
     75tidal = .true.
    7776
    78 # choix d'une grille reguliere                                         
    79   grireg=.true.
     77# Choix d'une grille reguliere
     78grireg = .true.
    8079
    81 # frequence (en pas) de l'ecriture du fichier diagfi               
    82  ecritphy=2000
     80# Frequence (en pas) de l'ecriture du fichier diagfi
     81ecritphy = 8000
    8382
    84 # longitude en degres du centre du zoom                                 
    85    clon=0.
     83# Longitude en degres du centre du zoom
     84clon = 0.
    8685
    87 # latitude en degres du centre du zoom                                 
    88    clat=0.
     86# Latitude en degres du centre du zoom
     87clat = 0.
    8988
    90 # facteur de grossissement du zoom,selon longitude                     
    91   grossismx=1.
     89# Facteur de grossissement du zoom, selon longitude
     90grossismx = 1.
    9291
    93 # facteur de grossissement du zoom ,selon latitude                     
    94  grossismy=1.
     92# Facteur de grossissement du zoom, selon latitude
     93grossismy = 1.
    9594
    96 #  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
    97   fxyhypb=.false.
     95# Fonction f(y) hyperbolique si = .true., sinon sinusoidale
     96fxyhypb = .false.
    9897
    99 # extension en longitude  de la zone du zoom  ( fraction de la zone totale)
    100    dzoomx= 0.
     98# Extension en longitude de la zone du zoom (fraction de la zone totale)
     99dzoomx = 0.
    101100
    102 # extension en latitude de la zone  du zoom  ( fraction de la zone totale)
    103    dzoomy=0.
     101# Extension en latitude de la zone du zoom (fraction de la zone totale)
     102dzoomy = 0.
    104103
    105 #  raideur du zoom en X
    106     taux=2.
     104# Raideur du zoom en X
     105taux = 2.
    107106
    108 #  raideur du zoom en Y
    109     tauy=2.
     107# Raideur du zoom en Y
     108tauy = 2.
    110109
    111 #  Fonction  f(y) avec y = Sin(latit.) si = .TRUE. ,  Sinon y = latit.
    112   ysinus= .false.
     110# Fonction f(y) avec y = Sin(latit.), si = .true., sinon y = latit.
     111ysinus = .false.
    113112
    114113# Avec sponge layer
    115   callsponge = .true.
     114callsponge = .true.
    116115
    117 # Sponge: mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
    118   mode_sponge= 2
     116# Sponge : mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
     117mode_sponge = 2
    119118
    120119# Sponge : nsponge
    121   nsponge = 5
     120nsponge = 4
    122121
    123 # Sponge: tetasponge (secondes)
    124   tetasponge = 10000
     122# Sponge : tetasponge (secondes)
     123tetasponge = 10000
    125124
    126 # some definitions for the physics, in file 'callphys.def'
     125# Some definitions for the physics, in file 'callphys.def'
    127126INCLUDEDEF=callphys.def

  • trunk/LMDZ.TITAN/deftank/run.def.12

    r1681 r3090  
    1 #
    2 #-----------------------------------------------------------------------
    3 #Parametres de controle du run:                                         
    4 #------------------------------                                         
     1# -----------------------------------------------------------------------
     2# Parametres de controle du run :
     3# -----------------------------------------------------------------------
    54
    65# Planete
    … …  
    87
    98# Nombre de jours d'integration                                         
    10      nday=12
     9nday = 12
    1110
    12 # nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )     
    13  day_step = 8000
     11# Nombre de pas par jour (multiple de iperiod) (ici pour dt = 1 min)
     12day_step = 8000
    1413
    15 # periode pour le pas Matsuno (en pas)                                 
    16   iperiod=5
     14# Periode pour le pas Matsuno (en pas)                       
     15iperiod = 5
    1716
    18 # periode de sortie des variables de controle (en pas)                 
    19   iconser=500
     17# Periode de sortie des variables de controle (en pas)
     18iconser = 1000
    2019
    21 # periode d'ecriture du fichier histoire (en jour)                     
    22     iecri=100
     20# Periode d'ecriture du fichier histoire (en jour)
     21iecri = 100
    2322
    24 # periode de stockage fichier histmoy (en jour)                         
    25  periodav=60.
     23# Periode de stockage fichier histmoy (en jour)
     24periodav = 60.
    2625
    27 # periode de la dissipation (en pas)                                   
    28   dissip_period=5
     26# Periode de la dissipation (en pas)
     27dissip_period = 5
    2928
    30 # choix de l'operateur de dissipation (star ou  non star )             
    31  lstardis=.true.
     29# Choix de l'operateur de dissipation (star ou non star)
     30lstardis = .true.
    3231
    33 # avec ou sans coordonnee hybrides
    34  hybrid=.true.
     32# Avec ou sans coordonnee hybrides
     33hybrid = .true.
    3534
    36 # nombre d'iterations de l'operateur de dissipation   gradiv           
    37 nitergdiv=1
     35# Nombre d'iterations de l'operateur de dissipation gradiv
     36nitergdiv = 1
    3837
    39 # nombre d'iterations de l'operateur de dissipation  nxgradrot         
    40 nitergrot=2
     38# Nombre d'iterations de l'operateur de dissipation nxgradrot
     39nitergrot = 2
    4140
    42 # nombre d'iterations de l'operateur de dissipation  divgrad           
    43    niterh=2
     41# Nombre d'iterations de l'operateur de dissipation divgrad
     42niterh = 2
    4443
    45 # temps de dissipation des plus petites long.d ondes pour u,v (gradiv) 
    46  tetagdiv=50000.
     44# Temps de dissipation des plus petites long. d'ondes pour u,v (gradiv)
     45tetagdiv = 25000.
    4746
    48 # temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
    49  tetagrot=50000.
     47# Temps de dissipation des plus petites long. d'ondes pour u,v (nxgradrot)
     48tetagrot = 25000.
    5049
    51 # temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
    52  tetatemp=50000.
     50# Temps de dissipation des plus petites long. d'ondes pour h (divgrad)
     51tetatemp = 25000.
    5352
    5453# Facteur d'amplification de la dissip en haut du modele
    55   dissip_fac_up = 30.
     54dissip_fac_up = 30.
    5655
    5756# Intervalle (km) pour le changement mid / up de la dissip
    58 #   dissip_deltaz = 50.
     57# dissip_deltaz = 30.
    5958
    6059# Pression (Pa) au bas la transition mid / up
    61    dissip_pupstart = 1.e1
     60dissip_pupstart = 1.e1
    6261
    63 # coefficient pour gamdissip                                           
    64   coefdis=0.
     62# Coefficient pour gamdissip
     63coefdis = 0.
    6564
    66 # choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
    67   purmats=.false.
     65# Choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
     66purmats = .false.
    6867
    69 # avec ou sans physique                                                 
    70    physic=.true.
     68# Avec ou sans physique
     69physic = .true.
    7170
    72 # periode de la physique (en pas)                                       
    73   iphysiq=40
     71# Periode de la physique (en pas)
     72iphysiq = 40
    7473
    75 # Titan: marees gravitationnelles
    76   tidal=.true.
     74# Titan : marees gravitationnelles
     75tidal = .true.
    7776
    78 # choix d'une grille reguliere                                         
    79   grireg=.true.
     77# Choix d'une grille reguliere
     78grireg = .true.
    8079
    81 # frequence (en pas) de l'ecriture du fichier diagfi               
    82  ecritphy=2000
     80# Frequence (en pas) de l'ecriture du fichier diagfi
     81ecritphy = 8000
    8382
    84 # longitude en degres du centre du zoom                                 
    85    clon=0.
     83# Longitude en degres du centre du zoom
     84clon = 0.
    8685
    87 # latitude en degres du centre du zoom                                 
    88    clat=0.
     86# Latitude en degres du centre du zoom
     87clat = 0.
    8988
    90 # facteur de grossissement du zoom,selon longitude                     
    91   grossismx=1.
     89# Facteur de grossissement du zoom, selon longitude
     90grossismx = 1.
    9291
    93 # facteur de grossissement du zoom ,selon latitude                     
    94  grossismy=1.
     92# Facteur de grossissement du zoom, selon latitude
     93grossismy = 1.
    9594
    96 #  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
    97   fxyhypb=.false.
     95# Fonction f(y) hyperbolique si = .true., sinon sinusoidale
     96fxyhypb = .false.
    9897
    99 # extension en longitude  de la zone du zoom  ( fraction de la zone totale)
    100    dzoomx= 0.
     98# Extension en longitude de la zone du zoom (fraction de la zone totale)
     99dzoomx = 0.
    101100
    102 # extension en latitude de la zone  du zoom  ( fraction de la zone totale)
    103    dzoomy=0.
     101# Extension en latitude de la zone du zoom (fraction de la zone totale)
     102dzoomy = 0.
    104103
    105 #  raideur du zoom en X
    106     taux=2.
     104# Raideur du zoom en X
     105taux = 2.
    107106
    108 #  raideur du zoom en Y
    109     tauy=2.
     107# Raideur du zoom en Y
     108tauy = 2.
    110109
    111 #  Fonction  f(y) avec y = Sin(latit.) si = .TRUE. ,  Sinon y = latit.
    112   ysinus= .false.
     110# Fonction f(y) avec y = Sin(latit.), si = .true., sinon y = latit.
     111ysinus = .false.
    113112
    114113# Avec sponge layer
    115   callsponge = .true.
     114callsponge = .true.
    116115
    117 # Sponge: mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
    118   mode_sponge= 2
     116# Sponge : mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
     117mode_sponge = 2
    119118
    120119# Sponge : nsponge
    121   nsponge = 5
     120nsponge = 4
    122121
    123 # Sponge: tetasponge (secondes)
    124   tetasponge = 10000
     122# Sponge : tetasponge (secondes)
     123tetasponge = 10000
    125124
    126 # some definitions for the physics, in file 'callphys.def'
     125# Some definitions for the physics, in file 'callphys.def'
    127126INCLUDEDEF=callphys.def

  • trunk/LMDZ.TITAN/deftank/z2sig.def

    r1681 r3090  
    1  40.00000      H: atmospheric scale height (km) (used as a reference only)
    2   0.040000000     Typical log(sigma) for 1st layer (z=H*log(sigma))
    3   0.240000000     Typical log(sigma) for 2nd layer, etc...
    4   0.640000000
    5   1.600000000
     1 40.00000       H: Atmospheric scale height [km] --> used as a reference only
     2  0.040000000   Typical log(sigma) for 1st layer (z = H*log(sigma)),
     3  0.240000000   Typical log(sigma) for 2nd layer,
     4  0.640000000   Typical log(sigma) for 3rd layer,
     5  1.600000000   etc...
    66  3.600000000
    77  7.200000000
Note: See TracChangeset for help on using the changeset viewer.