71 |
bsec = 2650. + 5000. * slat2 |
bsec = 2650. + 5000. * slat2 |
72 |
a = 4. * bsec**1.5 * ppm**1.5 * (1. + (bsec / ps)**1.5) & |
a = 4. * bsec**1.5 * ppm**1.5 * (1. + (bsec / ps)**1.5) & |
73 |
/ (bsec**1.5 + ppm**1.5)**2 |
/ (bsec**1.5 + ppm**1.5)**2 |
74 |
aprim = max(0., (2.666666 * (1.74E-5 - 7.5E-6 * slat2 & |
asec = max(0., (gms - max(0., (2.666666 * (1.74E-5 - 7.5E-6 * slat2 & |
75 |
- 1.7E-6 * cost * slat) * ppm - a * gms) / (1. - a)) |
- 1.7E-6 * cost * slat) * ppm - a * gms) / (1. - a))) & |
76 |
asec = max(0., (gms - aprim) * (1. + (bsec / ps)**1.5)) |
* (1. + (bsec / ps)**1.5)) |
77 |
aprim = gms - asec / (1. + (bsec / ps)**1.5) |
aprim = gms - asec / (1. + (bsec / ps)**1.5) |
78 |
|
|
79 |
DO k = 1, llm |
DO k = 1, llm |