!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/conf_phys.F90,v 1.3 2005/02/07 15:15:31 fairhead Exp $
!
!
!

  subroutine conf_phys(ok_journe, ok_mensuel, ok_instan, &
 &                     if_ebil)
   use init_print_control_mod, only: init_print_control
   use print_control_mod, only: lunout
   use IOIPSL, only: getin
   use age_of_air_mod, only: ok_aoa, reinit_aoa, lev_aoa

   implicit none

   include "YOMCST.h"
   include "clesphys.h"
   include "compbl.h"

! ok_journe:  sorties journalieres
! ok_mensuel: sorties mensuelles
! ok_instan:  sorties instantanees


! Sortie:
  logical,intent(out)  :: ok_journe, ok_mensuel, ok_instan        
  integer,intent(out)  :: if_ebil



!
! Configuration de la "physique" de LMDZ a l'aide de la fonction
! GETIN de IOIPSL
!
! LF 05/2001
!
!--- Ca lit le physiq.def ---

       ! initialize print_control module variables
       call init_print_control 

!******************* parametres anciennement lus dans gcm.def

       ! do we read a startphy.nc file? (default: .true.)
       startphy_file=.true.
       CALL getin("startphy_file",startphy_file)

!Config  Key  = cycle_diurne
!Config  Desc = Cycle diurne
!Config  Def  = y
!Config  Help = Cette option permet d'eteidre le cycle diurne.
!Config         Peut etre util pour accelerer le code !
       cycle_diurne = .TRUE.       
       call getin('cycle_diurne',cycle_diurne)

!Config  Key  = soil_model
!Config  Desc = Modele de sol
!Config  Def  = y
!Config  Help = Choix du modele de sol (Thermique ?)
!Config         Option qui pourait un string afin de pouvoir
!Config         plus de choix ! Ou meme une liste d'options !
       soil_model = .true.
       call getin('soil_model',soil_model)

!Config  Key  = ok_orodr
!Config  Desc = Oro drag
!Config  Def  = y
!Config  Help = GW drag orographie
!Config         
       ok_orodr = .false.
       call getin('ok_orodr',ok_orodr)

!Config  Key  =  ok_orolf
!Config  Desc = Oro lift
!Config  Def  = n
!Config  Help = GW lift orographie (pas utilise)
       ok_orolf = .false.
       call getin('ok_orolf', ok_orolf)

!Config  Key  = ok_gw_nonoro
!Config  Desc = Gravity waves parameterization
!Config  Def  = n
!Config  Help = GW drag non-orographique
       ok_gw_nonoro = .false.
       call getin('ok_gw_nonoro',ok_gw_nonoro)

!Config  Key  = nbapp_rad
!Config  Desc = Frequence d'appel au rayonnement
!Config  Def  = 12
!Config  Help = Nombre  d'appels des routines de rayonnements
!Config         par jour.
       nbapp_rad = 12
       call getin('nbapp_rad',nbapp_rad)
       print*,"nbapp_rad",nbapp_rad
!Config  Key  = nbapp_chem
!Config  Desc = Frequence d'appel a la chimie
!Config  Def  = 24000
!Config  Help = Nombre  d'appels des routines de chimie
!Config         par jour.
       nbapp_chem = 24000
       call getin('nbapp_chem',nbapp_chem)

!Config  Key  = iflag_con
!Config  Desc = Flag de convection
!Config  Def  = 0
!Config  Help = Flag  pour la convection les options suivantes existent :
!Config         0 : ajsec simple (VENUS, TITAN)
!Config         1 pour LMD,
!Config         2 pour Tiedtke,
!Config         3 pour CCM(NCAR)  
       iflag_con = 0
       call getin('iflag_con',iflag_con)

!******************* fin parametres anciennement lus dans gcm.def

!Config Key  = OK_journe
!Config Desc = Pour des sorties journalieres 
!Config Def  = .false.
!Config Help = Pour creer le fichier histday contenant les sorties
!              journalieres 
!
  ok_journe = .false.
  call getin('OK_journe', ok_journe)
!
!Config Key  = OK_mensuel
!Config Desc = Pour des sorties mensuelles 
!Config Def  = .false.
!Config Help = Pour creer le fichier histmth contenant les sorties
!              mensuelles 
!
  ok_mensuel = .false.
  call getin('OK_mensuel', ok_mensuel)
!
!Config Key  = OK_instan
!Config Desc = Pour des sorties instantanees 
!Config Def  = .false.
!Config Help = Pour creer le fichier histins contenant les sorties
!              instantanees 
!
  ok_instan = .false.
  call getin('OK_instan', ok_instan)
!
!Config  Key  = ecritphy
!Config  Desc = Frequence d'ecriture dans histmth et histins
!Config  Def  = 1
!Config  Help = frequence de l'ecriture du fichier histmth et histins
!Config         en jours.
!
       ecriphy = 1.
       call getin('ecritphy', ecriphy)
!
!
!Config Key  = if_ebil
!Config Desc = Niveau de sortie pour les diags bilan d'energie 
!Config Def  = 0
!Config Help = 
!               
!
  if_ebil = 0
  call getin('if_ebil', if_ebil)
!!
!! Parametres orbitaux
!!
!Config Key  = R_ecc
!Config Desc = Excentricite
!Config Def  = 0.006787
!Config Help = 
!               
! VENUS
! R_ecc = 0.006787
  R_ecc   = 0.0
  call getin('R_ecc', R_ecc)
!!
!Config Key  = R_peri
!Config Desc = Equinoxe
!Config Def  = 
!Config Help = 
!               
! VENUS
  R_peri = 0.
  call getin('R_peri', R_peri)
!!
!Config Key  = R_incl
!Config Desc = Inclinaison
!Config Def  = 
!Config Help = 
!               
! VENUS
  R_incl = 0.0
  call getin('R_incl', R_incl)
!
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! PARAMETER FOR THE PLANETARY BOUNDARY LAYER AND SOIL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!Config Key  = iflag_pbl
!Config Desc =
!Config Def  = 1
!Config Help =
!
! 2   = calculs Cd et K simples pour VENUS :
!       parametres = z0, lmixmin, ksta (en dur: umin2,ric,cepdu2,karman)
! 1   = calculs Cd et K issus LMDZ Terre
!       parametres = ksta, ok_kzmin (et plein d'autres en dur...)
! 6-9 = schema des thermiques Fred
  iflag_pbl = 1
  call getin('iflag_pbl',iflag_pbl)

!
!Config Key  = ksta
!Config Desc =
!Config Def  = 1.0e-7
!Config Help =
!
  ksta = 1.0e-7
  call getin('ksta',ksta)

!
!Config Key  = z0
!Config Desc =
!Config Def  = 1.0e-2
!Config Help =
!
  z0 = 1.0e-2
  call getin('z0',z0)

!
!Config Key  = lmixmin
!Config Desc =
!Config Def  = 35.
!Config Help =
!
  lmixmin = 35.
  call getin('lmixmin',lmixmin)

!
!Config Key  = ok_kzmin
!Config Desc =
!Config Def  = .false.
!Config Help =
!
  ok_kzmin = .false.
  call getin('ok_kzmin',ok_kzmin)

  ok_clmain = .true.
  call getin('ok_clmain',ok_clmain)

  physideal = .false.
  call getin('physideal',physideal)

!Config Key  = iflag_ajs
!Config Desc =
!Config Def  = 0
!Config Help =
!
  iflag_ajs = 0
  call getin('iflag_ajs',iflag_ajs)

!
!Config Key  = inertie
!Config Desc =
!Config Def  = 2000.
!Config Help =
!
  inertie = 2000.
  call getin('inertie',inertie)
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! PARAMETER FOR THE OUTPUT LEVELS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!Config Key  = lev_histins
!Config Desc =
!Config Def  = 0
!Config Help =
!
  lev_histins = 0
  call getin('lev_histins',lev_histins)

!
!Config Key  = lev_histday
!Config Desc =
!Config Def  = 1
!Config Help =
!
  lev_histday = 1
  call getin('lev_histday',lev_histday)

!
!Config Key  = lev_histmth
!Config Desc =
!Config Def  = 2
!Config Help =
!
  lev_histmth = 2
  call getin('lev_histmth',lev_histmth)


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! PARAMETER FOR THE TRACERS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!Config Key  = tr_scheme
!Config Desc =
!Config Def  = 0
!Config Help =
!
! 0   = Nothing is done (passive tracers)
! 1   = pseudo-chemistry with relaxation toward fixed profile
!       See Marcq&Lebonnois 2013
! 2   = surface emission
!       For the moment, inspired from Mars version
!       However, the variable 'source' could be used in physiq
!       so the call to phytrac_emiss could be to initialise it.
! 3   = Full chemistry and/or clouds => phytrac_chimie
!       Need ok_chem or ok_cloud
  tr_scheme = 0
  call getin('tr_scheme',tr_scheme)

!Config Key = ok_aoa
!Config Desc = Flag to run age of air subroutine
!Config Def = .false.
  ok_aoa=.false.
  call getin('ok_aoa',ok_aoa)

!Config Key = reinit_aoa
!Config Desc = Flag to reinitiate age of air array to 0
!Config Def = .true.
  reinit_aoa=.true.
  call getin('reinit_aoa',reinit_aoa)
  
!Config Key = lev_aoa
!Config Desc = Model level where tracer associated with age of air is released
!Config Def = 1 (surface)
  lev_aoa=1
  call getin('lev_aoa',lev_aoa)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!   PARAMETRES DE LA CHIMIE/NUAGE dans physiq.def
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!
!Config Key  = reinit_trac
!Config Desc =  
!Config Def  = .FALSE.
!Config Help = 
!
  reinit_trac = .FALSE.
  call getin('reinit_trac',reinit_trac)
  
!
!Config Key  = ok_cloud
!Config Desc =  
!Config Def  = .FALSE.
!Config Help = 
!
  ok_cloud = .false.
  call getin('ok_cloud',ok_cloud)

!
!Config Key  = cl_scheme
!Config Desc =
!Config Def  = 2
!Config Help =
!
! 1   = Simple microphysics (Aurelien Stolzenbach's PhD)
! 2   = Full microphysics (momentum scheme, Sabrina Guilbon's PhD)

  cl_scheme = 2
  call getin('cl_scheme',cl_scheme)

!
!Config Key  = ok_chem
!Config Desc =  
!Config Def  = .FALSE.
!Config Help = 
!
  ok_chem = .false.
  call getin('ok_chem',ok_chem)

  if (((tr_scheme.ne.3).and.(ok_chem.or.ok_cloud)).or. &
      ((tr_scheme.eq.3).and.(.not.ok_chem.and..not.ok_cloud))) then
    write(*,*) "Attention, incoherence :" 
    write(*,*) "tr_scheme=",tr_scheme," / ok_chem=",ok_chem, &
                                     " / ok_cloud=",ok_cloud
    write(*,*) "Verifier votre physiq.def"
    stop
  endif


! Check that chemistry and age of air flages are not both active
  if ((ok_chem) .and. (ok_aoa)) then
     write(*,*) "Chemistry and age of air cannot be run together"
     write(*,*) "ok_chem=",ok_chem," / ok_aoa=",ok_aoa
     stop
  endif



!
!Config Key  = ok_sedim
!Config Desc =  
!Config Def  = .FALSE.
!Config Help = 
!
  ok_sedim = .false.
  call getin('ok_sedim',ok_sedim)

!
!Config Key  = nb_mode
!Config Desc =  
!Config Def  = 0
!Config Help = 
!
  nb_mode = 0
  call getin('nb_mode',nb_mode)
  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! PARAMETER FOR SOLAR RADIATION
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!
!Config Key  = solarchoice
!Config Desc = 
!Config Def  = 1
!Config Help = 
!
! 1 = RH Tables
! 2 = Generic module
  solarchoice = 1
  call getin('solarchoice',solarchoice)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! PARAMETER FOR NLTE PHYSICS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!
!Config Key  = callnlte
!Config Desc = 
!Config Def  = .false.
!Config Help = 
!
  callnlte = .false.
  call getin('callnlte',callnlte)

!
!Config Key  = callnirco2
!Config Desc = 
!Config Def  = .false.
!
  callnirco2 = .false.
  call getin('callnirco2',callnirco2)

!
!Config Key  = nircorr
!Config Desc = 
!Config Def  = 0
!Config Help =
!
  nircorr = 0
  call getin('nircorr',nircorr)

!
!Config Key  = callthermos
!Config Desc = 
!Config Def  = .false.
!Config Help =
!
  callthermos = .false.
  call getin('callthermos',callthermos)

!
!Config Key  = nltemodel
!Config Desc = 
!Config Def  = 0
!Config Help =
!
  nltemodel = 0
  call getin('nltemodel',nltemodel)

!
!Config Key  = solvarmod
!Config Desc =
!Config Def  = 1
!Config Help =
!
  solvarmod = 1
  call getin('solvarmod',solvarmod)

!
!Config Key  = fixed_euv_value
!Config Desc = 
!Config Def  = 140.  ## Average solar cycle condition
!Config Help =
!
  fixed_euv_value =140. 
  call getin('fixed_euv_value',fixed_euv_value)

!
!Config Key  = euveff
!Config Desc =
!Config Def  = 0.22
!Config Help =
!
  euveff = 0.22
  call getin('euveff',euveff)

!Config Key  = tuneupperatm
!Config Desc =
!Config Def  = .false.
!Config Help = Activate thermosphere tuning for oxygen compensation
!
  tuneupperatm = .false.
  call getin('tuneupperatm',tuneupperatm)

!
!Config Key  = ok_jonline
!Config Desc =
!Config Def  = .false.
!Config Help =
!
  ok_jonline = .false.
  call getin('ok_jonline',ok_jonline)

!
!Config Key  = ok_ionchem
!Config Desc =
!Config Def  = .false.
!Config Help = Activate the ionosphere chemistry
!
  ok_ionchem = .false.
  call getin('ok_ionchem',ok_ionchem)

  if ((ok_jonline.eqv..false.).and.(.true..eqv.ok_ionchem)) then
    write(*,*) "Error incoherent flags :" 
    write(*,*) "ok_jonline=",ok_jonline," / ok_ionchem=",ok_ionchem
    write(*,*) "If you include ions, ok_jonline==.true."
    write(*,*) "If you do not include ions, ok_ionchem==.false."
    write(*,*) "Check physiq.def"
    stop
  endif

!
!Config Key  = ok_iondiff
!Config Desc =
!Config Def  = .false.
!Config Help = Activate the ambipolar ion diffusion in physic
!
  ok_iondiff = .false.
  call getin('ok_iondiff',ok_iondiff)

  if ((ok_ionchem.eqv..false.).and.(.true..eqv.ok_iondiff)) then
    write(*,*) "Error incoherent flags :"
    write(*,*) "ok_ionchem=",ok_ionchem," / ok_iondiff=",ok_iondiff
    write(*,*) "If you include ions diffusion,"
    write(*,*) "   you need also ok_ionchem==.true."
    write(*,*) "If you do not include ions diffusion,"
    write(*,*) "   you need ok_iondiff==.false."
    write(*,*) "Check physiq.def"
    stop
  endif

!
!
!Config Key  = 
!Config Desc =  
!Config Def  =
!Config Help = 
!
!   =
!  call getin('',)
!
!
!
!

  write(lunout,*)' ##############################################'
  write(lunout,*)' Configuration des parametres de la physique: '
  write(lunout,*)' cycle_diurne = ', cycle_diurne
  write(lunout,*)' soil_model = ', soil_model
  write(lunout,*)' ok_orodr = ', ok_orodr
  write(lunout,*)' ok_orolf = ', ok_orolf
  write(lunout,*)' ok_gw_nonoro = ', ok_gw_nonoro
  write(lunout,*)' nbapp_rad = ', nbapp_rad
  write(lunout,*)' nbapp_chem = ', nbapp_chem
  write(lunout,*)' iflag_con = ', iflag_con
  write(lunout,*)' Sortie journaliere = ', ok_journe
  write(lunout,*)' Sortie mensuelle = ', ok_mensuel
  write(lunout,*)' Sortie instantanee = ', ok_instan
  write(lunout,*)' frequence sorties = ', ecriphy  
  write(lunout,*)' Sortie bilan d''energie, if_ebil =', if_ebil
  write(lunout,*)' Excentricite = ',R_ecc
  write(lunout,*)' Equinoxe = ',R_peri
  write(lunout,*)' Inclinaison =',R_incl
  write(lunout,*)' tr_scheme = ', tr_scheme
  write(lunout,*)' ok_aoa = ', ok_aoa
  write(lunout,*)' reinit_aoa = ', reinit_aoa
  write(lunout,*)' lev_aoa = ', lev_aoa
  write(lunout,*)' iflag_pbl = ', iflag_pbl
  write(lunout,*)' z0 = ',z0 
  write(lunout,*)' lmixmin = ',lmixmin 
  write(lunout,*)' ksta = ',ksta 
  write(lunout,*)' ok_kzmin = ',ok_kzmin 
  write(lunout,*)' inertie = ', inertie 
  write(lunout,*)' ok_clmain = ',ok_clmain
  write(lunout,*)' physideal = ',physideal
  write(lunout,*)' iflag_ajs = ', iflag_ajs
  write(lunout,*)' lev_histins = ',lev_histins 
  write(lunout,*)' lev_histday = ',lev_histday 
  write(lunout,*)' lev_histmth = ',lev_histmth 
  write(lunout,*)' reinit_trac = ',reinit_trac
  write(lunout,*)' ok_cloud = ',ok_cloud
  write(lunout,*)' cl_scheme = ',cl_scheme
  write(lunout,*)' ok_chem = ',ok_chem
  write(lunout,*)' ok_sedim = ',ok_sedim
  write(lunout,*)' nb_mode = ',nb_mode
  write(lunout,*)' solarchoice = ',solarchoice
  write(lunout,*)' callnlte = ',callnlte
  write(lunout,*)' nltemodel = ',nltemodel
  write(lunout,*)' callnirco2 = ',callnirco2
  write(lunout,*)' nircorr = ',nircorr
  write(lunout,*)' callthermos = ',callthermos
  write(lunout,*)' solvarmod = ',solvarmod
  write(lunout,*)' fixed_euv_value = ',fixed_euv_value
  write(lunout,*)' euveff = ',euveff
  write(lunout,*)' tuneupperatm = ',tuneupperatm
  write(lunout,*)' ok_jonline = ',ok_jonline
  write(lunout,*)' ok_ionchem = ',ok_ionchem
  write(lunout,*)' ok_iondiff = ',ok_iondiff

  end subroutine conf_phys