1 | SUBROUTINE PYPSVFCJ(NPTS,DTT,TTHETA,SIG,GAM,SPH,PRFUNC) |
---|

2 | C DTT in degrees |
---|

3 | C TTHETA in degrees |
---|

4 | C SPH is S/L + H/L |
---|

5 | C RETURNS FUNCTION ONLY |
---|

6 | Cf2py intent(in) NPTS |
---|

7 | Cf2py intent(in) DTT |
---|

8 | cf2py depend(NPTS) DTT |
---|

9 | Cf2py intent(in) TTHETA |
---|

10 | Cf2py intent(in) SIG |
---|

11 | Cf2py intent(in) GAM |
---|

12 | Cf2py intent(in) SPH |
---|

13 | Cf2py intent(out) PRFUNC |
---|

14 | Cf2py depend(NPTS) PRFUNC |
---|

15 | |
---|

16 | REAL*4 DTT(0:NPTS-1),PRFUNC(0:NPTS-1) |
---|

17 | FW = (2.355*SQRT(SIG)+GAM)/100.0 |
---|

18 | FMIN = 10.0*(-FW-SPH*COSD(TTHETA)) |
---|

19 | FMAX = 15.0*FW |
---|

20 | DO I=0,NPTS-1 |
---|

21 | CALL PSVFCJ(DTT(I)*100.,TTHETA*100.,SIG,GAM,SPH, |
---|

22 | 1 PRFUNC(I),DPRDT,SIGPART,GAMPART,SLPART) |
---|

23 | END DO |
---|

24 | RETURN |
---|

25 | END |
---|

26 | |
---|

27 | SUBROUTINE PYDPSVFCJ(NPTS,DTT,TTHETA,SIG,GAM,SPH,PRFUNC, |
---|

28 | 1 DPRDT,SIGPART,GAMPART,SLPART) |
---|

29 | C DTT in degrees |
---|

30 | C TTHETA in degrees |
---|

31 | C SPH is S/L + H/L |
---|

32 | C RETURNS FUNCTION & DERIVATIVES |
---|

33 | Cf2py intent(in) NPTS |
---|

34 | Cf2py intent(in) DTT |
---|

35 | cf2py depend(NPTS) DTT |
---|

36 | Cf2py intent(in) TTHETA |
---|

37 | Cf2py intent(in) SIG |
---|

38 | Cf2py intent(in) GAM |
---|

39 | Cf2py intent(in) SPH |
---|

40 | Cf2py intent(out) PRFUNC |
---|

41 | Cf2py depend(NPTS) PRFUNC |
---|

42 | Cf2py intent(out) DPRDT |
---|

43 | Cf2py depend(NPTS) DPRDT |
---|

44 | Cf2py intent(out) SIGPART |
---|

45 | Cf2py depend(NPTS) SIGPART |
---|

46 | Cf2py intent(out) GAMPART |
---|

47 | Cf2py depend(NPTS) GAMPART |
---|

48 | Cf2py intent(out) SLPART |
---|

49 | Cf2py depend(NPTS) SLPART |
---|

50 | |
---|

51 | REAL*4 DTT(0:NPTS-1),DPRDT(0:NPTS-1),SIGPART(0:NPTS-1), |
---|

52 | 1 GAMPART(0:NPTS-1),SLPART(0:NPTS-1),PRFUNC(0:NPTS-1) |
---|

53 | FW = (2.355*SQRT(SIG)+GAM)/100.0 |
---|

54 | FMIN = 10.0*(-FW-SPH*COSD(TTHETA)) |
---|

55 | FMAX = 15.0*FW |
---|

56 | DO I=0,NPTS-1 |
---|

57 | CALL PSVFCJ(DTT(I)*100.,TTHETA*100.,SIG,GAM,SPH, |
---|

58 | 1 PRFUNC(I),DPRDT(I),SIGPART(I),GAMPART(I),SLPART(I)) |
---|

59 | DPRDT(I) = DPRDT(I)*100. |
---|

60 | END DO |
---|

61 | RETURN |
---|

62 | END |
---|

63 | |
---|