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

Diff of /trunk/phylmd/clesphys.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 328 by guez, Thu Jun 13 14:40:06 2019 UTC revision 345 by guez, Tue Nov 26 20:47:46 2019 UTC
# Line 4  module clesphys Line 4  module clesphys
4    
5    implicit none    implicit none
6    
7    REAL, protected:: co2_ppm = 348.    real, protected:: solaire = 1365. ! AMIP II, constante solaire en W / m2
8    ! AMIP II, concentration du gaz carbonique en ppmv    double precision, save, protected:: RCO2 ! mass mixing ratio of CO2
9        double precision, save, protected:: RCH4 ! mass mixing ratio of CH4
10    real, protected:: solaire = 1365. ! AMIP II, constante solaire en W/m2    double precision, save, protected:: RN2O ! mass mixing ratio of N2O
11    double precision, save, protected:: RCO2 ! Concentration du CO2    double precision, save, protected:: RCFC11 ! mass mixing ratio of CFC11
12    double precision, save, protected:: RCH4 ! Concentration du CH4    double precision, save, protected:: RCFC12 ! mass mixing ratio of CFC12
   double precision, save, protected:: RN2O ! Concentration du N2O  
   double precision, save, protected:: RCFC11 ! Concentration du CFC11  
   double precision, save, protected:: RCFC12 ! Concentration du CFC12  
   double precision, protected:: CH4_ppb = 1650., N2O_ppb = 306.  
   double precision, protected:: CFC11_ppt = 280., CFC12_ppt = 484.  
13    REAL, protected:: cdmmax = 1.3E-3, cdhmax = 1.1E-3 ! seuils cdrm, cdrh    REAL, protected:: cdmmax = 1.3E-3, cdhmax = 1.1E-3 ! seuils cdrm, cdrh
14        
15    REAL, protected:: ksta = 1e-10, ksta_ter = 1e-10    REAL, protected:: ksta = 1e-10, ksta_ter = 1e-10
# Line 23  module clesphys Line 18  module clesphys
18    LOGICAL, protected:: ok_kzmin  = .true.    LOGICAL, protected:: ok_kzmin  = .true.
19    ! calcul de Kzmin dans la couche limite de surface    ! calcul de Kzmin dans la couche limite de surface
20        
21    INTEGER, save, protected:: ecrit_ins = 1    INTEGER, protected:: ecrit_ins = 1
22    ! number of time steps of physics between outputs    ! number of time steps of physics between outputs
23    
24    logical, protected:: ok_instan = .false.    logical, protected:: ok_instan = .false.
# Line 38  contains Line 33  contains
33    
34      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
35    
36        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      namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &      namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &
41           CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &           CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &
42           ok_instan, f_cdrag_ter, f_cdrag_oce           ok_instan, f_cdrag_ter, f_cdrag_oce
# Line 48  contains Line 47  contains
47      read(unit=*, nml=clesphys_nml)      read(unit=*, nml=clesphys_nml)
48      write(unit_nml, nml=clesphys_nml)      write(unit_nml, nml=clesphys_nml)
49    
50      RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97      RCO2 = co2_ppm * 1e-06 * 44.011 / 28.97
51      RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97      RCH4 = CH4_ppb * 1E-09 * 16.043 / 28.97
52      RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97      RN2O = N2O_ppb * 1E-09 * 44.013 / 28.97
53      RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97      RCFC11=CFC11_ppt* 1E-12 * 137.3686 / 28.97
54      RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97      RCFC12 = CFC12_ppt * 1E-12 * 120.9140 / 28.97
55    
56      print *, ' RCO2 = ', RCO2      print *, ' RCO2 = ', RCO2
57      print *, ' RCH4 = ', RCH4      print *, ' RCH4 = ', RCH4
Legend:
Removed from v.328  
changed lines
  Added in v.345
  ViewVC Help
Powered by ViewVC 1.1.21