]>
Commit | Line | Data |
---|---|---|
ef94df36 | 1 | *CMZ : 23/08/93 13.30.10 by Jonathan Butterworth |
2 | *-- Author : | |
3 | SUBROUTINE TSTEXT (N,XEXTR,XLOW,XUP,EDGFAC,FRACT,NCUT,NCDIM,ICUT,D | |
4 | 1ELPLS,DELNEG) | |
5 | INTEGER N, NCDIM, NCUT | |
6 | INTEGER ICUT(NCDIM) | |
7 | DOUBLE PRECISION EDGFAC, FRACT | |
8 | DOUBLE PRECISION XEXTR(N), XLOW(N), XUP(N), DELPLS(N), DELNEG(N) | |
9 | INTEGER I | |
10 | DOUBLE PRECISION DIFX | |
11 | SAVE | |
12 | NCUT=0 | |
13 | DO 10 I=1,N | |
14 | DIFX=XUP(I)-XEXTR(I) | |
15 | DELPLS(I)=DIFX | |
16 | * write(*,*) 'delpls(i),i,xup(i),xextr(i)',delpls(i),i,xup(i) | |
17 | * & ,xextr(i) | |
18 | IF(DIFX.LT.EDGFAC*(XUP(I)-XLOW(I))) GOTO 10 | |
19 | DELPLS(I)=FRACT*DIFX | |
20 | * write(*,*) 'delpls(i),i,fract',delpls(i),i,fract | |
21 | NCUT=NCUT+1 | |
22 | ICUT(NCUT)=I | |
23 | 10 CONTINUE | |
24 | DO 20 I=1,N | |
25 | DIFX=XEXTR(I)-XLOW(I) | |
26 | DELNEG(I)=DIFX | |
27 | * write(*,*) 'delneg(i),i,xup(i),xextr(i)',delneg(i),i,xup(i) | |
28 | * & ,xextr(i) | |
29 | IF(DIFX.LT.EDGFAC*(XUP(I)-XLOW(I))) GOTO 20 | |
30 | DELNEG(I)=FRACT*DIFX | |
31 | * write(*,*) 'delpls(i),i,fract',delpls(i),i,fract | |
32 | NCUT=NCUT+1 | |
33 | ICUT(NCUT)=-I | |
34 | 20 CONTINUE | |
35 | RETURN | |
36 | END | |
37 | ||
38 | ||
39 | ||
40 | ||
41 | ||
42 | ||
43 | ||
44 |