libf/phylmd/conf_phys.f90

module conf_phys_m

  implicit none

  integer iflag_pbl

contains

  subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &
       fact_cldcon, facttemps, ok_newmicro, iflag_cldcon, ratqsbas, &
       ratqshaut, if_ebil, ok_ade, ok_aie, bl95_b0, bl95_b1, iflag_thermals, &
       nsplit_thermals)

    ! From phylmd/conf_phys.F90, v 1.7 2005/07/05 07:21:23

    use getincom, only: getin
    use clesphys, only: solaire, co2_ppm, rco2, ch4_ppb, rch4, n2o_ppb, rn2o, &
         cfc11_ppt, rcfc11, cfc12_ppt, rcfc12, top_height, overlap, cdmmax, &
         cdhmax, ksta, ksta_ter, ok_kzmin, lev_histhf, lev_histday, &
         lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, &
         latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, &
         ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone
    use yomcst
    use conema3_m
    use comfisrtilp
    use nuagecom

    !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12

    ! Configuration de la "physique" de LMDZ a l'aide de la fonction
    ! GETIN de IOIPSL

    ! ocean: type d'ocean (force, slab, couple)
    ! ok_veget: type de modele de vegetation
    ! ok_journe: sorties journalieres
    ! ok_mensuel: sorties mensuelles
    ! ok_instan: sorties instantanees
    ! 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, intent(out):: ok_ade ! apply aerosol direct effect
    LOGICAL, intent(out):: ok_aie ! ! apply aerosol indirect effect
    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)

    ok_ade = .false.
    call getin('ok_ade', ok_ade)

    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

    iflag_pbl = 1
    call getin('iflag_pbl', iflag_pbl)

    !Config Key = iflag_thermals
    !Config Desc =
    !Config Def = 0
    !Config Help =

    iflag_thermals = 0
    call getin('iflag_thermals', iflag_thermals)

    !Config Key = nsplit_thermals
    !Config Desc =
    !Config Def = 1
    !Config Help =

    nsplit_thermals = 1
    call getin('nsplit_thermals', nsplit_thermals)

    !Config Key = lev_histhf
    !Config Desc =
    !Config Def = 0
    !Config Help =

    lev_histhf = 0
    call getin('lev_histhf', lev_histhf)

    !Config Key = lev_histday
    !Config Desc =
    !Config Def = 1
    !Config Help =

    lev_histday = 1
    call getin('lev_histday', lev_histday)

    !Config Key = lev_histmth
    !Config Desc =
    !Config Def = 2
    !Config Help =

    lev_histmth = 2
    call getin('lev_histmth', lev_histmth)

    !Config Key = type_run
    !Config Desc =
    !Config Def = 'AMIP' ou 'ENSP'
    !Config Help =

    type_run = 'AMIP'
    call getin('type_run', type_run)

    !Config Key = ok_isccp
    !Config Desc =
    !Config Def = .true.
    !Config Help =

    ! ok_isccp = .true.
    ok_isccp = .false.
    call getin('ok_isccp', ok_isccp)

    !Config Key = ok_regdyn
    !Config Desc =
    !Config Def = 'AMIP'
    !Config Help =

    ! ok_regdyn = .true.
    ok_regdyn = .false.
    call getin('ok_regdyn', ok_regdyn)

    ! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone 
    ! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"

    !Config Key = lonmin_ins
    !Config Desc = 100. 
    !Config Def = longitude minimale sorties "bilKP_ins"
    !Config Help = 

    lonmin_ins = 100.
    call getin('lonmin_ins', lonmin_ins)

    !Config Key = lonmax_ins
    !Config Desc = 130. 
    !Config Def = longitude maximale sorties "bilKP_ins"
    !Config Help =

    lonmax_ins = 130.
    call getin('lonmax_ins', lonmax_ins)

    !Config Key = latmin_ins
    !Config Desc = -20. 
    !Config Def = latitude minimale sorties "bilKP_ins"
    !Config Help = 

    latmin_ins = -20.
    call getin('latmin_ins', latmin_ins)

    !Config Key = latmax_ins
    !Config Desc = 20. 
    !Config Def = latitude maximale sorties "bilKP_ins"
    !Config Help =

    latmax_ins = 20.
    call getin('latmax_ins', latmax_ins)

    !Config Key = ecrit_ins
    !Config Desc =
    !Config Def = NINT(86400./dtime/48.) !a chaque pas de temps physique
    !Config Help =

    ! ecrit_ins = NINT(86400./dtime/48.)
    ecrit_ins = NINT(86400./48.)
    call getin('ecrit_ins', ecrit_ins)

    !Config Key = ecrit_hf
    !Config Desc =
    !Config Def = NINT(86400./dtime *0.25) !toutes les 6h
    !Config Help =

    ! ecrit_hf = NINT(86400./dtime *0.25)
    ecrit_hf = NINT(86400. *0.25)
    call getin('ecrit_hf', ecrit_hf)

    !Config Key = ecrit_hf2mth
    !Config Desc =
    !Config Def = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h
    !Config Help =

    ecrit_hf2mth = 4*30
    call getin('ecrit_hf2mth', ecrit_hf2mth)

    !Config Key = ecrit_day
    ecrit_day = 86400
    call getin('ecrit_day', ecrit_day)

    ecrit_mth = 86400
    call getin('ecrit_mth', ecrit_mth)

    ecrit_tra = 1
    call getin('ecrit_tra', ecrit_tra)

    !Config Key = ecrit_reg
    !Config Desc =
    !Config Def = NINT(86400./dtime *0.25) !4 fois par jour
    !Config Help =

    ! ecrit_reg = NINT(86400./dtime *0.25) !4 fois par jour
    ecrit_reg = NINT(86400. *0.25) !4 fois par jour
    call getin('ecrit_reg', ecrit_reg)

    !Config Key = bug_ozone
    !Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true
    !Config Def = false
    !Config Help =

    bug_ozone = .false.
    call getin('bug_ozone', bug_ozone)

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

  end subroutine conf_phys

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