1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768 |
- SUBROUTINE PRFROAD (PRFR,TC,USCT,UCET,USET,FSC,FCE,AA,
- 1 PSC,PCE,BB,TAUSC,TAUCE,TAUSE)
- IMPLICIT REAL*8 (A-H,O-Z)
- DIMENSION Z(14), X(3), XL(3), XU(3)
- DIMENSION W(3000), IACT(51)
- INCLUDE 'param.inc'
- COMMON /PB/ NPROB
- COMMON /VARS/ Z
- ZERO = 0.0D0
- ONE = 1.0D0
- TWO = 2.0D0
-
- NPROB=3
- Z(7) = PSC
- Z(8) = PCE
- Z(14) = BB
- Z(10) = FSC
- Z(11) = FCE
- Z(13) = AA
-
- C INITIAL CONDITION
- X(1)=POPSC/TWO
- X(2)=POPCE/TWO
- X(3)=POPSE/TWO
- C SET BOUNDS FOR VARIABLES
- XL(1)=ZERO
- XL(2)=ZERO
- XL(3)=ZERO
- XU(1)=POPSC
- XU(2)=POPCE
- XU(3)=POPSE
- C SET PARAMETERS FOR BOBYQA
- N=3
- NPT=N+2
- RHOBEG=5.D0
- RHOEND=1.D-6
- IPRINT=0
- MAXFUN=2000
- CALL BOBYQA (N,NPT,X,XL,XU,RHOBEG,RHOEND,IPRINT,MAXFUN,W,IACT)
- USCT = X(1)
- UCET = X(2)
- USET = X(3)
- TAUSC = Z(4)
- TAUCE = Z(5)
- TAUSE = Z(6)
- CALL CALFUN(N,X,F)
- PRFR = -F
- NPROB = 1
- CALL CALFUN(N,X,F)
- TC = F
-
- RETURN
- END
-
|