| 189 |
! difference equation for downdraft theta and mass |
! difference equation for downdraft theta and mass |
| 190 |
! flux from two simultaneous differential equations |
! flux from two simultaneous differential equations |
| 191 |
|
|
| 192 |
amfac = sigd * sigd * 70. * ph(il, i) & |
amfac = sigd**2 * 70. * ph(il, i) * (p(il, i - 1) - p(il, i)) & |
|
* (p(il, i - 1) - p(il, i)) & |
|
| 193 |
* (th(il, i) - th(il, i - 1)) / (tv(il, i) * th(il, i)) |
* (th(il, i) - th(il, i - 1)) / (tv(il, i) * th(il, i)) |
| 194 |
amp2 = abs(mp(il, i + 1) * mp(il, i + 1) - mp(il, i) & |
amp2 = abs(mp(il, i + 1)**2 - mp(il, i)**2) |
|
* mp(il, i)) |
|
| 195 |
|
|
| 196 |
if (amp2 > 0.1 * amfac) then |
if (amp2 > 0.1 * amfac) then |
| 197 |
xf = 100. * sigd * sigd * sigd * (ph(il, i) - ph(il, i + 1)) |
xf = 100. * sigd**3 * (ph(il, i) - ph(il, i + 1)) |
| 198 |
tf = b(il, i) - 5. * (th(il, i) - th(il, i - 1)) & |
tf = b(il, i) - 5. * (th(il, i) - th(il, i - 1)) & |
| 199 |
* t(il, i) / (lvcp(il, i) * sigd * th(il, i)) |
* t(il, i) / (lvcp(il, i) * sigd * th(il, i)) |
| 200 |
af = xf * tf + mp(il, i + 1) * mp(il, i + 1) * tinv |
af = xf * tf + mp(il, i + 1)**2 * tinv |
| 201 |
bf = 2. * (tinv * mp(il, i + 1))**3 + tinv & |
bf = 2. * (tinv * mp(il, i + 1))**3 + tinv & |
| 202 |
* mp(il, i + 1) * xf * tf + 50. * (p(il, i - 1) & |
* mp(il, i + 1) * xf * tf + 50. * (p(il, i - 1) & |
| 203 |
- p(il, i)) * xf * tevap |
- p(il, i)) * xf * tevap |
Parent Directory
| 
Revision Log
| 
Patch