Odds and Ends
Creating Analytic Boundaries
Analytic Boundaries can only be created outside Grid3D. The coordinates of
the boundary are calculated in a fortran program which is called by Grid
or FeatFlow to get the coordinates for a given parameter value. The fortran
program has to use exactly the method signatures from the code sample:
C calculates and returns the x coordinate based on the
C given parameters:
C T parameter value, must be DOUBLE PRECISION TO WORK!
C must be in range [0,TMAX(IBCT)]
C IBCT Number of boundary
C I1,I2 utility integers (use for whatever you deem
C neccessary)
C D1,D2 utility doubles (use for whatever you deem
C neccessary)
DOUBLE PRECISION FUNCTION PARX(T,IBCT,I1,I2,D1,D2)
IMPLICIT NONE
DOUBLE PRECISION T
INTEGER IBCT,I1,I2
DOUBLE PRECISION D1,D2
99998 END
C
C
C
C calculates and returns the y coordinate based on the
C given parameters:
C T parameter value, must be DOUBLE PRECISION TO WORK!
C must be in range [0,TMAX(IBCT)]
C IBCT Number of boundary
C I1,I2 utility integers (use for whatever you deem
C neccessary)
C D1,D2 utility doubles (use for whatever you deem
C neccessary)
DOUBLE PRECISION FUNCTION PARY(T,IBCT,I1,I2,D1,D2)
IMPLICIT NONE
DOUBLE PRECISION T
INTEGER IBCT,I1,I2
DOUBLE PRECISION D1,D2
PARY=T;
99999 END
C
C
C
C Returns the maximum parameter value for boundary number
C IBCT. Must be DOUBLE PRECISIION TO WORK
C
DOUBLE PRECISION FUNCTION TMAX(IBCT)
IMPLICIT NONE
INTEGER IBCT
TMAX=1D0;
END
Note that the parameter interval is fixed at [0, TMAX(IBCT)] for each
boundary. You can create several boundaries within the same file. To add this
file to Grid use the browser.