Sources/phylmd/conf_phys.f

module conf_phys_m

  ! This module is clean: no C preprocessor directive, no include line.

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