source: LMDZ6/branches/IPSLCM6.0.13/libf/phylmd/rrtm/surdi15.F90 @ 3012

Last change on this file since 3012 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM


Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 3.5 KB
Line 
1!OPTIONS XOPT(NOEVAL)
2SUBROUTINE SURDI15
3
4!**** *SURDI15*   - INITIALIZE COMMON YOMRDI15 CONTROLLING RADINT
5!****               FROZEN VERSION (CYCLE 15) OF SURDI
6
7!     PURPOSE.
8!     --------
9!           INITIALIZE YOMRDI15, THE COMMON THAT CONTROLS THE
10!           RADIATION INTERFACE
11
12!**   INTERFACE.
13!     ----------
14!        CALL *SURDI15* FROM *SUECRAD*
15!              -------        -------
16
17!        EXPLICIT ARGUMENTS :
18!        --------------------
19!        NONE
20
21!        IMPLICIT ARGUMENTS :
22!        --------------------
23!        COMMON YOMRDI15
24
25!     METHOD.
26!     -------
27!        SEE DOCUMENTATION
28
29!     EXTERNALS.
30!     ----------
31!        NONE
32
33!     REFERENCE.
34!     ----------
35!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
36!     "IN CORE MODEL"
37
38!     AUTHOR.
39!     -------
40!        96-11: Ph. Dandin. Meteo-France
41!        ORIGINAL : 88-12-15 BY JEAN-JACQUES MORCRETTE  *ECMWF*
42
43!     MODIFICATIONS.
44!     --------------
45!        M.Hamrud      01-Oct-2003 CY28 Cleaning
46!        A.Alias       05-12-2005 greenhouse gases variables (M.Deque)
47!        A.Alias       13-06-2006 RI0 value can be changed via namscen.h
48
49!     ------------------------------------------------------------------
50
51USE PARKIND1        ,ONLY : JPIM     ,JPRB
52USE YOMHOOK         ,ONLY : LHOOK,   DR_HOOK
53! Ce qui concerne NULNAM commente par MPL le 15.04.09
54!USE YOMLUN          , ONLY : NULNAM
55USE YOMLUN_IFSAUX   , ONLY : NULOUT
56
57USE YOMCST          , ONLY : RI0
58USE YOMRDI15        , ONLY : RSDTSN15 ,RRAE15   ,RMU0015  ,RALBICE15,&
59 & RALBSEA15,RALBSNM15,RALBSNO15,RCARDI15 ,REMISS15 ,&
60 & RSNOWAL15,RVLBDC15 ,RCH415   ,RN2O15   ,RCFC1115 ,&
61 & RCFC1215 ,REPALB15 ,REPCLC15 ,REPH2O15 
62USE YOMRDU15        , ONLY : REPSEC15
63
64IMPLICIT NONE
65
66REAL(KIND=JPRB) :: XCARDI, XCFC11, XCFC12, XCH4, XN2O
67REAL(KIND=JPRB) :: ZAIRMWG, ZC11MWG, ZC12MWG, ZCH4MWG, ZCO2MWG, ZN2OMWG, ZSUPSAT
68REAL(KIND=JPRB) :: ZHOOK_HANDLE
69
70#include "namscen.h"
71#include "posnam.intfb.h"
72!      ----------------------------------------------------------------
73
74!*       1.    SET DEFAULT VALUES.
75!              -------------------
76
77IF (LHOOK) CALL DR_HOOK('SURDI15',0,ZHOOK_HANDLE)
78RSDTSN15  = 5.0_JPRB
79
80RRAE15 = 0.1277E-02_JPRB
81RMU0015 = RRAE15/SQRT(RRAE15*(RRAE15+2.0_JPRB))
82
83RALBICE15 = 0.55_JPRB
84RALBSEA15 = 0.07_JPRB
85RALBSNO15 = 0.80_JPRB
86RALBSNM15 = 0.40_JPRB
87RSNOWAL15 = 0.01_JPRB
88!*  Concentration of the various trace gases (IPCC/SACC values for 1990)
89!        CO2         CH4        N2O        CFC11       CFC12
90!      353ppmv     1.72ppmv   310ppbv     280pptv     484pptv
91
92XCARDI  = 353.E-06_JPRB
93XCH4    = 1.72E-06_JPRB
94XN2O    = 310.E-09_JPRB
95XCFC11  = 280.E-12_JPRB
96XCFC12  = 484.E-12_JPRB
97
98ZAIRMWG = 28.970_JPRB
99ZCO2MWG = 44.011_JPRB
100ZCH4MWG = 16.043_JPRB
101ZN2OMWG = 44.013_JPRB
102ZC11MWG = 137.3686_JPRB
103ZC12MWG = 120.9140_JPRB
104
105! Ce qui concerne NAMSCEN commente par MPL le 15.04.09
106!CALL POSNAM(NULNAM,'NAMSCEN')
107!READ       (NULNAM, NAMSCEN)
108
109WRITE(NULOUT,'( &
110 & '' CO2   = '',E14.7,'' CH4   = '',E14.7 &
111 & ,'' N2O   = '',E14.7,'' CFC11 = '',E14.7 &
112 & ,'' CFC12 = '',E14.7,'' RI0 = '',E14.7 &
113 & )') XCARDI,XCH4,XN2O,XCFC11,XCFC12,RI0 
114
115RCARDI15  = XCARDI*ZCO2MWG/ZAIRMWG
116RCH415    = XCH4*ZCH4MWG/ZAIRMWG
117RN2O15    = XN2O*ZN2OMWG/ZAIRMWG
118RCFC1115  = XCFC11*ZC11MWG/ZAIRMWG
119RCFC1215  = XCFC12*ZC12MWG/ZAIRMWG
120REMISS15  = 0.996_JPRB
121!ZSUPSAT = 0.01_JPRB
122RVLBDC15  = 0.5_JPRB
123
124REPSEC15=1.E-12_JPRB
125REPCLC15=1.E-12_JPRB
126REPH2O15=1.E-12_JPRB
127REPALB15=1.E-12_JPRB
128
129!     -----------------------------------------------------------------
130
131IF (LHOOK) CALL DR_HOOK('SURDI15',1,ZHOOK_HANDLE)
132END SUBROUTINE SURDI15
Note: See TracBrowser for help on using the repository browser.