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
    ! 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	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			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	guez	30	use getincom, only: getin
19	guez	3	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	guez	38	use yomcst
26	guez	3	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	guez	12	character(len=*), intent(out):: ocean
45	guez	3	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	guez	17	real, intent(out):: fact_cldcon
50			real, intent(out):: facttemps
51			real ratqsbas,ratqshaut
52	guez	3	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	guez	12	!Config Help = Type d'ocean utilise: force ou slab mais pas couple
67	guez	3	!
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	guez	17	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	guez	3	lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
757	guez	17	print *, 'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, &
758	guez	3	ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
759	guez	17	print *, ' bug_ozone = ', bug_ozone
760	guez	3
761			end subroutine conf_phys
762
763			end module conf_phys_m