CTITLEMODELM -- MATRIX MULTIPLY WITH CONSTRAINTS (RNMO VERSION) 00000010 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR D. THOMPSON 00000020 CA DESIGNER D. THOMPSON 00000030 CA LANGUAGE S/370 FORTRAN H 00000040 CA WRITTEN 2/15/77 00000050 CA REVISED 00000060 CA 00000070 CA CALL MODELM (IND, DSQ, NT, LD, LR, LS, X, Y) 00000080 CA 00000100 CA INPUT IND = SUBSCRIPTS OF UNITY ELEMENTS IN B, ALL I400000110 CA OTHER ELEMENTS ARE ZERO. THIS MATRIX IS 00000120 CA THE GEOMETRY OF THE LINE IN A COMPRESSED FORM. 00000130 CA INPUT DSQ = NORMALIZED OFFSET DISTANCES SQUARED. R400000131 CA INPUT NT = NUMBER OF TRACES I400000140 CA INPUT LD = LAST BIASED DEPTH POINT INDEX I400000150 CA INPUT LR = LAST BIASED REC. INDEX I400000160 CA INPUT LS = LAST BIASED SHOT INDEX I400000170 CA INPUT X = TRIAL SOLUTION VECTOR R800000180 CA OUTPUT Y = TRIAL OBSERVED VALUE VECTOR R800000190 CA 00000270 CA SPECIAL VERSION WITHOUT CDP AND RNMO COMPONENTS. 00000271 CA 00000272 CA THIS ROUTINE IS CALLED BY MCGRDM TO COMPUTE X FOR A GIVEN Y 00000280 CA IT COMPUTES X=BY WHERE B IS THE MODELM MATRIX DEFINED BY IND. 00000290 CAEND 00000295 C 00000300 SUBROUTINE MODELM (IND,DSQ,NT,LD,LR,LS,X,Y,ISAWL,ISAWH, 00000310 1 BIAS,SCALE,KSURF,QRA,QRB) 00000320 C 00000330 IMPLICIT REAL*8 (A-H,O-Z) 00000340 DIMENSION X(1) 00000341 DIMENSION Y(1) 00000342 INTEGER IND(3,1) 00000350 REAL DSQ(1) 00000351 C 00000530 C SCAN ALL TRACES AND SUM CDP, REC AND SHOT COMPONENTS 00000540 C OF X TO FORM MODELM PREDICTED Y(I) VALUE. 00000550 C 00000560 DO 2000000570 * I = 1, NT 00000571 K=IND(1,I) 00000580 20 Y(I)= X(IND(2,I)) + X(IND(3,I)) 00000590 RETURN 00001160 END 00001170