source: LMDZ6/branches/LMDZ-COSP/libf/phylmd/rrtm/suphy3.F90 @ 3847

Last change on this file since 3847 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: 13.8 KB
Line 
1!OPTIONS XOPT(NOEVAL)
2SUBROUTINE 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
53USE PARKIND1  ,ONLY : JPIM     ,JPRB
54USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
55
56! Ce qui concerne NULNAM commente par MPL le 15.04.09
57!USE YOMLUN   , ONLY : NULNAM
58USE YOMCT0B  , ONLY : LECMWF
59USE YOMPHY   , ONLY : LCLSATUR
60USE 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
74IMPLICIT NONE
75
76INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT
77REAL(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
89IF (LHOOK) CALL DR_HOOK('SUPHY3',0,ZHOOK_HANDLE)
90BSFSA=.3490_JPRB
91BSFSI=.3238_JPRB
92BSFSN=.3300_JPRB
93BSFTA=.3471_JPRB
94BSFTI=.3231_JPRB
95BSFTN=.3584_JPRB
96EARRT=0.001324_JPRB
97EOASA=.6727E-01_JPRB
98EOASI=.1692E-01_JPRB
99EOASN=.3353E-01_JPRB
100EOATA=.5037E-01_JPRB
101EOATI=.3455E+01_JPRB
102EOATN=.9077E+01_JPRB
103EODSA=.3665E+00_JPRB
104EODSI=.1406E+01_JPRB
105EODSN=.6558E+01_JPRB
106EODTA=.1766E-01_JPRB
107EODTI=.8435E+00_JPRB
108EODTN=.4319E+01_JPRB
109EORAY=.8606E-06_JPRB
110GCA(1)=.8041E-01_JPRB
111GCA(2)=.1456E+00_JPRB
112GCA(3)=.4787E+01_JPRB
113GCA(4)=.2102E+04_JPRB
114GCA(5)=.1334E+01_JPRB
115GCA(6)=.1551E-04_JPRB
116GCB(1)=.8968E+07_JPRB
117GCB(2)=.2413E+10_JPRB
118GCB(3)=.3548E+05_JPRB
119GCB(4)=.6370E+10_JPRB
120GCB(5)=.8499E+11_JPRB
121GCB(6)=.1012E+06_JPRB
122GCC(1)=.5925E-10_JPRB
123GCC(2)=.1842E-10_JPRB
124GCC(3)=.2532E-07_JPRB
125GCC(4)=.1953E+07_JPRB
126GCC(5)=.1734E-11_JPRB
127GCC(6)=.1225E-16_JPRB
128GCD4=.3608E-69_JPRB
129GCE4=.7563E+04_JPRB
130GIREC1=0.0_JPRB
131GIREC2=0.0_JPRB
132GIREC3=0.0_JPRB
133GIREC4=0.0_JPRB
134QCO2=.5366E-03_JPRB
135QLIMI=15000._JPRB
136QLIP0=8.E+06_JPRB
137USAA=-.3020_JPRB
138USAI=-.3524_JPRB
139USAN=-.3400_JPRB
140USBA=0._JPRB
141USBI=0._JPRB
142USBN=0._JPRB
143VDP(1,1)=.21868E+02_JPRB
144VDP(2,1)=.17453E+02_JPRB
145VDP(3,1)=.68918E+00_JPRB
146VDP(4,1)=.23456E-03_JPRB
147VDP(5,1)=.22317E-09_JPRB
148VDP(1,2)=.48401E+02_JPRB
149VDP(2,2)=.18648E+02_JPRB
150VDP(3,2)=.35199E-01_JPRB
151VDP(4,2)=.63691E-06_JPRB
152VDP(5,2)=.86395E-13_JPRB
153VDP(1,3)=.76948E+02_JPRB
154VDP(2,3)=.41056E+01_JPRB
155VDP(3,3)=.42667E-03_JPRB
156VDP(4,3)=0._JPRB
157VDP(5,3)=0._JPRB
158VDP(1,4)=.27241E+03_JPRB
159VDP(2,4)=.57091E+04_JPRB
160VDP(3,4)=.14393E+05_JPRB
161VDP(4,4)=.29879E+04_JPRB
162VDP(5,4)=.25382E+02_JPRB
163VDP(1,5)=.86942E+02_JPRB
164VDP(2,5)=.32186E+03_JPRB
165VDP(3,5)=.10775E+03_JPRB
166VDP(4,5)=.21261E+01_JPRB
167VDP(5,5)=.40003E-02_JPRB
168VDP(1,6)=.24408E+02_JPRB
169VDP(2,6)=.81919E+01_JPRB
170VDP(3,6)=.72193E+01_JPRB
171VDP(4,6)=.56230E+00_JPRB
172VDP(5,6)=.25384E-02_JPRB
173VNP(1,1)=.69926E+01_JPRB
174VNP(2,1)=.63915E+00_JPRB
175VNP(3,1)=.28896E-03_JPRB
176VNP(4,1)=.10050E-08_JPRB
177VNP(5,1)=.99037E-16_JPRB
178VNP(1,2)=.31105E+01_JPRB
179VNP(2,2)=.14225E-01_JPRB
180VNP(3,2)=.69355E-06_JPRB
181VNP(4,2)=.36087E-12_JPRB
182VNP(5,2)=.44113E-20_JPRB
183VNP(1,3)=.23659E+01_JPRB
184VNP(2,3)=.11139E-02_JPRB
185VNP(3,3)=.10618E-05_JPRB
186VNP(4,3)=0._JPRB
187VNP(5,3)=0._JPRB
188VNP(1,4)=.19810E+03_JPRB
189VNP(2,4)=.46954E+04_JPRB
190VNP(3,4)=.22512E+04_JPRB
191VNP(4,4)=.52461E+02_JPRB
192VNP(5,4)=.11645E+00_JPRB
193VNP(1,5)=.46348E+02_JPRB
194VNP(2,5)=.35630E+02_JPRB
195VNP(3,5)=.33005E+01_JPRB
196VNP(4,5)=.18045E-01_JPRB
197VNP(5,5)=.88667E-05_JPRB
198VNP(1,6)=.47413E+01_JPRB
199VNP(2,6)=.16334E+01_JPRB
200VNP(3,6)=.48164E+00_JPRB
201VNP(4,6)=.56140E-02_JPRB
202VNP(5,6)=.67790E-04_JPRB
203
204! Default values for cloud model - first index denotes spectral band:
205!   1 - solar
206!   2 - thermal
207FCM_DEL_A(1)  =  9.55496e-01_JPRB
208FCM_DEL_A(2)  =  1.79152e+02_JPRB
209FCM_DEL_D(1)  =  5.72959e+01_JPRB
210FCM_DEL_D(2)  =  2.00159e+02_JPRB
211FCM_MU_A(1)   =  8.66200e-01_JPRB
212FCM_MU_A(2)   =  6.65532e-01_JPRB
213FCM_MU_D(1)   =  1.00000e+00_JPRB
214FCM_MU_D(2)   =  6.70517e-01_JPRB
215FCM_N_I       =  1.00000e-01_JPRB
216FCM_N_L       =  1.00000e-01_JPRB
217FCM_P_AI(1,0) =  1.48920e+01_JPRB
218FCM_P_AI(1,1) = -3.71379e+01_JPRB
219FCM_P_AI(1,2) =  0.00000e+00_JPRB
220FCM_P_AI(1,3) =  0.00000e+00_JPRB
221FCM_P_AI(2,0) =  9.76920e+00_JPRB
222FCM_P_AI(2,1) = -7.54840e+00_JPRB
223FCM_P_AI(2,2) =  0.00000e+00_JPRB
224FCM_P_AI(2,3) =  0.00000e+00_JPRB
225FCM_P_AL(1,0) =  1.42062e+00_JPRB
226FCM_P_AL(1,1) = -3.54037e+00_JPRB
227FCM_P_AL(1,2) =  1.73535e+00_JPRB
228FCM_P_AL(1,3) =  0.00000e+00_JPRB
229FCM_P_AL(2,0) =  5.08522e+00_JPRB
230FCM_P_AL(2,1) = -1.08144e+01_JPRB
231FCM_P_AL(2,2) =  5.16331e+00_JPRB
232FCM_P_AL(2,3) =  0.00000e+00_JPRB
233FCM_P_DI(1,0) =  1.30487e+01_JPRB
234FCM_P_DI(1,1) =  1.98670e+01_JPRB
235FCM_P_DI(1,2) = -5.08578e+01_JPRB
236FCM_P_DI(1,3) =  0.00000e+00_JPRB
237FCM_P_DI(2,0) =  1.23173e+01_JPRB
238FCM_P_DI(2,1) = -3.32808e+01_JPRB
239FCM_P_DI(2,2) =  5.77295e+01_JPRB
240FCM_P_DI(2,3) = -7.35728e+01_JPRB
241FCM_P_DL(1,0) =  4.68646e+00_JPRB
242FCM_P_DL(1,1) = -8.77652e+00_JPRB
243FCM_P_DL(1,2) =  2.12789e+00_JPRB
244FCM_P_DL(1,3) =  0.00000e+00_JPRB
245FCM_P_DL(2,0) =  4.52292e+00_JPRB
246FCM_P_DL(2,1) = -1.43327e+01_JPRB
247FCM_P_DL(2,2) =  1.16650e+01_JPRB
248FCM_P_DL(2,3) =  0.00000e+00_JPRB
249FCM_P_GI(1,0) = -2.01472e+00_JPRB
250FCM_P_GI(1,1) =  9.16284e+01_JPRB
251FCM_P_GI(1,2) =  7.82275e+01_JPRB
252FCM_P_GI(1,3) =  0.00000e+00_JPRB
253FCM_P_GI(2,0) = -7.55124e-01_JPRB
254FCM_P_GI(2,1) = -4.68837e+00_JPRB
255FCM_P_GI(2,2) =  6.72877e+01_JPRB
256FCM_P_GI(2,3) =  0.00000e+00_JPRB
257FCM_P_GL(1,0) =  8.82686e-01_JPRB
258FCM_P_GL(1,1) = -2.02445e+00_JPRB
259FCM_P_GL(1,2) =  1.45259e+00_JPRB
260FCM_P_GL(1,3) =  0.00000e+00_JPRB
261FCM_P_GL(2,0) =  1.93784e-01_JPRB
262FCM_P_GL(2,1) = -6.07676e-01_JPRB
263FCM_P_GL(2,2) =  5.37734e-01_JPRB
264FCM_P_GL(2,3) =  0.00000e+00_JPRB
265FCM_Q_AI(1,1) =  1.86588e+01_JPRB
266FCM_Q_AI(1,2) =  0.00000e+00_JPRB
267FCM_Q_AI(1,3) =  0.00000e+00_JPRB
268FCM_Q_AI(2,1) =  2.78063e+00_JPRB
269FCM_Q_AI(2,2) =  0.00000e+00_JPRB
270FCM_Q_AI(2,3) =  0.00000e+00_JPRB
271FCM_Q_AL(1,1) = -2.81210e+00_JPRB
272FCM_Q_AL(1,2) =  2.35758e+00_JPRB
273FCM_Q_AL(1,3) =  0.00000e+00_JPRB
274FCM_Q_AL(2,1) = -2.31252e+00_JPRB
275FCM_Q_AL(2,2) =  1.56848e+00_JPRB
276FCM_Q_AL(2,3) =  0.00000e+00_JPRB
277FCM_Q_DI(1,1) =  1.04582e+01_JPRB
278FCM_Q_DI(1,2) =  6.15641e+00_JPRB
279FCM_Q_DI(1,3) =  0.00000e+00_JPRB
280FCM_Q_DI(2,1) =  1.18370e+00_JPRB
281FCM_Q_DI(2,2) =  2.23187e+00_JPRB
282FCM_Q_DI(2,3) =  7.00424e+00_JPRB
283FCM_Q_DL(1,1) = -2.18205e+00_JPRB
284FCM_Q_DL(1,2) =  1.34831e+00_JPRB
285FCM_Q_DL(1,3) =  0.00000e+00_JPRB
286FCM_Q_DL(2,1) = -3.12310e+00_JPRB
287FCM_Q_DL(2,2) =  2.57959e+00_JPRB
288FCM_Q_DL(2,3) =  0.00000e+00_JPRB
289FCM_Q_GI(1,1) =  7.89999e+01_JPRB
290FCM_Q_GI(1,2) =  5.24929e+01_JPRB
291FCM_Q_GI(1,3) =  0.00000e+00_JPRB
292FCM_Q_GI(2,1) =  4.93801e+00_JPRB
293FCM_Q_GI(2,2) =  2.38096e+01_JPRB
294FCM_Q_GI(2,3) =  0.00000e+00_JPRB
295FCM_Q_GL(1,1) = -2.20830e+00_JPRB
296FCM_Q_GL(1,2) =  1.39832e+00_JPRB
297FCM_Q_GL(1,3) =  0.00000e+00_JPRB
298FCM_Q_GL(2,1) = -3.27770e+00_JPRB
299FCM_Q_GL(2,2) =  2.78446e+00_JPRB
300FCM_Q_GL(2,3) =  0.00000e+00_JPRB
301REXP_NEB      = 8.0_JPRB
302
303!        1.2 Modify default values according to LECMWF
304
305IF (LECMWF) THEN
306ELSE
307ENDIF
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
324WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY3 '')')
325WRITE(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 & BSFSA,BSFSI,BSFSN,BSFTA,BSFTI,BSFTN,EARRT,EOASA,EOASI,EOASN,&
345 & EOATA,EOATI,EOATN,EODSA,EODSI,EODSN,EODTA,EODTI,EODTN,EORAY,&
346 & GCA,GCB,GCC,GCD4,GCE4,QCO2,QLIP0,USAA,USAI,USAN,USBA,USBI,USBN,&
347 & VDP,VNP,QLIMI,GIREC1,GIREC2,GIREC3,GIREC4 
348
349IF (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)
400ENDIF
401!     ------------------------------------------------------------------
402
403IF (LHOOK) CALL DR_HOOK('SUPHY3',1,ZHOOK_HANDLE)
404END SUBROUTINE SUPHY3
Note: See TracBrowser for help on using the repository browser.