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

Diff of /trunk/Sources/phylmd/conf_phys.f

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

trunk/libf/phylmd/conf_phys.f90 revision 68 by guez, Wed Nov 14 16:59:30 2012 UTC trunk/Sources/phylmd/conf_phys.f revision 191 by guez, Mon May 9 19:56:28 2016 UTC
# Line 2  module conf_phys_m Line 2  module conf_phys_m
2    
3    implicit none    implicit none
4    
5    integer iflag_pbl ! parameter for the planetary boundary layer    integer:: iflag_pbl = 1 ! for the planetary boundary layer
6      REAL:: rad_chau1 = 13., rad_chau2 = 9.
7    
8  contains  contains
9    
# Line 12  contains Line 13  contains
13    
14      ! Configuration de la "physique" de LMDZ.      ! Configuration de la "physique" de LMDZ.
15    
16      USE clesphys, ONLY: bug_ozone, cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, &      USE clesphys, ONLY: read_clesphys
17           ch4_ppb, co2_ppm, ecrit_day, ecrit_hf, ecrit_hf2mth, ecrit_ins, &      use clesphys2, only: read_clesphys2
          ecrit_mth, ecrit_reg, ecrit_tra, ksta, ksta_ter, latmax_ins, &  
          latmin_ins, lev_histday, lev_histhf, lev_histmth, lonmax_ins, &  
          lonmin_ins, n2o_ppb, ok_isccp, ok_kzmin, ok_regdyn, overlap, rcfc11, &  
          rcfc12, rch4, rco2, rn2o, solaire, top_height, type_run  
18      USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &      USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
19           cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice           cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice
20      USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema      USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema
     USE nuagecom, ONLY: rad_chau1, rad_chau2, rad_froid  
21      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
22      USE yomcst, ONLY: r_ecc, r_incl, r_peri      USE yomcst, ONLY: read_YOMCST
23    
24      namelist /conf_phys_nml/ R_ecc, R_peri, R_incl, solaire, co2_ppm, &      namelist /conf_phys_nml/ epmax, ok_adj_ema, iflag_clw, cld_lc_lsc, &
25           CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, epmax, ok_adj_ema, &           cld_lc_con, cld_tau_lsc, cld_tau_con, ffallv_lsc, ffallv_con, &
26           iflag_clw, cld_lc_lsc, cld_lc_con, cld_tau_lsc, cld_tau_con, &           coef_eva, reevap_ice, iflag_pdf, iflag_pbl
27           ffallv_lsc, ffallv_con, coef_eva, reevap_ice, iflag_pdf, &  
28           rad_froid, rad_chau1, rad_chau2, top_height, overlap, cdmmax, &      namelist /nuagecom/ rad_chau1, rad_chau2
          cdhmax, ksta, ksta_ter, ok_kzmin, iflag_pbl, 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  
29    
30      !-----------------------------------------------------------      !-----------------------------------------------------------
31    
32      print *, "Call sequence information: conf_phys"      print *, "Call sequence information: conf_phys"
33        call read_clesphys2
34        call read_YOMCST
35    
     R_ecc = 0.016715 ! AMIP II  
     R_peri = 102.7 ! AMIP II  
     R_incl = 23.441 ! AMIP II  
     solaire = 1365. ! AMIP II  
     co2_ppm = 348. ! AMIP II  
     CH4_ppb = 1650.  
     N2O_ppb = 306.  
     CFC11_ppt = 280.  
     CFC12_ppt = 484.  
36      epmax = .993      epmax = .993
37      ok_adj_ema = .false.      ok_adj_ema = .false.
38      iflag_clw = 0      iflag_clw = 0
# Line 61  contains Line 45  contains
45      coef_eva = 2.e-5      coef_eva = 2.e-5
46      reevap_ice = .false.      reevap_ice = .false.
47      iflag_pdf = 0      iflag_pdf = 0
     rad_froid = 35.0  
     rad_chau1 = 13.0  
     rad_chau2 = 9.0  
     top_height = 3  
     overlap = 3  
     cdmmax = 1.3E-3  
     cdhmax = 1.1E-3  
     ksta = 1.0e-10  
     ksta_ter = 1.0e-10  
     ok_kzmin = .true.  
     iflag_pbl = 1  
     lev_histhf = 0  
     lev_histday = 1  
     lev_histmth = 2  
     type_run = 'AMIP'  
     ok_isccp = .false.  
     ok_regdyn = .false.  
     lonmin_ins = 100.  
     lonmax_ins = 130.  
     latmin_ins = -20.  
     latmax_ins = 20.  
     ecrit_ins = NINT(86400./48.)  
     ecrit_hf = NINT(86400. *0.25)  
     ecrit_hf2mth = 4*30 ! ecriture mens. a partir de val. inst. toutes les 6h  
     ecrit_day = 86400  
     ecrit_mth = 86400  
     ecrit_tra = 1  
     ecrit_reg = NINT(86400. *0.25) ! 4 fois par jour  
     bug_ozone = .false.  
48    
49      print *, "Enter namelist 'conf_phys_nml'."      print *, "Enter namelist 'conf_phys_nml'."
50      read(unit=*, nml=conf_phys_nml)      read(unit=*, nml=conf_phys_nml)
51      write(unit_nml, nml=conf_phys_nml)      write(unit_nml, nml=conf_phys_nml)
52    
53      RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97      call read_clesphys
54      RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97  
55      RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97      print *, "Enter namelist 'nuagecom'."
56      RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97      read(unit=*, nml=nuagecom)
57      RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97      write(unit_nml, nml=nuagecom)
   
     print *, ' RCO2 = ', RCO2  
     print *, ' RCH4 = ', RCH4  
     print *, ' RN2O = ', RN2O  
     print *, ' RCFC11 = ', RCFC11  
     print *, ' RCFC12 = ', RCFC12  
58    
59    end subroutine conf_phys    end subroutine conf_phys
60    
Legend:
Removed from v.68  
changed lines
  Added in v.191
  ViewVC Help
Powered by ViewVC 1.1.21