trunk/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 IOIPSL, 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 = 6)  :: ocean
    logical              :: ok_veget, ok_newmicro
    logical              :: ok_journe, ok_mensuel, ok_instan        
    LOGICAL              :: ok_ade, ok_aie
    REAL                 :: bl95_b0, bl95_b1
    real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
    integer              :: iflag_cldcon, if_ebil

    ! Local
    integer              :: numout = 6
    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, slab,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)


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

  end subroutine conf_phys

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