5 SUBROUTINE f_pt1 ( fres )
11 REAL,
INTENT(OUT) :: fres
13 REAL :: i1_dft, i2_dft, order1
16 order1 = mseg(1)*mseg(1)*parame(1,2)**3 *parame(1,3) / t
18 CALL f_dft (i1_dft,i2_dft)
20 fres = 2.0 * pi * rho * i1_dft * order1
30 SUBROUTINE phipt1 ( my_pt1, p_pt1 )
36 REAL,
INTENT(OUT) :: my_pt1
37 REAL,
INTENT(IN OUT) :: p_pt1
40 REAL :: fres, pgesdz, pgesd2, pgesd3
48 CALL pressure_pt1 ( p_pt1, pgesdz, pgesd2, pgesd3 )
53 zres = (p_pt1 * 1.e-30) / (kbol*t*rho)
70 SUBROUTINE pressure_pt1 ( pges, pgesdz, pgesd2, pgesd3 )
76 REAL,
INTENT(OUT) :: pges
77 REAL,
INTENT(OUT) :: pgesdz
78 REAL,
INTENT(OUT) :: pgesd2
79 REAL,
INTENT(OUT) :: pgesd3
82 REAL :: dzetdv, dicht, dist, fact, z3t
83 REAL :: fres1, fres2, fres3, fres4, fres5, fres
84 REAL :: df_dr, df_drdr
85 REAL :: tfr_1, tfr_2, tfr_3, tfr_4, tfr_5
93 ELSE IF (eta <= 0.1.AND.eta > 0.01)
THEN 98 dist = eta*3.e-3 *fact
101 eta = dicht - 2.0*dist
116 eta = dicht + 2.0*dist
144 df_dr = (-fres4+8.0*fres3-8.0*fres2+fres1)/(12.0*dist)
145 df_drdr = (-fres4+16.0*fres3-3.e1*fres5+16.0*fres2-fres1) &
148 pges = (-fres4+8.0*fres3-8.0*fres2+fres1) &
149 /(12.0*dist) *dzetdv*(kbol*t)/1.e-30
151 pgesdz = (-fres4+16.0*fres3-3.e1*fres5+16.0*fres2-fres1) &
152 /(12.0*(dist**2 ))* dzetdv*(kbol*t)/1.e-30 &
153 + (-fres4+8.0*fres3-8.0*fres2+fres1) /(12.0*dist) * 2.0 *rho &
156 pgesd2 = (fres4-2.0*fres3+2.0*fres2-fres1) /(2.0*dist**3 ) &
157 * dzetdv*(kbol*t)/1.e-30 &
158 + (-fres4+16.0*fres3-3.e1*fres5+16.0*fres2-fres1) /(12.0*(dist**2 )) &
159 * 4.0 *rho *(kbol*t)/1.e-30 + (-fres4+8.0*fres3-8.0*fres2+fres1) &
160 /(12.0*dist) * 2.0 /z3t *(kbol*t)/1.e-30
161 pgesd3 = (fres4-4.0*fres3+6.0*fres5-4.0*fres2+fres1) /(dist**4) &
162 * dzetdv*(kbol*t)/1.e-30 + (fres4-2.0*fres3+2.0*fres2-fres1) &
163 /(2.0*dist**3 ) * 6.0 *rho *(kbol*t)/1.e-30 &
164 + (-fres4+16.0*fres3-3.e1*fres5+16.0*fres2-fres1) &
165 /(12.0*dist**2 )* 6.0 /z3t *(kbol*t)/1.e-30
168 END SUBROUTINE pressure_pt1
WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW This module contains paramete...