Ignore:
Timestamp:
Aug 25, 2011 12:14:25 PM (11 years ago)
Author:
vondreele
Message:

done with peak fitting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/fsource/pypowder.for

    r353 r355  
    6262      END
    6363
     64      SUBROUTINE PYPSVFCJO(NPTS,DTT,TTHETA,SIG,GAM,SPH,PRFUNC)
     65C DTT in degrees
     66C TTHETA in degrees
     67C SPH is S/L + H/L
     68C RETURNS FUNCTION ONLY
     69Cf2py intent(in) NPTS
     70Cf2py intent(in) DTT
     71cf2py depend(NPTS) DTT
     72Cf2py intent(in) TTHETA
     73Cf2py intent(in) SIG
     74Cf2py intent(in) GAM
     75Cf2py intent(in) SPH
     76Cf2py intent(out) PRFUNC
     77Cf2py depend(NPTS) PRFUNC
     78
     79      REAL*4 DTT(0:NPTS-1),PRFUNC(0:NPTS-1)
     80      FW = (2.355*SQRT(SIG)+GAM)/100.0
     81      FMIN = 10.0*(-FW-SPH*COSD(TTHETA))
     82      FMAX = 15.0*FW
     83      DO I=0,NPTS-1
     84        CALL PSVFCJO(DTT(I)*100.,TTHETA*100.,SIG,GAM,SPH/2.0,SPH/2.0,
     85     1    PRFUNC(I),DPRDT,SIGPART,GAMPART,SLPART,HLPART)
     86      END DO
     87      RETURN
     88      END
     89
     90      SUBROUTINE PYDPSVFCJO(NPTS,DTT,TTHETA,SIG,GAM,SHL,PRFUNC,
     91     1  DPRDT,SIGPART,GAMPART,SLPART)
     92C DTT in degrees
     93C TTHETA in degrees
     94C SPH is S/L + H/L
     95C RETURNS FUNCTION & DERIVATIVES
     96Cf2py intent(in) NPTS
     97Cf2py intent(in) DTT
     98cf2py depend(NPTS) DTT
     99Cf2py intent(in) TTHETA
     100Cf2py intent(in) SIG
     101Cf2py intent(in) GAM
     102Cf2py intent(in) SHL
     103Cf2py intent(out) PRFUNC
     104Cf2py depend(NPTS) PRFUNC
     105Cf2py intent(out) DPRDT
     106Cf2py depend(NPTS) DPRDT
     107Cf2py intent(out) SIGPART
     108Cf2py depend(NPTS) SIGPART
     109Cf2py intent(out) GAMPART
     110Cf2py depend(NPTS) GAMPART
     111Cf2py intent(out) SLPART
     112Cf2py depend(NPTS) SLPART
     113
     114      REAL*4 DTT(0:NPTS-1),DPRDT(0:NPTS-1),SIGPART(0:NPTS-1),
     115     1  GAMPART(0:NPTS-1),SLPART(0:NPTS-1),PRFUNC(0:NPTS-1)
     116      FW = (2.355*SQRT(SIG)+GAM)/100.0
     117      FMIN = 10.0*(-FW-SPH*COSD(TTHETA))
     118      FMAX = 15.0*FW
     119      DO I=0,NPTS-1
     120        CALL PSVFCJO(DTT(I)*100.,TTHETA*100.,SIG,GAM,SHL/2.,SHL/2.,
     121     1    PRFUNC(I),DPRDT(I),SIGPART(I),GAMPART(I),SPART,HPART)
     122          SLPART(I) = SPART
     123        DPRDT(I) = DPRDT(I)*100.
     124      END DO
     125      RETURN
     126      END
     127
Note: See TracChangeset for help on using the changeset viewer.