/[lmdze]/trunk/libf/phylmd/fisrtilp.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/fisrtilp.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC revision 36 by guez, Thu Dec 2 17:11:04 2010 UTC
# Line 15  SUBROUTINE fisrtilp(dtime,paprs,pplay,t, Line 15  SUBROUTINE fisrtilp(dtime,paprs,pplay,t,
15    USE yoethf    USE yoethf
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    
# Line 338  SUBROUTINE fisrtilp(dtime,paprs,pplay,t, Line 338  SUBROUTINE fisrtilp(dtime,paprs,pplay,t,
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)
Legend:
Removed from v.22  
changed lines
  Added in v.36
  ViewVC Help
Powered by ViewVC 1.1.21