MoDeNa  1.0
Software framework facilitating sequential multi-scale modelling
eos_const.f90
1 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
2 ! SUBROUTINE eos_const
3 !
4 ! This subroutine provides the constants of the PC-SAFT EOS.
5 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
6 
7 SUBROUTINE eos_const ( ap, bp )
8 
9  IMPLICIT NONE
10 
11  !-----------------------------------------------------------------------------
12  REAL, INTENT(OUT) :: ap(0:6,3)
13  REAL, INTENT(OUT) :: bp(0:6,3)
14  !-----------------------------------------------------------------------------
15 
16 
17  !-----------------------------------------------------------------------------
18  ! dispersion term constants
19  !-----------------------------------------------------------------------------
20 
21  ap(0,1) = 0.91056314451539
22  ap(0,2) = -0.30840169182720
23  ap(0,3) = -0.09061483509767
24  ap(1,1) = 0.63612814494991
25  ap(1,2) = 0.18605311591713
26  ap(1,3) = 0.45278428063920
27  ap(2,1) = 2.68613478913903
28  ap(2,2) = -2.50300472586548
29  ap(2,3) = 0.59627007280101
30  ap(3,1) = -26.5473624914884
31  ap(3,2) = 21.4197936296668
32  ap(3,3) = -1.72418291311787
33  ap(4,1) = 97.7592087835073
34  ap(4,2) = -65.2558853303492
35  ap(4,3) = -4.13021125311661
36  ap(5,1) = -159.591540865600
37  ap(5,2) = 83.3186804808856
38  ap(5,3) = 13.7766318697211
39  ap(6,1) = 91.2977740839123
40  ap(6,2) = -33.7469229297323
41  ap(6,3) = -8.67284703679646
42 
43  bp(0,1) = 0.72409469413165
44  bp(0,2) = -0.57554980753450
45  bp(0,3) = 0.09768831158356
46  bp(1,1) = 1.11913959304690 *2.0
47  bp(1,2) = 0.34975477607218 *2.0
48  bp(1,3) = -0.12787874908050 *2.0
49  bp(2,1) = -1.33419498282114 *3.0
50  bp(2,2) = 1.29752244631769 *3.0
51  bp(2,3) = -3.05195205099107 *3.0
52  bp(3,1) = -5.25089420371162 *4.0
53  bp(3,2) = -4.30386791194303 *4.0
54  bp(3,3) = 5.16051899359931 *4.0
55  bp(4,1) = 5.37112827253230 *5.0
56  bp(4,2) = 38.5344528930499 *5.0
57  bp(4,3) = -7.76088601041257 *5.0
58  bp(5,1) = 34.4252230677698 *6.0
59  bp(5,2) = -26.9710769414608 *6.0
60  bp(5,3) = 15.6044623461691 *6.0
61  bp(6,1) = -50.8003365888685 *7.0
62  bp(6,2) = -23.6010990650801 *7.0
63  bp(6,3) = -4.23812936930675 *7.0
64 
65 
66  !-----------------------------------------------------------------------------
67  ! square-well fluid
68  !-----------------------------------------------------------------------------
69 
70  ! ap(1,1)= 0.79152347258784
71  ! ap(1,2)= -0.62269805320654
72  ! ap(1,3)= -0.06798823934067
73  ! ap(2,1)= 1.07120982251709
74  ! ap(2,2)= 0.48628215731716
75  ! ap(2,3)= 0.02837828512515
76  ! ap(3,1)= 0.92084839459226
77  ! ap(3,2)= 1.11652038059747
78  ! ap(3,3)= 0.09713202077943
79  ! ap(4,1)= -7.84708350369249
80  ! ap(4,2)= -2.04200599876547
81  ! ap(4,3)= 0.06475764015088
82  ! ap(5,1)= 25.90284137818050
83  ! ap(5,2)= 9.27791640100603
84  ! ap(5,3)= 0.07729792971827
85  ! ap(6,1)= -57.1528726997640
86  ! ap(6,2)= -16.8377999920957
87  ! ap(6,3)= 0.24883598436184
88  ! ap(7,1)= 42.02314637860930
89  ! ap(7,2)= 7.62432635016420
90  ! ap(7,3)= -0.72472024688888
91 
92  ! bp(1,1)= 0.79152347258784
93  ! bp(1,2)= -0.62269805320654
94  ! bp(1,3)= -0.06798823934067
95  ! bp(2,1)= 1.07120982251709 *2.0
96  ! bp(2,2)= 0.48628215731716 *2.0
97  ! bp(2,3)= 0.02837828512515 *2.0
98  ! bp(3,1)= 0.92084839459226 *3.0
99  ! bp(3,2)= 1.11652038059747 *3.0
100  ! bp(3,3)= 0.09713202077943 *3.0
101  ! bp(4,1)= -7.84708350369249 *4.0
102  ! bp(4,2)= -2.04200599876547 *4.0
103  ! bp(4,3)= 0.06475764015088 *4.0
104  ! bp(5,1)= 25.90284137818050 *5.0
105  ! bp(5,2)= 9.27791640100603 *5.0
106  ! bp(5,3)= 0.07729792971827 *5.0
107  ! bp(6,1)= -57.1528726997640 *6.0
108  ! bp(6,2)= -16.8377999920957 *6.0
109  ! bp(6,3)= 0.24883598436184 *6.0
110  ! bp(7,1)= 42.02314637860930 *7.0
111  ! bp(7,2)= 7.62432635016420 *7.0
112  ! bp(7,3)= -0.72472024688888 *7.0
113 
114 
115 END SUBROUTINE eos_const
116 
117 
118 
119 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
120 ! SUBROUTINE dq_const
121 !
122 ! This subr. provides the constants of the dipole-quadrupole term.
123 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
124 
125 SUBROUTINE dq_const ( dqp2, dqp3, dqp4 )
126 
127  USE parameters, ONLY: nc
128  USE eos_variables, ONLY: ncomp, parame
129  IMPLICIT NONE
130 
131  !-----------------------------------------------------------------------------
132  REAL, INTENT(OUT) :: dqp2(nc,nc,0:4)
133  REAL, INTENT(OUT) :: dqp3(nc,nc,nc,0:4)
134  REAL, INTENT(OUT) :: dqp4(nc,nc,0:4)
135 
136  !-----------------------------------------------------------------------------
137  INTEGER :: i, j, k
138  REAL :: mdq(nc)
139  REAL :: mf1, mf2, msegij
140  !-----------------------------------------------------------------------------
141 
142  DO i=1,ncomp
143  mdq(i) = parame(i,1)
144  IF (mdq(i) > 2.0) mdq(i) = 2.0
145  END DO
146 
147 
148  DO i=1,ncomp
149  DO j=1,ncomp
150 
151  msegij = ( mdq(i) * mdq(j) )**0.5
152  mf1 = ( msegij - 1.0 ) / msegij
153  mf2 = mf1 * ( msegij - 2.0 ) / msegij
154 
155  dqp2(i,j,0) = 0.697094963 + mf1*(-0.673459279) + mf2*0.670340770
156  dqp2(i,j,1) = -0.633554144 + mf1*(-1.425899106) + mf2*(-4.338471826)
157  dqp2(i,j,2) = 2.945509028 + mf1 * 4.19441392 + mf2*7.234168360
158  dqp2(i,j,3) = -1.467027314 + mf1 * 1.0266216
159  dqp2(i,j,4) = 0.0
160 
161  dqp4(i,j,0) = -0.484038322 + mf1 * 0.67651011 + mf2*(-1.167560146)
162  dqp4(i,j,1) = 1.970405465 + mf1*(-3.013867512) + mf2*2.13488432
163  dqp4(i,j,2) = -2.118572671 + mf1 * 0.46742656
164  dqp4(i,j,3) = 0.0
165  dqp4(i,j,4) = 0.0
166 
167 
168  DO k=1,ncomp
169  msegij = ( mdq(i) * mdq(j) * mdq(k) )**(1.0/3.0)
170  mf1 = ( msegij - 1.0 ) / msegij
171  mf2 = ( msegij - 2.0 ) / msegij
172  dqp3(i,j,k,0) = 0.795009692 + mf1*(-2.099579397)
173  dqp3(i,j,k,1) = 3.386863396 + mf1*(-5.941376392)
174  dqp3(i,j,k,2) = 0.475106328 + mf1*(-0.178820384)
175  dqp3(i,j,k,3) = 0.0
176  dqp3(i,j,k,4) = 0.0
177  END DO
178 
179  END DO
180  END DO
181 
182 END SUBROUTINE dq_const
183 
184 
185 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
186 ! SUBROUTINE dd_const
187 !
188 ! This subroutine provides the constants of the dipole-term.
189 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
190 
191 SUBROUTINE dd_const ( ddp2, ddp3, ddp4 )
192 
193  USE parameters, ONLY: nc, pi
194  USE eos_variables, ONLY: ncomp, parame
195  IMPLICIT NONE
196 
197  !-----------------------------------------------------------------------------
198  REAL, INTENT(OUT) :: ddp2(nc,nc,0:4)
199  REAL, INTENT(OUT) :: ddp3(nc,nc,nc,0:4)
200  REAL, INTENT(OUT) :: ddp4(nc,nc,0:4)
201 
202  !-----------------------------------------------------------------------------
203  INTEGER :: i, j, k
204  REAL :: pardd(nc)
205  REAL :: mf1,mf2,msegij,sin2t
206  !-----------------------------------------------------------------------------
207 
208  sin2t = sin( 0.0 * pi / 180.0 )
209  sin2t = sin2t*sin2t
210 
211  DO i = 1, ncomp
212  pardd(i) = parame(i,1)
213  IF (pardd(i) > 2.0) pardd(i) = 2.0
214  END DO
215 
216  DO i=1,ncomp
217  DO j=1,ncomp
218  ! IF (parame(i,6) /= 0.0.AND.parame(j,6) /= 0.0) THEN
219 
220  msegij = ( pardd(i) * pardd(j) )**0.5
221  mf1 = ( msegij - 1.0 ) / msegij
222  mf2 = mf1 * ( msegij - 2.0 ) / msegij
223 
224  ddp2(i,j,0) = 0.30435038064 + mf1*(0.95346405973+0.201436*sin2t) &
225  + mf2*(-1.16100802773-1.74114*sin2t)
226  ddp2(i,j,1) = -0.13585877707 + mf1*(-1.83963831920+1.31649*sin2t) + mf2*4.52586067320
227  ddp2(i,j,2) = 1.44933285154 + mf1 * 2.01311801180 + mf2*0.97512223853
228  ddp2(i,j,3) = 0.35569769252 + mf1*(-7.37249576667) + mf2*(-12.2810377713)
229  ddp2(i,j,4) = -2.06533084541 + mf1 * 8.23741345333 + mf2*5.93975747420
230 
231  ddp4(i,j,0) = 0.21879385627 + mf1*(-0.58731641193) + mf2*3.48695755800
232  ddp4(i,j,1) = -1.18964307357 + mf1 * 1.24891317047 + mf2*(-14.9159739347)
233  ddp4(i,j,2) = 1.16268885692 + mf1*(-0.50852797392) + mf2*15.3720218600
234  ddp4(i,j,3) = 0.0
235  ddp4(i,j,4) = 0.0
236 
237  DO k=1,ncomp
238  ! IF (parame(k,6) /= 0.0) THEN
239  msegij = ( pardd(i) * pardd(j) * pardd(k) )**(1.0/3.0)
240  mf1 = ( msegij - 1.0 ) / msegij
241  mf2 = mf1 * ( msegij - 2.0 ) / msegij
242  ddp3(i,j,k,0) = -0.06467735252 + mf1*(-0.95208758351+0.28503*sin2t) &
243  + mf2*(-0.62609792333+2.2195*sin2t)
244  ddp3(i,j,k,1) = 0.19758818347 + mf1 * 2.99242575222 + mf2*1.29246858189
245  ddp3(i,j,k,2) = -0.80875619458 + mf1*(-2.38026356489) + mf2*1.65427830900
246  ddp3(i,j,k,3) = 0.69028490492 + mf1*(-0.27012609786) + mf2*(-3.43967436378)
247  ddp3(i,j,k,4) = 0.0
248 
249  ! ENDIF
250  END DO
251 
252  ! ENDIF
253  END DO
254  END DO
255 
256 END SUBROUTINE dd_const
257 
258 
259 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
260 ! SUBROUTINE qq_const
261 !
262 ! This subroutine provides the constants of the quadrupole-term.
263 !WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
264 
265 SUBROUTINE qq_const ( qqp2,qqp3,qqp4 )
266 
267  USE parameters, ONLY: nc
268  USE eos_variables, ONLY: ncomp, parame
269  IMPLICIT NONE
270 
271  !-----------------------------------------------------------------------------
272  REAL, INTENT(OUT) :: qqp2(nc,nc,0:4)
273  REAL, INTENT(OUT) :: qqp3(nc,nc,nc,0:4)
274  REAL, INTENT(OUT) :: qqp4(nc,nc,0:4)
275 
276  !-----------------------------------------------------------------------------
277  INTEGER :: i, j, k
278  REAL :: mqq(nc)
279  REAL :: mf1, mf2, msegij
280  !-----------------------------------------------------------------------------
281 
282  DO i = 1,ncomp
283  mqq(i) = parame(i,1)
284  IF (mqq(i) > 2.0) mqq(i) = 2.0
285  END DO
286 
287  DO i = 1,ncomp
288  DO j = 1,ncomp
289  IF (parame(i,7) /= 0.0 .AND. parame(j,7) /= 0.0) THEN
290 
291  msegij = ( mqq(i) * mqq(j) )**0.5
292  ! msegij = ( parame(i,1) * parame(j,1) )**0.50
293  mf1 = ( msegij - 1.0 ) / msegij
294  mf2 = mf1 * ( msegij - 2.0 ) / msegij
295 
296  qqp2(i,j,0) = 1.237830788 + mf1 * 1.285410878 + mf2*1.794295401
297  qqp2(i,j,1) = 2.435503144 + mf1*(-11.46561451) + mf2*0.769510293
298  qqp2(i,j,2) = 1.633090469 + mf1 *22.08689285 + mf2*7.264792255
299  qqp2(i,j,3) = -1.611815241 + mf1 * 7.46913832 + mf2*94.48669892
300  qqp2(i,j,4) = 6.977118504 + mf1*(-17.19777208) + mf2*(-77.1484579)
301 
302  qqp4(i,j,0) = 0.454271755 + mf1*(-0.813734006) + mf2*6.868267516
303  qqp4(i,j,1) = -4.501626435 + mf1 * 10.06402986 + mf2*(-5.173223765)
304  qqp4(i,j,2) = 3.585886783 + mf1*(-10.87663092) + mf2*(-17.2402066)
305  qqp4(i,j,3) = 0.0
306  qqp4(i,j,4) = 0.0
307 
308  DO k = 1,ncomp
309  IF (parame(k,7) /= 0.0) THEN
310  msegij = ( mqq(i) * mqq(j) * mqq(k) )**(1.0/3.0)
311  ! msegij = ( parame(i,1)*parame(j,1)*parame(k,1) )**(1.0/3.0)
312  mf1 = ( msegij - 1.0 ) / msegij
313  mf2 = mf1 * ( msegij - 2.0 ) / msegij
314  qqp3(i,j,k,0) = -0.500043713 + mf1 * 2.000209381 + mf2*3.135827145
315  qqp3(i,j,k,1) = 6.531869153 + mf1*(-6.78386584) + mf2*7.247588801
316  qqp3(i,j,k,2) = -16.01477983 + mf1 * 20.38324603 + mf2*3.075947834
317  qqp3(i,j,k,3) = 14.42597018 + mf1*(-10.89598394)
318  qqp3(i,j,k,4) = 0.0
319  END IF
320  END DO
321 
322  END IF
323  END DO
324  END DO
325 
326 END SUBROUTINE qq_const
327 
328 
329 
330 !!$!WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
331 !!$!
332 !!$!WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW
333 !!$
334 !!$SUBROUTINE dd_fischer (c_dd,n_dd,m_dd,k_dd,o_dd)
335 !!$
336 !!$ IMPLICIT NONE
337 !!$
338 !!$ !-----------------------------------------------------------------------------
339 !!$ REAL, INTENT(OUT) :: c_dd(28)
340 !!$ REAL, INTENT(OUT) :: n_dd(28)
341 !!$ REAL, INTENT(OUT) :: m_dd(28)
342 !!$ REAL, INTENT(OUT) :: k_dd(28)
343 !!$ REAL, INTENT(OUT) :: o_dd(28)
344 !!$
345 !!$ !-----------------------------------------------------------------------------
346 !!$ REAL, DIMENSION(17) :: c_qq, n_qq, m_qq, k_qq, o_qq
347 !!$ !-----------------------------------------------------------------------------
348 !!$
349 !!$ c_dd(1) = - 0.423652173318E-1
350 !!$ c_dd(2) = 0.204459397242E-1
351 !!$ c_dd(3) = 0.664266837321E-1
352 !!$ c_dd(4) = - 0.324168341478E-1
353 !!$ c_dd(5) = - 0.741263275720E-2
354 !!$ c_dd(6) = - 0.160855507113E-1
355 !!$ c_dd(7) = 0.435623305093E-2
356 !!$ c_dd(8) = - 0.105933370736E-3
357 !!$ c_dd(9) = - 0.132000046519E-5
358 !!$ c_dd(10) = 0.838157718194E-5
359 !!$ c_dd(11) = 0.109144074057E-1
360 !!$ c_dd(12) = 0.257960188278E-1
361 !!$ c_dd(13) = - 0.544140085185E-3
362 !!$ c_dd(14) = 0.349568484468E-2
363 !!$ c_dd(15) = - 0.421407562467E-1
364 !!$ c_dd(16) = - 0.745992658113E-2
365 !!$ c_dd(17) = 0.146102252152E-3
366 !!$ c_dd(18) = 0.566611094911E-3
367 !!$ c_dd(19) = - 0.378643890614E-2
368 !!$ c_dd(20) = - 0.365824539450E-1
369 !!$ c_dd(21) = 0.169287932475E-1
370 !!$ c_dd(22) = 0.663866480778E-2
371 !!$ c_dd(23) = 0.294409406715E-1
372 !!$ c_dd(24) = - 0.112110434947E-1
373 !!$ c_dd(25) = - 0.182144939032E-5
374 !!$ c_dd(26) = 0.758594753989E-7
375 !!$ c_dd(27) = - 0.216942306418E-4
376 !!$ c_dd(28) = - 0.274025042954E-5
377 !!$
378 !!$ n_dd(1) = -5.0
379 !!$ n_dd(2) = -8.0
380 !!$ n_dd(3) = -4.0
381 !!$ n_dd(4) = -3.0
382 !!$ n_dd(5) = -10.0
383 !!$ n_dd(6) = -7.0
384 !!$ n_dd(7) = -10.0
385 !!$ n_dd(8) = -11.0
386 !!$ n_dd(9) = -15.0
387 !!$ n_dd(10) = -10.0
388 !!$ n_dd(11) = -2.0
389 !!$ n_dd(12) = -2.0
390 !!$ n_dd(13) = -1.0
391 !!$ n_dd(14) = -5.0
392 !!$ n_dd(15) = -3.0
393 !!$ n_dd(16) = -1.0
394 !!$ n_dd(17) = 1.0
395 !!$ n_dd(18) = -9.0
396 !!$ n_dd(19) = -7.0
397 !!$ n_dd(20) = -2.0
398 !!$ n_dd(21) = -1.0
399 !!$ n_dd(22) = -5.0
400 !!$ n_dd(23) = -2.0
401 !!$ n_dd(24) = -1.0
402 !!$ n_dd(25) = -8.0
403 !!$ n_dd(26) = -5.0
404 !!$ n_dd(27) = 1.0
405 !!$ n_dd(28) = -4.0
406 !!$
407 !!$ m_dd(1) = 2.0
408 !!$ m_dd(2) = 2.0
409 !!$ m_dd(3) = 2.0
410 !!$ m_dd(4) = 2.0
411 !!$ m_dd(5) = 2.0
412 !!$ m_dd(6) = 2.0
413 !!$ m_dd(7) = 2.0
414 !!$ m_dd(8) = 2.0
415 !!$ m_dd(9) = 2.0
416 !!$ m_dd(10) = 3.0
417 !!$ m_dd(11) = 2.0
418 !!$ m_dd(12) = 3.0
419 !!$ m_dd(13) = 6.0
420 !!$ m_dd(14) = 2.0
421 !!$ m_dd(15) = 3.0
422 !!$ m_dd(16) = 3.0
423 !!$ m_dd(17) = 6.0
424 !!$ m_dd(18) = 2.0
425 !!$ m_dd(19) = 2.0
426 !!$ m_dd(20) = 2.0
427 !!$ m_dd(21) = 2.0
428 !!$ m_dd(22) = 3.0
429 !!$ m_dd(23) = 3.0
430 !!$ m_dd(24) = 3.0
431 !!$ m_dd(25) = 10.0
432 !!$ m_dd(26) = 16.0
433 !!$ m_dd(27) = 4.0
434 !!$ m_dd(28) = 9.0
435 !!$
436 !!$ k_dd(1) = 5.0
437 !!$ k_dd(2) = 6.0
438 !!$ k_dd(3) = 7.0
439 !!$ k_dd(4) = 7.0
440 !!$ k_dd(5) = 9.0
441 !!$ k_dd(6) = 9.0
442 !!$ k_dd(7) = 11.0
443 !!$ k_dd(8) = 15.0
444 !!$ k_dd(9) = 18.0
445 !!$ k_dd(10) = 18.0
446 !!$ k_dd(11) = 5.0
447 !!$ k_dd(12) = 5.0
448 !!$ k_dd(13) = 5.0
449 !!$ k_dd(14) = 6.0
450 !!$ k_dd(15) = 6.0
451 !!$ k_dd(16) = 6.0
452 !!$ k_dd(17) = 6.0
453 !!$ k_dd(18) = 7.0
454 !!$ k_dd(19) = 7.0
455 !!$ k_dd(20) = 7.0
456 !!$ k_dd(21) = 7.0
457 !!$ k_dd(22) = 7.0
458 !!$ k_dd(23) = 7.0
459 !!$ k_dd(24) = 7.0
460 !!$ k_dd(25) = 7.0
461 !!$ k_dd(26) = 7.0
462 !!$ k_dd(27) = 8.0
463 !!$ k_dd(28) = 10.0
464 !!$
465 !!$ o_dd(1) = 1.0
466 !!$ o_dd(2) = 1.0
467 !!$ o_dd(3) = 1.0
468 !!$ o_dd(4) = 1.0
469 !!$ o_dd(5) = 1.0
470 !!$ o_dd(6) = 1.0
471 !!$ o_dd(7) = 1.0
472 !!$ o_dd(8) = 1.0
473 !!$ o_dd(9) = 1.0
474 !!$ o_dd(10) = 1.0
475 !!$ o_dd(11) = 0.0
476 !!$ o_dd(12) = 0.0
477 !!$ o_dd(13) = 0.0
478 !!$ o_dd(14) = 0.0
479 !!$ o_dd(15) = 0.0
480 !!$ o_dd(16) = 0.0
481 !!$ o_dd(17) = 0.0
482 !!$ o_dd(18) = 0.0
483 !!$ o_dd(19) = 0.0
484 !!$ o_dd(20) = 0.0
485 !!$ o_dd(21) = 0.0
486 !!$ o_dd(22) = 0.0
487 !!$ o_dd(23) = 0.0
488 !!$ o_dd(24) = 0.0
489 !!$ o_dd(25) = 0.0
490 !!$ o_dd(26) = 0.0
491 !!$ o_dd(27) = 0.0
492 !!$ o_dd(28) = 0.0
493 !!$
494 !!$ c_qq(1) = - 0.412154280896E-2
495 !!$ c_qq(2) = 0.355780441736E-2
496 !!$ c_qq(3) = - 0.888093798389E-3
497 !!$ c_qq(4) = 0.973791559609E-4
498 !!$ c_qq(5) = - 0.604233719326E-7
499 !!$ c_qq(6) = - 0.304478633146E-4
500 !!$ c_qq(7) = - 0.378930196337E-3
501 !!$ c_qq(8) = - 0.275388267352E-1
502 !!$ c_qq(9) = 0.118301888420E-1
503 !!$ c_qq(10) = - 0.283451230562E-2
504 !!$ c_qq(11) = - 0.567703873212E-4
505 !!$ c_qq(12) = 0.314708573212E-2
506 !!$ c_qq(13) = 0.963786052569E-3
507 !!$ c_qq(14) = - 0.127591002424E-2
508 !!$ c_qq(15) = 0.363746463238E-3
509 !!$ c_qq(16) = 0.301067943096E-4
510 !!$ c_qq(17) = 0.291778231128E-6
511 !!$
512 !!$ n_qq(1) = -8.0
513 !!$ n_qq(2) = -6.0
514 !!$ n_qq(3) = -4.0
515 !!$ n_qq(4) = -10.0
516 !!$ n_qq(5) = -20.0
517 !!$ n_qq(6) = -8.0
518 !!$ n_qq(7) = -3.0
519 !!$ n_qq(8) = -3.0
520 !!$ n_qq(9) = -2.0
521 !!$ n_qq(10) = 0.0
522 !!$ n_qq(11) = -5.0
523 !!$ n_qq(12) = -1.0
524 !!$ n_qq(13) = -3.0
525 !!$ n_qq(14) = -1.0
526 !!$ n_qq(15) = 0.0
527 !!$ n_qq(16) = 0.0
528 !!$ n_qq(17) = -10.0
529 !!$
530 !!$ m_qq(1) = 2.0
531 !!$ m_qq(2) = 2.0
532 !!$ m_qq(3) = 2.0
533 !!$ m_qq(4) = 2.0
534 !!$ m_qq(5) = 2.0
535 !!$ m_qq(6) = 2.0
536 !!$ m_qq(7) = 8.0
537 !!$ m_qq(8) = 2.0
538 !!$ m_qq(9) = 2.0
539 !!$ m_qq(10) = 2.0
540 !!$ m_qq(11) = 8.0
541 !!$ m_qq(12) = 2.0
542 !!$ m_qq(13) = 5.0
543 !!$ m_qq(14) = 5.0
544 !!$ m_qq(15) = 5.0
545 !!$ m_qq(16) = 8.0
546 !!$ m_qq(17) = 7.0
547 !!$
548 !!$ k_qq(1) = 11.0
549 !!$ k_qq(2) = 12.0
550 !!$ k_qq(3) = 13.0
551 !!$ k_qq(4) = 16.0
552 !!$ k_qq(5) = 19.0
553 !!$ k_qq(6) = 20.0
554 !!$ k_qq(7) = 7.0
555 !!$ k_qq(8) = 8.0
556 !!$ k_qq(9) = 8.0
557 !!$ k_qq(10) = 8.0
558 !!$ k_qq(11) = 8.0
559 !!$ k_qq(12) = 9.0
560 !!$ k_qq(13) = 10.0
561 !!$ k_qq(14) = 10.0
562 !!$ k_qq(15) = 10.0
563 !!$ k_qq(16) = 10.0
564 !!$ k_qq(17) = 18.0
565 !!$
566 !!$ o_qq(1) = 1.0
567 !!$ o_qq(2) = 1.0
568 !!$ o_qq(3) = 1.0
569 !!$ o_qq(4) = 1.0
570 !!$ o_qq(5) = 1.0
571 !!$ o_qq(6) = 1.0
572 !!$ o_qq(7) = 0.0
573 !!$ o_qq(8) = 0.0
574 !!$ o_qq(9) = 0.0
575 !!$ o_qq(10) = 0.0
576 !!$ o_qq(11) = 0.0
577 !!$ o_qq(12) = 0.0
578 !!$ o_qq(13) = 0.0
579 !!$ o_qq(14) = 0.0
580 !!$ o_qq(15) = 0.0
581 !!$ o_qq(16) = 0.0
582 !!$ o_qq(17) = 0.0
583 !!$
584 !!$END SUBROUTINE dd_fischer
585 
parameters
WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW This module contains constant...
Definition: modules.f90:6
eos_variables
WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW This module contains paramete...
Definition: modules.f90:120