1 |
|
|
!OPTIONS XOPT(NOEVAL) |
2 |
|
1 |
SUBROUTINE SUPHY3(KULOUT) |
3 |
|
|
|
4 |
|
|
!**** *SUPHY3* - Initialize common YOMPHY3 physics radiative |
5 |
|
|
! constants |
6 |
|
|
|
7 |
|
|
! Purpose. |
8 |
|
|
! -------- |
9 |
|
|
! Initialize YOMPHY3, the common that contains the parameters |
10 |
|
|
! for the radiation part of the physics of the model. |
11 |
|
|
|
12 |
|
|
!** Interface. |
13 |
|
|
! ---------- |
14 |
|
|
! *CALL* *SUPHY3(KULOUT) |
15 |
|
|
|
16 |
|
|
! Explicit arguments : |
17 |
|
|
! -------------------- |
18 |
|
|
! KULOUT : Logical unit for the output |
19 |
|
|
|
20 |
|
|
! Implicit arguments : |
21 |
|
|
! -------------------- |
22 |
|
|
! COMMON YOMPHY3 |
23 |
|
|
|
24 |
|
|
! Method. |
25 |
|
|
! ------- |
26 |
|
|
! See documentation |
27 |
|
|
|
28 |
|
|
! Externals. |
29 |
|
|
! ---------- |
30 |
|
|
|
31 |
|
|
! Reference. |
32 |
|
|
! ---------- |
33 |
|
|
! Documentation ARPEGE |
34 |
|
|
|
35 |
|
|
! Author. |
36 |
|
|
! ------- |
37 |
|
|
! J.-F. Geleyn . |
38 |
|
|
|
39 |
|
|
! Modifications. |
40 |
|
|
! -------------- |
41 |
|
|
! Original : 90-9-1 |
42 |
|
|
! Ajout de GCE4 et |
43 |
|
|
! modification de toutes les constantes optiques : 91-2-1 |
44 |
|
|
! Ajout des nuages de glace : 92-04, L. Labbe. |
45 |
|
|
! Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF |
46 |
|
|
! Modified 97-04-17 J.M. Piriou: radiative coefficients default values, and new QSUSX default. |
47 |
|
|
! Modified 2000-08 R. Randriamampianina, J.F. Geleyn et J.M. Piriou: change infra-red interaction between layers (GIREC*). |
48 |
|
|
! Modified 2000-12-12 by E. Bazile: CYCORA's default value. |
49 |
|
|
! M.Hamrud 01-Oct-2003 CY28 Cleaning |
50 |
|
|
! R.Brozkova 24-Sep-2006 ALARO-0 (new cloud model) |
51 |
|
|
! ------------------------------------------------------------------ |
52 |
|
|
|
53 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
54 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
55 |
|
|
|
56 |
|
|
! Ce qui concerne NULNAM commente par MPL le 15.04.09 |
57 |
|
|
!USE YOMLUN , ONLY : NULNAM |
58 |
|
|
USE YOMCT0B , ONLY : LECMWF |
59 |
|
|
USE YOMPHY , ONLY : LCLSATUR |
60 |
|
|
USE YOMPHY3 , ONLY : GCA ,GCB ,GCC ,VDP ,& |
61 |
|
|
& VNP ,BSFSA ,BSFSI ,BSFSN ,BSFTA ,& |
62 |
|
|
& BSFTI ,BSFTN ,EARRT ,EOASA ,EOASI ,& |
63 |
|
|
& EOASN ,EOATA ,EOATI ,EOATN ,EODSA ,& |
64 |
|
|
& EODSI ,EODSN ,EODTA ,EODTI ,EODTN ,& |
65 |
|
|
& EORAY ,GCD4 ,GCE4 ,QCO2 ,QLIMI ,& |
66 |
|
|
& QLIP0 ,RII0 ,USAA ,USAI ,& |
67 |
|
|
& USAN ,USBA ,USBI ,USBN ,& |
68 |
|
|
& GIREC1 ,GIREC2 ,GIREC3 ,GIREC4 ,& |
69 |
|
|
& FCM_DEL_A,FCM_DEL_D,FCM_MU_A ,FCM_MU_D ,FCM_N_I ,& |
70 |
|
|
& FCM_N_L ,FCM_P_AI ,FCM_P_AL ,FCM_P_DI ,FCM_P_DL ,& |
71 |
|
|
& FCM_P_GI ,FCM_P_GL ,FCM_Q_AI ,FCM_Q_AL ,FCM_Q_DI ,& |
72 |
|
|
& FCM_Q_DL ,FCM_Q_GI ,FCM_Q_GL ,N_SPBAND ,REXP_NEB |
73 |
|
|
|
74 |
|
|
IMPLICIT NONE |
75 |
|
|
|
76 |
|
|
INTEGER(KIND=JPIM),INTENT(IN) :: KULOUT |
77 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
78 |
|
|
|
79 |
|
|
#include "posnam.intfb.h" |
80 |
|
|
|
81 |
|
|
#include "namphy3.h" |
82 |
|
|
! ------------------------------------------------------------------ |
83 |
|
|
|
84 |
|
|
!* 1. Set default values. |
85 |
|
|
! ------------------- |
86 |
|
|
|
87 |
|
|
! 1.1 Set implicit default values |
88 |
|
|
|
89 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SUPHY3',0,ZHOOK_HANDLE) |
90 |
|
1 |
BSFSA=.3490_JPRB |
91 |
|
1 |
BSFSI=.3238_JPRB |
92 |
|
1 |
BSFSN=.3300_JPRB |
93 |
|
1 |
BSFTA=.3471_JPRB |
94 |
|
1 |
BSFTI=.3231_JPRB |
95 |
|
1 |
BSFTN=.3584_JPRB |
96 |
|
1 |
EARRT=0.001324_JPRB |
97 |
|
1 |
EOASA=.6727E-01_JPRB |
98 |
|
1 |
EOASI=.1692E-01_JPRB |
99 |
|
1 |
EOASN=.3353E-01_JPRB |
100 |
|
1 |
EOATA=.5037E-01_JPRB |
101 |
|
1 |
EOATI=.3455E+01_JPRB |
102 |
|
1 |
EOATN=.9077E+01_JPRB |
103 |
|
1 |
EODSA=.3665E+00_JPRB |
104 |
|
1 |
EODSI=.1406E+01_JPRB |
105 |
|
1 |
EODSN=.6558E+01_JPRB |
106 |
|
1 |
EODTA=.1766E-01_JPRB |
107 |
|
1 |
EODTI=.8435E+00_JPRB |
108 |
|
1 |
EODTN=.4319E+01_JPRB |
109 |
|
1 |
EORAY=.8606E-06_JPRB |
110 |
|
1 |
GCA(1)=.8041E-01_JPRB |
111 |
|
1 |
GCA(2)=.1456E+00_JPRB |
112 |
|
1 |
GCA(3)=.4787E+01_JPRB |
113 |
|
1 |
GCA(4)=.2102E+04_JPRB |
114 |
|
1 |
GCA(5)=.1334E+01_JPRB |
115 |
|
1 |
GCA(6)=.1551E-04_JPRB |
116 |
|
1 |
GCB(1)=.8968E+07_JPRB |
117 |
|
1 |
GCB(2)=.2413E+10_JPRB |
118 |
|
1 |
GCB(3)=.3548E+05_JPRB |
119 |
|
1 |
GCB(4)=.6370E+10_JPRB |
120 |
|
1 |
GCB(5)=.8499E+11_JPRB |
121 |
|
1 |
GCB(6)=.1012E+06_JPRB |
122 |
|
1 |
GCC(1)=.5925E-10_JPRB |
123 |
|
1 |
GCC(2)=.1842E-10_JPRB |
124 |
|
1 |
GCC(3)=.2532E-07_JPRB |
125 |
|
1 |
GCC(4)=.1953E+07_JPRB |
126 |
|
1 |
GCC(5)=.1734E-11_JPRB |
127 |
|
1 |
GCC(6)=.1225E-16_JPRB |
128 |
|
1 |
GCD4=.3608E-69_JPRB |
129 |
|
1 |
GCE4=.7563E+04_JPRB |
130 |
|
1 |
GIREC1=0.0_JPRB |
131 |
|
1 |
GIREC2=0.0_JPRB |
132 |
|
1 |
GIREC3=0.0_JPRB |
133 |
|
1 |
GIREC4=0.0_JPRB |
134 |
|
1 |
QCO2=.5366E-03_JPRB |
135 |
|
1 |
QLIMI=15000._JPRB |
136 |
|
1 |
QLIP0=8.E+06_JPRB |
137 |
|
1 |
USAA=-.3020_JPRB |
138 |
|
1 |
USAI=-.3524_JPRB |
139 |
|
1 |
USAN=-.3400_JPRB |
140 |
|
1 |
USBA=0._JPRB |
141 |
|
1 |
USBI=0._JPRB |
142 |
|
1 |
USBN=0._JPRB |
143 |
|
1 |
VDP(1,1)=.21868E+02_JPRB |
144 |
|
1 |
VDP(2,1)=.17453E+02_JPRB |
145 |
|
1 |
VDP(3,1)=.68918E+00_JPRB |
146 |
|
1 |
VDP(4,1)=.23456E-03_JPRB |
147 |
|
1 |
VDP(5,1)=.22317E-09_JPRB |
148 |
|
1 |
VDP(1,2)=.48401E+02_JPRB |
149 |
|
1 |
VDP(2,2)=.18648E+02_JPRB |
150 |
|
1 |
VDP(3,2)=.35199E-01_JPRB |
151 |
|
1 |
VDP(4,2)=.63691E-06_JPRB |
152 |
|
1 |
VDP(5,2)=.86395E-13_JPRB |
153 |
|
1 |
VDP(1,3)=.76948E+02_JPRB |
154 |
|
1 |
VDP(2,3)=.41056E+01_JPRB |
155 |
|
1 |
VDP(3,3)=.42667E-03_JPRB |
156 |
|
1 |
VDP(4,3)=0._JPRB |
157 |
|
1 |
VDP(5,3)=0._JPRB |
158 |
|
1 |
VDP(1,4)=.27241E+03_JPRB |
159 |
|
1 |
VDP(2,4)=.57091E+04_JPRB |
160 |
|
1 |
VDP(3,4)=.14393E+05_JPRB |
161 |
|
1 |
VDP(4,4)=.29879E+04_JPRB |
162 |
|
1 |
VDP(5,4)=.25382E+02_JPRB |
163 |
|
1 |
VDP(1,5)=.86942E+02_JPRB |
164 |
|
1 |
VDP(2,5)=.32186E+03_JPRB |
165 |
|
1 |
VDP(3,5)=.10775E+03_JPRB |
166 |
|
1 |
VDP(4,5)=.21261E+01_JPRB |
167 |
|
1 |
VDP(5,5)=.40003E-02_JPRB |
168 |
|
1 |
VDP(1,6)=.24408E+02_JPRB |
169 |
|
1 |
VDP(2,6)=.81919E+01_JPRB |
170 |
|
1 |
VDP(3,6)=.72193E+01_JPRB |
171 |
|
1 |
VDP(4,6)=.56230E+00_JPRB |
172 |
|
1 |
VDP(5,6)=.25384E-02_JPRB |
173 |
|
1 |
VNP(1,1)=.69926E+01_JPRB |
174 |
|
1 |
VNP(2,1)=.63915E+00_JPRB |
175 |
|
1 |
VNP(3,1)=.28896E-03_JPRB |
176 |
|
1 |
VNP(4,1)=.10050E-08_JPRB |
177 |
|
1 |
VNP(5,1)=.99037E-16_JPRB |
178 |
|
1 |
VNP(1,2)=.31105E+01_JPRB |
179 |
|
1 |
VNP(2,2)=.14225E-01_JPRB |
180 |
|
1 |
VNP(3,2)=.69355E-06_JPRB |
181 |
|
1 |
VNP(4,2)=.36087E-12_JPRB |
182 |
|
1 |
VNP(5,2)=.44113E-20_JPRB |
183 |
|
1 |
VNP(1,3)=.23659E+01_JPRB |
184 |
|
1 |
VNP(2,3)=.11139E-02_JPRB |
185 |
|
1 |
VNP(3,3)=.10618E-05_JPRB |
186 |
|
1 |
VNP(4,3)=0._JPRB |
187 |
|
1 |
VNP(5,3)=0._JPRB |
188 |
|
1 |
VNP(1,4)=.19810E+03_JPRB |
189 |
|
1 |
VNP(2,4)=.46954E+04_JPRB |
190 |
|
1 |
VNP(3,4)=.22512E+04_JPRB |
191 |
|
1 |
VNP(4,4)=.52461E+02_JPRB |
192 |
|
1 |
VNP(5,4)=.11645E+00_JPRB |
193 |
|
1 |
VNP(1,5)=.46348E+02_JPRB |
194 |
|
1 |
VNP(2,5)=.35630E+02_JPRB |
195 |
|
1 |
VNP(3,5)=.33005E+01_JPRB |
196 |
|
1 |
VNP(4,5)=.18045E-01_JPRB |
197 |
|
1 |
VNP(5,5)=.88667E-05_JPRB |
198 |
|
1 |
VNP(1,6)=.47413E+01_JPRB |
199 |
|
1 |
VNP(2,6)=.16334E+01_JPRB |
200 |
|
1 |
VNP(3,6)=.48164E+00_JPRB |
201 |
|
1 |
VNP(4,6)=.56140E-02_JPRB |
202 |
|
1 |
VNP(5,6)=.67790E-04_JPRB |
203 |
|
|
|
204 |
|
|
! Default values for cloud model - first index denotes spectral band: |
205 |
|
|
! 1 - solar |
206 |
|
|
! 2 - thermal |
207 |
|
1 |
FCM_DEL_A(1) = 9.55496e-01_JPRB |
208 |
|
1 |
FCM_DEL_A(2) = 1.79152e+02_JPRB |
209 |
|
1 |
FCM_DEL_D(1) = 5.72959e+01_JPRB |
210 |
|
1 |
FCM_DEL_D(2) = 2.00159e+02_JPRB |
211 |
|
1 |
FCM_MU_A(1) = 8.66200e-01_JPRB |
212 |
|
1 |
FCM_MU_A(2) = 6.65532e-01_JPRB |
213 |
|
1 |
FCM_MU_D(1) = 1.00000e+00_JPRB |
214 |
|
1 |
FCM_MU_D(2) = 6.70517e-01_JPRB |
215 |
|
1 |
FCM_N_I = 1.00000e-01_JPRB |
216 |
|
1 |
FCM_N_L = 1.00000e-01_JPRB |
217 |
|
1 |
FCM_P_AI(1,0) = 1.48920e+01_JPRB |
218 |
|
1 |
FCM_P_AI(1,1) = -3.71379e+01_JPRB |
219 |
|
1 |
FCM_P_AI(1,2) = 0.00000e+00_JPRB |
220 |
|
1 |
FCM_P_AI(1,3) = 0.00000e+00_JPRB |
221 |
|
1 |
FCM_P_AI(2,0) = 9.76920e+00_JPRB |
222 |
|
1 |
FCM_P_AI(2,1) = -7.54840e+00_JPRB |
223 |
|
1 |
FCM_P_AI(2,2) = 0.00000e+00_JPRB |
224 |
|
1 |
FCM_P_AI(2,3) = 0.00000e+00_JPRB |
225 |
|
1 |
FCM_P_AL(1,0) = 1.42062e+00_JPRB |
226 |
|
1 |
FCM_P_AL(1,1) = -3.54037e+00_JPRB |
227 |
|
1 |
FCM_P_AL(1,2) = 1.73535e+00_JPRB |
228 |
|
1 |
FCM_P_AL(1,3) = 0.00000e+00_JPRB |
229 |
|
1 |
FCM_P_AL(2,0) = 5.08522e+00_JPRB |
230 |
|
1 |
FCM_P_AL(2,1) = -1.08144e+01_JPRB |
231 |
|
1 |
FCM_P_AL(2,2) = 5.16331e+00_JPRB |
232 |
|
1 |
FCM_P_AL(2,3) = 0.00000e+00_JPRB |
233 |
|
1 |
FCM_P_DI(1,0) = 1.30487e+01_JPRB |
234 |
|
1 |
FCM_P_DI(1,1) = 1.98670e+01_JPRB |
235 |
|
1 |
FCM_P_DI(1,2) = -5.08578e+01_JPRB |
236 |
|
1 |
FCM_P_DI(1,3) = 0.00000e+00_JPRB |
237 |
|
1 |
FCM_P_DI(2,0) = 1.23173e+01_JPRB |
238 |
|
1 |
FCM_P_DI(2,1) = -3.32808e+01_JPRB |
239 |
|
1 |
FCM_P_DI(2,2) = 5.77295e+01_JPRB |
240 |
|
1 |
FCM_P_DI(2,3) = -7.35728e+01_JPRB |
241 |
|
1 |
FCM_P_DL(1,0) = 4.68646e+00_JPRB |
242 |
|
1 |
FCM_P_DL(1,1) = -8.77652e+00_JPRB |
243 |
|
1 |
FCM_P_DL(1,2) = 2.12789e+00_JPRB |
244 |
|
1 |
FCM_P_DL(1,3) = 0.00000e+00_JPRB |
245 |
|
1 |
FCM_P_DL(2,0) = 4.52292e+00_JPRB |
246 |
|
1 |
FCM_P_DL(2,1) = -1.43327e+01_JPRB |
247 |
|
1 |
FCM_P_DL(2,2) = 1.16650e+01_JPRB |
248 |
|
1 |
FCM_P_DL(2,3) = 0.00000e+00_JPRB |
249 |
|
1 |
FCM_P_GI(1,0) = -2.01472e+00_JPRB |
250 |
|
1 |
FCM_P_GI(1,1) = 9.16284e+01_JPRB |
251 |
|
1 |
FCM_P_GI(1,2) = 7.82275e+01_JPRB |
252 |
|
1 |
FCM_P_GI(1,3) = 0.00000e+00_JPRB |
253 |
|
1 |
FCM_P_GI(2,0) = -7.55124e-01_JPRB |
254 |
|
1 |
FCM_P_GI(2,1) = -4.68837e+00_JPRB |
255 |
|
1 |
FCM_P_GI(2,2) = 6.72877e+01_JPRB |
256 |
|
1 |
FCM_P_GI(2,3) = 0.00000e+00_JPRB |
257 |
|
1 |
FCM_P_GL(1,0) = 8.82686e-01_JPRB |
258 |
|
1 |
FCM_P_GL(1,1) = -2.02445e+00_JPRB |
259 |
|
1 |
FCM_P_GL(1,2) = 1.45259e+00_JPRB |
260 |
|
1 |
FCM_P_GL(1,3) = 0.00000e+00_JPRB |
261 |
|
1 |
FCM_P_GL(2,0) = 1.93784e-01_JPRB |
262 |
|
1 |
FCM_P_GL(2,1) = -6.07676e-01_JPRB |
263 |
|
1 |
FCM_P_GL(2,2) = 5.37734e-01_JPRB |
264 |
|
1 |
FCM_P_GL(2,3) = 0.00000e+00_JPRB |
265 |
|
1 |
FCM_Q_AI(1,1) = 1.86588e+01_JPRB |
266 |
|
1 |
FCM_Q_AI(1,2) = 0.00000e+00_JPRB |
267 |
|
1 |
FCM_Q_AI(1,3) = 0.00000e+00_JPRB |
268 |
|
1 |
FCM_Q_AI(2,1) = 2.78063e+00_JPRB |
269 |
|
1 |
FCM_Q_AI(2,2) = 0.00000e+00_JPRB |
270 |
|
1 |
FCM_Q_AI(2,3) = 0.00000e+00_JPRB |
271 |
|
1 |
FCM_Q_AL(1,1) = -2.81210e+00_JPRB |
272 |
|
1 |
FCM_Q_AL(1,2) = 2.35758e+00_JPRB |
273 |
|
1 |
FCM_Q_AL(1,3) = 0.00000e+00_JPRB |
274 |
|
1 |
FCM_Q_AL(2,1) = -2.31252e+00_JPRB |
275 |
|
1 |
FCM_Q_AL(2,2) = 1.56848e+00_JPRB |
276 |
|
1 |
FCM_Q_AL(2,3) = 0.00000e+00_JPRB |
277 |
|
1 |
FCM_Q_DI(1,1) = 1.04582e+01_JPRB |
278 |
|
1 |
FCM_Q_DI(1,2) = 6.15641e+00_JPRB |
279 |
|
1 |
FCM_Q_DI(1,3) = 0.00000e+00_JPRB |
280 |
|
1 |
FCM_Q_DI(2,1) = 1.18370e+00_JPRB |
281 |
|
1 |
FCM_Q_DI(2,2) = 2.23187e+00_JPRB |
282 |
|
1 |
FCM_Q_DI(2,3) = 7.00424e+00_JPRB |
283 |
|
1 |
FCM_Q_DL(1,1) = -2.18205e+00_JPRB |
284 |
|
1 |
FCM_Q_DL(1,2) = 1.34831e+00_JPRB |
285 |
|
1 |
FCM_Q_DL(1,3) = 0.00000e+00_JPRB |
286 |
|
1 |
FCM_Q_DL(2,1) = -3.12310e+00_JPRB |
287 |
|
1 |
FCM_Q_DL(2,2) = 2.57959e+00_JPRB |
288 |
|
1 |
FCM_Q_DL(2,3) = 0.00000e+00_JPRB |
289 |
|
1 |
FCM_Q_GI(1,1) = 7.89999e+01_JPRB |
290 |
|
1 |
FCM_Q_GI(1,2) = 5.24929e+01_JPRB |
291 |
|
1 |
FCM_Q_GI(1,3) = 0.00000e+00_JPRB |
292 |
|
1 |
FCM_Q_GI(2,1) = 4.93801e+00_JPRB |
293 |
|
1 |
FCM_Q_GI(2,2) = 2.38096e+01_JPRB |
294 |
|
1 |
FCM_Q_GI(2,3) = 0.00000e+00_JPRB |
295 |
|
1 |
FCM_Q_GL(1,1) = -2.20830e+00_JPRB |
296 |
|
1 |
FCM_Q_GL(1,2) = 1.39832e+00_JPRB |
297 |
|
1 |
FCM_Q_GL(1,3) = 0.00000e+00_JPRB |
298 |
|
1 |
FCM_Q_GL(2,1) = -3.27770e+00_JPRB |
299 |
|
1 |
FCM_Q_GL(2,2) = 2.78446e+00_JPRB |
300 |
|
1 |
FCM_Q_GL(2,3) = 0.00000e+00_JPRB |
301 |
|
1 |
REXP_NEB = 8.0_JPRB |
302 |
|
|
|
303 |
|
|
! 1.2 Modify default values according to LECMWF |
304 |
|
|
|
305 |
|
|
IF (LECMWF) THEN |
306 |
|
|
ELSE |
307 |
|
|
ENDIF |
308 |
|
|
|
309 |
|
|
! Remark : the values or RII0 is calculated and not set up |
310 |
|
|
|
311 |
|
|
! ------------------------------------------------------------------ |
312 |
|
|
|
313 |
|
|
!* 2. Modify default values. |
314 |
|
|
! ---------------------- |
315 |
|
|
|
316 |
|
|
! Ce qui concerne NAMPHY3 commente par MPL le 15.04.09 |
317 |
|
|
!CALL POSNAM(NULNAM,'NAMPHY3') |
318 |
|
|
!READ(NULNAM,NAMPHY3) |
319 |
|
|
! ------------------------------------------------------------------ |
320 |
|
|
|
321 |
|
|
!* 3. Print final values. |
322 |
|
|
! ------------------- |
323 |
|
|
|
324 |
|
1 |
WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY3 '')') |
325 |
|
|
WRITE(UNIT=KULOUT,FMT='('' BSFSA = '',E10.4,'' BSFSI = '',E10.4 & |
326 |
|
|
& ,'' BSFSN = '',E10.4,'' BSFTA = '',E10.4,'' BSFTI = '',E10.4,/& |
327 |
|
|
& ,'' BSFTN = '',E10.4,'' EARRT = '',E10.4,'' EOASA = ''& |
328 |
|
|
& ,E10.4,'' EOASI = '',E10.4,'' EOASN = '',E10.4 & |
329 |
|
|
& ,'' EOATA = '',E10.4,/,'' EOATI = '',E10.4,'' EOATN = '',E10.4 & |
330 |
|
|
& ,'' EODSA = '',E10.4,'' EODSI = '',E10.4,'' EODSN = '',E10.4 & |
331 |
|
|
& ,'' EODTA = '',E10.4,/,'' EODTI = '',E10.4 & |
332 |
|
|
& ,'' EODTN = '',E10.4,'' EORAY = '',E10.4,/& |
333 |
|
|
& ,'' GCA(6) = '',6E11.4,/,'' GCB(6) = '',6E11.4,/,'' GCC(6) = ''& |
334 |
|
|
& ,6E11.4,/,'' GCD4 = '',E11.4,'' GCE4 = '',E11.4,/& |
335 |
|
|
& ,'' QCO2 = '',E10.4,'' QLIP0 = '',E10.4 & |
336 |
|
|
& ,'' USAA = '',E10.4,'' USAI = '',E10.4,'' USAN = '',E10.4 & |
337 |
|
|
& ,'' USBA = '',E10.4,/,'' USBI = '',E10.4,'' USBN = '',E10.4,/& |
338 |
|
|
& ,'' VDP(5,6) = '',5E12.5,/,5(12X,5E12.5,/)& |
339 |
|
|
& ,'' VNP(5,6) = '',5E12.5,/,5(12X,5E12.5,/)& |
340 |
|
|
& ,'' QLIMI = '',E10.4 & |
341 |
|
|
& ,'' GIREC1 = '',E10.4,'' GIREC2 = '',E10.4,'' GIREC3 = '',E10.4 & |
342 |
|
|
& ,'' GIREC4 = '',E10.4 & |
343 |
|
|
& )')& |
344 |
|
1 |
& BSFSA,BSFSI,BSFSN,BSFTA,BSFTI,BSFTN,EARRT,EOASA,EOASI,EOASN,& |
345 |
|
1 |
& EOATA,EOATI,EOATN,EODSA,EODSI,EODSN,EODTA,EODTI,EODTN,EORAY,& |
346 |
|
1 |
& GCA,GCB,GCC,GCD4,GCE4,QCO2,QLIP0,USAA,USAI,USAN,USBA,USBI,USBN,& |
347 |
|
2 |
& VDP,VNP,QLIMI,GIREC1,GIREC2,GIREC3,GIREC4 |
348 |
|
|
|
349 |
✗✓ |
1 |
IF (LCLSATUR) THEN |
350 |
|
|
WRITE(KULOUT,'('' N_SPBAND = '',I2,'' REXP_NEB = '',F5.2)') & |
351 |
|
|
& N_SPBAND,REXP_NEB |
352 |
|
|
WRITE(KULOUT,'('' FCM_DEL_A(:) ='',2(1X,ES12.5))') FCM_DEL_A(:) |
353 |
|
|
WRITE(KULOUT,'('' FCM_DEL_D(:) ='',2(1X,ES12.5))') FCM_DEL_D(:) |
354 |
|
|
WRITE(KULOUT,'('' FCM_MU_A (:) ='',2(1X,ES12.5))') FCM_MU_A (:) |
355 |
|
|
WRITE(KULOUT,'('' FCM_MU_D (:) ='',2(1X,ES12.5))') FCM_MU_D (:) |
356 |
|
|
WRITE(KULOUT,'('' FCM_N_I ='',1X,ES12.5)') FCM_N_I |
357 |
|
|
WRITE(KULOUT,'('' FCM_N_L ='',1X,ES12.5)') FCM_N_L |
358 |
|
|
WRITE(KULOUT,'('' FCM_P_AI (:,0) ='',2(1X,ES12.5))') FCM_P_AI (:,0) |
359 |
|
|
WRITE(KULOUT,'('' FCM_P_AI (:,1) ='',2(1X,ES12.5))') FCM_P_AI (:,1) |
360 |
|
|
WRITE(KULOUT,'('' FCM_P_AI (:,2) ='',2(1X,ES12.5))') FCM_P_AI (:,2) |
361 |
|
|
WRITE(KULOUT,'('' FCM_P_AI (:,3) ='',2(1X,ES12.5))') FCM_P_AI (:,3) |
362 |
|
|
WRITE(KULOUT,'('' FCM_P_AL (:,0) ='',2(1X,ES12.5))') FCM_P_AL (:,0) |
363 |
|
|
WRITE(KULOUT,'('' FCM_P_AL (:,1) ='',2(1X,ES12.5))') FCM_P_AL (:,1) |
364 |
|
|
WRITE(KULOUT,'('' FCM_P_AL (:,2) ='',2(1X,ES12.5))') FCM_P_AL (:,2) |
365 |
|
|
WRITE(KULOUT,'('' FCM_P_AL (:,3) ='',2(1X,ES12.5))') FCM_P_AL (:,3) |
366 |
|
|
WRITE(KULOUT,'('' FCM_P_DI (:,0) ='',2(1X,ES12.5))') FCM_P_DI (:,0) |
367 |
|
|
WRITE(KULOUT,'('' FCM_P_DI (:,1) ='',2(1X,ES12.5))') FCM_P_DI (:,1) |
368 |
|
|
WRITE(KULOUT,'('' FCM_P_DI (:,2) ='',2(1X,ES12.5))') FCM_P_DI (:,2) |
369 |
|
|
WRITE(KULOUT,'('' FCM_P_DI (:,3) ='',2(1X,ES12.5))') FCM_P_DI (:,3) |
370 |
|
|
WRITE(KULOUT,'('' FCM_P_DL (:,0) ='',2(1X,ES12.5))') FCM_P_DL (:,0) |
371 |
|
|
WRITE(KULOUT,'('' FCM_P_DL (:,1) ='',2(1X,ES12.5))') FCM_P_DL (:,1) |
372 |
|
|
WRITE(KULOUT,'('' FCM_P_DL (:,2) ='',2(1X,ES12.5))') FCM_P_DL (:,2) |
373 |
|
|
WRITE(KULOUT,'('' FCM_P_DL (:,3) ='',2(1X,ES12.5))') FCM_P_DL (:,3) |
374 |
|
|
WRITE(KULOUT,'('' FCM_P_GI (:,0) ='',2(1X,ES12.5))') FCM_P_GI (:,0) |
375 |
|
|
WRITE(KULOUT,'('' FCM_P_GI (:,1) ='',2(1X,ES12.5))') FCM_P_GI (:,1) |
376 |
|
|
WRITE(KULOUT,'('' FCM_P_GI (:,2) ='',2(1X,ES12.5))') FCM_P_GI (:,2) |
377 |
|
|
WRITE(KULOUT,'('' FCM_P_GI (:,3) ='',2(1X,ES12.5))') FCM_P_GI (:,3) |
378 |
|
|
WRITE(KULOUT,'('' FCM_P_GL (:,0) ='',2(1X,ES12.5))') FCM_P_GL (:,0) |
379 |
|
|
WRITE(KULOUT,'('' FCM_P_GL (:,1) ='',2(1X,ES12.5))') FCM_P_GL (:,1) |
380 |
|
|
WRITE(KULOUT,'('' FCM_P_GL (:,2) ='',2(1X,ES12.5))') FCM_P_GL (:,2) |
381 |
|
|
WRITE(KULOUT,'('' FCM_P_GL (:,3) ='',2(1X,ES12.5))') FCM_P_GL (:,3) |
382 |
|
|
WRITE(KULOUT,'('' FCM_Q_AI (:,1) ='',2(1X,ES12.5))') FCM_Q_AI (:,1) |
383 |
|
|
WRITE(KULOUT,'('' FCM_Q_AI (:,2) ='',2(1X,ES12.5))') FCM_Q_AI (:,2) |
384 |
|
|
WRITE(KULOUT,'('' FCM_Q_AI (:,3) ='',2(1X,ES12.5))') FCM_Q_AI (:,3) |
385 |
|
|
WRITE(KULOUT,'('' FCM_Q_AL (:,1) ='',2(1X,ES12.5))') FCM_Q_AL (:,1) |
386 |
|
|
WRITE(KULOUT,'('' FCM_Q_AL (:,2) ='',2(1X,ES12.5))') FCM_Q_AL (:,2) |
387 |
|
|
WRITE(KULOUT,'('' FCM_Q_AL (:,3) ='',2(1X,ES12.5))') FCM_Q_AL (:,3) |
388 |
|
|
WRITE(KULOUT,'('' FCM_Q_DI (:,1) ='',2(1X,ES12.5))') FCM_Q_DI (:,1) |
389 |
|
|
WRITE(KULOUT,'('' FCM_Q_DI (:,2) ='',2(1X,ES12.5))') FCM_Q_DI (:,2) |
390 |
|
|
WRITE(KULOUT,'('' FCM_Q_DI (:,3) ='',2(1X,ES12.5))') FCM_Q_DI (:,3) |
391 |
|
|
WRITE(KULOUT,'('' FCM_Q_DL (:,1) ='',2(1X,ES12.5))') FCM_Q_DL (:,1) |
392 |
|
|
WRITE(KULOUT,'('' FCM_Q_DL (:,2) ='',2(1X,ES12.5))') FCM_Q_DL (:,2) |
393 |
|
|
WRITE(KULOUT,'('' FCM_Q_DL (:,3) ='',2(1X,ES12.5))') FCM_Q_DL (:,3) |
394 |
|
|
WRITE(KULOUT,'('' FCM_Q_GI (:,1) ='',2(1X,ES12.5))') FCM_Q_GI (:,1) |
395 |
|
|
WRITE(KULOUT,'('' FCM_Q_GI (:,2) ='',2(1X,ES12.5))') FCM_Q_GI (:,2) |
396 |
|
|
WRITE(KULOUT,'('' FCM_Q_GI (:,3) ='',2(1X,ES12.5))') FCM_Q_GI (:,3) |
397 |
|
|
WRITE(KULOUT,'('' FCM_Q_GL (:,1) ='',2(1X,ES12.5))') FCM_Q_GL (:,1) |
398 |
|
|
WRITE(KULOUT,'('' FCM_Q_GL (:,2) ='',2(1X,ES12.5))') FCM_Q_GL (:,2) |
399 |
|
|
WRITE(KULOUT,'('' FCM_Q_GL (:,3) ='',2(1X,ES12.5))') FCM_Q_GL (:,3) |
400 |
|
|
ENDIF |
401 |
|
|
! ------------------------------------------------------------------ |
402 |
|
|
|
403 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SUPHY3',1,ZHOOK_HANDLE) |
404 |
|
1 |
END SUBROUTINE SUPHY3 |