B.6. Source Programs
349
C
C********************************************************************
C
C
Subroutine FOB calculates solutions at observation locations
C
C********************************************************************
C
C
C
C
C
C
SUBROUTINE FOB (X,Y,IJM,HO,CO,HOB,COB,OBWX,OBWY,LEOB,
L,NP,NE,NOBW,NOBT,IHC)
IMPLICIT REAL*8 (A-H, O-Z)
DIMENSION HO(NP) ,CO (NP) ,HOB (NOBW, NOBT) ,COB(NOBW,NOBT),
1
X(NP) ,Y(NP),IJM(3,NE) ,OBWX(NOBW) ,OBWY(NOBW) ,
2
LEOB(NOBW)
DO J=l,NOBW
MI=IJM(l, LEOB(J))
MJ=IJM(2, LEOB(J))
MK=IJM(3, LEOB(J))
AI=X(MJ) *Y(MK)-X(MK) *Y(MJ)
AJ=X(MK)*Y(MI)-X(MI)*Y(MK)
AK=X(MI) *Y(MJ)-X(MJ) *Y(MI)
BI=Y(MJ)-Y(MK)
BJ=Y(MK) -Y(MI)
BK=Y(MI)-Y(MJ)
CI=X (MK) -x (MJ)
CJ=X(MI)-X(MK)
CK=X(MJ) -X(MI)
TRS=ABS(CJ*BI-BJ*CI)
WI=(AI+BI*OBWX(J)+CI*OBWY(J))/TRS
WJ=(AJ+BJ*OBWX(J)+CJ*OBWY(J))/TRS
WK=(AK+BK*OBWX(J)+CK*OBWY(J))/TRS
WI=ABS(WI)
WJ=ABS(WJ)
WK=ABS(WK)
WW=WI+WJ+WK
WI=WI/WW
WJ=WJ/WW
WK=WK/WW
HOB(J,L)=HO(MI)*WI+HO(MJ) *WJ+HO(MK) *WK
IF (IHC .EQ. 2) THEN
COB(J,L)=CO(MI)*WI+CO(MJ) *WJ+CO(MK) *WK
END IF
END DO
RETURN
END