1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162 |
- SUBROUTINE EQUILIBRIUM (TC,USCT,UCET,USET,FSC,FCE,FSE,AA,
- 1 PSC,PCE,BB,TAUSC,TAUCE,TAUSE,ID)
- C ID = 2, Equilibrium without direct train SE
- C ID = 11, Equilibrium with direct train SE
- 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=ID
- Z(4) = TAUSC
- Z(5) = TAUCE
- Z(6) = TAUSE
- Z(7) = PSC
- Z(8) = PCE
- Z(14) = BB
- Z(10) = FSC
- Z(11) = FCE
- Z(12) = FSE
- Z(13) = AA
-
- N=3
- NPT=N+2
- 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
- RHOBEG=5.0D-1
- RHOEND=1.0D-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)
- NPROB = ID - 1
- CALL CALFUN(N,X,F)
- TC = F
- RETURN
- END
-
|