CTITLEMSCINM -- MATRIX MULTIPLY A'*Y=X WITH CONSTRAINTS (RNMO VER) 00000010 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR GARY CREWS R&D 00000020 CA DESIGNER GARY CREWS R&D 00000030 CA LANGUAGE S/370 FORTRAN H 00000040 CA WRITTEN FEB 76 00000050 CA REVISED D.D. THOMPSON 2/15/77 00000060 CA 00000070 CA CALL MSCINM (IND, DSQ, NT, LD, LR, LS, X, Y) 00000080 CA 00000090 CA INPUT IND = BIASED STACKING TABLES I4 00000100 CA INPUT DSQ = NORMALIZED OFFSET DISTANCES SQUARED R4 00000101 CA INPUT NT = NUMBER OF TRACES I4 00000110 CA INPUT LD = INDEX OF LAST CDP SOLUTION I4 00000120 CA INPUT LR = INDEX OF LAST REC. SOLUTION I4 00000130 CA INPUT LS = INDEX OF LAST SHOT SOLUTION I4 00000140 CA OUTPUT X = OUTPUT ARRAY (SUMMED TO INPUT) R8 00000150 CA OUTPUT Y = INPUT ARRAY TO MULTIPLY R8 00000160 CA BIAS AND SCALE ARE DEFINED AS FOLLOWS: 00000170 CA IF THE SHOT FOR TRACE I IS AT THE SAME GROUND 00000180 CA POSITION AS THE RECEIVER FOR TRACE J THEN 00000190 CA SCALE*IND(3,I)+BIAS=IND(2,J) 00000200 CA 00000201 CA 00000202 CA SPECIAL VERSION WITHOUT CDP AND RNMO COMPONENTS. 00000203 CAEND 00000206 C 00000210 SUBROUTINE MSCINM (IND, DSQ, NT, LD, LR, LS, X, Y) 00000220 C 00000230 IMPLICIT REAL*8 (A-H,O-Z) 00000240 DIMENSION X(1) 00000250 DIMENSION Y(1) 00000260 INTEGER IND(3,1) 00000270 REAL DSQ(1) 00000280 C 00000380 C THIS LOOP SCANS ALL TRACES AND INCREMENTS APPROPRIATE 00000390 C CDP, REC. AND SHOT LOCATIONS IN X BY Y(I). 00000400 C 00000410 DO 20 00000420 * I = 1, NT 00000430 KR=IND(2,I) 00000450 KS=IND(3,I) 00000451 T=Y(I) 00000470 X(KR)=X(KR)+T 00000490 X(KS)=X(KS)+T 00000500 20 CONTINUE 00000510 RETURN 00000520 END 00000530