/[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 228 by guez, Fri Nov 3 12:38:47 2017 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      real:: epmax = 0.993 ! \'efficacit\'e de pr\'ecipitation
8      integer:: iflag_clw = 0
9    
10  contains  contains
11    
# Line 12  contains Line 15  contains
15    
16      ! Configuration de la "physique" de LMDZ.      ! Configuration de la "physique" de LMDZ.
17    
18      USE clesphys, ONLY: bug_ozone, cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, &      USE clesphys, ONLY: read_clesphys
19           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  
20      USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &      USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
21           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
22      USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema      use nr_util, only: assert
     USE nuagecom, ONLY: rad_chau1, rad_chau2, rad_froid  
23      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
24      USE yomcst, ONLY: r_ecc, r_incl, r_peri      USE yomcst, ONLY: read_YOMCST
25    
26      namelist /conf_phys_nml/ R_ecc, R_peri, R_incl, solaire, co2_ppm, &      namelist /conf_phys_nml/ epmax, iflag_clw, cld_lc_lsc, cld_lc_con, &
27           CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, epmax, ok_adj_ema, &           cld_tau_lsc, cld_tau_con, ffallv_lsc, ffallv_con, coef_eva, &
28           iflag_clw, cld_lc_lsc, cld_lc_con, cld_tau_lsc, cld_tau_con, &           reevap_ice, iflag_pdf, iflag_pbl
29           ffallv_lsc, ffallv_con, coef_eva, reevap_ice, iflag_pdf, &  
30           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  
31    
32      !-----------------------------------------------------------      !-----------------------------------------------------------
33    
34      print *, "Call sequence information: conf_phys"      print *, "Call sequence information: conf_phys"
35        call read_clesphys2
36        call read_YOMCST
37    
     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.  
     epmax = .993  
     ok_adj_ema = .false.  
     iflag_clw = 0  
38      cld_lc_lsc = 2.6e-4      cld_lc_lsc = 2.6e-4
39      cld_lc_con = 2.6e-4      cld_lc_con = 2.6e-4
40      cld_tau_lsc = 3600.      cld_tau_lsc = 3600.
# Line 61  contains Line 44  contains
44      coef_eva = 2.e-5      coef_eva = 2.e-5
45      reevap_ice = .false.      reevap_ice = .false.
46      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.  
47    
48      print *, "Enter namelist 'conf_phys_nml'."      print *, "Enter namelist 'conf_phys_nml'."
49      read(unit=*, nml=conf_phys_nml)      read(unit=*, nml=conf_phys_nml)
50      write(unit_nml, nml=conf_phys_nml)      write(unit_nml, nml=conf_phys_nml)
51    
52      RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97      call assert(any(iflag_pbl == [0, 1, 6, 8, 9]), &
53      RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97           "conf_phys: bad value for iflag_pbl")
54      RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97      call read_clesphys
55      RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97  
56      RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97      print *, "Enter namelist 'nuagecom'."
57        read(unit=*, nml=nuagecom)
58      print *, ' RCO2 = ', RCO2      write(unit_nml, nml=nuagecom)
     print *, ' RCH4 = ', RCH4  
     print *, ' RN2O = ', RN2O  
     print *, ' RCFC11 = ', RCFC11  
     print *, ' RCFC12 = ', RCFC12  
59    
60    end subroutine conf_phys    end subroutine conf_phys
61    
Legend:
Removed from v.68  
changed lines
  Added in v.228
  ViewVC Help
Powered by ViewVC 1.1.21