11 |
USE dimens_m |
USE dimens_m |
12 |
USE dimphy |
USE dimphy |
13 |
USE tracstoke |
USE tracstoke |
14 |
USE yomcst |
USE suphec_m |
15 |
USE yoethf |
USE yoethf_m |
16 |
USE fcttre |
USE fcttre |
17 |
USE comfisrtilp |
USE comfisrtilp |
18 |
use numer_rec, only: erf |
use numer_rec, only: nr_erf |
19 |
|
|
20 |
IMPLICIT NONE |
IMPLICIT NONE |
21 |
|
|
338 |
zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.)) |
zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.)) |
339 |
zpdf_e1(i) = zpdf_a(i) - zpdf_b(i) |
zpdf_e1(i) = zpdf_a(i) - zpdf_b(i) |
340 |
zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i)) |
zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i)) |
341 |
zpdf_e1(i) = 1. - erf(zpdf_e1(i)) |
zpdf_e1(i) = 1. - nr_erf(zpdf_e1(i)) |
342 |
zpdf_e2(i) = zpdf_a(i) + zpdf_b(i) |
zpdf_e2(i) = zpdf_a(i) + zpdf_b(i) |
343 |
zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i)) |
zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i)) |
344 |
zpdf_e2(i) = 1. - erf(zpdf_e2(i)) |
zpdf_e2(i) = 1. - nr_erf(zpdf_e2(i)) |
345 |
IF (zpdf_e1(i)<1.E-10) THEN |
IF (zpdf_e1(i)<1.E-10) THEN |
346 |
rneb(i,k) = 0. |
rneb(i,k) = 0. |
347 |
zqn(i) = zqs(i) |
zqn(i) = zqs(i) |