Ignore:
Timestamp:
Mar 22, 2011, 5:25:44 PM (14 years ago)
Author:
slebonnois
Message:

Serie de modifs SL pour homogeneisation des phytitan et phyvenus
Ca touche aussi aux liens phy/dyn (surtout a propos de clesphy0),
a verifier avec les autres, donc...

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/libf/phytitan/conf_phys.F90

    r3 r97  
    1010   implicit none
    1111
    12 #include "YOMCST.inc"
    13 #include "clesphys.inc"
     12#include "YOMCST.h"
     13#include "clesphys.h"
    1414#include "compbl.h"
    15 #include "comorbit.inc"
     15#include "comorbit.h"
     16
     17! ok_journe:  sorties journalieres
     18! ok_mensuel: sorties mensuelles
     19! ok_instan:  sorties instantanees
     20
     21
     22! Sortie:
     23  logical              :: ok_journe, ok_mensuel, ok_instan       
     24  integer              :: if_ebil
     25
     26! Local
     27  integer              :: numout = 6
     28
    1629!
    1730! Configuration de la "physique" de LMDZ a l'aide de la fonction
    … …  
    2235!--- Ca lit le physiq.def ---
    2336
    24 !
    25 ! ok_mensuel: sorties mensuelles
    26 ! ok_journe:  sorties journalieres
    27 ! ok_instan:  sorties instantanees
    28 
    29 
    30 ! Sortie:
    31   logical              :: ok_mensuel, ok_journe, ok_instan
    32   integer              :: if_ebil
    33 
    34 ! Local
    35   integer              :: numout = 6
    36   real*4               :: zzz
     37!******************* parametres anciennement lus dans gcm.def
     38
     39!Config  Key  = cycle_diurne
     40!Config  Desc = Cycle ddiurne
     41!Config  Def  = y
     42!Config  Help = Cette option permet d'eteidre le cycle diurne.
     43!Config         Peut etre util pour accelerer le code !
     44       cycle_diurne = .TRUE.
     45       CALL getin('cycle_diurne',cycle_diurne)
     46
     47!Config  Key  = soil_model
     48!Config  Desc = Modele de sol
     49!Config  Def  = y
     50!Config  Help = Choix du modele de sol (Thermique ?)
     51!Config         Option qui pourait un string afin de pouvoir
     52!Config         plus de choix ! Ou meme une liste d'options !
     53       soil_model = .TRUE.
     54       CALL getin('soil_model',soil_model)
     55
     56!Config  Key  = ok_orodr
     57!Config  Desc = Oro drag
     58!Config  Def  = y
     59!Config  Help = GW drag orographie
     60!Config         
     61       ok_orodr = .TRUE.
     62       CALL getin('ok_orodr',ok_orodr)
     63
     64!Config  Key  =  ok_orolf
     65!Config  Desc = Oro lift
     66!Config  Def  = n
     67!Config  Help = GW lift orographie (pas utilise)
     68       ok_orolf = .TRUE.
     69       CALL getin('ok_orolf', ok_orolf)
     70
     71!Config  Key  = ok_gw_nonoro
     72!Config  Desc = Gravity waves parameterization
     73!Config  Def  = n
     74!Config  Help = GW drag non-orographique
     75       ok_gw_nonoro = .FALSE.
     76       CALL getin('ok_gw_nonoro',ok_gw_nonoro)
     77
     78!Config  Key  = nbapp_rad
     79!Config  Desc = Frequence d'appel au rayonnement
     80!Config  Def  = 12
     81!Config  Help = Nombre  d'appels des routines de rayonnements
     82!Config         par jour.
     83       nbapp_rad = 12
     84       CALL getin('nbapp_rad',nbapp_rad)
     85
     86!Config  Key  = nbapp_chim
     87!Config  Desc = Frequence d'appel a la chimie
     88!Config  Def  = 1
     89!Config  Help = Nombre  d'appels des routines de chimie
     90!Config         par jour.
     91       nbapp_chim = 1
     92       CALL getin('nbapp_chim',nbapp_chim)
     93
     94!Config  Key  = iflag_con
     95!Config  Desc = Flag de convection
     96!Config  Def  = 0
     97!Config  Help = Flag  pour la convection les options suivantes existent :
     98!Config         0 : ajsec simple (VENUS, TITAN)
     99!Config         1 pour LMD,
     100!Config         2 pour Tiedtke,
     101!Config         3 pour CCM(NCAR) 
     102       iflag_con = 0
     103       CALL getin('iflag_con',iflag_con)
     104
     105!******************* fin parametres anciennement lus dans gcm.def
    37106
    38107!Config Key  = OK_mensuel
    … …  
    64133!
    65134!
     135!Config  Key  = ecritphy
     136!Config  Desc = Frequence d'ecriture dans histins
     137!Config  Def  = 1
     138!Config  Help = frequence de l'ecriture du fichier histins
     139!Config         en jours.
     140!
     141       ecritphy = 1.
     142       CALL getin('ecritphy', ecritphy)
     143!
    66144!Config Key  = if_ebil
    67145!Config Desc = Niveau de sortie pour les diags bilan d'energie
    … …  
    85163!Config Help =
    86164!               
    87   zzz = 673.
    88   call getin('year_day', zzz)
    89   year_day = zzz
     165  year_day = 673.
     166  call getin('year_day', year_day)
    90167!
    91168!Config Key  = peri_day
    … …  
    94171!Config Help =
    95172!               
    96   zzz = 533.
    97   call getin('peri_day', zzz)
    98   peri_day = zzz
     173  peri_day = 533.
     174  call getin('peri_day', peri_day)
    99175!
    100176!Config Key  = periheli
    … …  
    103179!Config Help =
    104180!               
    105   zzz = 1354.5
    106   call getin('periheli', zzz)
    107   periheli = zzz
     181  periheli = 1354.5
     182  call getin('periheli', periheli)
    108183!!
    109184!Config Key  = aphelie
    … …  
    112187!Config Help =
    113188!               
    114   zzz = 1506.0
    115   call getin('aphelie', zzz)
    116   aphelie = zzz
     189  aphelie = 1506.0
     190  call getin('aphelie', aphelie)
    117191!!
    118192!Config Key  = obliquit
    … …  
    121195!Config Help =
    122196!               
    123   zzz = 26.7
    124   call getin('obliquit', zzz)
    125   obliquit = zzz
     197  obliquit = 26.7
     198  call getin('obliquit', obliquit)
    126199!
    127200!Config Key  = solaire
    … …  
    131204!Config Help =
    132205!
    133   zzz = 2620.
    134     call getin('solaire', zzz)
    135   solaire = zzz
     206  solaire = 2620.
     207    call getin('solaire', solaire)
    136208!
    137209!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    158230!Config Help =
    159231!
    160   zzz = 1.0e-7
    161   call getin('ksta',zzz)
    162   ksta = zzz
     232  ksta = 1.0e-7
     233  call getin('ksta',ksta)
    163234
    164235!
    … …  
    168239!Config Help =
    169240!
    170   zzz = 1.0e-2
    171   call getin('z0',zzz)
    172   z0 = zzz
     241  z0 = 1.0e-2
     242  call getin('z0',z0)
    173243
    174244!
    … …  
    178248!Config Help =
    179249!
    180   zzz = 35.
    181   call getin('lmixmin',zzz)
    182   lmixmin = zzz
    183 
    184 !
     250  lmixmin = 35.
     251  call getin('lmixmin',lmixmin)
     252
    185253!
    186254!Config Key  = ok_kzmin
    … …  
    204272!Config Key  = inertie
    205273!Config Desc =
    206 !Config Def  = 2000.
    207 !Config Help =
    208 !
    209   zzz = 2000.
    210   call getin('inertie',zzz)
    211   inertie = zzz
     274!Config Def  = 340.
     275!Config Help =
     276!
     277  inertie = 340.
     278  call getin('inertie',inertie)
    212279!
    213280!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    287354!Config Help =
    288355!
    289   zzz = 3.5
    290   call getin('tx',zzz)
    291   tx = zzz
     356  tx = 3.5
     357  call getin('tx',tx)
    292358
    293359!
    … …  
    297363!Config Help =
    298364!
    299   zzz = 1.
    300   call getin('tcorrect',zzz)
    301   tcorrect = zzz
     365  tcorrect = 1.
     366  call getin('tcorrect',tcorrect)
    302367
    303368!
    … …  
    349414  write(numout,*)' ##############################################'
    350415  write(numout,*)' Configuration des parametres de la physique: '
     416  write(numout,*)' cycle_diurne = ', cycle_diurne
     417  write(numout,*)' soil_model = ', soil_model
     418  write(numout,*)' ok_orodr = ', ok_orodr
     419  write(numout,*)' ok_orolf = ', ok_orolf
     420  write(numout,*)' ok_gw_nonoro = ', ok_gw_nonoro
     421  write(numout,*)' nbapp_rad = ', nbapp_rad
     422  write(numout,*)' nbapp_chim = ', nbapp_chim
     423  write(numout,*)' iflag_con = ', iflag_con
    351424  write(numout,*)' Sortie mensuelle = ', ok_mensuel
    352425  write(numout,*)' Sortie journaliere = ', ok_journe
    353426  write(numout,*)' Sortie instantanee = ', ok_instan
     427  write(numout,*)' frequence sorties = ', ecritphy 
    354428  write(numout,*)' Sortie bilan d energie, if_ebil =', if_ebil
    355429  write(numout,*)' Duree de l annee = ',year_day
Note: See TracChangeset for help on using the changeset viewer.