libf/phylmd/conf_phys.f90

module conf_phys_m

  implicit none

  integer iflag_pbl

contains

  subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
       fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
       ratqsbas,ratqshaut,if_ebil, &
       ok_ade, ok_aie, &
       bl95_b0, bl95_b1,&
       iflag_thermals,nsplit_thermals)

    ! From phylmd/conf_phys.F90,v 1.7 2005/07/05 07:21:23

    use getincom, only: getin
    use clesphys, only: solaire, co2_ppm, rco2, ch4_ppb, rch4, n2o_ppb, rn2o, &
         cfc11_ppt, rcfc11, cfc12_ppt, rcfc12, top_height, overlap, cdmmax, &
         cdhmax, ksta, ksta_ter, ok_kzmin, lev_histhf, lev_histday, &
         lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, &
         latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, &
         ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone
    use yomcst
    use conema3_m
    use comfisrtilp
    use nuagecom

    !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12

    ! Configuration de la "physique" de LMDZ a l'aide de la fonction
    ! GETIN de IOIPSL

    ! ocean:      type d'ocean (force, slab, couple)
    ! ok_veget:   type de modele de vegetation
    ! ok_journe:  sorties journalieres
    ! ok_mensuel: sorties mensuelles
    ! ok_instan:  sorties instantanees
    ! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
    ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc 

    ! Sortie:
    character(len=*), intent(out):: ocean
    logical              :: ok_veget, ok_newmicro
    logical              :: ok_journe, ok_mensuel, ok_instan        
    LOGICAL              :: ok_ade, ok_aie
    REAL                 :: bl95_b0, bl95_b1
    real, intent(out):: fact_cldcon
    real, intent(out):: facttemps
    real ratqsbas,ratqshaut
    integer              :: iflag_cldcon, if_ebil

    ! Local
    real                 :: zzz

    integer :: iflag_thermals,nsplit_thermals

    !-----------------------------------------------------------

    print *, "Call sequence information: conf_phys"

    !Config Key  = OCEAN 
    !Config Desc = Type d'ocean
    !Config Def  = force
    !Config Help = Type d'ocean utilise: force ou slab mais pas couple
    !
    ocean = 'force '
    call getin('OCEAN', ocean)
    !
    !Config Key  = VEGET 
    !Config Desc = Type de modele de vegetation
    !Config Def  = .false.
    !Config Help = Type de modele de vegetation utilise
    !
    ok_veget = .false.
    call getin('VEGET', ok_veget)
    !
    !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  = .true.
    !Config Help = Pour creer le fichier histmth contenant les sorties
    !              mensuelles 
    !
    ok_mensuel = .true.
    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  = ok_ade
    !Config Desc = Aerosol direct effect or not?
    !Config Def  = .false.
    !Config Help = Used in radlwsw.F
    !
    ok_ade = .false.
    call getin('ok_ade', ok_ade)

    !
    !Config Key  = ok_aie
    !Config Desc = Aerosol indirect effect or not?
    !Config Def  = .false.
    !Config Help = Used in nuage.F and radlwsw.F
    !
    ok_aie = .false.
    call getin('ok_aie', ok_aie)

    !
    !Config Key  = bl95_b0
    !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
    !Config Def  = .false.
    !Config Help = Used in nuage.F
    !
    bl95_b0 = 2.
    call getin('bl95_b0', bl95_b0)

    !Config Key  = bl95_b1
    !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
    !Config Def  = .false.
    !Config Help = Used in nuage.F
    !
    bl95_b1 = 0.2
    call getin('bl95_b1', bl95_b1)

    !
    !
    !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)
    !!
    !! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
    !!
    !Config Key  = R_ecc
    !Config Desc = Excentricite
    !Config Def  = 0.016715
    !Config Help = 
    !               
    !valeur AMIP II
    R_ecc = 0.016715
    call getin('R_ecc', R_ecc)
    !!
    !Config Key  = R_peri
    !Config Desc = Equinoxe
    !Config Def  = 
    !Config Help = 
    !               
    !
    !valeur AMIP II
    R_peri = 102.7
    call getin('R_peri', R_peri)
    !!
    !Config Key  = R_incl
    !Config Desc = Inclinaison
    !Config Def  = 
    !Config Help = 
    !               
    !
    !valeur AMIP II
    R_incl = 23.441
    call getin('R_incl', R_incl)
    !!
    !Config Key  = solaire
    !Config Desc = Constante solaire en W/m2
    !Config Def  = 1365.
    !Config Help = 
    !               
    !
    !valeur AMIP II
    solaire = 1365.
    call getin('solaire', solaire)
    !!
    !Config Key  = co2_ppm
    !Config Desc = concentration du gaz carbonique en ppmv
    !Config Def  = 348.
    !Config Help = 
    !               
    !
    !valeur AMIP II
    co2_ppm = 348.
    call getin('co2_ppm', co2_ppm)
    !!
    !Config Key  = RCO2
    !Config Desc = Concentration du CO2
    !Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
    !Config Def  = 348. * 1.0e-06  * 44.011/28.97
    !Config Help = 
    !               
    ! RCO2 = 5.286789092164308E-04
    !ancienne valeur
    RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.

    !!  call getin('RCO2', RCO2)
    !!
    !Config Key  = RCH4
    !Config Desc = Concentration du CH4
    !Config Def  = 1.65E-06* 16.043/28.97
    !Config Help = 
    !               
    !
    !valeur AMIP II
    !OK  RCH4 = 1.65E-06* 16.043/28.97
    ! RCH4 = 9.137366240938903E-07
    !
    !ancienne valeur
    ! RCH4 = 1.72E-06* 16.043/28.97
    !OK call getin('RCH4', RCH4)
    zzz = 1650.
    call getin('CH4_ppb', zzz)
    CH4_ppb = zzz
    RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
    !!
    !Config Key  = RN2O
    !Config Desc = Concentration du N2O
    !Config Def  = 306.E-09* 44.013/28.97
    !Config Help = 
    !               
    !
    !valeur AMIP II
    !OK  RN2O = 306.E-09* 44.013/28.97
    ! RN2O = 4.648939592682085E-07
    !
    !ancienne valeur
    ! RN2O = 310.E-09* 44.013/28.97
    !OK  call getin('RN2O', RN2O)
    zzz=306.
    call getin('N2O_ppb', zzz)
    N2O_ppb = zzz
    RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
    !!
    !Config Key  = RCFC11
    !Config Desc = Concentration du CFC11
    !Config Def  = 280.E-12* 137.3686/28.97
    !Config Help = 
    !               
    !
    !OK RCFC11 = 280.E-12* 137.3686/28.97
    zzz = 280.
    call getin('CFC11_ppt',zzz)
    CFC11_ppt = zzz
    RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
    ! RCFC11 = 1.327690990680013E-09
    !OK call getin('RCFC11', RCFC11)
    !!
    !Config Key  = RCFC12
    !Config Desc = Concentration du CFC12
    !Config Def  = 484.E-12* 120.9140/28.97
    !Config Help = 
    !               
    !
    !OK RCFC12 = 484.E-12* 120.9140/28.97
    zzz = 484.
    call getin('CFC12_ppt',zzz)
    CFC12_ppt = zzz
    RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
    ! RCFC12 = 2.020102726958923E-09
    !OK call getin('RCFC12', RCFC12)
    !!
    !! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
    !!
    !! KE
    !
    !Config Key  = epmax
    !Config Desc = Efficacite precip
    !Config Def  = 0.993
    !Config Help = 
    !
    epmax = .993
    call getin('epmax', epmax)
    !
    !Config Key  = ok_adj_ema
    !Config Desc =  
    !Config Def  = false
    !Config Help = 
    !
    ok_adj_ema = .false.
    call getin('ok_adj_ema',ok_adj_ema)
    !
    !Config Key  = iflag_clw
    !Config Desc =  
    !Config Def  = 0
    !Config Help = 
    !
    iflag_clw = 0
    call getin('iflag_clw',iflag_clw)
    !
    !Config Key  = cld_lc_lsc 
    !Config Desc =  
    !Config Def  = 2.6e-4
    !Config Help = 
    !
    cld_lc_lsc = 2.6e-4
    call getin('cld_lc_lsc',cld_lc_lsc)
    !
    !Config Key  = cld_lc_con
    !Config Desc =  
    !Config Def  = 2.6e-4
    !Config Help = 
    !
    cld_lc_con = 2.6e-4
    call getin('cld_lc_con',cld_lc_con)
    !
    !Config Key  = cld_tau_lsc
    !Config Desc =  
    !Config Def  = 3600.
    !Config Help = 
    !
    cld_tau_lsc = 3600.
    call getin('cld_tau_lsc',cld_tau_lsc)
    !
    !Config Key  = cld_tau_con
    !Config Desc =  
    !Config Def  = 3600.
    !Config Help = 
    !
    cld_tau_con = 3600.
    call getin('cld_tau_con',cld_tau_con)
    !
    !Config Key  = ffallv_lsc
    !Config Desc =  
    !Config Def  = 1.
    !Config Help = 
    !
    ffallv_lsc = 1.
    call getin('ffallv_lsc',ffallv_lsc)
    !
    !Config Key  = ffallv_con
    !Config Desc =  
    !Config Def  = 1.
    !Config Help = 
    !
    ffallv_con = 1.
    call getin('ffallv_con',ffallv_con)
    !
    !Config Key  = coef_eva
    !Config Desc =  
    !Config Def  = 2.e-5
    !Config Help = 
    !
    coef_eva = 2.e-5
    call getin('coef_eva',coef_eva)
    !
    !Config Key  = reevap_ice
    !Config Desc =  
    !Config Def  = .false.
    !Config Help = 
    !
    reevap_ice = .false.
    call getin('reevap_ice',reevap_ice)
    !
    !Config Key  = iflag_cldcon 
    !Config Desc =  
    !Config Def  = 1
    !Config Help = 
    !
    iflag_cldcon = 1
    call getin('iflag_cldcon',iflag_cldcon)

    !
    !Config Key  = iflag_pdf 
    !Config Desc =  
    !Config Def  = 0
    !Config Help = 
    !
    iflag_pdf = 0
    call getin('iflag_pdf',iflag_pdf)
    !
    !Config Key  = fact_cldcon
    !Config Desc =  
    !Config Def  = 0.375
    !Config Help = 
    !
    fact_cldcon = 0.375
    call getin('fact_cldcon',fact_cldcon)

    !
    !Config Key  = facttemps
    !Config Desc =  
    !Config Def  = 1.e-4
    !Config Help = 
    !
    facttemps = 1.e-4
    call getin('facttemps',facttemps)

    !
    !Config Key  = ok_newmicro
    !Config Desc =  
    !Config Def  = .true.
    !Config Help = 
    !
    ok_newmicro = .true.
    call getin('ok_newmicro',ok_newmicro)
    !
    !Config Key  = ratqsbas
    !Config Desc =  
    !Config Def  = 0.01
    !Config Help = 
    !
    ratqsbas = 0.01
    call getin('ratqsbas',ratqsbas)
    !
    !Config Key  = ratqshaut
    !Config Desc =  
    !Config Def  = 0.3
    !Config Help = 
    !
    ratqshaut = 0.3
    call getin('ratqshaut',ratqshaut)

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

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

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

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

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


    !
    !
    !Config Key  = cdmmax
    !Config Desc =
    !Config Def  = 1.3E-3
    !Config Help =
    !
    cdmmax = 1.3E-3
    call getin('cdmmax',cdmmax)

    !
    !Config Key  = cdhmax
    !Config Desc =
    !Config Def  = 1.1E-3
    !Config Help =
    !
    cdhmax = 1.1E-3
    call getin('cdhmax',cdhmax)

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

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

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

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    ! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    !Config Key  = iflag_pbl
    !Config Desc =
    !Config Def  = 1
    !Config Help =
    !
    iflag_pbl = 1
    call getin('iflag_pbl',iflag_pbl)
    !
    !Config Key  = iflag_thermals
    !Config Desc =
    !Config Def  = 0
    !Config Help =
    !
    iflag_thermals = 0
    call getin('iflag_thermals',iflag_thermals)
    !
    !
    !Config Key  = nsplit_thermals
    !Config Desc =
    !Config Def  = 1
    !Config Help =
    !
    nsplit_thermals = 1
    call getin('nsplit_thermals',nsplit_thermals)


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

    !
    !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)

    !
    !Config Key  = type_run
    !Config Desc =
    !Config Def  = 'AMIP' ou 'ENSP'
    !Config Help =
    !
    type_run = 'AMIP'
    call getin('type_run',type_run)

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

    !
    !
    !Config Key  = ok_regdyn
    !Config Desc =
    !Config Def  = 'AMIP'
    !Config Help =
    !
    ! ok_regdyn = .true.
    ok_regdyn = .false.
    call getin('ok_regdyn',ok_regdyn)
    !
    ! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone 
    ! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
    !
    !Config Key  = lonmin_ins
    !Config Desc = 100.  
    !Config Def  = longitude minimale sorties "bilKP_ins"
    !Config Help = 
    !
    lonmin_ins = 100.
    call getin('lonmin_ins',lonmin_ins)
    !
    !Config Key  = lonmax_ins
    !Config Desc = 130. 
    !Config Def  = longitude maximale sorties "bilKP_ins"
    !Config Help =
    !
    lonmax_ins = 130.
    call getin('lonmax_ins',lonmax_ins)
    !
    !Config Key  = latmin_ins
    !Config Desc = -20.  
    !Config Def  = latitude minimale sorties "bilKP_ins"
    !Config Help = 
    !
    latmin_ins = -20.
    call getin('latmin_ins',latmin_ins)
    !
    !Config Key  = latmax_ins
    !Config Desc = 20. 
    !Config Def  = latitude maximale sorties "bilKP_ins"
    !Config Help =
    !
    latmax_ins = 20.
    call getin('latmax_ins',latmax_ins)
    !
    !Config Key  = ecrit_ins
    !Config Desc =
    !Config Def  = NINT(86400./dtime/48.) !a chaque pas de temps physique
    !Config Help =
    !
    ! ecrit_ins = NINT(86400./dtime/48.)
    ecrit_ins = NINT(86400./48.)
    call getin('ecrit_ins',ecrit_ins)
    !
    !Config Key  = ecrit_hf
    !Config Desc =
    !Config Def  = NINT(86400./dtime *0.25) !toutes les 6h
    !Config Help =
    !
    ! ecrit_hf = NINT(86400./dtime *0.25)
    ecrit_hf = NINT(86400. *0.25)
    call getin('ecrit_hf',ecrit_hf)
    !
    !Config Key  = ecrit_hf2mth
    !Config Desc =
    !Config Def  = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h
    !Config Help =
    !
    ecrit_hf2mth = 4*30
    call getin('ecrit_hf2mth',ecrit_hf2mth)
    !
    !Config Key  = ecrit_day
    ecrit_day = 86400
    call getin('ecrit_day',ecrit_day)
    !
    ecrit_mth = 86400
    call getin('ecrit_mth',ecrit_mth)
    !
    ecrit_tra = 1
    call getin('ecrit_tra',ecrit_tra)
    !
    !Config Key  = ecrit_reg
    !Config Desc =
    !Config Def  = NINT(86400./dtime *0.25)  !4 fois par jour
    !Config Help =
    !
    ! ecrit_reg = NINT(86400./dtime *0.25)  !4 fois par jour
    ecrit_reg = NINT(86400. *0.25)  !4 fois par jour
    call getin('ecrit_reg',ecrit_reg)
    !
    !
    !
    !Config Key  = bug_ozone
    !Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true
    !Config Def  = false
    !Config Help =
    !
    bug_ozone = .false.
    call getin('bug_ozone',bug_ozone)

    print *, ' *********'
    print *, ' Configuration des parametres de la physique: '
    print *, ' Config ocean = ', ocean
    print *, ' Config veget = ', ok_veget
    print *, ' Sortie journaliere = ', ok_journe
    print *, ' Sortie mensuelle = ', ok_mensuel
    print *, ' Sortie instantanee = ', ok_instan
    print *, ' Sortie bilan d''energie, if_ebil =', if_ebil
    print *, ' Excentricite = ',R_ecc
    print *, ' Equinoxe = ',R_peri
    print *, ' Inclinaison =',R_incl
    print *, ' Constante solaire =',solaire
    print *, ' co2_ppm =',co2_ppm
    print *, ' RCO2 = ',RCO2
    print *, ' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4
    print *, ' N2O_ppb =',N2O_ppb,' RN2O =  ',RN2O
    print *, ' CFC11_ppt=',CFC11_ppt,' RCFC11 =  ',RCFC11
    print *, ' CFC12_ppt=',CFC12_ppt,' RCFC12 =  ',RCFC12
    print *, ' epmax = ', epmax
    print *, ' ok_adj_ema = ', ok_adj_ema
    print *, ' iflag_clw = ', iflag_clw
    print *, ' cld_lc_lsc = ', cld_lc_lsc
    print *, ' cld_lc_con = ', cld_lc_con
    print *, ' cld_tau_lsc = ', cld_tau_lsc
    print *, ' cld_tau_con = ', cld_tau_con
    print *, ' ffallv_lsc = ', ffallv_lsc
    print *, ' ffallv_con = ', ffallv_con
    print *, ' coef_eva = ', coef_eva
    print *, ' reevap_ice = ', reevap_ice
    print *, ' iflag_pdf = ', iflag_pdf
    print *, ' iflag_cldcon = ', iflag_cldcon
    print *, ' fact_cldcon = ', fact_cldcon
    print *, ' facttemps = ', facttemps
    print *, ' ok_newmicro = ',ok_newmicro 
    print *, ' ratqsbas = ',ratqsbas 
    print *, ' ratqshaut = ',ratqshaut 
    print *, ' top_height = ',top_height 
    print *, ' overlap = ',overlap 
    print *, ' cdmmax = ',cdmmax 
    print *, ' cdhmax = ',cdhmax 
    print *, ' ksta = ',ksta 
    print *, ' ksta_ter = ',ksta_ter 
    print *, ' ok_kzmin = ',ok_kzmin 
    print *, ' ok_ade = ',ok_ade
    print *, ' ok_aie = ',ok_aie
    print *, ' bl95_b0 = ',bl95_b0
    print *, ' bl95_b1 = ',bl95_b1
    print *, ' lev_histhf = ',lev_histhf 
    print *, ' lev_histday = ',lev_histday 
    print *, ' lev_histmth = ',lev_histmth 
    print *, ' iflag_pbl = ', iflag_pbl
    print *, ' iflag_thermals = ', iflag_thermals
    print *, ' type_run = ',type_run 
    print *, ' ok_isccp = ',ok_isccp 
    print *, ' ok_regdyn = ',ok_regdyn
    print *, ' lonmin lonmax latmin latmax bilKP_ins =',&
         lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
    print *,  'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, &
         ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
    print *, ' bug_ozone = ', bug_ozone

  end subroutine conf_phys

end module conf_phys_m
1	module conf_phys_m
2
3	implicit none
4
5	integer iflag_pbl
6
7	contains
8
9	subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
10	fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, &
11	ratqsbas,ratqshaut,if_ebil, &
12	ok_ade, ok_aie, &
13	bl95_b0, bl95_b1,&
14	iflag_thermals,nsplit_thermals)
15
16	! From phylmd/conf_phys.F90,v 1.7 2005/07/05 07:21:23
17
18	use getincom, only: getin
19	use clesphys, only: solaire, co2_ppm, rco2, ch4_ppb, rch4, n2o_ppb, rn2o, &
20	cfc11_ppt, rcfc11, cfc12_ppt, rcfc12, top_height, overlap, cdmmax, &
21	cdhmax, ksta, ksta_ter, ok_kzmin, lev_histhf, lev_histday, &
22	lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, &
23	latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, &
24	ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone
25	use yomcst
26	use conema3_m
27	use comfisrtilp
28	use nuagecom
29
30	!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
31
32	! Configuration de la "physique" de LMDZ a l'aide de la fonction
33	! GETIN de IOIPSL
34
35	! ocean: type d'ocean (force, slab, couple)
36	! ok_veget: type de modele de vegetation
37	! ok_journe: sorties journalieres
38	! ok_mensuel: sorties mensuelles
39	! ok_instan: sorties instantanees
40	! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
41	! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
42
43	! Sortie:
44	character(len=*), intent(out):: ocean
45	logical :: ok_veget, ok_newmicro
46	logical :: ok_journe, ok_mensuel, ok_instan
47	LOGICAL :: ok_ade, ok_aie
48	REAL :: bl95_b0, bl95_b1
49	real, intent(out):: fact_cldcon
50	real, intent(out):: facttemps
51	real ratqsbas,ratqshaut
52	integer :: iflag_cldcon, if_ebil
53
54	! Local
55	real :: zzz
56
57	integer :: iflag_thermals,nsplit_thermals
58
59	!-----------------------------------------------------------
60
61	print *, "Call sequence information: conf_phys"
62
63	!Config Key = OCEAN
64	!Config Desc = Type d'ocean
65	!Config Def = force
66	!Config Help = Type d'ocean utilise: force ou slab mais pas couple
67	!
68	ocean = 'force '
69	call getin('OCEAN', ocean)
70	!
71	!Config Key = VEGET
72	!Config Desc = Type de modele de vegetation
73	!Config Def = .false.
74	!Config Help = Type de modele de vegetation utilise
75	!
76	ok_veget = .false.
77	call getin('VEGET', ok_veget)
78	!
79	!Config Key = OK_journe
80	!Config Desc = Pour des sorties journalieres
81	!Config Def = .false.
82	!Config Help = Pour creer le fichier histday contenant les sorties
83	! journalieres
84	!
85	ok_journe = .false.
86	call getin('OK_journe', ok_journe)
87	!
88	!Config Key = OK_mensuel
89	!Config Desc = Pour des sorties mensuelles
90	!Config Def = .true.
91	!Config Help = Pour creer le fichier histmth contenant les sorties
92	! mensuelles
93	!
94	ok_mensuel = .true.
95	call getin('OK_mensuel', ok_mensuel)
96	!
97	!Config Key = OK_instan
98	!Config Desc = Pour des sorties instantanees
99	!Config Def = .false.
100	!Config Help = Pour creer le fichier histins contenant les sorties
101	! instantanees
102	!
103	ok_instan = .false.
104	call getin('OK_instan', ok_instan)
105	!
106	!Config Key = ok_ade
107	!Config Desc = Aerosol direct effect or not?
108	!Config Def = .false.
109	!Config Help = Used in radlwsw.F
110	!
111	ok_ade = .false.
112	call getin('ok_ade', ok_ade)
113
114	!
115	!Config Key = ok_aie
116	!Config Desc = Aerosol indirect effect or not?
117	!Config Def = .false.
118	!Config Help = Used in nuage.F and radlwsw.F
119	!
120	ok_aie = .false.
121	call getin('ok_aie', ok_aie)
122
123	!
124	!Config Key = bl95_b0
125	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
126	!Config Def = .false.
127	!Config Help = Used in nuage.F
128	!
129	bl95_b0 = 2.
130	call getin('bl95_b0', bl95_b0)
131
132	!Config Key = bl95_b1
133	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
134	!Config Def = .false.
135	!Config Help = Used in nuage.F
136	!
137	bl95_b1 = 0.2
138	call getin('bl95_b1', bl95_b1)
139
140	!
141	!
142	!Config Key = if_ebil
143	!Config Desc = Niveau de sortie pour les diags bilan d'energie
144	!Config Def = 0
145	!Config Help =
146	!
147	!
148	if_ebil = 0
149	call getin('if_ebil', if_ebil)
150	!!
151	!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
152	!!
153	!Config Key = R_ecc
154	!Config Desc = Excentricite
155	!Config Def = 0.016715
156	!Config Help =
157	!
158	!valeur AMIP II
159	R_ecc = 0.016715
160	call getin('R_ecc', R_ecc)
161	!!
162	!Config Key = R_peri
163	!Config Desc = Equinoxe
164	!Config Def =
165	!Config Help =
166	!
167	!
168	!valeur AMIP II
169	R_peri = 102.7
170	call getin('R_peri', R_peri)
171	!!
172	!Config Key = R_incl
173	!Config Desc = Inclinaison
174	!Config Def =
175	!Config Help =
176	!
177	!
178	!valeur AMIP II
179	R_incl = 23.441
180	call getin('R_incl', R_incl)
181	!!
182	!Config Key = solaire
183	!Config Desc = Constante solaire en W/m2
184	!Config Def = 1365.
185	!Config Help =
186	!
187	!
188	!valeur AMIP II
189	solaire = 1365.
190	call getin('solaire', solaire)
191	!!
192	!Config Key = co2_ppm
193	!Config Desc = concentration du gaz carbonique en ppmv
194	!Config Def = 348.
195	!Config Help =
196	!
197	!
198	!valeur AMIP II
199	co2_ppm = 348.
200	call getin('co2_ppm', co2_ppm)
201	!!
202	!Config Key = RCO2
203	!Config Desc = Concentration du CO2
204	!Config Def = co2_ppm * 1.0e-06 * 44.011/28.97
205	!Config Def = 348. * 1.0e-06 * 44.011/28.97
206	!Config Help =
207	!
208	! RCO2 = 5.286789092164308E-04
209	!ancienne valeur
210	RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97 ! pour co2_ppm=348.
211
212	!! call getin('RCO2', RCO2)
213	!!
214	!Config Key = RCH4
215	!Config Desc = Concentration du CH4
216	!Config Def = 1.65E-06* 16.043/28.97
217	!Config Help =
218	!
219	!
220	!valeur AMIP II
221	!OK RCH4 = 1.65E-06* 16.043/28.97
222	! RCH4 = 9.137366240938903E-07
223	!
224	!ancienne valeur
225	! RCH4 = 1.72E-06* 16.043/28.97
226	!OK call getin('RCH4', RCH4)
227	zzz = 1650.
228	call getin('CH4_ppb', zzz)
229	CH4_ppb = zzz
230	RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
231	!!
232	!Config Key = RN2O
233	!Config Desc = Concentration du N2O
234	!Config Def = 306.E-09* 44.013/28.97
235	!Config Help =
236	!
237	!
238	!valeur AMIP II
239	!OK RN2O = 306.E-09* 44.013/28.97
240	! RN2O = 4.648939592682085E-07
241	!
242	!ancienne valeur
243	! RN2O = 310.E-09* 44.013/28.97
244	!OK call getin('RN2O', RN2O)
245	zzz=306.
246	call getin('N2O_ppb', zzz)
247	N2O_ppb = zzz
248	RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
249	!!
250	!Config Key = RCFC11
251	!Config Desc = Concentration du CFC11
252	!Config Def = 280.E-12* 137.3686/28.97
253	!Config Help =
254	!
255	!
256	!OK RCFC11 = 280.E-12* 137.3686/28.97
257	zzz = 280.
258	call getin('CFC11_ppt',zzz)
259	CFC11_ppt = zzz
260	RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
261	! RCFC11 = 1.327690990680013E-09
262	!OK call getin('RCFC11', RCFC11)
263	!!
264	!Config Key = RCFC12
265	!Config Desc = Concentration du CFC12
266	!Config Def = 484.E-12* 120.9140/28.97
267	!Config Help =
268	!
269	!
270	!OK RCFC12 = 484.E-12* 120.9140/28.97
271	zzz = 484.
272	call getin('CFC12_ppt',zzz)
273	CFC12_ppt = zzz
274	RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
275	! RCFC12 = 2.020102726958923E-09
276	!OK call getin('RCFC12', RCFC12)
277	!!
278	!! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
279	!!
280	!! KE
281	!
282	!Config Key = epmax
283	!Config Desc = Efficacite precip
284	!Config Def = 0.993
285	!Config Help =
286	!
287	epmax = .993
288	call getin('epmax', epmax)
289	!
290	!Config Key = ok_adj_ema
291	!Config Desc =
292	!Config Def = false
293	!Config Help =
294	!
295	ok_adj_ema = .false.
296	call getin('ok_adj_ema',ok_adj_ema)
297	!
298	!Config Key = iflag_clw
299	!Config Desc =
300	!Config Def = 0
301	!Config Help =
302	!
303	iflag_clw = 0
304	call getin('iflag_clw',iflag_clw)
305	!
306	!Config Key = cld_lc_lsc
307	!Config Desc =
308	!Config Def = 2.6e-4
309	!Config Help =
310	!
311	cld_lc_lsc = 2.6e-4
312	call getin('cld_lc_lsc',cld_lc_lsc)
313	!
314	!Config Key = cld_lc_con
315	!Config Desc =
316	!Config Def = 2.6e-4
317	!Config Help =
318	!
319	cld_lc_con = 2.6e-4
320	call getin('cld_lc_con',cld_lc_con)
321	!
322	!Config Key = cld_tau_lsc
323	!Config Desc =
324	!Config Def = 3600.
325	!Config Help =
326	!
327	cld_tau_lsc = 3600.
328	call getin('cld_tau_lsc',cld_tau_lsc)
329	!
330	!Config Key = cld_tau_con
331	!Config Desc =
332	!Config Def = 3600.
333	!Config Help =
334	!
335	cld_tau_con = 3600.
336	call getin('cld_tau_con',cld_tau_con)
337	!
338	!Config Key = ffallv_lsc
339	!Config Desc =
340	!Config Def = 1.
341	!Config Help =
342	!
343	ffallv_lsc = 1.
344	call getin('ffallv_lsc',ffallv_lsc)
345	!
346	!Config Key = ffallv_con
347	!Config Desc =
348	!Config Def = 1.
349	!Config Help =
350	!
351	ffallv_con = 1.
352	call getin('ffallv_con',ffallv_con)
353	!
354	!Config Key = coef_eva
355	!Config Desc =
356	!Config Def = 2.e-5
357	!Config Help =
358	!
359	coef_eva = 2.e-5
360	call getin('coef_eva',coef_eva)
361	!
362	!Config Key = reevap_ice
363	!Config Desc =
364	!Config Def = .false.
365	!Config Help =
366	!
367	reevap_ice = .false.
368	call getin('reevap_ice',reevap_ice)
369	!
370	!Config Key = iflag_cldcon
371	!Config Desc =
372	!Config Def = 1
373	!Config Help =
374	!
375	iflag_cldcon = 1
376	call getin('iflag_cldcon',iflag_cldcon)
377
378	!
379	!Config Key = iflag_pdf
380	!Config Desc =
381	!Config Def = 0
382	!Config Help =
383	!
384	iflag_pdf = 0
385	call getin('iflag_pdf',iflag_pdf)
386	!
387	!Config Key = fact_cldcon
388	!Config Desc =
389	!Config Def = 0.375
390	!Config Help =
391	!
392	fact_cldcon = 0.375
393	call getin('fact_cldcon',fact_cldcon)
394
395	!
396	!Config Key = facttemps
397	!Config Desc =
398	!Config Def = 1.e-4
399	!Config Help =
400	!
401	facttemps = 1.e-4
402	call getin('facttemps',facttemps)
403
404	!
405	!Config Key = ok_newmicro
406	!Config Desc =
407	!Config Def = .true.
408	!Config Help =
409	!
410	ok_newmicro = .true.
411	call getin('ok_newmicro',ok_newmicro)
412	!
413	!Config Key = ratqsbas
414	!Config Desc =
415	!Config Def = 0.01
416	!Config Help =
417	!
418	ratqsbas = 0.01
419	call getin('ratqsbas',ratqsbas)
420	!
421	!Config Key = ratqshaut
422	!Config Desc =
423	!Config Def = 0.3
424	!Config Help =
425	!
426	ratqshaut = 0.3
427	call getin('ratqshaut',ratqshaut)
428
429	!
430	!Config Key = rad_froid
431	!Config Desc =
432	!Config Def = 35.0
433	!Config Help =
434	!
435	rad_froid = 35.0
436	call getin('rad_froid',rad_froid)
437
438	!
439	!Config Key = rad_chau1
440	!Config Desc =
441	!Config Def = 13.0
442	!Config Help =
443	!
444	rad_chau1 = 13.0
445	call getin('rad_chau1',rad_chau1)
446
447	!
448	!Config Key = rad_chau2
449	!Config Desc =
450	!Config Def = 9.0
451	!Config Help =
452	!
453	rad_chau2 = 9.0
454	call getin('rad_chau2',rad_chau2)
455
456	!
457	!Config Key = top_height
458	!Config Desc =
459	!Config Def = 3
460	!Config Help =
461	!
462	top_height = 3
463	call getin('top_height',top_height)
464
465	!
466	!Config Key = overlap
467	!Config Desc =
468	!Config Def = 3
469	!Config Help =
470	!
471	overlap = 3
472	call getin('overlap',overlap)
473
474
475	!
476	!
477	!Config Key = cdmmax
478	!Config Desc =
479	!Config Def = 1.3E-3
480	!Config Help =
481	!
482	cdmmax = 1.3E-3
483	call getin('cdmmax',cdmmax)
484
485	!
486	!Config Key = cdhmax
487	!Config Desc =
488	!Config Def = 1.1E-3
489	!Config Help =
490	!
491	cdhmax = 1.1E-3
492	call getin('cdhmax',cdhmax)
493
494	!261103
495	!
496	!Config Key = ksta
497	!Config Desc =
498	!Config Def = 1.0e-10
499	!Config Help =
500	!
501	ksta = 1.0e-10
502	call getin('ksta',ksta)
503
504	!
505	!Config Key = ksta_ter
506	!Config Desc =
507	!Config Def = 1.0e-10
508	!Config Help =
509	!
510	ksta_ter = 1.0e-10
511	call getin('ksta_ter',ksta_ter)
512
513	!
514	!Config Key = ok_kzmin
515	!Config Desc =
516	!Config Def = .true.
517	!Config Help =
518	!
519	ok_kzmin = .true.
520	call getin('ok_kzmin',ok_kzmin)
521
522	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
523	! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
524	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
525
526	!Config Key = iflag_pbl
527	!Config Desc =
528	!Config Def = 1
529	!Config Help =
530	!
531	iflag_pbl = 1
532	call getin('iflag_pbl',iflag_pbl)
533	!
534	!Config Key = iflag_thermals
535	!Config Desc =
536	!Config Def = 0
537	!Config Help =
538	!
539	iflag_thermals = 0
540	call getin('iflag_thermals',iflag_thermals)
541	!
542	!
543	!Config Key = nsplit_thermals
544	!Config Desc =
545	!Config Def = 1
546	!Config Help =
547	!
548	nsplit_thermals = 1
549	call getin('nsplit_thermals',nsplit_thermals)
550
551
552
553	!
554	!Config Key = lev_histhf
555	!Config Desc =
556	!Config Def = 0
557	!Config Help =
558	!
559	lev_histhf = 0
560	call getin('lev_histhf',lev_histhf)
561
562	!
563	!Config Key = lev_histday
564	!Config Desc =
565	!Config Def = 1
566	!Config Help =
567	!
568	lev_histday = 1
569	call getin('lev_histday',lev_histday)
570
571	!
572	!Config Key = lev_histmth
573	!Config Desc =
574	!Config Def = 2
575	!Config Help =
576	!
577	lev_histmth = 2
578	call getin('lev_histmth',lev_histmth)
579
580	!
581	!Config Key = type_run
582	!Config Desc =
583	!Config Def = 'AMIP' ou 'ENSP'
584	!Config Help =
585	!
586	type_run = 'AMIP'
587	call getin('type_run',type_run)
588
589	!
590	!Config Key = ok_isccp
591	!Config Desc =
592	!Config Def = .true.
593	!Config Help =
594	!
595	! ok_isccp = .true.
596	ok_isccp = .false.
597	call getin('ok_isccp',ok_isccp)
598
599	!
600	!
601	!Config Key = ok_regdyn
602	!Config Desc =
603	!Config Def = 'AMIP'
604	!Config Help =
605	!
606	! ok_regdyn = .true.
607	ok_regdyn = .false.
608	call getin('ok_regdyn',ok_regdyn)
609	!
610	! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone
611	! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
612	!
613	!Config Key = lonmin_ins
614	!Config Desc = 100.
615	!Config Def = longitude minimale sorties "bilKP_ins"
616	!Config Help =
617	!
618	lonmin_ins = 100.
619	call getin('lonmin_ins',lonmin_ins)
620	!
621	!Config Key = lonmax_ins
622	!Config Desc = 130.
623	!Config Def = longitude maximale sorties "bilKP_ins"
624	!Config Help =
625	!
626	lonmax_ins = 130.
627	call getin('lonmax_ins',lonmax_ins)
628	!
629	!Config Key = latmin_ins
630	!Config Desc = -20.
631	!Config Def = latitude minimale sorties "bilKP_ins"
632	!Config Help =
633	!
634	latmin_ins = -20.
635	call getin('latmin_ins',latmin_ins)
636	!
637	!Config Key = latmax_ins
638	!Config Desc = 20.
639	!Config Def = latitude maximale sorties "bilKP_ins"
640	!Config Help =
641	!
642	latmax_ins = 20.
643	call getin('latmax_ins',latmax_ins)
644	!
645	!Config Key = ecrit_ins
646	!Config Desc =
647	!Config Def = NINT(86400./dtime/48.) !a chaque pas de temps physique
648	!Config Help =
649	!
650	! ecrit_ins = NINT(86400./dtime/48.)
651	ecrit_ins = NINT(86400./48.)
652	call getin('ecrit_ins',ecrit_ins)
653	!
654	!Config Key = ecrit_hf
655	!Config Desc =
656	!Config Def = NINT(86400./dtime *0.25) !toutes les 6h
657	!Config Help =
658	!
659	! ecrit_hf = NINT(86400./dtime *0.25)
660	ecrit_hf = NINT(86400. *0.25)
661	call getin('ecrit_hf',ecrit_hf)
662	!
663	!Config Key = ecrit_hf2mth
664	!Config Desc =
665	!Config Def = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h
666	!Config Help =
667	!
668	ecrit_hf2mth = 4*30
669	call getin('ecrit_hf2mth',ecrit_hf2mth)
670	!
671	!Config Key = ecrit_day
672	ecrit_day = 86400
673	call getin('ecrit_day',ecrit_day)
674	!
675	ecrit_mth = 86400
676	call getin('ecrit_mth',ecrit_mth)
677	!
678	ecrit_tra = 1
679	call getin('ecrit_tra',ecrit_tra)
680	!
681	!Config Key = ecrit_reg
682	!Config Desc =
683	!Config Def = NINT(86400./dtime *0.25) !4 fois par jour
684	!Config Help =
685	!
686	! ecrit_reg = NINT(86400./dtime *0.25) !4 fois par jour
687	ecrit_reg = NINT(86400. *0.25) !4 fois par jour
688	call getin('ecrit_reg',ecrit_reg)
689	!
690	!
691	!
692	!Config Key = bug_ozone
693	!Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true
694	!Config Def = false
695	!Config Help =
696	!
697	bug_ozone = .false.
698	call getin('bug_ozone',bug_ozone)
699
700	print , ' ********'
701	print *, ' Configuration des parametres de la physique: '
702	print *, ' Config ocean = ', ocean
703	print *, ' Config veget = ', ok_veget
704	print *, ' Sortie journaliere = ', ok_journe
705	print *, ' Sortie mensuelle = ', ok_mensuel
706	print *, ' Sortie instantanee = ', ok_instan
707	print *, ' Sortie bilan d''energie, if_ebil =', if_ebil
708	print *, ' Excentricite = ',R_ecc
709	print *, ' Equinoxe = ',R_peri
710	print *, ' Inclinaison =',R_incl
711	print *, ' Constante solaire =',solaire
712	print *, ' co2_ppm =',co2_ppm
713	print *, ' RCO2 = ',RCO2
714	print *, ' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4
715	print *, ' N2O_ppb =',N2O_ppb,' RN2O = ',RN2O
716	print *, ' CFC11_ppt=',CFC11_ppt,' RCFC11 = ',RCFC11
717	print *, ' CFC12_ppt=',CFC12_ppt,' RCFC12 = ',RCFC12
718	print *, ' epmax = ', epmax
719	print *, ' ok_adj_ema = ', ok_adj_ema
720	print *, ' iflag_clw = ', iflag_clw
721	print *, ' cld_lc_lsc = ', cld_lc_lsc
722	print *, ' cld_lc_con = ', cld_lc_con
723	print *, ' cld_tau_lsc = ', cld_tau_lsc
724	print *, ' cld_tau_con = ', cld_tau_con
725	print *, ' ffallv_lsc = ', ffallv_lsc
726	print *, ' ffallv_con = ', ffallv_con
727	print *, ' coef_eva = ', coef_eva
728	print *, ' reevap_ice = ', reevap_ice
729	print *, ' iflag_pdf = ', iflag_pdf
730	print *, ' iflag_cldcon = ', iflag_cldcon
731	print *, ' fact_cldcon = ', fact_cldcon
732	print *, ' facttemps = ', facttemps
733	print *, ' ok_newmicro = ',ok_newmicro
734	print *, ' ratqsbas = ',ratqsbas
735	print *, ' ratqshaut = ',ratqshaut
736	print *, ' top_height = ',top_height
737	print *, ' overlap = ',overlap
738	print *, ' cdmmax = ',cdmmax
739	print *, ' cdhmax = ',cdhmax
740	print *, ' ksta = ',ksta
741	print *, ' ksta_ter = ',ksta_ter
742	print *, ' ok_kzmin = ',ok_kzmin
743	print *, ' ok_ade = ',ok_ade
744	print *, ' ok_aie = ',ok_aie
745	print *, ' bl95_b0 = ',bl95_b0
746	print *, ' bl95_b1 = ',bl95_b1
747	print *, ' lev_histhf = ',lev_histhf
748	print *, ' lev_histday = ',lev_histday
749	print *, ' lev_histmth = ',lev_histmth
750	print *, ' iflag_pbl = ', iflag_pbl
751	print *, ' iflag_thermals = ', iflag_thermals
752	print *, ' type_run = ',type_run
753	print *, ' ok_isccp = ',ok_isccp
754	print *, ' ok_regdyn = ',ok_regdyn
755	print *, ' lonmin lonmax latmin latmax bilKP_ins =',&
756	lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
757	print *, 'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, &
758	ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
759	print *, ' bug_ozone = ', bug_ozone
760
761	end subroutine conf_phys
762
763	end module conf_phys_m