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

    ! Local
    real                 :: zzz

    integer :: iflag_thermals,nsplit_thermals

    !-----------------------------------------------------------

    print *, "Call sequence information: conf_phys"

    !Config Key  = OCEAN 
    !Config Desc = Type d'ocean
    !Config Def  = force
    !Config Help = Type d'ocean utilise: force ou slab mais pas couple
    !
    ocean = 'force '
    call getin('OCEAN', ocean)
    !
    !Config Key  = VEGET 
    !Config Desc = Type de modele de vegetation
    !Config Def  = .false.
    !Config Help = Type de modele de vegetation utilise
    !
    ok_veget = .false.
    call getin('VEGET', ok_veget)
    !
    !Config Key  = OK_journe
    !Config Desc = Pour des sorties journalieres 
    !Config Def  = .false.
    !Config Help = Pour creer le fichier histday contenant les sorties
    !              journalieres 
    !
    ok_journe = .false.
    call getin('OK_journe', ok_journe)
    !
    !Config Key  = OK_mensuel
    !Config Desc = Pour des sorties mensuelles 
    !Config Def  = .true.
    !Config Help = Pour creer le fichier histmth contenant les sorties
    !              mensuelles 
    !
    ok_mensuel = .true.
    call getin('OK_mensuel', ok_mensuel)
    !
    !Config Key  = OK_instan
    !Config Desc = Pour des sorties instantanees 
    !Config Def  = .false.
    !Config Help = Pour creer le fichier histins contenant les sorties
    !              instantanees 
    !
    ok_instan = .false.
    call getin('OK_instan', ok_instan)
    !
    !Config Key  = ok_ade
    !Config Desc = Aerosol direct effect or not?
    !Config Def  = .false.
    !Config Help = Used in radlwsw.F
    !
    ok_ade = .false.
    call getin('ok_ade', ok_ade)

    !
    !Config Key  = ok_aie
    !Config Desc = Aerosol indirect effect or not?
    !Config Def  = .false.
    !Config Help = Used in nuage.F and radlwsw.F
    !
    ok_aie = .false.
    call getin('ok_aie', ok_aie)

    !
    !Config Key  = bl95_b0
    !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
    !Config Def  = .false.
    !Config Help = Used in nuage.F
    !
    bl95_b0 = 2.
    call getin('bl95_b0', bl95_b0)

    !Config Key  = bl95_b1
    !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
    !Config Def  = .false.
    !Config Help = Used in nuage.F
    !
    bl95_b1 = 0.2
    call getin('bl95_b1', bl95_b1)

    !
    !
    !Config Key  = if_ebil
    !Config Desc = Niveau de sortie pour les diags bilan d'energie 
    !Config Def  = 0
    !Config Help = 
    !               
    !
    if_ebil = 0
    call getin('if_ebil', if_ebil)
    !!
    !! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
    !!
    !Config Key  = R_ecc
    !Config Desc = Excentricite
    !Config Def  = 0.016715
    !Config Help = 
    !               
    !valeur AMIP II
    R_ecc = 0.016715
    call getin('R_ecc', R_ecc)
    !!
    !Config Key  = R_peri
    !Config Desc = Equinoxe
    !Config Def  = 
    !Config Help = 
    !               
    !
    !valeur AMIP II
    R_peri = 102.7
    call getin('R_peri', R_peri)
    !!
    !Config Key  = R_incl
    !Config Desc = Inclinaison
    !Config Def  = 
    !Config Help = 
    !               
    !
    !valeur AMIP II
    R_incl = 23.441
    call getin('R_incl', R_incl)
    !!
    !Config Key  = solaire
    !Config Desc = Constante solaire en W/m2
    !Config Def  = 1365.
    !Config Help = 
    !               
    !
    !valeur AMIP II
    solaire = 1365.
    call getin('solaire', solaire)
    !!
    !Config Key  = co2_ppm
    !Config Desc = concentration du gaz carbonique en ppmv
    !Config Def  = 348.
    !Config Help = 
    !               
    !
    !valeur AMIP II
    co2_ppm = 348.
    call getin('co2_ppm', co2_ppm)
    !!
    !Config Key  = RCO2
    !Config Desc = Concentration du CO2
    !Config Def  = co2_ppm * 1.0e-06  * 44.011/28.97
    !Config Def  = 348. * 1.0e-06  * 44.011/28.97
    !Config Help = 
    !               
    ! RCO2 = 5.286789092164308E-04
    !ancienne valeur
    RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 ! pour co2_ppm=348.

    !!  call getin('RCO2', RCO2)
    !!
    !Config Key  = RCH4
    !Config Desc = Concentration du CH4
    !Config Def  = 1.65E-06* 16.043/28.97
    !Config Help = 
    !               
    !
    !valeur AMIP II
    !OK  RCH4 = 1.65E-06* 16.043/28.97
    ! RCH4 = 9.137366240938903E-07
    !
    !ancienne valeur
    ! RCH4 = 1.72E-06* 16.043/28.97
    !OK call getin('RCH4', RCH4)
    zzz = 1650.
    call getin('CH4_ppb', zzz)
    CH4_ppb = zzz
    RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
    !!
    !Config Key  = RN2O
    !Config Desc = Concentration du N2O
    !Config Def  = 306.E-09* 44.013/28.97
    !Config Help = 
    !               
    !
    !valeur AMIP II
    !OK  RN2O = 306.E-09* 44.013/28.97
    ! RN2O = 4.648939592682085E-07
    !
    !ancienne valeur
    ! RN2O = 310.E-09* 44.013/28.97
    !OK  call getin('RN2O', RN2O)
    zzz=306.
    call getin('N2O_ppb', zzz)
    N2O_ppb = zzz
    RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
    !!
    !Config Key  = RCFC11
    !Config Desc = Concentration du CFC11
    !Config Def  = 280.E-12* 137.3686/28.97
    !Config Help = 
    !               
    !
    !OK RCFC11 = 280.E-12* 137.3686/28.97
    zzz = 280.
    call getin('CFC11_ppt',zzz)
    CFC11_ppt = zzz
    RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
    ! RCFC11 = 1.327690990680013E-09
    !OK call getin('RCFC11', RCFC11)
    !!
    !Config Key  = RCFC12
    !Config Desc = Concentration du CFC12
    !Config Def  = 484.E-12* 120.9140/28.97
    !Config Help = 
    !               
    !
    !OK RCFC12 = 484.E-12* 120.9140/28.97
    zzz = 484.
    call getin('CFC12_ppt',zzz)
    CFC12_ppt = zzz
    RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
    ! RCFC12 = 2.020102726958923E-09
    !OK call getin('RCFC12', RCFC12)
    !!
    !! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
    !!
    !! KE
    !
    !Config Key  = epmax
    !Config Desc = Efficacite precip
    !Config Def  = 0.993
    !Config Help = 
    !
    epmax = .993
    call getin('epmax', epmax)
    !
    !Config Key  = ok_adj_ema
    !Config Desc =  
    !Config Def  = false
    !Config Help = 
    !
    ok_adj_ema = .false.
    call getin('ok_adj_ema',ok_adj_ema)
    !
    !Config Key  = iflag_clw
    !Config Desc =  
    !Config Def  = 0
    !Config Help = 
    !
    iflag_clw = 0
    call getin('iflag_clw',iflag_clw)
    !
    !Config Key  = cld_lc_lsc 
    !Config Desc =  
    !Config Def  = 2.6e-4
    !Config Help = 
    !
    cld_lc_lsc = 2.6e-4
    call getin('cld_lc_lsc',cld_lc_lsc)
    !
    !Config Key  = cld_lc_con
    !Config Desc =  
    !Config Def  = 2.6e-4
    !Config Help = 
    !
    cld_lc_con = 2.6e-4
    call getin('cld_lc_con',cld_lc_con)
    !
    !Config Key  = cld_tau_lsc
    !Config Desc =  
    !Config Def  = 3600.
    !Config Help = 
    !
    cld_tau_lsc = 3600.
    call getin('cld_tau_lsc',cld_tau_lsc)
    !
    !Config Key  = cld_tau_con
    !Config Desc =  
    !Config Def  = 3600.
    !Config Help = 
    !
    cld_tau_con = 3600.
    call getin('cld_tau_con',cld_tau_con)
    !
    !Config Key  = ffallv_lsc
    !Config Desc =  
    !Config Def  = 1.
    !Config Help = 
    !
    ffallv_lsc = 1.
    call getin('ffallv_lsc',ffallv_lsc)
    !
    !Config Key  = ffallv_con
    !Config Desc =  
    !Config Def  = 1.
    !Config Help = 
    !
    ffallv_con = 1.
    call getin('ffallv_con',ffallv_con)
    !
    !Config Key  = coef_eva
    !Config Desc =  
    !Config Def  = 2.e-5
    !Config Help = 
    !
    coef_eva = 2.e-5
    call getin('coef_eva',coef_eva)
    !
    !Config Key  = reevap_ice
    !Config Desc =  
    !Config Def  = .false.
    !Config Help = 
    !
    reevap_ice = .false.
    call getin('reevap_ice',reevap_ice)
    !
    !Config Key  = iflag_cldcon 
    !Config Desc =  
    !Config Def  = 1
    !Config Help = 
    !
    iflag_cldcon = 1
    call getin('iflag_cldcon',iflag_cldcon)

    !
    !Config Key  = iflag_pdf 
    !Config Desc =  
    !Config Def  = 0
    !Config Help = 
    !
    iflag_pdf = 0
    call getin('iflag_pdf',iflag_pdf)
    !
    !Config Key  = fact_cldcon
    !Config Desc =  
    !Config Def  = 0.375
    !Config Help = 
    !
    fact_cldcon = 0.375
    call getin('fact_cldcon',fact_cldcon)

    !
    !Config Key  = facttemps
    !Config Desc =  
    !Config Def  = 1.e-4
    !Config Help = 
    !
    facttemps = 1.e-4
    call getin('facttemps',facttemps)

    !
    !Config Key  = ok_newmicro
    !Config Desc =  
    !Config Def  = .true.
    !Config Help = 
    !
    ok_newmicro = .true.
    call getin('ok_newmicro',ok_newmicro)
    !
    !Config Key  = ratqsbas
    !Config Desc =  
    !Config Def  = 0.01
    !Config Help = 
    !
    ratqsbas = 0.01
    call getin('ratqsbas',ratqsbas)
    !
    !Config Key  = ratqshaut
    !Config Desc =  
    !Config Def  = 0.3
    !Config Help = 
    !
    ratqshaut = 0.3
    call getin('ratqshaut',ratqshaut)

    !
    !Config Key  = rad_froid
    !Config Desc =  
    !Config Def  = 35.0
    !Config Help = 
    !
    rad_froid = 35.0
    call getin('rad_froid',rad_froid)

    !
    !Config Key  = rad_chau1
    !Config Desc =  
    !Config Def  = 13.0
    !Config Help = 
    !
    rad_chau1 = 13.0
    call getin('rad_chau1',rad_chau1)

    !
    !Config Key  = rad_chau2
    !Config Desc =  
    !Config Def  = 9.0
    !Config Help = 
    !
    rad_chau2 = 9.0
    call getin('rad_chau2',rad_chau2)

    !
    !Config Key  = top_height
    !Config Desc =
    !Config Def  = 3
    !Config Help =
    !
    top_height = 3
    call getin('top_height',top_height)

    !
    !Config Key  = overlap
    !Config Desc =
    !Config Def  = 3
    !Config Help =
    !
    overlap = 3
    call getin('overlap',overlap)


    !
    !
    !Config Key  = cdmmax
    !Config Desc =
    !Config Def  = 1.3E-3
    !Config Help =
    !
    cdmmax = 1.3E-3
    call getin('cdmmax',cdmmax)

    !
    !Config Key  = cdhmax
    !Config Desc =
    !Config Def  = 1.1E-3
    !Config Help =
    !
    cdhmax = 1.1E-3
    call getin('cdhmax',cdhmax)

    !261103
    !
    !Config Key  = ksta
    !Config Desc =
    !Config Def  = 1.0e-10
    !Config Help =
    !
    ksta = 1.0e-10
    call getin('ksta',ksta)

    !
    !Config Key  = ksta_ter
    !Config Desc =
    !Config Def  = 1.0e-10
    !Config Help =
    !
    ksta_ter = 1.0e-10
    call getin('ksta_ter',ksta_ter)

    !
    !Config Key  = ok_kzmin
    !Config Desc =
    !Config Def  = .true.
    !Config Help =
    !
    ok_kzmin = .true.
    call getin('ok_kzmin',ok_kzmin)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    ! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    !Config Key  = iflag_pbl
    !Config Desc =
    !Config Def  = 1
    !Config Help =
    !
    iflag_pbl = 1
    call getin('iflag_pbl',iflag_pbl)
    !
    !Config Key  = iflag_thermals
    !Config Desc =
    !Config Def  = 0
    !Config Help =
    !
    iflag_thermals = 0
    call getin('iflag_thermals',iflag_thermals)
    !
    !
    !Config Key  = nsplit_thermals
    !Config Desc =
    !Config Def  = 1
    !Config Help =
    !
    nsplit_thermals = 1
    call getin('nsplit_thermals',nsplit_thermals)


    !
    !Config Key  = lev_histhf
    !Config Desc =
    !Config Def  = 0
    !Config Help =
    !
    lev_histhf = 0
    call getin('lev_histhf',lev_histhf)

    !
    !Config Key  = lev_histday
    !Config Desc =
    !Config Def  = 1
    !Config Help =
    !
    lev_histday = 1
    call getin('lev_histday',lev_histday)

    !
    !Config Key  = lev_histmth
    !Config Desc =
    !Config Def  = 2
    !Config Help =
    !
    lev_histmth = 2
    call getin('lev_histmth',lev_histmth)

    !
    !Config Key  = type_run
    !Config Desc =
    !Config Def  = 'AMIP' ou 'ENSP'
    !Config Help =
    !
    type_run = 'AMIP'
    call getin('type_run',type_run)

    !
    !Config Key  = ok_isccp
    !Config Desc =
    !Config Def  = .true.
    !Config Help =
    !
    ! ok_isccp = .true.
    ok_isccp = .false.
    call getin('ok_isccp',ok_isccp)

    !
    !
    !Config Key  = ok_regdyn
    !Config Desc =
    !Config Def  = 'AMIP'
    !Config Help =
    !
    ! ok_regdyn = .true.
    ok_regdyn = .false.
    call getin('ok_regdyn',ok_regdyn)
    !
    ! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone 
    ! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
    !
    !Config Key  = lonmin_ins
    !Config Desc = 100.  
    !Config Def  = longitude minimale sorties "bilKP_ins"
    !Config Help = 
    !
    lonmin_ins = 100.
    call getin('lonmin_ins',lonmin_ins)
    !
    !Config Key  = lonmax_ins
    !Config Desc = 130. 
    !Config Def  = longitude maximale sorties "bilKP_ins"
    !Config Help =
    !
    lonmax_ins = 130.
    call getin('lonmax_ins',lonmax_ins)
    !
    !Config Key  = latmin_ins
    !Config Desc = -20.  
    !Config Def  = latitude minimale sorties "bilKP_ins"
    !Config Help = 
    !
    latmin_ins = -20.
    call getin('latmin_ins',latmin_ins)
    !
    !Config Key  = latmax_ins
    !Config Desc = 20. 
    !Config Def  = latitude maximale sorties "bilKP_ins"
    !Config Help =
    !
    latmax_ins = 20.
    call getin('latmax_ins',latmax_ins)
    !
    !Config Key  = ecrit_ins
    !Config Desc =
    !Config Def  = NINT(86400./dtime/48.) !a chaque pas de temps physique
    !Config Help =
    !
    ! ecrit_ins = NINT(86400./dtime/48.)
    ecrit_ins = NINT(86400./48.)
    call getin('ecrit_ins',ecrit_ins)
    !
    !Config Key  = ecrit_hf
    !Config Desc =
    !Config Def  = NINT(86400./dtime *0.25) !toutes les 6h
    !Config Help =
    !
    ! ecrit_hf = NINT(86400./dtime *0.25)
    ecrit_hf = NINT(86400. *0.25)
    call getin('ecrit_hf',ecrit_hf)
    !
    !Config Key  = ecrit_hf2mth
    !Config Desc =
    !Config Def  = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h
    !Config Help =
    !
    ecrit_hf2mth = 4*30
    call getin('ecrit_hf2mth',ecrit_hf2mth)
    !
    !Config Key  = ecrit_day
    ecrit_day = 86400
    call getin('ecrit_day',ecrit_day)
    !
    ecrit_mth = 86400
    call getin('ecrit_mth',ecrit_mth)
    !
    ecrit_tra = 1
    call getin('ecrit_tra',ecrit_tra)
    !
    !Config Key  = ecrit_reg
    !Config Desc =
    !Config Def  = NINT(86400./dtime *0.25)  !4 fois par jour
    !Config Help =
    !
    ! ecrit_reg = NINT(86400./dtime *0.25)  !4 fois par jour
    ecrit_reg = NINT(86400. *0.25)  !4 fois par jour
    call getin('ecrit_reg',ecrit_reg)
    !
    !
    !
    !Config Key  = bug_ozone
    !Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true
    !Config Def  = false
    !Config Help =
    !
    bug_ozone = .false.
    call getin('bug_ozone',bug_ozone)

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

  end subroutine conf_phys

end module conf_phys_m
1	guez	3	module conf_phys_m
2
3			! 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	guez	12	character(len=*), intent(out):: ocean
47	guez	3	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	guez	17	real, intent(out):: fact_cldcon
52			real, intent(out):: facttemps
53			real ratqsbas,ratqshaut
54	guez	3	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	guez	12	!Config Help = Type d'ocean utilise: force ou slab mais pas couple
69	guez	3	!
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	guez	17	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	guez	3	lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
759	guez	17	print *, 'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, &
760	guez	3	ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra
761	guez	17	print *, ' bug_ozone = ', bug_ozone
762	guez	3
763			end subroutine conf_phys
764
765			end module conf_phys_m