Encyclopaedia Index

The required GROUND coding is:

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