Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/misc/pchfe.F
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/misc/pchfe.F
r5082 r5086 145 145 IF ( N<2 ) GO TO 5001 146 146 IF ( INCFD<1 ) GO TO 5002 147 DO 1I = 2, N147 DO I = 2, N 148 148 IF ( X(I)<=X(I-1) ) GO TO 5003 149 1 CONTINUE149 END DO 150 150 C 151 151 C FUNCTION DEFINITION IS OK, GO ON. … … 168 168 C LOCATE ALL POINTS IN INTERVAL. 169 169 C 170 DO 20J = JFIRST, NE170 DO J = JFIRST, NE 171 171 IF (XE(J) >= X(IR)) GO TO 30 172 20 CONTINUE172 END DO 173 173 J = NE + 1 174 174 GO TO 40 … … 228 228 C 229 229 C FIRST, LOCATE FIRST POINT TO LEFT OF X(IR-1). 230 DO 44I = JFIRST, J-1230 DO I = JFIRST, J-1 231 231 IF (XE(I) < X(IR-1)) GO TO 45 232 44 CONTINUE232 END DO 233 233 C NOTE-- CANNOT DROP THROUGH HERE UNLESS THERE IS AN ERROR 234 234 C IN CHFEV. … … 240 240 C 241 241 C NOW FIND OUT HOW FAR TO BACK UP IN THE X-ARRAY. 242 DO 46I = 1, IR-1242 DO I = 1, IR-1 243 243 IF (XE(J) < X(I)) GO TO 47 244 46 CONTINUE244 END DO 245 245 C NB-- CAN NEVER DROP THROUGH HERE, SINCE XE(J).LT.X(IR-1). 246 246 C
Note: See TracChangeset
for help on using the changeset viewer.