Changeset 352


Ignore:
Timestamp:
Aug 20, 2011 12:11:35 PM (12 years ago)
Author:
vondreele
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/fsource/pypowder.for

    r351 r352  
    2424      RETURN
    2525      END
    26 
    27       SUBROUTINE PYPSVFCJ2(NPTS,DTT,TTHETA,SIG,GAM,SPH,PRFUNC)
    28 C DTT in degrees
    29 C TTHETA in degrees
    30 C SPH is S/L + H/L
    31 Cf2py intent(in) NPTS
    32 Cf2py intent(in) DTT
    33 cf2py depend(NPTS) DTT
    34 Cf2py intent(in) TTHETA
    35 Cf2py intent(in) SIG
    36 Cf2py intent(in) GAM
    37 Cf2py intent(in) SPH
    38 Cf2py intent(out) PRFUNC
    39 Cf2py depend(NPTS) PRFUNC
    40 
    41       REAL*4 DTT(0:NPTS-1),PRFUNC(0:NPTS-1)
    42       SL = SPH/2.0
    43       FW = (2.355*SQRT(SIG)+GAM)/100.0
    44       FMIN = 10.0*(-FW-SPH*COSD(TTHETA))
    45       FMAX = 15.0*FW
    46       DO I=0,NPTS-1
    47         CALL PSVFCJ2(DTT(I)*100.,TTHETA*100.,SL,SL,SIG,GAM,
    48      1    PRFUNC(I))
    49       END DO
    50       RETURN
    51       END
Note: See TracChangeset for help on using the changeset viewer.