3 REAL(KIND=8),
parameter ::
zpdh2o = 0.8
4 REAL(KIND=8),
parameter ::
zpdumg = 0.75
5 REAL(KIND=8),
parameter ::
zprh2o = 30000.0
6 REAL(KIND=8),
parameter ::
zprumg = 30000.0
7 REAL(KIND=8),
parameter ::
rtdh2o = 0.40
8 REAL(KIND=8),
parameter ::
rtdumg = 0.375
9 REAL(KIND=8),
parameter ::
rth2o = 240.0
10 REAL(KIND=8),
parameter ::
rtumg = 240.0
12 REAL(KIND=8),
dimension(2),
parameter ::
wg1 = (/1.0, 1.0/)
13 REAL(KIND=8),
dimension(11),
parameter ::
tintp = (/ 187.5, 200., 212.5, 225., 237.5, 250., 262.5, 275., 287.5, 300., 312.5 /)
15 real(kind=8),
dimension(11,16,3),
parameter ::
ga = reshape ( (/ &
16 0.63499072e-02, 0.65566348e-02, 0.67849730e-02, 0.70481947e-02, 0.73585943e-02, 0.77242818e-02, 0.81472693e-02, 0.86227527e-02,&
17 0.91396814e-02, 0.96825438e-02, 0.10233955e-01, 0.77266491e-02, 0.81323287e-02, 0.86507620e-02, 0.92776391e-02, 0.99806312e-02,&
18 0.10709803e-01, 0.11414739e-01, 0.12058772e-01, 0.12623992e-01, 0.13108146e-01, 0.13518390e-01, 0.11644593e+01, 0.11747203e+01,&
19 0.11837872e+01, 0.11918561e+01, 0.11990757e+01, 0.12055643e+01, 0.12114186e+01, 0.12167192e+01, 0.12215344e+01, 0.12259226e+01,&
20 0.12299344e+01, 0.12006968e+01, 0.12108196e+01, 0.12196717e+01, 0.12274493e+01, 0.12343189e+01, 0.12404147e+01, 0.12458431e+01,&
21 0.12506907e+01, 0.12550299e+01, 0.12589256e+01, 0.12624402e+01, 0.15750172e+00, 0.16174076e+00, 0.16548628e+00, 0.16881124e+00,&
22 0.17177839e+00, 0.17443933e+00, 0.17683622e+00, 0.17900375e+00, 0.18097099e+00, 0.18276283e+00, 0.18440117e+00, 0.17770551e+00,&
23 0.18176757e+00, 0.18527967e+00, 0.18833348e+00, 0.19100108e+00, 0.19334122e+00, 0.19540288e+00, 0.19722732e+00, 0.19884918e+00,&
24 0.20029696e+00, 0.20159300e+00, 0.10192131e+02, 0.97258602e+01, 0.92992890e+01, 0.89154021e+01, 0.85730084e+01, 0.82685838e+01,&
25 0.79978921e+01, 0.77568055e+01, 0.75416266e+01, 0.73491694e+01, 0.71767400e+01, 0.92439050e+01, 0.87567422e+01, 0.83270144e+01,&
26 0.79528337e+01, 0.76286839e+01, 0.73477879e+01, 0.71035818e+01, 0.68903312e+01, 0.67032875e+01, 0.65386461e+01, 0.63934377e+01,&
27 0.24870635e+02, 0.24725591e+02, 0.24600320e+02, 0.24487300e+02, 0.24384935e+02, 0.24292341e+02, 0.24208572e+02, 0.24132642e+02,&
28 0.24063614e+02, 0.24000649e+02, 0.23943021e+02, 0.24586283e+02, 0.24441465e+02, 0.24311657e+02, 0.24196167e+02, 0.24093406e+02,&
29 0.24001597e+02, 0.23919098e+02, 0.23844511e+02, 0.23776708e+02, 0.23714816e+02, 0.23658197e+02, 0.11990218e+02, 0.10904073e+02,&
30 0.89126838e+01, 0.85622405e+01, 0.94892164e+01, 0.13580937e+02,-0.32050918e+03,-0.37133165e+01, 0.18890836e+00, 0.14209226e+01,&
31 0.19817679e+01, 0.79709806e+01, 0.75400737e+01, 0.81804377e+01, 0.10564339e+02, 0.46896789e+02,-0.30926524e+01, 0.85742941e+00,&
32 0.19164038e+01, 0.23513199e+01, 0.25566644e+01, 0.26555181e+01, 0.87668459e-01, 0.83754276e-01, 0.80460283e-01, 0.77659686e-01,&
33 0.75257056e-01, 0.73179175e-01, 0.71369063e-01, 0.69781812e-01, 0.68381606e-01, 0.67139539e-01, 0.66032012e-01, 0.74878820e-01,&
34 0.71650966e-01, 0.68979615e-01, 0.66745345e-01, 0.64857571e-01, 0.63248495e-01, 0.61866970e-01, 0.60673632e-01, 0.59637277e-01,&
35 0.58732178e-01, 0.57936092e-01, 0.13230067e+02, 0.13213564e+02, 0.13209140e+02, 0.13213894e+02, 0.13225963e+02, 0.13243806e+02,&
36 0.13266104e+02, 0.13291782e+02, 0.13319961e+02, 0.13349927e+02, 0.13381108e+02, 0.13183816e+02, 0.13189991e+02, 0.13209485e+02,&
37 0.13238789e+02, 0.13275017e+02, 0.13316096e+02, 0.13360555e+02, 0.13407324e+02, 0.13455544e+02, 0.13504450e+02, 0.13553282e+02,&
38 -0.99506586e-03,-0.10184169e-02,-0.10404730e-02,-0.10621792e-02,-0.10847662e-02,-0.11094726e-02,-0.11372949e-02,-0.11687683e-02,&
39 -0.12038314e-02,-0.12418367e-02,-0.12817135e-02,-0.11661515e-02,-0.11886130e-02,-0.12139929e-02,-0.12445811e-02,-0.12807672e-02,&
40 -0.13208251e-02,-0.13619034e-02,-0.14014165e-02,-0.14378639e-02,-0.14708488e-02,-0.15006791e-02, 0.41243390e+00, 0.43407282e+00,&
41 0.45331413e+00, 0.47048604e+00, 0.48586286e+00, 0.49968044e+00, 0.51214132e+00, 0.52341830e+00, 0.53365803e+00, 0.54298448e+00,&
42 0.55150227e+00, 0.48318936e+00, 0.50501827e+00, 0.52409502e+00, 0.54085277e+00, 0.55565422e+00, 0.56878618e+00, 0.58047395e+00,&
43 0.59089894e+00, 0.60021475e+00, 0.60856112e+00, 0.61607594e+00,-0.22159303e-01,-0.22748917e-01,-0.23269898e-01,-0.23732392e-01,&
44 -0.24145123e-01,-0.24515269e-01,-0.24848690e-01,-0.25150210e-01,-0.25423873e-01,-0.25673139e-01,-0.25901055e-01,-0.24972399e-01,&
45 -0.25537247e-01,-0.26025624e-01,-0.26450280e-01,-0.26821236e-01,-0.27146657e-01,-0.27433354e-01,-0.27687065e-01,-0.27912608e-01,&
46 -0.28113944e-01,-0.28294180e-01, 0.80737799e+01, 0.79171158e+01, 0.77609605e+01, 0.76087371e+01, 0.74627112e+01, 0.73239981e+01,&
47 0.71929934e+01, 0.70697065e+01, 0.69539626e+01, 0.68455144e+01, 0.67441020e+01, 0.77425778e+01, 0.75443460e+01, 0.73526151e+01,&
48 0.71711188e+01, 0.70015571e+01, 0.68442532e+01, 0.66987996e+01, 0.65644820e+01, 0.64405267e+01, 0.63262376e+01, 0.62210701e+01,&
49 0.10542131e+02, 0.10515895e+02, 0.10492949e+02, 0.10472049e+02, 0.10452961e+02, 0.10435562e+02, 0.10419710e+02, 0.10405247e+02,&
50 0.10392022e+02, 0.10379892e+02, 0.10368736e+02, 0.10490353e+02, 0.10463512e+02, 0.10439183e+02, 0.10417324e+02, 0.10397704e+02,&
51 0.10380038e+02, 0.10364052e+02, 0.10349509e+02, 0.10336215e+02, 0.10324018e+02, 0.10312808e+02,-0.12823142e+01,-0.10571588e+01,&
52 -0.74864953e+00,-0.58705980e+00,-0.49305772e+00,-0.51461431e+00, 0.12373350e+02, 0.44809588e+00, 0.46548918e+00, 0.59121475e+00,&
53 0.74676119e+00,-0.74805226e+00,-0.56252739e+00,-0.46188072e+00,-0.40712065e+00,-0.15295996e+01, 0.43555255e+00, 0.50380874e+00,&
54 0.68537352e+00, 0.89437630e+00, 0.11127003e+01, 0.13329782e+01, 0.13845511e+01, 0.13187042e+01, 0.12644396e+01, 0.12191543e+01,&
55 0.11809511e+01, 0.11484154e+01, 0.11204723e+01, 0.10962918e+01, 0.10752229e+01, 0.10567474e+01, 0.10404465e+01, 0.11718758e+01,&
56 0.11216131e+01, 0.10809473e+01, 0.10476396e+01, 0.10200373e+01, 0.99692726e+00, 0.97740923e+00, 0.96080188e+00, 0.94657562e+00,&
57 0.93430511e+00, 0.92363528e+00, 0.22042132e+02, 0.22107298e+02, 0.22180915e+02, 0.22259478e+02, 0.22341039e+02, 0.22424247e+02,&
58 0.22508089e+02, 0.22591771e+02, 0.22674661e+02, 0.22756246e+02, 0.22836093e+02, 0.22169501e+02, 0.22270075e+02, 0.22379193e+02,&
59 0.22492992e+02, 0.22608508e+02, 0.22723843e+02, 0.22837837e+02, 0.22949751e+02, 0.23059032e+02, 0.23165146e+02, 0.23267456e+02,&
60 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
61 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
62 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
63 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
64 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
65 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
66 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
67 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
68 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
69 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
70 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
71 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
72 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
73 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
74 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
75 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
76 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
77 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
78 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
79 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
80 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,&
81 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00 /) &
84 real(kind=8),
dimension(11,16,3),
parameter ::
gb = reshape ( (/ &
85 0.63499072e-02, 0.65566348e-02, 0.67849730e-02, 0.70481947e-02, 0.73585943e-02, 0.77242818e-02, 0.81472693e-02, 0.86227527e-02,&
86 0.91396814e-02, 0.96825438e-02, 0.10233955e-01, 0.77266491e-02, 0.81323287e-02, 0.86507620e-02, 0.92776391e-02, 0.99806312e-02,&
87 0.10709803e-01, 0.11414739e-01, 0.12058772e-01, 0.12623992e-01, 0.13108146e-01, 0.13518390e-01, 0.11644593e+01, 0.11747203e+01,&
88 0.11837872e+01, 0.11918561e+01, 0.11990757e+01, 0.12055643e+01, 0.12114186e+01, 0.12167192e+01, 0.12215344e+01, 0.12259226e+01,&
89 0.12299344e+01, 0.12006968e+01, 0.12108196e+01, 0.12196717e+01, 0.12274493e+01, 0.12343189e+01, 0.12404147e+01, 0.12458431e+01,&
90 0.12506907e+01, 0.12550299e+01, 0.12589256e+01, 0.12624402e+01, 0.15750172e+00, 0.16174076e+00, 0.16548628e+00, 0.16881124e+00,&
91 0.17177839e+00, 0.17443933e+00, 0.17683622e+00, 0.17900375e+00, 0.18097099e+00, 0.18276283e+00, 0.18440117e+00, 0.17770551e+00,&
92 0.18176757e+00, 0.18527967e+00, 0.18833348e+00, 0.19100108e+00, 0.19334122e+00, 0.19540288e+00, 0.19722732e+00, 0.19884918e+00,&
93 0.20029696e+00, 0.20159300e+00, 0.10192131e+02, 0.97258602e+01, 0.92992890e+01, 0.89154021e+01, 0.85730084e+01, 0.82685838e+01,&
94 0.79978921e+01, 0.77568055e+01, 0.75416266e+01, 0.73491694e+01, 0.71767400e+01, 0.92439050e+01, 0.87567422e+01, 0.83270144e+01,&
95 0.79528337e+01, 0.76286839e+01, 0.73477879e+01, 0.71035818e+01, 0.68903312e+01, 0.67032875e+01, 0.65386461e+01, 0.63934377e+01,&
96 0.24870635e+02, 0.24725591e+02, 0.24600320e+02, 0.24487300e+02, 0.24384935e+02, 0.24292341e+02, 0.24208572e+02, 0.24132642e+02,&
97 0.24063614e+02, 0.24000649e+02, 0.23943021e+02, 0.24586283e+02, 0.24441465e+02, 0.24311657e+02, 0.24196167e+02, 0.24093406e+02,&
98 0.24001597e+02, 0.23919098e+02, 0.23844511e+02, 0.23776708e+02, 0.23714816e+02, 0.23658197e+02, 0.11990218e+02, 0.10904073e+02,&
99 0.89126838e+01, 0.85622405e+01, 0.94892164e+01, 0.13580937e+02,-0.32050918e+03,-0.37133165e+01, 0.18890836e+00, 0.14209226e+01,&
100 0.19817679e+01, 0.79709806e+01, 0.75400737e+01, 0.81804377e+01, 0.10564339e+02, 0.46896789e+02,-0.30926524e+01, 0.85742941e+00,&
101 0.19164038e+01, 0.23513199e+01, 0.25566644e+01, 0.26555181e+01, 0.87668459e-01, 0.83754276e-01, 0.80460283e-01, 0.77659686e-01,&
102 0.75257056e-01, 0.73179175e-01, 0.71369063e-01, 0.69781812e-01, 0.68381606e-01, 0.67139539e-01, 0.66032012e-01, 0.74878820e-01,&
103 0.71650966e-01, 0.68979615e-01, 0.66745345e-01, 0.64857571e-01, 0.63248495e-01, 0.61866970e-01, 0.60673632e-01, 0.59637277e-01,&
104 0.58732178e-01, 0.57936092e-01, 0.13230067e+02, 0.13213564e+02, 0.13209140e+02, 0.13213894e+02, 0.13225963e+02, 0.13243806e+02,&
105 0.13266104e+02, 0.13291782e+02, 0.13319961e+02, 0.13349927e+02, 0.13381108e+02, 0.13183816e+02, 0.13189991e+02, 0.13209485e+02,&
106 0.13238789e+02, 0.13275017e+02, 0.13316096e+02, 0.13360555e+02, 0.13407324e+02, 0.13455544e+02, 0.13504450e+02, 0.13553282e+02,&
107 0.97222852e-01, 0.98862238e-01, 0.10061504e+00, 0.10256222e+00, 0.10475952e+00, 0.10720986e+00, 0.10985370e+00, 0.11257633e+00,&
108 0.11522980e+00, 0.11766343e+00, 0.11975320e+00, 0.10681591e+00, 0.10921298e+00, 0.11198225e+00, 0.11487826e+00, 0.11751113e+00,&
109 0.11951535e+00, 0.12069945e+00, 0.12108524e+00, 0.12084229e+00, 0.12019005e+00, 0.11932684e+00, 0.10346097e+01, 0.10433655e+01,&
110 0.10511933e+01, 0.10582150e+01, 0.10645317e+01, 0.10702313e+01, 0.10753907e+01, 0.10800762e+01, 0.10843446e+01, 0.10882439e+01,&
111 0.10918144e+01, 0.10626130e+01, 0.10716026e+01, 0.10795108e+01, 0.10865006e+01, 0.10927103e+01, 0.10982489e+01, 0.11032019e+01,&
112 0.11076379e+01, 0.11116160e+01, 0.11151910e+01, 0.11184188e+01, 0.38103212e+00, 0.38913800e+00, 0.39613651e+00, 0.40222421e+00,&
113 0.40756010e+00, 0.41226954e+00, 0.41645142e+00, 0.42018474e+00, 0.42353379e+00, 0.42655211e+00, 0.42928533e+00, 0.41646579e+00,&
114 0.42345095e+00, 0.42937476e+00, 0.43444062e+00, 0.43880316e+00, 0.44258354e+00, 0.44587882e+00, 0.44876776e+00, 0.45131451e+00,&
115 0.45357095e+00, 0.45557797e+00, 0.82623280e+01, 0.81072291e+01, 0.79523834e+01, 0.78012527e+01, 0.76561458e+01, 0.75182174e+01,&
116 0.73878952e+01, 0.72652133e+01, 0.71500151e+01, 0.70420667e+01, 0.69411177e+01, 0.79342219e+01, 0.77373458e+01, 0.75467334e+01,&
117 0.73661786e+01, 0.71974319e+01, 0.70408543e+01, 0.68960649e+01, 0.67623672e+01, 0.66389989e+01, 0.65252707e+01, 0.64206412e+01,&
118 0.10656640e+02, 0.10630910e+02, 0.10608399e+02, 0.10587891e+02, 0.10569156e+02, 0.10552075e+02, 0.10536510e+02, 0.10522307e+02,&
119 0.10509317e+02, 0.10497402e+02, 0.10486443e+02, 0.10605856e+02, 0.10579514e+02, 0.10555632e+02, 0.10534169e+02, 0.10514900e+02,&
120 0.10497547e+02, 0.10481842e+02, 0.10467553e+02, 0.10454488e+02, 0.10442501e+02, 0.10431483e+02, 0.26681588e+02, 0.24728346e+02,&
121 0.20551342e+02, 0.19955244e+02, 0.22227100e+02, 0.31770288e+02,-0.74061287e+03,-0.81329826e+01, 0.90279822e+00, 0.37532746e+01,&
122 0.50437916e+01, 0.18377807e+02, 0.17643148e+02, 0.19296161e+02, 0.24951120e+02, 0.10957372e+03,-0.67432659e+01, 0.24550746e+01,&
123 0.49089917e+01, 0.59008712e+01, 0.63532616e+01, 0.65558627e+01, 0.23203798e+01, 0.22288925e+01, 0.21515593e+01, 0.20855896e+01,&
124 0.20288489e+01, 0.19796791e+01, 0.19367778e+01, 0.18991112e+01, 0.18658501e+01, 0.18363226e+01, 0.18099779e+01, 0.20206726e+01,&
125 0.19441824e+01, 0.18807257e+01, 0.18275618e+01, 0.17825910e+01, 0.17442308e+01, 0.17112809e+01, 0.16828137e+01, 0.16580908e+01,&
126 0.16365014e+01, 0.16175164e+01, 0.22051750e+02, 0.22116850e+02, 0.22190410e+02, 0.22268925e+02, 0.22350445e+02, 0.22433617e+02,&
127 0.22517429e+02, 0.22601086e+02, 0.22683956e+02, 0.22765522e+02, 0.22845354e+02, 0.22178972e+02, 0.22279484e+02, 0.22388551e+02,&
128 0.22502309e+02, 0.22617792e+02, 0.22733099e+02, 0.22847071e+02, 0.22958967e+02, 0.23068234e+02, 0.23174336e+02, 0.23276638e+02,&
129 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
130 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
131 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
132 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
133 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
134 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
135 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
136 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
137 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
138 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
139 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
140 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
141 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
142 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
143 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
144 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
145 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
146 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
147 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
148 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
149 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01,&
150 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01, 0.10000000e+01 /) &
153 real(kind=8),
dimension(6,6),
parameter ::
xp = reshape ( (/ &
154 0.46430621e+02, 0.12928299e+03, 0.20732648e+03, 0.31398411e+03, 0.18373177e+03, -0.11412303e+03, &
155 0.73604774e+02, 0.27887914e+03, 0.27076947e+03, -0.57322111e+02, -0.64742459e+02, 0.87238280e+02, &
156 0.37050866e+02, 0.20498759e+03, 0.37558029e+03, 0.17401171e+03, -0.13350302e+03, -0.37651795e+02, &
157 0.14930141e+02, 0.89161160e+02, 0.17793062e+03, 0.93433860e+02, -0.70646020e+02, -0.26373150e+02, &
158 0.40386780e+02, 0.10855270e+03, 0.50755010e+02, -0.31496190e+02, 0.12791300e+00, 0.18017770e+01, &
159 0.90811926e+01, 0.75073923e+02, 0.24654438e+03, 0.39332612e+03, 0.29385281e+03, 0.89107921e+02 /) , (/ 6,6 /) )
162 REAL(KIND=8),
dimension(2) ::
rsun = (/ 0.441676 , 0.558324 /)
163 REAL(KIND=8),
dimension(2,6),
parameter ::
rray = reshape ( &
164 (/ .428937e-01, .697200e-02,&
165 .890743e+00, .173297e-01,&
166 -.288555e+01,-.850903e-01,&
167 .522744e+01, .248261e+00,&
168 -.469173e+01,-.302031e+00,&
169 .161645e+01, .129662e+00 /) , (/2,6/) )
171 REAL(KIND=8),
dimension(2,5),
parameter ::
taua = reshape ( &
172 (/ 0.730719, 0.730719, 0.912819, 0.912819, 0.725059, &
173 0.725059, 0.745405, 0.745405, 0.682188, 0.682188 /),(/2,5/) )
174 REAL(KIND=8),
dimension(2,5),
parameter ::
rpiza = reshape ( &
175 (/ 0.872212, 0.872212, 0.982545, 0.982545, 0.623143, &
176 0.623143, 0.944887, 0.944887, 0.997975, 0.997975 /),(/2,5/) )
177 REAL(KIND=8),
dimension(2,5),
parameter ::
rcga = reshape ( &
178 (/ 0.647596, 0.647596, 0.739002, 0.739002, 0.580845, &
179 0.580845, 0.662657, 0.662657, 0.624246, 0.624246 /),(/2,5/) )
181 REAL(KIND=8),
dimension(2,3,7),
parameter ::
apad = reshape ( &
182 (/ 0.912418292e+05, 0.376655383e-08, 0.000000000e-00,&
183 0.739646016e-08, 0.925887084e-04, 0.410177786e+03,&
184 0.723613782e+05, 0.978576773e-04, 0.000000000e-00,&
185 0.131849595e-03, 0.129353723e-01, 0.672595424e+02,&
186 0.596037057e+04, 0.387714006e+00, 0.000000000e-00,&
187 0.437772681e+00, 0.800821928e+00, 0.000000000e-00,&
188 0.000000000e-00, 0.118461660e+03, 0.000000000e-00,&
189 0.151345118e+03, 0.242715973e+02, 0.000000000e-00,&
190 0.000000000e-00, 0.119079797e+04, 0.000000000e-00,&
191 0.233628890e+04, 0.878331486e+02, 0.000000000e-00,&
192 0.000000000e-00, 0.293353397e+03, 0.000000000e-00,&
193 0.797219934e+03, 0.191559725e+02, 0.000000000e-00,&
194 0.000000000e-00, 0.000000000e+00, 0.000000000e-00,&
195 0.000000000e+00, 0.000000000e+00, 0.000000000e+00 /) , (/2,3,7/) )
196 REAL(KIND=8),
dimension(2,3,7),
parameter ::
bpad = reshape ( &
197 (/ 0.912418292e+05, 0.376655383e-08, 0.000000000e-00,&
198 0.739646016e-08, 0.925887084e-04, 0.410177786e+03,&
199 0.724555318e+05, 0.979023421e-04, 0.000000000e-00,&
200 0.131861712e-03, 0.131812683e-01, 0.731185438e+02,&
201 0.602593328e+04, 0.388611139e+00, 0.000000000e-00,&
202 0.437949001e+00, 0.812706117e+00, 0.100000000e+01,&
203 0.100000000e+01, 0.120291383e+03, 0.000000000e-00,&
204 0.151692730e+03, 0.249863591e+02, 0.000000000e+00,&
205 0.000000000e-00, 0.130531005e+04, 0.000000000e-00,&
206 0.237071130e+04, 0.931071925e+02, 0.000000000e+00,&
207 0.000000000e-00, 0.415049409e+03, 0.000000000e-00,&
208 0.867914360e+03, 0.252233437e+02, 0.000000000e+00,&
209 0.000000000e-00, 0.100000000e+01, 0.000000000e-00,&
210 0.100000000e+01, 0.100000000e+01, 0.000000000e+00 /) , (/2,3,7/) )
211 REAL(KIND=8),
dimension(2,3),
parameter ::
d = reshape ( &
212 (/ 0.0, 0.0, 0.0, 0.0, 0.0, 0.8 /) , (/2,3/) )
214 REAL(KIND=8),
parameter ::
tref = 250.0
215 REAL(KIND=8),
dimension(2),
parameter ::
rt1 = (/ -0.577350269, +0.577350269 /)
216 REAL(KIND=8),
dimension(5,5),
parameter ::
raer= reshape ( &
217 (/ .038520, .037196, .040532, .054934, .038520 &
218 , .12613 , .18313 , .10357 , .064106, .126130 &
219 , .012579, .013649, .018652, .025181, .012579 &
220 , .011890, .016142, .021105, .028908, .011890 &
221 , .013792, .026810, .052203, .066338, .013792 /) , (/5,5/) )
223 REAL(KIND=8),
dimension(8,3),
parameter :: at= reshape ( &
224 (/ 0.298199e-02,0.143676e-01,0.197861e-01,0.289560e-01,&
225 0.103800e-01,0.868859e-02,0.250073e-03,0.307423e-01,&
226 -.394023e-03,0.366501e-02,0.315541e-02,-.208807e-02,&
227 0.436296e-02,-.972752e-03,0.455875e-03,0.110879e-02,&
228 0.319566e-04,-.160822e-02,-.174547e-02,-.121943e-02,&
229 -.161431e-02,0.000000e-00,0.109242e-03,-.322172e-03 /) , (/8,3/) )
231 REAL(KIND=8),
dimension(8,3),
parameter ::
bt= reshape ( &
232 (/ -0.106432e-04,-0.553979e-04,-0.877012e-04,-0.165960e-03,&
233 -0.276744e-04,-0.278412e-04, 0.199846e-05,-0.108482e-03,&
234 0.660324e-06,-0.101701e-04, 0.513302e-04, 0.157704e-03,&
235 -0.327381e-04,-0.713940e-06,-0.216313e-05, 0.258096e-05,&
236 0.174356e-06, 0.920868e-05, 0.523138e-06,-0.146427e-04,&
237 0.127646e-04 ,0.117469e-05, 0.175991e-06,-0.814575e-06 /) , (/8,3/) )
239 REAL(KIND=8),
dimension(4),
parameter ::
oct = (/ -.326e-03, -.102e-05, .137e-02, -.535e-05 /)
real(kind=8), dimension(5, 5), parameter raer
real(kind=8), dimension(11, 16, 3), parameter gb
real(kind=8), dimension(11), parameter tintp
real(kind=8), parameter rtumg
real(kind=8), dimension(11, 16, 3), parameter ga
real(kind=8), parameter zprumg
real(kind=8), dimension(4), parameter oct
real(kind=8), dimension(2, 5), parameter taua
real(kind=8), dimension(2) rsun
real(kind=8), parameter zpdh2o
real(kind=8), dimension(2, 3, 7), parameter apad
real(kind=8), dimension(2, 6), parameter rray
real(kind=8), dimension(2), parameter rt1
real(kind=8), dimension(2, 5), parameter rcga
real(kind=8), parameter zprh2o
real(kind=8), dimension(2, 3, 7), parameter bpad
real(kind=8), parameter rtdh2o
real(kind=8), dimension(2), parameter wg1
real(kind=8), parameter tref
real(kind=8), dimension(6, 6), parameter xp
real(kind=8), dimension(2, 5), parameter rpiza
real(kind=8), parameter rtdumg
real(kind=8), parameter zpdumg
real(kind=8), parameter rth2o
real(kind=8), dimension(2, 3), parameter d
real(kind=8), dimension(8, 3), parameter bt