472 lines
13 KiB
Fortran
472 lines
13 KiB
Fortran
C
|
|
C Complex example FORTRAN program to write a
|
|
C binary data file for Tecplot. This example
|
|
C does the following:
|
|
C
|
|
C 1. Open a data file called "field.plt."
|
|
C 2. Open a data file called "line.plt."
|
|
C 3. Assign values for X, Y and P. These will be used
|
|
C in both the ordered and FE data files.
|
|
C 4. Write out an ordered zone dimensioned 4 x 5 to "field.plt."
|
|
C 5. Assign values for XL and YL arrays.
|
|
C 6. Write out data for line plot to "line.plt." Make the data
|
|
C use double precision.
|
|
C 7. Write out a finite element zone to "field.plt."
|
|
C 8. Write out a text record to "field.plt."
|
|
C 9. Write out a geometry (circle) record to "field.plt."
|
|
C 10. Close file 1.
|
|
C 11. Close file 2.
|
|
C
|
|
Program ComplexTest
|
|
|
|
C This code was written for FORTRAN77 compilers. If you are using a
|
|
C FORTRAN90/95 compiler, we recommend using the .f90 source code files.
|
|
C Otherwise, you will need to alter the include statement below to include
|
|
C tecio.for instead of tecio.inc.
|
|
|
|
Include "tecio.inc"
|
|
|
|
REAL*4 X(4,5), Y(4,5), P(4,5)
|
|
REAL*8 XL(50), YL(50)
|
|
REAL*4 XLDummy(1), YLDummy(1)
|
|
EQUIVALENCE (XLDummy(1), XL(1))
|
|
EQUIVALENCE (YLDummy(1), YL(1))
|
|
REAL*8 SolTime
|
|
INTEGER*4 Debug,I,J,K,L,III,NPts,NElm,DIsDouble,VIsDouble
|
|
INTEGER*4 IMax,JMax,KMax,NM(4,12),FileType
|
|
INTEGER*4 StrandID,ParentZn
|
|
INTEGER*4 SharingZone(3)
|
|
REAL*8 XP, YP, ZP, FH, LineSpacing, PatternLength
|
|
REAL*8 BoxMargin, BoxLineThickness, TextAngle
|
|
INTEGER*4 AttachToZone, Zone, Scope, PositionCoordSys
|
|
INTEGER*4 Clipping
|
|
INTEGER*4 FontType, HeightUnits, Anchor, BoxType
|
|
INTEGER*4 IsFilled, GeomType, LinePattern, NumEllipsePts
|
|
INTEGER*4 BoxColor, BoxFillColor, TextColor, Color, FillColor
|
|
INTEGER*4 ArrowheadStyle, ArrowheadAttachment, NumSegments
|
|
INTEGER*4 NumSegPts(1)
|
|
REAL*8 LineThickness, ArrowheadSize, ArrowheadAngle
|
|
REAL*4 XGeomData(1), YGeomData(1), ZGeomData(1)
|
|
CHARACTER*1 NULCHAR
|
|
INTEGER*4 Zero
|
|
POINTER (NullPtr,Null)
|
|
INTEGER*4 Null(*)
|
|
|
|
Debug = 2
|
|
VIsDouble = 0
|
|
FileType = 0
|
|
DIsDouble = 0
|
|
NULCHAR = CHAR(0)
|
|
Zero = 0
|
|
NullPtr = 0
|
|
C
|
|
C Open field.plt and write the header information.
|
|
C
|
|
I = TECINI112('DATASET WITH 1 ORDERED ZONE, '//
|
|
& '1 QUAD ZONE OVER 2 TIME STEPS'//NULCHAR,
|
|
& 'X Y P'//NULCHAR,
|
|
& 'field.plt'//NULCHAR,
|
|
& '.'//NULCHAR,
|
|
& FileType,
|
|
& Debug,
|
|
& VIsDouble)
|
|
C
|
|
C Open line.plt and write the header information.
|
|
C
|
|
VIsDouble = 1
|
|
I = TECINI112('DATASET WITH ONE I-ORDERED ZONE'//NULCHAR,
|
|
& 'X Y'//NULCHAR,
|
|
& 'line.plt'//NULCHAR,
|
|
& '.'//NULCHAR,
|
|
& FileType,
|
|
& Debug,
|
|
& VIsDouble)
|
|
|
|
C
|
|
C Calculate values for the field variables.
|
|
C
|
|
Do 10 J = 1,5
|
|
Do 10 I = 1,4
|
|
X(I,J) = I
|
|
Y(I,J) = J
|
|
P(I,J) = I*J
|
|
10 Continue
|
|
|
|
C
|
|
C Make sure writing to file #1.
|
|
C
|
|
III = 1
|
|
I = TECFIL112(III)
|
|
|
|
C
|
|
C Write the zone header information for the ordered zone.
|
|
C
|
|
IMax = 4
|
|
JMax = 5
|
|
KMax = 1
|
|
SolTime = 10.0
|
|
StrandID = 1
|
|
ParentZn = 0
|
|
I = TECZNE112('Ordered Zone 1'//NULCHAR,
|
|
& 0, ! ZONETYPE
|
|
& IMax,
|
|
& JMax,
|
|
& KMax,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& SolTime,
|
|
& StrandID,
|
|
& ParentZn,
|
|
& 1, ! ISBLOCK
|
|
& 0, ! NumFaceConnections
|
|
& 0, ! FaceNeighborMode
|
|
& 0, ! TotalNumFaceNodes
|
|
& 0, ! NumConnectedBoundaryFaces
|
|
& 0, ! TotalNumBoundaryConnections
|
|
& Null, ! PassiveVarList
|
|
& Null, ! ValueLocation
|
|
& Null, ! ShareVarFromZone
|
|
& 0) ! ShareConnectivityFromZone)
|
|
|
|
C
|
|
C Write out the field data for the ordered zone.
|
|
C
|
|
III = IMax*JMax
|
|
I = TECDAT112(III,X,DIsDouble)
|
|
I = TECDAT112(III,Y,DIsDouble)
|
|
I = TECDAT112(III,P,DIsDouble)
|
|
|
|
C
|
|
C Calculate values for the I-ordered zone.
|
|
C
|
|
|
|
Do 20 I = 1,50
|
|
XL(I) = I
|
|
YL(I) = sin(I/20.0)
|
|
20 Continue
|
|
C
|
|
C Switch to the 'line.plt' file (file number 2)
|
|
C and write out the line plot data.
|
|
C
|
|
III = 2
|
|
I = TECFIL112(III)
|
|
C
|
|
C Write the zone header information for the XY-data.
|
|
C
|
|
IMax = 50
|
|
JMax = 1
|
|
KMax = 1
|
|
SolTime = 0.0
|
|
StrandID = 0
|
|
I = TECZNE112('XY Line plot'//NULCHAR,
|
|
& 0,
|
|
& IMax,
|
|
& JMax,
|
|
& KMax,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& SolTime,
|
|
& StrandID,
|
|
& ParentZn,
|
|
& 1,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& Null,
|
|
& Null,
|
|
& Null,
|
|
& 0)
|
|
C
|
|
C Write out the line plot.
|
|
C
|
|
DIsDouble = 1
|
|
III = IMax
|
|
I = TECDAT112(III,XLDummy,DIsDouble)
|
|
I = TECDAT112(III,YLDummy,DIsDouble)
|
|
|
|
C
|
|
C Switch back to the field plot file and write out
|
|
C the finite-element zone.
|
|
C
|
|
III = 1
|
|
I = TECFIL112(III)
|
|
C
|
|
C Move the coordinates so this zone's not on top of the other
|
|
C
|
|
Do 30 J = 1,5
|
|
Do 30 I = 1,4
|
|
X(I,J) = I+5
|
|
Y(I,J) = J
|
|
P(I,J) = I*J
|
|
30 Continue
|
|
C
|
|
C Write the zone header information for the finite-element zone.
|
|
C
|
|
NPts = 20
|
|
NElm = 12
|
|
KMax = 1
|
|
SolTime = 10.0
|
|
StrandID = 2
|
|
I = TECZNE112('Finite Zone 1'//NULCHAR,
|
|
& 3, ! FEQUADRILATERAL
|
|
& NPts,
|
|
& NElm,
|
|
& KMax,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& SolTime,
|
|
& StrandID,
|
|
& ParentZn,
|
|
& 1,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& Null,
|
|
& Null,
|
|
& Null,
|
|
& 0)
|
|
C
|
|
C Write out the field data for the finite-element zone.
|
|
C
|
|
IMax = 4
|
|
JMax = 5
|
|
III = IMax*JMax
|
|
DIsDouble = 0
|
|
I = TECDAT112(III,X,DIsDouble)
|
|
I = TECDAT112(III,Y,DIsDouble)
|
|
I = TECDAT112(III,P,DIsDouble)
|
|
|
|
C
|
|
C Calculate and then write out the connectivity list.
|
|
C Note: The NM array references cells starting with
|
|
C offset of 1.
|
|
C
|
|
|
|
Do 40 I = 1,IMax-1
|
|
Do 40 J = 1,JMax-1
|
|
K = I+(J-1)*(IMax-1)
|
|
L = I+(J-1)*IMax
|
|
NM(1,K) = L
|
|
NM(2,K) = L+1
|
|
NM(3,K) = L+IMax+1
|
|
NM(4,K) = L+IMax
|
|
40 Continue
|
|
|
|
I = TECNOD112(NM)
|
|
C
|
|
C Calculate vlues for the new solution variable.
|
|
C
|
|
Do 50 J = 1,5
|
|
Do 50 I = 1,4
|
|
P(I,J) = 2*I*J
|
|
50 Continue
|
|
C
|
|
C Write the zone header information for time step 2
|
|
C
|
|
IMax = 4
|
|
JMax = 5
|
|
KMax = 1
|
|
SolTime = 20.0
|
|
StrandID = 1
|
|
SharingZone(1) = 1
|
|
SharingZone(2) = 1
|
|
SharingZone(3) = 0
|
|
I = TECZNE112('Ordered Zone 2'//NULCHAR,
|
|
& 0, ! ORDERED
|
|
& IMax,
|
|
& JMax,
|
|
& KMax,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& SolTime,
|
|
& StrandID,
|
|
& ParentZn,
|
|
& 1,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& Null,
|
|
& Null,
|
|
& SharingZone,
|
|
& 0)
|
|
C
|
|
C Write out the solution variable the grid variables are shared.
|
|
C
|
|
IMax = 4
|
|
JMax = 5
|
|
III = IMax*JMax
|
|
DIsDouble = 0
|
|
I = TECDAT112(III,P,DIsDouble)
|
|
C
|
|
C Calculate values for the new solution variable.
|
|
C
|
|
Do 60 J = 1,5
|
|
Do 60 I = 1,4
|
|
P(I,J) = 3*I*J
|
|
60 Continue
|
|
C
|
|
C Write another time step for the FEZone and share from the first
|
|
C
|
|
SolTime = 20.0
|
|
StrandID = 2
|
|
KMax = 0
|
|
SharingZone(1) = 2
|
|
SharingZone(2) = 2
|
|
SharingZone(3) = 0
|
|
I = TECZNE112('Finite Zone 2'//NULCHAR,
|
|
& 3, ! FEQUADRILATERAL
|
|
& NPts,
|
|
& NElm,
|
|
& KMax,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& SolTime,
|
|
& StrandID,
|
|
& ParentZn,
|
|
& 1,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& 0,
|
|
& Null,
|
|
& Null,
|
|
& SharingZone,
|
|
& 2)
|
|
C
|
|
C Write out the solution variable the grid variables are shared.
|
|
C
|
|
IMax = 4
|
|
JMax = 5
|
|
III = IMax*JMax
|
|
DIsDouble = 0
|
|
I = TECDAT112(III,P,DIsDouble)
|
|
|
|
C
|
|
C Prepare to write out text record. Text is positioned
|
|
C at 50, 50 in frame units and has a height 5 frame units.
|
|
C
|
|
XP = 50
|
|
YP = 50
|
|
FH = 5
|
|
Scope = 1
|
|
Clipping = 0
|
|
PositionCoordSys = 1
|
|
FontType = 1
|
|
HeightUnits = 1
|
|
AttachToZone = 0
|
|
Zone = 0
|
|
BoxType = 0
|
|
BoxMargin = 5.0
|
|
BoxLineThickness = 0.5
|
|
BoxColor = 3
|
|
BoxFillColor = 7
|
|
TextAngle = 0.0
|
|
Anchor = 0
|
|
LineSpacing = 1.5
|
|
TextColor = 0
|
|
|
|
III = TECTXT112(XP,
|
|
& YP,
|
|
& 0.0d0,
|
|
& PositionCoordSys,
|
|
& AttachToZone,
|
|
& Zone,
|
|
& FontType,
|
|
& HeightUnits,
|
|
& FH,
|
|
& BoxType,
|
|
& BoxMargin,
|
|
& BoxLineThickness,
|
|
& BoxColor,
|
|
& BoxFillColor,
|
|
& TextAngle,
|
|
& Anchor,
|
|
& LineSpacing,
|
|
& TextColor,
|
|
& Scope,
|
|
& Clipping,
|
|
& 'Hi Mom'//NULCHAR,
|
|
& NULCHAR)
|
|
|
|
C
|
|
C Prepare to write out geometry record (circle). Circle is
|
|
C positioned at 25, 25 in frame units and has a radius of 30.
|
|
C Circle is drawn using a dashed line pattern.
|
|
C
|
|
|
|
|
|
XP = 25
|
|
YP = 25
|
|
ZP = 0.0
|
|
IsFilled = 0
|
|
Color = 0
|
|
FillColor = 7
|
|
GeomType = 2
|
|
LinePattern = 1
|
|
LineThickness = 0.3
|
|
PatternLength = 1
|
|
NumEllipsePts = 72
|
|
ArrowheadStyle = 0
|
|
ArrowheadAttachment = 0
|
|
ArrowheadSize = 0.0
|
|
ArrowheadAngle = 15.0
|
|
NumSegments = 1
|
|
NumSegPts(1) = 1
|
|
|
|
XGeomData(1) = 30
|
|
YGeomData(1) = 0.0
|
|
ZGeomData(1) = 0.0
|
|
|
|
|
|
III = TECGEO112(XP,
|
|
& YP,
|
|
& ZP,
|
|
& PositionCoordSys,
|
|
& AttachToZone,
|
|
& Zone,
|
|
& Color,
|
|
& FillColor,
|
|
& IsFilled,
|
|
& GeomType,
|
|
& LinePattern,
|
|
& PatternLength,
|
|
& LineThickness,
|
|
& NumEllipsePts,
|
|
& ArrowheadStyle,
|
|
& ArrowheadAttachment,
|
|
& ArrowheadSize,
|
|
& ArrowheadAngle,
|
|
& Scope,
|
|
& Clipping,
|
|
& NumSegments,
|
|
& NumSegPts,
|
|
& XGeomData,
|
|
& YGeomData,
|
|
& ZGeomData,
|
|
& NULCHAR)
|
|
|
|
C
|
|
C Close out file 1.
|
|
C
|
|
I = TECEND112()
|
|
|
|
C
|
|
C Close out file 2.
|
|
C
|
|
III = 2
|
|
I = TECFIL112(III)
|
|
I = TECEND112()
|
|
STOP
|
|
END
|