source:
trunk/procs/alpha.pro
@
12
| Last change on this file since 12 was 2, checked in by post_it, 17 years ago | |
|---|---|
| File size: 208 bytes | |
| Line | |
|---|---|
| 1 | FUNCTION alpha, t, s, p |
| 2 | ; |
| 3 | ; compute alpha=-1/rho (d rho / d T) |
| 4 | ; |
| 5 | dt = 0.05 |
| 6 | siga = eos(t, s)-1000. |
| 7 | sigb = eos(t+dt, s)-1000. |
| 8 | |
| 9 | al = -0.001*(sigb-siga)/dt/(1.+1.e-3*siga) |
| 10 | |
| 11 | return, al |
| 12 | |
| 13 | END |
Note: See TracBrowser
for help on using the repository browser.
