/[lmdze]/trunk/phylmd/clesphys.f90
ViewVC logotype

Annotation of /trunk/phylmd/clesphys.f90

Parent Directory Parent Directory | Revision Log Revision Log

Revision 346 - (hide annotations)
Mon Dec 9 20:15:29 2019 UTC (4 years, 5 months ago) by guez
File size: 2094 byte(s) 
Rename block to `my_block` in procedure `CLOUDS_GNO` because block is
a Fortran keyword.

Remove computation of palpbla in procedure sw. It was not used nor
output. (Not used nor output either in LMDZ.)

In procedure physiq, define `d_[uv]_con` and add them to `[uv]_seri`
only if `conv_Emanuel`. Thus, we do not need to initialize
`d_[uv]_con` to 0, we do not have to save them and we do not add 0 to
`[uv]_seri`.

In procedure physiq, no need to initialize rnebcon to 0, it is defined
by phyetat0 afterwards.

Check that `iflag_cldcon` is between - 2 and 3.
1 guez 3 module clesphys
2    
3 guez 51 ! From version 1.3 2005/06/06 13:16:33
4 guez 3
5     implicit none
6    
7 guez 346 real, protected:: solaire = 1365. ! constante solaire en W / m2
8 guez 345 double precision, save, protected:: RCO2 ! mass mixing ratio of CO2
9     double precision, save, protected:: RCH4 ! mass mixing ratio of CH4
10     double precision, save, protected:: RN2O ! mass mixing ratio of N2O
11     double precision, save, protected:: RCFC11 ! mass mixing ratio of CFC11
12     double precision, save, protected:: RCFC12 ! mass mixing ratio of CFC12
13 guez 288 REAL, protected:: cdmmax = 1.3E-3, cdhmax = 1.1E-3 ! seuils cdrm, cdrh
14    
15     REAL, protected:: ksta = 1e-10, ksta_ter = 1e-10
16 guez 251 ! param\`etres de stabilit\'e sur terre et en dehors
17 guez 68
18 guez 288 LOGICAL, protected:: ok_kzmin = .true.
19     ! calcul de Kzmin dans la couche limite de surface
20    
21 guez 345 INTEGER, protected:: ecrit_ins = 1
22 guez 288 ! number of time steps of physics between outputs
23 guez 3
24 guez 288 logical, protected:: ok_instan = .false.
25     ! sorties instantanees dans le fichier histins
26    
27     real, protected:: f_cdrag_ter = 1., f_cdrag_oce = 0.8
28     ! surface drag coefficients
29    
30 guez 191 contains
31    
32     subroutine read_clesphys
33    
34     use unit_nml_m, only: unit_nml
35    
36 guez 345 REAL:: co2_ppm = 348. ! number mixing ratio, in ppm
37     real:: CH4_ppb = 1650., N2O_ppb = 306.
38     real:: CFC11_ppt = 280., CFC12_ppt = 484.
39    
40 guez 191 namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &
41 guez 212 CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &
42 guez 273 ok_instan, f_cdrag_ter, f_cdrag_oce
43 guez 191
44     !---------------------------------------------------------
45    
46     print *, "Enter namelist 'clesphys_nml'."
47     read(unit=*, nml=clesphys_nml)
48     write(unit_nml, nml=clesphys_nml)
49    
50 guez 345 RCO2 = co2_ppm * 1e-06 * 44.011 / 28.97
51     RCH4 = CH4_ppb * 1E-09 * 16.043 / 28.97
52     RN2O = N2O_ppb * 1E-09 * 44.013 / 28.97
53     RCFC11=CFC11_ppt* 1E-12 * 137.3686 / 28.97
54     RCFC12 = CFC12_ppt * 1E-12 * 120.9140 / 28.97
55 guez 191
56     print *, ' RCO2 = ', RCO2
57     print *, ' RCH4 = ', RCH4
58     print *, ' RN2O = ', RN2O
59     print *, ' RCFC11 = ', RCFC11
60     print *, ' RCFC12 = ', RCFC12
61    
62     end subroutine read_clesphys
63    
64 guez 3 end module clesphys
  ViewVC Help
Powered by ViewVC 1.1.21