source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/suclop.F90 @ 5308

Last change on this file since 5308 was 3331, checked in by acozic, 6 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 2.5 KB
Line 
1SUBROUTINE SUCLOP
2
3!**** *SUCLOP*  - INITIALIZE COMMON YOECLOP
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOMCLOP, WITH CLOUD OPTICAL PARAMETERS
8
9!**   INTERFACE.
10!     ----------
11!        *CALL*  SUCLOP
12!     FROM *SUECRAD*
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        COMMON YOECLOP
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28!        NONE
29
30!     REFERENCE.
31!     ----------
32!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
33!     "INTEGRATED FORECASTING SYSTEM"
34
35!     AUTHOR.
36!     -------
37!        JEAN-JACQUES MORCRETTE  *ECMWF*
38
39!     MODIFICATIONS.
40!     --------------
41!        ORIGINAL : 92-02-29
42!        M.Hamrud      01-Oct-2003 CY28 Cleaning
43!     ------------------------------------------------------------------
44
45USE PARKIND1  ,ONLY : JPIM     ,JPRB
46USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
47
48USE YOECLOP  , ONLY : RYFWCA   ,RYFWCB   ,RYFWCC   ,RYFWCD   ,&
49 & RYFWCE   ,RYFWCF   ,REBCUA   ,REBCUB   ,REBCUC   ,&
50 & REBCUD   ,REBCUE   ,REBCUF   ,REBCUG   ,REBCUH   ,&
51 & REFFIA   ,REFFIB   ,RTIW     ,RRIW 
52
53IMPLICIT NONE
54
55!      ----------------------------------------------------------------
56
57!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992
58
59! SW : 2 spectral intervals
60
61REAL(KIND=JPRB) :: ZHOOK_HANDLE
62IF (LHOOK) CALL DR_HOOK('SUCLOP',0,ZHOOK_HANDLE)
63REBCUA(1)= 3.448E-03_JPRB
64REBCUA(2)= 3.448E-03_JPRB
65REBCUB(1)= 2.431_JPRB
66REBCUB(2)= 2.431_JPRB
67REBCUC(1)= 0.99999_JPRB
68REBCUC(2)= 0.975634_JPRB
69REBCUD(1)= 0._JPRB
70REBCUD(2)= 2.487E-04_JPRB
71REBCUE(1)= 0.7661_JPRB
72REBCUE(2)= 0.7866_JPRB
73REBCUF(1)= 5.851E-04_JPRB
74REBCUF(2)= 5.937E-04_JPRB
75
76! LW : spectrally averaged with reference Planck function at 257 K
77
78REBCUG= 1.07677_JPRB
79REBCUH= 0.00267_JPRB
80
81! Ice particle Effective Radius as a function of LWC
82
83REFFIA= 40._JPRB
84REFFIB=  0._JPRB
85
86!* Water cloud properties - from Fouquart (1987)
87
88! SW : 2 spectral intervals: parameters as a function of Reff
89
90RYFWCA(1)= 0._JPRB
91RYFWCA(2)= 0._JPRB
92RYFWCB(1)= 1.5_JPRB
93RYFWCB(2)= 1.5_JPRB
94RYFWCC(1)= 0.9999_JPRB
95RYFWCC(2)= 0.9988_JPRB
96RYFWCD(1)= 5.000E-04_JPRB
97RYFWCD(2)= 2.500E-03_JPRB
98RYFWCE(1)= 0.5_JPRB
99RYFWCE(2)= 0.05_JPRB
100RYFWCF(1)= 0.865_JPRB
101RYFWCF(2)= 0.910_JPRB
102
103!* Liquid/Solid water transition
104
105RTIW= 263._JPRB
106RRIW= 20._JPRB
107IF (LHOOK) CALL DR_HOOK('SUCLOP',1,ZHOOK_HANDLE)
108
109!     ------------------------------------------------------------------
110
111END SUBROUTINE SUCLOP
Note: See TracBrowser for help on using the repository browser.