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

Contents of /trunk/Sources/phylmd/clesphys.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 212 - (show annotations)
Thu Jan 12 12:31:31 2017 UTC (7 years, 4 months ago) by guez
File size: 1858 byte(s) 
Moved variables from module com_io_dyn to module inithist_m, where
they are defined.

Split grid_atob.f into grille_m.f and dist_sphe.f. Extracted ASCCI art
to documentation. In grille_m, use automatic arrays instead of maximum
size. In grille_m, instead of printing data for every problematic
point, print a single diagnostic message.

Removed variables top_height, overlap, lev_histhf, lev_histday,
lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins,
latmin_ins, latmax_ins of module clesphys, not used.

Removed variable itap of module histwrite_phy_m, not used. There is a
variable itap in module time_phylmdz.

Added output of tro3.

In physiq, no need to compute wo at every time-step, since we only use
it in radlwsw.
1 module clesphys
2
3 ! From version 1.3 2005/06/06 13:16:33
4
5 implicit none
6
7 REAL:: co2_ppm = 348. ! AMIP II, concentration du gaz carbonique en ppmv
8 real:: solaire = 1365. ! AMIP II, constante solaire en W/m2
9 double precision, save:: RCO2 ! Concentration du CO2
10 double precision, save:: RCH4 ! Concentration du CH4
11 double precision, save:: RN2O ! Concentration du N2O
12 double precision, save:: RCFC11 ! Concentration du CFC11
13 double precision, save:: RCFC12 ! Concentration du CFC12
14 double precision:: CH4_ppb = 1650., N2O_ppb = 306.
15 double precision:: CFC11_ppt = 280., CFC12_ppt = 484.
16 REAL:: cdmmax = 1.3E-3, cdhmax = 1.1E-3 ! seuils cdrm, cdrh
17
18 REAL:: ksta = 1.0e-10, ksta_ter = 1.0e-10
19 ! param\`etres stabilit\'e sur terre et en dehors
20
21 LOGICAL:: ok_kzmin = .true.! calcul Kzmin dans la couche limite de surface
22 INTEGER, save:: ecrit_ins
23 logical:: ok_instan = .false. ! sorties instantanees dans le fichier histins
24
25 contains
26
27 subroutine read_clesphys
28
29 use unit_nml_m, only: unit_nml
30
31 namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &
32 CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &
33 ok_instan
34
35 !---------------------------------------------------------
36
37 ecrit_ins = NINT(86400./48.)
38
39 print *, "Enter namelist 'clesphys_nml'."
40 read(unit=*, nml=clesphys_nml)
41 write(unit_nml, nml=clesphys_nml)
42
43 RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97
44 RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
45 RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
46 RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
47 RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
48
49 print *, ' RCO2 = ', RCO2
50 print *, ' RCH4 = ', RCH4
51 print *, ' RN2O = ', RN2O
52 print *, ' RCFC11 = ', RCFC11
53 print *, ' RCFC12 = ', RCFC12
54
55 end subroutine read_clesphys
56
57 end module clesphys
  ViewVC Help
Powered by ViewVC 1.1.21