/[lmdze]/trunk/phylmd/conf_phys.f
ViewVC logotype

Diff of /trunk/phylmd/conf_phys.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 62 by guez, Thu Jul 26 14:37:37 2012 UTC revision 68 by guez, Wed Nov 14 16:59:30 2012 UTC
# Line 2  module conf_phys_m Line 2  module conf_phys_m
2    
3    implicit none    implicit none
4    
5    integer iflag_pbl    integer iflag_pbl ! parameter for the planetary boundary layer
6    
7  contains  contains
8    
9    subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, &    subroutine conf_phys
        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  
10    
11      ! Local      ! From phylmd/conf_phys.F90, version 1.7 2005/07/05 07:21:23
     real :: zzz  
12    
13      integer :: iflag_thermals, nsplit_thermals      ! Configuration de la "physique" de LMDZ.
14    
15        USE clesphys, ONLY: bug_ozone, cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, &
16             ch4_ppb, co2_ppm, ecrit_day, ecrit_hf, ecrit_hf2mth, ecrit_ins, &
17             ecrit_mth, ecrit_reg, ecrit_tra, ksta, ksta_ter, latmax_ins, &
18             latmin_ins, lev_histday, lev_histhf, lev_histmth, lonmax_ins, &
19             lonmin_ins, n2o_ppb, ok_isccp, ok_kzmin, ok_regdyn, overlap, rcfc11, &
20             rcfc12, rch4, rco2, rn2o, solaire, top_height, type_run
21        USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
22             cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice
23        USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema
24        USE nuagecom, ONLY: rad_chau1, rad_chau2, rad_froid
25        use unit_nml_m, only: unit_nml
26        USE yomcst, ONLY: r_ecc, r_incl, r_peri
27    
28        namelist /conf_phys_nml/ R_ecc, R_peri, R_incl, solaire, co2_ppm, &
29             CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, epmax, ok_adj_ema, &
30             iflag_clw, cld_lc_lsc, cld_lc_con, cld_tau_lsc, cld_tau_con, &
31             ffallv_lsc, ffallv_con, coef_eva, reevap_ice, iflag_pdf, &
32             rad_froid, rad_chau1, rad_chau2, top_height, overlap, cdmmax, &
33             cdhmax, ksta, ksta_ter, ok_kzmin, iflag_pbl, lev_histhf, &
34             lev_histday, lev_histmth, type_run, ok_isccp, ok_regdyn, &
35             lonmin_ins, lonmax_ins, latmin_ins, latmax_ins, ecrit_ins, &
36             ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_tra, &
37             ecrit_reg, bug_ozone
38    
39      !-----------------------------------------------------------      !-----------------------------------------------------------
40    
41      print *, "Call sequence information: conf_phys"      print *, "Call sequence information: conf_phys"
42    
43      !Config Key = OCEAN      R_ecc = 0.016715 ! AMIP II
44      !Config Desc = Type d'ocean      R_peri = 102.7 ! AMIP II
45      !Config Def = force      R_incl = 23.441 ! AMIP II
46      !Config Help = Type d'ocean utilise: force ou slab mais pas couple      solaire = 1365. ! AMIP II
47        co2_ppm = 348. ! AMIP II
48      ocean = 'force '      CH4_ppb = 1650.
49      call getin('OCEAN', ocean)      N2O_ppb = 306.
50        CFC11_ppt = 280.
51      !Config Key = VEGET      CFC12_ppt = 484.
     !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 =  
   
52      epmax = .993      epmax = .993
     call getin('epmax', epmax)  
   
     !Config Key = ok_adj_ema  
     !Config Desc =  
     !Config Def = false  
     !Config Help =  
   
53      ok_adj_ema = .false.      ok_adj_ema = .false.
     call getin('ok_adj_ema', ok_adj_ema)  
   
     !Config Key = iflag_clw  
     !Config Desc =  
     !Config Def = 0  
     !Config Help =  
   
54      iflag_clw = 0      iflag_clw = 0
     call getin('iflag_clw', iflag_clw)  
   
     !Config Key = cld_lc_lsc  
     !Config Desc =  
     !Config Def = 2.6e-4  
     !Config Help =  
   
55      cld_lc_lsc = 2.6e-4      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 =  
   
56      cld_lc_con = 2.6e-4      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 =  
   
57      cld_tau_lsc = 3600.      cld_tau_lsc = 3600.
     call getin('cld_tau_lsc', cld_tau_lsc)  
   
     !Config Key = cld_tau_con  
     !Config Desc =  
     !Config Def = 3600.  
     !Config Help =  
   
58      cld_tau_con = 3600.      cld_tau_con = 3600.
     call getin('cld_tau_con', cld_tau_con)  
   
     !Config Key = ffallv_lsc  
     !Config Desc =  
     !Config Def = 1.  
     !Config Help =  
   
59      ffallv_lsc = 1.      ffallv_lsc = 1.
     call getin('ffallv_lsc', ffallv_lsc)  
   
     !Config Key = ffallv_con  
     !Config Desc =  
     !Config Def = 1.  
     !Config Help =  
   
60      ffallv_con = 1.      ffallv_con = 1.
     call getin('ffallv_con', ffallv_con)  
   
     !Config Key = coef_eva  
     !Config Desc =  
     !Config Def = 2.e-5  
     !Config Help =  
   
61      coef_eva = 2.e-5      coef_eva = 2.e-5
     call getin('coef_eva', coef_eva)  
   
     !Config Key = reevap_ice  
     !Config Desc =  
     !Config Def = .false.  
     !Config Help =  
   
62      reevap_ice = .false.      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 =  
   
63      iflag_pdf = 0      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 =  
   
64      rad_froid = 35.0      rad_froid = 35.0
     call getin('rad_froid', rad_froid)  
   
     !Config Key = rad_chau1  
     !Config Desc =  
     !Config Def = 13.0  
     !Config Help =  
   
65      rad_chau1 = 13.0      rad_chau1 = 13.0
     call getin('rad_chau1', rad_chau1)  
   
     !Config Key = rad_chau2  
     !Config Desc =  
     !Config Def = 9.0  
     !Config Help =  
   
66      rad_chau2 = 9.0      rad_chau2 = 9.0
     call getin('rad_chau2', rad_chau2)  
   
     !Config Key = top_height  
     !Config Desc =  
     !Config Def = 3  
     !Config Help =  
   
67      top_height = 3      top_height = 3
     call getin('top_height', top_height)  
   
     !Config Key = overlap  
     !Config Desc =  
     !Config Def = 3  
     !Config Help =  
   
68      overlap = 3      overlap = 3
     call getin('overlap', overlap)  
   
     !Config Key = cdmmax  
     !Config Desc =  
     !Config Def = 1.3E-3  
     !Config Help =  
   
69      cdmmax = 1.3E-3      cdmmax = 1.3E-3
     call getin('cdmmax', cdmmax)  
   
     !Config Key = cdhmax  
     !Config Desc =  
     !Config Def = 1.1E-3  
     !Config Help =  
   
70      cdhmax = 1.1E-3      cdhmax = 1.1E-3
     call getin('cdhmax', cdhmax)  
   
     !261103  
   
     !Config Key = ksta  
     !Config Desc =  
     !Config Def = 1.0e-10  
     !Config Help =  
   
71      ksta = 1.0e-10      ksta = 1.0e-10
     call getin('ksta', ksta)  
   
     !Config Key = ksta_ter  
     !Config Desc =  
     !Config Def = 1.0e-10  
     !Config Help =  
   
72      ksta_ter = 1.0e-10      ksta_ter = 1.0e-10
     call getin('ksta_ter', ksta_ter)  
   
     !Config Key = ok_kzmin  
     !Config Desc =  
     !Config Def = .true.  
     !Config Help =  
   
73      ok_kzmin = .true.      ok_kzmin = .true.
     call getin('ok_kzmin', ok_kzmin)  
   
     ! PARAMETER FOR THE PLANETARY BOUNDARY LAYER  
   
74      iflag_pbl = 1      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 =  
   
75      lev_histhf = 0      lev_histhf = 0
     call getin('lev_histhf', lev_histhf)  
   
     !Config Key = lev_histday  
     !Config Desc =  
     !Config Def = 1  
     !Config Help =  
   
76      lev_histday = 1      lev_histday = 1
     call getin('lev_histday', lev_histday)  
   
     !Config Key = lev_histmth  
     !Config Desc =  
     !Config Def = 2  
     !Config Help =  
   
77      lev_histmth = 2      lev_histmth = 2
     call getin('lev_histmth', lev_histmth)  
   
     !Config Key = type_run  
     !Config Desc =  
     !Config Def = 'AMIP' ou 'ENSP'  
     !Config Help =  
   
78      type_run = 'AMIP'      type_run = 'AMIP'
     call getin('type_run', type_run)  
   
     !Config Key = ok_isccp  
     !Config Desc =  
     !Config Def = .true.  
     !Config Help =  
   
     ! ok_isccp = .true.  
79      ok_isccp = .false.      ok_isccp = .false.
     call getin('ok_isccp', ok_isccp)  
   
     !Config Key = ok_regdyn  
     !Config Desc =  
     !Config Def = 'AMIP'  
     !Config Help =  
   
     ! ok_regdyn = .true.  
80      ok_regdyn = .false.      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 =  
   
81      lonmin_ins = 100.      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 =  
   
82      lonmax_ins = 130.      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 =  
   
83      latmin_ins = -20.      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 =  
   
84      latmax_ins = 20.      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.)  
85      ecrit_ins = NINT(86400./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)  
86      ecrit_hf = NINT(86400. *0.25)      ecrit_hf = NINT(86400. *0.25)
87      call getin('ecrit_hf', ecrit_hf)      ecrit_hf2mth = 4*30 ! ecriture mens. a partir de val. inst. toutes les 6h
   
     !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  
88      ecrit_day = 86400      ecrit_day = 86400
     call getin('ecrit_day', ecrit_day)  
   
89      ecrit_mth = 86400      ecrit_mth = 86400
     call getin('ecrit_mth', ecrit_mth)  
   
90      ecrit_tra = 1      ecrit_tra = 1
91      call getin('ecrit_tra', ecrit_tra)      ecrit_reg = NINT(86400. *0.25) ! 4 fois par jour
92        bug_ozone = .false.
93    
94      !Config Key = ecrit_reg      print *, "Enter namelist 'conf_phys_nml'."
95      !Config Desc =      read(unit=*, nml=conf_phys_nml)
96      !Config Def = NINT(86400./dtime *0.25) !4 fois par jour      write(unit_nml, nml=conf_phys_nml)
     !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 =  
97    
98      bug_ozone = .false.      RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97
99      call getin('bug_ozone', bug_ozone)      RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
100        RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
101        RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
102        RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
103    
     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  
104      print *, ' RCO2 = ', RCO2      print *, ' RCO2 = ', RCO2
105      print *, ' CH4_ppb =', CH4_ppb, ' RCH4 = ', RCH4      print *, ' RCH4 = ', RCH4
106      print *, ' N2O_ppb =', N2O_ppb, ' RN2O = ', RN2O      print *, ' RN2O = ', RN2O
107      print *, ' CFC11_ppt=', CFC11_ppt, ' RCFC11 = ', RCFC11      print *, ' RCFC11 = ', RCFC11
108      print *, ' CFC12_ppt=', CFC12_ppt, ' RCFC12 = ', RCFC12      print *, ' 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  
109    
110    end subroutine conf_phys    end subroutine conf_phys
111    
Legend:
Removed from v.62  
changed lines
  Added in v.68
  ViewVC Help
Powered by ViewVC 1.1.21