1 |
|
|
!OPTIONS XOPT(NOEVAL) |
2 |
|
|
SUBROUTINE SUCLOP15 |
3 |
|
|
|
4 |
|
|
!**** *SUCLOP15* - INITIALIZE COMMON YOMCLOP15 |
5 |
|
|
!**** FROZEN VERSION (CYCLE 15) OF SUCLOP |
6 |
|
|
|
7 |
|
|
! PURPOSE. |
8 |
|
|
! -------- |
9 |
|
|
! INITIALIZE YOMCLOP15, WITH CLOUD OPTICAL PARAMETERS |
10 |
|
|
|
11 |
|
|
!** INTERFACE. |
12 |
|
|
! ---------- |
13 |
|
|
! *CALL* SUCLOP15 |
14 |
|
|
! FROM *SUPHEC* |
15 |
|
|
|
16 |
|
|
! EXPLICIT ARGUMENTS : |
17 |
|
|
! -------------------- |
18 |
|
|
! NONE |
19 |
|
|
|
20 |
|
|
! IMPLICIT ARGUMENTS : |
21 |
|
|
! -------------------- |
22 |
|
|
! COMMON YOMCLOP15 |
23 |
|
|
|
24 |
|
|
! METHOD. |
25 |
|
|
! ------- |
26 |
|
|
! SEE DOCUMENTATION |
27 |
|
|
|
28 |
|
|
! EXTERNALS. |
29 |
|
|
! ---------- |
30 |
|
|
! NONE |
31 |
|
|
|
32 |
|
|
! REFERENCE. |
33 |
|
|
! ---------- |
34 |
|
|
! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE |
35 |
|
|
! "INTEGRATED FORECASTING SYSTEM" |
36 |
|
|
|
37 |
|
|
! AUTHOR. |
38 |
|
|
! ------- |
39 |
|
|
! 96-11: Ph. Dandin. Meteo-France |
40 |
|
|
! ORIGINAL : J.-J. MORCRETTE *ECMWF* |
41 |
|
|
|
42 |
|
|
! MODIFICATIONS. |
43 |
|
|
! -------------- |
44 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
45 |
|
|
! P. Marquet 14-Feb-2006 REFFWIA15 + NAMCLOP15 introduced |
46 |
|
|
|
47 |
|
|
! ------------------------------------------------------------------ |
48 |
|
|
|
49 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
50 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
51 |
|
|
|
52 |
|
|
! Ce qui concerne NULNAM commente par MPL le 15.04.09 |
53 |
|
|
!USE YOMLUN , ONLY : NULNAM, NULOUT |
54 |
|
|
USE YOMLUN , ONLY : NULOUT |
55 |
|
|
USE YOMCLOP15, ONLY : RYFWCA15 ,RYFWCB15 ,RYFWCC15 ,RYFWCD15 ,& |
56 |
|
|
& RYFWCE15 ,RYFWCF15 ,REBCUA15 ,REBCUB15 ,REBCUC15 ,& |
57 |
|
|
& REBCUD15 ,REBCUE15 ,REBCUF15 ,REBCUG15 ,REBCUH15 ,& |
58 |
|
|
& REFFIA15 ,REFFIB15 ,RTIW15 ,RRIW15 ,REFFWIA15 |
59 |
|
|
|
60 |
|
|
IMPLICIT NONE |
61 |
|
|
|
62 |
|
|
#include "posnam.intfb.h" |
63 |
|
|
#include "namclop15.h" |
64 |
|
|
|
65 |
|
|
|
66 |
|
|
!* 1. Set default values. |
67 |
|
|
! ------------------- |
68 |
|
|
|
69 |
|
|
!* Ice cloud properties - crystal: adapted from Ebert and Curry, 1992 |
70 |
|
|
|
71 |
|
|
! SW : 2 spectral intervals |
72 |
|
|
|
73 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
74 |
|
|
IF (LHOOK) CALL DR_HOOK('SUCLOP15',0,ZHOOK_HANDLE) |
75 |
|
|
REBCUA15(1)= 3.448E-03_JPRB |
76 |
|
|
REBCUA15(2)= 3.448E-03_JPRB |
77 |
|
|
REBCUB15(1)= 2.431_JPRB |
78 |
|
|
REBCUB15(2)= 2.431_JPRB |
79 |
|
|
REBCUC15(1)= 0.99999_JPRB |
80 |
|
|
REBCUC15(2)= 0.975634_JPRB |
81 |
|
|
REBCUD15(1)= 0._JPRB |
82 |
|
|
REBCUD15(2)= 2.487E-04_JPRB |
83 |
|
|
REBCUE15(1)= 0.7661_JPRB |
84 |
|
|
REBCUE15(2)= 0.7866_JPRB |
85 |
|
|
REBCUF15(1)= 5.851E-04_JPRB |
86 |
|
|
REBCUF15(2)= 5.937E-04_JPRB |
87 |
|
|
|
88 |
|
|
! LW : spectrally averaged with reference Planck function at 257 K |
89 |
|
|
|
90 |
|
|
REBCUG15= 1.07677_JPRB |
91 |
|
|
REBCUH15= 0.00267_JPRB |
92 |
|
|
|
93 |
|
|
! Ice particle Effective Radius as a function of LWC |
94 |
|
|
|
95 |
|
|
REFFIA15= 40._JPRB |
96 |
|
|
REFFIB15= 0._JPRB |
97 |
|
|
|
98 |
|
|
! Water and Ice particle Effective Radius in the RADLSW15 formulae |
99 |
|
|
|
100 |
|
|
REFFWIA15= 10._JPRB |
101 |
|
|
|
102 |
|
|
!* Water cloud properties - from Fouquart (1987) |
103 |
|
|
|
104 |
|
|
! SW : 2 spectral intervals: parameters as a function of Reff |
105 |
|
|
|
106 |
|
|
RYFWCA15(1)= 0._JPRB |
107 |
|
|
RYFWCA15(2)= 0._JPRB |
108 |
|
|
RYFWCB15(1)= 1.5_JPRB |
109 |
|
|
RYFWCB15(2)= 1.5_JPRB |
110 |
|
|
RYFWCC15(1)= 0.9999_JPRB |
111 |
|
|
RYFWCC15(2)= 0.9988_JPRB |
112 |
|
|
RYFWCD15(1)= 5.000E-04_JPRB |
113 |
|
|
RYFWCD15(2)= 2.500E-03_JPRB |
114 |
|
|
RYFWCE15(1)= 0.5_JPRB |
115 |
|
|
RYFWCE15(2)= 0.05_JPRB |
116 |
|
|
RYFWCF15(1)= 0.865_JPRB |
117 |
|
|
RYFWCF15(2)= 0.910_JPRB |
118 |
|
|
|
119 |
|
|
!* Liquid/Solid water transition |
120 |
|
|
|
121 |
|
|
RTIW15= 263._JPRB |
122 |
|
|
RRIW15= 20._JPRB |
123 |
|
|
|
124 |
|
|
|
125 |
|
|
!* 2. Modify default values. |
126 |
|
|
! ---------------------- |
127 |
|
|
|
128 |
|
|
! Ce qui concerne NAMCLOP15 commente par MPL le 15.04.09 |
129 |
|
|
!CALL POSNAM(NULNAM,'NAMCLOP15') |
130 |
|
|
!READ (NULNAM, NAMCLOP15) |
131 |
|
|
|
132 |
|
|
|
133 |
|
|
!* 3. Print final values. |
134 |
|
|
! ------------------- |
135 |
|
|
|
136 |
|
|
WRITE(UNIT=NULOUT,FMT='('' - - - - - - - - -'')') |
137 |
|
|
WRITE(UNIT=NULOUT,FMT='('' COMMON YOMCLOP15 '')') |
138 |
|
|
WRITE(UNIT=NULOUT,FMT='('' - - - - - - - - -'')') |
139 |
|
|
WRITE(UNIT=NULOUT,FMT='( '' REFFWIA15 = '',E11.4 )') & |
140 |
|
|
& REFFWIA15 |
141 |
|
|
|
142 |
|
|
IF (LHOOK) CALL DR_HOOK('SUCLOP15',1,ZHOOK_HANDLE) |
143 |
|
|
|
144 |
|
|
! ------------------------------------------------------------------ |
145 |
|
|
|
146 |
|
|
END SUBROUTINE SUCLOP15 |