trunk/phylmd/clesphys.f

module clesphys

  ! From version 1.3 2005/06/06 13:16:33

  implicit none

  REAL:: co2_ppm = 348. ! AMIP II, concentration du gaz carbonique en ppmv
  real:: solaire = 1365. ! AMIP II, constante solaire en W/m2
  double precision, save:: RCO2 ! Concentration du CO2
  double precision, save:: RCH4 ! Concentration du CH4
  double precision, save:: RN2O ! Concentration du N2O
  double precision, save:: RCFC11 ! Concentration du CFC11
  double precision, save:: RCFC12 ! Concentration du CFC12
  double precision:: CH4_ppb = 1650., N2O_ppb = 306.
  double precision:: CFC11_ppt = 280., CFC12_ppt = 484.
  REAL:: cdmmax = 1.3E-3, cdhmax = 1.1E-3 ! seuils cdrm, cdrh
  
  REAL:: ksta = 1.0e-10, ksta_ter = 1.0e-10
  ! param\`etres de stabilit\'e sur terre et en dehors

  LOGICAL:: ok_kzmin  = .true.! calcul de Kzmin dans la couche limite de surface
  INTEGER, save:: ecrit_ins
  logical:: ok_instan = .false. ! sorties instantanees dans le fichier histins
  real:: f_cdrag_ter = 1., f_cdrag_oce = 0.8 ! surface drag coefficients

contains

  subroutine read_clesphys

    use unit_nml_m, only: unit_nml

    namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &
         CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &
         ok_instan, f_cdrag_ter, f_cdrag_oce

    !---------------------------------------------------------

    ecrit_ins = NINT(86400./48.)

    print *, "Enter namelist 'clesphys_nml'."
    read(unit=*, nml=clesphys_nml)
    write(unit_nml, nml=clesphys_nml)

    RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97
    RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
    RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
    RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
    RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97

    print *, ' RCO2 = ', RCO2
    print *, ' RCH4 = ', RCH4
    print *, ' RN2O = ', RN2O
    print *, ' RCFC11 = ', RCFC11
    print *, ' RCFC12 = ', RCFC12

  end subroutine read_clesphys

end module clesphys
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 de stabilit\'e sur terre et en dehors
20
21	LOGICAL:: ok_kzmin = .true.! calcul de Kzmin dans la couche limite de surface
22	INTEGER, save:: ecrit_ins
23	logical:: ok_instan = .false. ! sorties instantanees dans le fichier histins
24	real:: f_cdrag_ter = 1., f_cdrag_oce = 0.8 ! surface drag coefficients
25
26	contains
27
28	subroutine read_clesphys
29
30	use unit_nml_m, only: unit_nml
31
32	namelist /clesphys_nml/ solaire, co2_ppm, CH4_ppb, N2O_ppb, CFC11_ppt, &
33	CFC12_ppt, cdmmax, cdhmax, ksta, ksta_ter, ok_kzmin, ecrit_ins, &
34	ok_instan, f_cdrag_ter, f_cdrag_oce
35
36	!---------------------------------------------------------
37
38	ecrit_ins = NINT(86400./48.)
39
40	print *, "Enter namelist 'clesphys_nml'."
41	read(unit=*, nml=clesphys_nml)
42	write(unit_nml, nml=clesphys_nml)
43
44	RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97
45	RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
46	RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
47	RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
48	RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
49
50	print *, ' RCO2 = ', RCO2
51	print *, ' RCH4 = ', RCH4
52	print *, ' RN2O = ', RN2O
53	print *, ' RCFC11 = ', RCFC11
54	print *, ' RCFC12 = ', RCFC12
55
56	end subroutine read_clesphys
57
58	end module clesphys