--- trunk/libf/phylmd/conf_phys.f90 2012/11/14 16:59:30 68 +++ trunk/Sources/phylmd/conf_phys.f 2015/04/29 15:47:56 134 @@ -2,7 +2,8 @@ implicit none - integer iflag_pbl ! parameter for the planetary boundary layer + integer:: iflag_pbl = 1 ! for the planetary boundary layer + REAL:: rad_chau1 = 13., rad_chau2 = 9. contains @@ -12,37 +13,36 @@ ! Configuration de la "physique" de LMDZ. - USE clesphys, ONLY: bug_ozone, cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, & - ch4_ppb, co2_ppm, ecrit_day, ecrit_hf, ecrit_hf2mth, ecrit_ins, & - 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 + USE clesphys, ONLY: cdhmax, cdmmax, cfc11_ppt, cfc12_ppt, ch4_ppb, & + co2_ppm, ecrit_day, ecrit_hf, ecrit_hf2mth, ecrit_ins, 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 + use clesphys2, only: read_clesphys2 USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, & cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice USE conema3_m, ONLY: epmax, iflag_clw, ok_adj_ema - USE nuagecom, ONLY: rad_chau1, rad_chau2, rad_froid use unit_nml_m, only: unit_nml - USE yomcst, ONLY: r_ecc, r_incl, r_peri + USE yomcst, ONLY: read_YOMCST - namelist /conf_phys_nml/ R_ecc, R_peri, R_incl, solaire, co2_ppm, & - CH4_ppb, N2O_ppb, CFC11_ppt, CFC12_ppt, epmax, ok_adj_ema, & - iflag_clw, cld_lc_lsc, cld_lc_con, cld_tau_lsc, cld_tau_con, & - ffallv_lsc, ffallv_con, coef_eva, reevap_ice, iflag_pdf, & - rad_froid, rad_chau1, rad_chau2, top_height, overlap, cdmmax, & - 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 + namelist /conf_phys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, & + CFC12_ppt, epmax, ok_adj_ema, iflag_clw, cld_lc_lsc, cld_lc_con, & + cld_tau_lsc, cld_tau_con, ffallv_lsc, ffallv_con, coef_eva, & + reevap_ice, iflag_pdf, top_height, overlap, cdmmax, 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 + + namelist /nuagecom/ rad_chau1, rad_chau2 !----------------------------------------------------------- print *, "Call sequence information: conf_phys" + call read_clesphys2 + call read_YOMCST - 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. @@ -61,9 +61,6 @@ coef_eva = 2.e-5 reevap_ice = .false. iflag_pdf = 0 - rad_froid = 35.0 - rad_chau1 = 13.0 - rad_chau2 = 9.0 top_height = 3 overlap = 3 cdmmax = 1.3E-3 @@ -71,7 +68,6 @@ ksta = 1.0e-10 ksta_ter = 1.0e-10 ok_kzmin = .true. - iflag_pbl = 1 lev_histhf = 0 lev_histday = 1 lev_histmth = 2 @@ -89,7 +85,6 @@ ecrit_mth = 86400 ecrit_tra = 1 ecrit_reg = NINT(86400. *0.25) ! 4 fois par jour - bug_ozone = .false. print *, "Enter namelist 'conf_phys_nml'." read(unit=*, nml=conf_phys_nml) @@ -107,6 +102,10 @@ print *, ' RCFC11 = ', RCFC11 print *, ' RCFC12 = ', RCFC12 + print *, "Enter namelist 'nuagecom'." + read(unit=*, nml=nuagecom) + write(unit_nml, nml=nuagecom) + end subroutine conf_phys end module conf_phys_m