trunk/phylmd/conf_phys.f

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
    ! 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, intent(out):: ok_ade ! apply aerosol direct effect
    LOGICAL, intent(out):: ok_aie ! ! apply aerosol indirect effect
    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)

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

    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

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