(t) Creation of new FN-subroutines

It may happen that GROUND users recognise the power and economy of using FN-type subroutines, but need one which does not appear in the library. In this case, they may care to create the subroutine which they need by copying, with appropriate changes, one of the examples provided in PHENC entry: FNxxx

Inspection of these listings will reveal that they differ in the numbers of LB's and real numbers in their arguments. Thus FN1, having one LB, calls L0F1 to get values of I and IADD; and so on.

Suppose that a GROUND user wished to create a function which would take the square root of the product of a second and a third arrayed variable and put it into the storage allocated to the first; the subroutine would probably appear as follows:

SUBROUTINE FNNEW(K1,K2,K3)
INCLUDE '/phoenics/d_includ/farray'
COMMON/IGE/IXF,IXL,IYF,IYL,IGFILL(21)
CALL L0F3(K1,K2,K3,I,I2M1,I3M1,IADD,'FNNEW ')
DO 1 IX=IXF,IXL
I=I+IADD
DO 1 IY=IYF,IYL
I=I+1
1 F(I)=SQRT(F(I2M1+I)*F(I3M1+I))
END

CALL FNNEW(K1,K2,K3)

wbs