-
Notifications
You must be signed in to change notification settings - Fork 0
/
WaterData.f
116 lines (116 loc) · 7.7 KB
/
WaterData.f
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
subroutine Water(l_lnk,n_lnk,k_lnk,nlam)
IMPLICIT NONE
integer nlam,j
real*8 l_lnk(*),n_lnk(*),k_lnk(*)
real*8 l0( 169),n0( 169),k0( 169)
data (l0(j),j=1, 169) /
& 0.20000E+00, 0.22500E+00, 0.25000E+00, 0.27500E+00, 0.30000E+00,
& 0.32500E+00, 0.35000E+00, 0.37500E+00, 0.40000E+00, 0.42500E+00,
& 0.45000E+00, 0.47500E+00, 0.50000E+00, 0.52500E+00, 0.55000E+00,
& 0.57500E+00, 0.60000E+00, 0.62500E+00, 0.65000E+00, 0.67500E+00,
& 0.70000E+00, 0.72500E+00, 0.75000E+00, 0.77500E+00, 0.80000E+00,
& 0.82500E+00, 0.85000E+00, 0.87500E+00, 0.90000E+00, 0.92500E+00,
& 0.95000E+00, 0.97500E+00, 0.10000E+01, 0.12000E+01, 0.14000E+01,
& 0.16000E+01, 0.18000E+01, 0.20000E+01, 0.22000E+01, 0.24000E+01,
& 0.26000E+01, 0.26500E+01, 0.27000E+01, 0.27500E+01, 0.28000E+01,
& 0.28500E+01, 0.29000E+01, 0.29500E+01, 0.30000E+01, 0.30500E+01,
& 0.31000E+01, 0.31500E+01, 0.32000E+01, 0.32500E+01, 0.33000E+01,
& 0.33500E+01, 0.34000E+01, 0.34500E+01, 0.35000E+01, 0.36000E+01,
& 0.37000E+01, 0.38000E+01, 0.39000E+01, 0.40000E+01, 0.41000E+01,
& 0.42000E+01, 0.43000E+01, 0.44000E+01, 0.45000E+01, 0.46000E+01,
& 0.47000E+01, 0.48000E+01, 0.49000E+01, 0.50000E+01, 0.51000E+01,
& 0.52000E+01, 0.53000E+01, 0.54000E+01, 0.55000E+01, 0.56000E+01,
& 0.57000E+01, 0.58000E+01, 0.59000E+01, 0.60000E+01, 0.61000E+01,
& 0.62000E+01, 0.63000E+01, 0.64000E+01, 0.65000E+01, 0.66000E+01,
& 0.67000E+01, 0.68000E+01, 0.69000E+01, 0.70000E+01, 0.71000E+01,
& 0.72000E+01, 0.73000E+01, 0.74000E+01, 0.75000E+01, 0.76000E+01,
& 0.77000E+01, 0.78000E+01, 0.79000E+01, 0.80000E+01, 0.82000E+01,
& 0.84000E+01, 0.86000E+01, 0.88000E+01, 0.90000E+01, 0.92000E+01,
& 0.94000E+01, 0.96000E+01, 0.98000E+01, 0.10000E+02, 0.10500E+02,
& 0.11000E+02, 0.11500E+02, 0.12000E+02, 0.12500E+02, 0.13000E+02,
& 0.13500E+02, 0.14000E+02, 0.14500E+02, 0.15000E+02, 0.15500E+02,
& 0.16000E+02, 0.16500E+02, 0.17000E+02, 0.17500E+02, 0.18000E+02,
& 0.18500E+02, 0.19000E+02, 0.19500E+02, 0.20000E+02, 0.21000E+02,
& 0.22000E+02, 0.23000E+02, 0.24000E+02, 0.25000E+02, 0.26000E+02,
& 0.27000E+02, 0.28000E+02, 0.29000E+02, 0.30000E+02, 0.32000E+02,
& 0.34000E+02, 0.36000E+02, 0.38000E+02, 0.40000E+02, 0.42000E+02,
& 0.44000E+02, 0.46000E+02, 0.48000E+02, 0.50000E+02, 0.60000E+02,
& 0.70000E+02, 0.80000E+02, 0.90000E+02, 0.10000E+03, 0.11000E+03,
& 0.12000E+03, 0.13000E+03, 0.14000E+03, 0.15000E+03, 0.16000E+03,
& 0.17000E+03, 0.18000E+03, 0.19000E+03, 0.20000E+03 /
data (n0(j),j=1, 169) /
& 0.13960E+01, 0.13730E+01, 0.13620E+01, 0.13540E+01, 0.13490E+01,
& 0.13460E+01, 0.13430E+01, 0.13410E+01, 0.13390E+01, 0.13380E+01,
& 0.13370E+01, 0.13360E+01, 0.13350E+01, 0.13340E+01, 0.13330E+01,
& 0.13330E+01, 0.13320E+01, 0.13320E+01, 0.13310E+01, 0.13310E+01,
& 0.13310E+01, 0.13300E+01, 0.13300E+01, 0.13300E+01, 0.13290E+01,
& 0.13290E+01, 0.13290E+01, 0.13280E+01, 0.13280E+01, 0.13280E+01,
& 0.13270E+01, 0.13270E+01, 0.13270E+01, 0.13240E+01, 0.13210E+01,
& 0.13170E+01, 0.13120E+01, 0.13060E+01, 0.12960E+01, 0.12790E+01,
& 0.12420E+01, 0.12190E+01, 0.11880E+01, 0.11570E+01, 0.11420E+01,
& 0.11490E+01, 0.12010E+01, 0.12920E+01, 0.13710E+01, 0.14260E+01,
& 0.14670E+01, 0.14830E+01, 0.14780E+01, 0.14670E+01, 0.14500E+01,
& 0.14320E+01, 0.14200E+01, 0.14100E+01, 0.14000E+01, 0.13850E+01,
& 0.13740E+01, 0.13640E+01, 0.13570E+01, 0.13510E+01, 0.13460E+01,
& 0.13420E+01, 0.13380E+01, 0.13340E+01, 0.13320E+01, 0.13300E+01,
& 0.13300E+01, 0.13300E+01, 0.13280E+01, 0.13250E+01, 0.13220E+01,
& 0.13170E+01, 0.13120E+01, 0.13050E+01, 0.12980E+01, 0.12890E+01,
& 0.12770E+01, 0.12620E+01, 0.12480E+01, 0.12650E+01, 0.13190E+01,
& 0.13630E+01, 0.13570E+01, 0.13470E+01, 0.13390E+01, 0.13340E+01,
& 0.13290E+01, 0.13240E+01, 0.13210E+01, 0.13170E+01, 0.13140E+01,
& 0.13120E+01, 0.13090E+01, 0.13070E+01, 0.13040E+01, 0.13020E+01,
& 0.12990E+01, 0.12970E+01, 0.12940E+01, 0.12910E+01, 0.12860E+01,
& 0.12810E+01, 0.12750E+01, 0.12690E+01, 0.12620E+01, 0.12550E+01,
& 0.12470E+01, 0.12390E+01, 0.12290E+01, 0.12180E+01, 0.11850E+01,
& 0.11530E+01, 0.11260E+01, 0.11110E+01, 0.11230E+01, 0.11460E+01,
& 0.11770E+01, 0.12100E+01, 0.12410E+01, 0.12700E+01, 0.12970E+01,
& 0.13250E+01, 0.13510E+01, 0.13760E+01, 0.14010E+01, 0.14230E+01,
& 0.14430E+01, 0.14610E+01, 0.14760E+01, 0.14800E+01, 0.14870E+01,
& 0.15000E+01, 0.15110E+01, 0.15210E+01, 0.15310E+01, 0.15390E+01,
& 0.15450E+01, 0.15490E+01, 0.15510E+01, 0.15510E+01, 0.15460E+01,
& 0.15360E+01, 0.15270E+01, 0.15220E+01, 0.15190E+01, 0.15220E+01,
& 0.15300E+01, 0.15410E+01, 0.15550E+01, 0.15870E+01, 0.17030E+01,
& 0.18210E+01, 0.18860E+01, 0.19240E+01, 0.19570E+01, 0.19660E+01,
& 0.20040E+01, 0.20360E+01, 0.20560E+01, 0.20690E+01, 0.20810E+01,
& 0.20940E+01, 0.21070E+01, 0.21190E+01, 0.21300E+01 /
data (k0(j),j=1, 169) /
& 0.11000E-06, 0.49000E-07, 0.33500E-07, 0.23500E-07, 0.16000E-07,
& 0.10800E-07, 0.65000E-08, 0.35000E-08, 0.18600E-08, 0.13000E-08,
& 0.10200E-08, 0.93500E-09, 0.10000E-08, 0.13200E-08, 0.19600E-08,
& 0.36000E-08, 0.10900E-07, 0.13900E-07, 0.16400E-07, 0.22300E-07,
& 0.33500E-07, 0.91500E-07, 0.15600E-06, 0.14800E-06, 0.12500E-06,
& 0.18200E-06, 0.29300E-06, 0.39100E-06, 0.48600E-06, 0.10600E-05,
& 0.29300E-05, 0.34800E-05, 0.28900E-05, 0.98900E-05, 0.13800E-03,
& 0.85500E-04, 0.11500E-03, 0.11000E-02, 0.28900E-03, 0.95600E-03,
& 0.31700E-02, 0.67000E-02, 0.19000E-01, 0.59000E-01, 0.11500E+00,
& 0.18500E+00, 0.26800E+00, 0.29800E+00, 0.27200E+00, 0.24000E+00,
& 0.19200E+00, 0.13500E+00, 0.92400E-01, 0.61000E-01, 0.36800E-01,
& 0.26100E-01, 0.19500E-01, 0.13200E-01, 0.94000E-02, 0.51500E-02,
& 0.36000E-02, 0.34000E-02, 0.38000E-02, 0.46000E-02, 0.56200E-02,
& 0.68800E-02, 0.84500E-02, 0.10300E-01, 0.13400E-01, 0.14700E-01,
& 0.15700E-01, 0.15000E-01, 0.13700E-01, 0.12400E-01, 0.11100E-01,
& 0.10100E-01, 0.98000E-02, 0.10300E-01, 0.11600E-01, 0.14200E-01,
& 0.20300E-01, 0.33000E-01, 0.62200E-01, 0.10700E+00, 0.13100E+00,
& 0.88000E-01, 0.57000E-01, 0.44900E-01, 0.39200E-01, 0.35600E-01,
& 0.33700E-01, 0.32700E-01, 0.32200E-01, 0.32000E-01, 0.32000E-01,
& 0.32100E-01, 0.32200E-01, 0.32400E-01, 0.32600E-01, 0.32800E-01,
& 0.33100E-01, 0.33500E-01, 0.33900E-01, 0.34300E-01, 0.35100E-01,
& 0.36100E-01, 0.37200E-01, 0.38500E-01, 0.39900E-01, 0.41500E-01,
& 0.43300E-01, 0.45400E-01, 0.47900E-01, 0.50800E-01, 0.66200E-01,
& 0.96800E-01, 0.14200E+00, 0.19900E+00, 0.25900E+00, 0.30500E+00,
& 0.34300E+00, 0.37000E+00, 0.38800E+00, 0.40200E+00, 0.41400E+00,
& 0.42200E+00, 0.42800E+00, 0.42900E+00, 0.42900E+00, 0.42600E+00,
& 0.42100E+00, 0.41400E+00, 0.40400E+00, 0.39300E+00, 0.38200E+00,
& 0.37300E+00, 0.36700E+00, 0.36100E+00, 0.35600E+00, 0.35000E+00,
& 0.34400E+00, 0.33800E+00, 0.33300E+00, 0.32800E+00, 0.32400E+00,
& 0.32900E+00, 0.34300E+00, 0.36100E+00, 0.38500E+00, 0.40900E+00,
& 0.43600E+00, 0.46200E+00, 0.48800E+00, 0.51400E+00, 0.58700E+00,
& 0.57600E+00, 0.54700E+00, 0.53600E+00, 0.53200E+00, 0.53100E+00,
& 0.52600E+00, 0.51400E+00, 0.50000E+00, 0.49500E+00, 0.49600E+00,
& 0.49700E+00, 0.49900E+00, 0.50100E+00, 0.50400E+00 /
nlam= 169
l_lnk(1:nlam)=l0(1:nlam)
n_lnk(1:nlam)=n0(1:nlam)
k_lnk(1:nlam)=k0(1:nlam)
return
end