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