1 |
|
2 |
#include "ctrparam.h" |
3 |
|
4 |
! ========================================================== |
5 |
! |
6 |
! VFUNCT.F: "unknown" |
7 |
! |
8 |
! ---------------------------------------------------------- |
9 |
! |
10 |
! Revision History: |
11 |
! |
12 |
! When Who What |
13 |
! ----- ---------- ------- |
14 |
! 080100 Chien Wang repack based on M24x11, add cpp |
15 |
! |
16 |
! ========================================================== |
17 |
|
18 |
subroutine vfunct(v,z) |
19 |
|
20 |
zm=11000. |
21 |
d=2.71364 |
22 |
a=1.27 |
23 |
c=50900. |
24 |
zz=abs(z-zm) |
25 |
v=d*exp(-(zz**a/c)) |
26 |
|
27 |
return |
28 |
end |