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

Annotation of /trunk/phylmd/clesphys.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 288 - (hide annotations)
Tue Jul 24 16:27:12 2018 UTC (5 years, 10 months ago) by guez
Original Path: trunk/phylmd/clesphys.f
File size: 2163 byte(s)
Remove tests on richum, tvirtu and opt_ec in procedure coefkz (not
used in LMDZ either).

Change the meaning of variable ecrit_ins of module clesphys from
number of seconds (integer, weird), which was modified in physiq to
make a multiple of the time step of physics, to number of physics time
steps.

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

  ViewVC Help
Powered by ViewVC 1.1.21