Existing code is bold.
91 CONTINUE
C * ------------------- SECTION 1 ---------------------------
C For RHO1.LE.GRND--- density for phase 1 Index DEN1
IF(RHO1.EQ.GRND) THEN
C... User's own density laws
IF(IG(1).EQ.1) THEN
C... Density law 1
C FN15 -> Y = (A + B * X1) / X2
CALL FN15(DEN1, P1, TEMP1, PRESS0/RG(1), 1.0/RG(1))
ELSE IF(IG(1).EQ.2) THEN
C... Density Law 2
L0D=L0F(DEN1)
L0P=L0F(P1)
L0T=L0F(TEMP1)
DO I=1,NX*NY
F(L0D+I)=(F(L0P+I)+PRESS0)/(RG(1)*F(L0T+I))
ENDDO
ELSE
C... Error message for unknown density law
CALL WRIT40('Unknown density option ')
STOP
ENDIF
ENDIF
RETURN