applications/utilities/postProcessing/dataConversion/foamToTecplot360/tecio/examples/comtest/comtest.f

   1 C
   2 C  Complex example FORTRAN program to write a
   3 C  binary data file for Tecplot. This example
   4 C  does the following:
   5 C
   6 C    1.  Open a data file called "field.plt."
   7 C    2.  Open a data file called "line.plt."
   8 C    3.  Assign values for X, Y and P. These will be used
   9 C        in both the ordered and FE data files.
  10 C    4.  Write out an ordered zone dimensioned 4 x 5 to "field.plt."
  11 C    5.  Assign values for XL and YL arrays.
  12 C    6.  Write out data for line plot to "line.plt."  Make the data
  13 C        use double precision.
  14 C    7.  Write out a finite element zone to "field.plt."
  15 C    8.  Write out a text record to "field.plt."
  16 C    9.  Write out a geometry (circle) record to "field.plt."
  17 C   10.  Close file 1.
  18 C   11.  Close file 2.
  19 C
  20       Program ComplexTest
  21
  22       Include "tecio.inc"
  23
  24       REAL*4      X(4,5),  Y(4,5), P(4,5)
  25       REAL*8      XL(50), YL(50)
  26       REAL*4      XLDummy(1), YLDummy(1)
  27       EQUIVALENCE (XLDummy(1), XL(1))
  28       EQUIVALENCE (YLDummy(1), YL(1))
  29       REAL*8      SolTime
  30       INTEGER*4   Debug,I,J,K,L,III,NPts,NElm,DIsDouble,VIsDouble
  31       INTEGER*4   IMax,JMax,KMax,NM(4,12),FileType
  32       INTEGER*4   StrandID,ParentZn
  33       INTEGER*4   SharingZone(3)
  34       REAL*8      XP, YP, ZP, FH, LineSpacing, PatternLength
  35       REAL*8      BoxMargin, BoxLineThickness, TextAngle
  36       INTEGER*4   AttachToZone, Zone, Scope, PositionCoordSys
  37       INTEGER*4   Clipping
  38       INTEGER*4   FontType, HeightUnits, Anchor, BoxType
  39       INTEGER*4   IsFilled, GeomType, LinePattern, NumEllipsePts
  40       INTEGER*4   BoxColor, BoxFillColor, TextColor, Color, FillColor
  41       INTEGER*4   ArrowheadStyle, ArrowheadAttachment, NumSegments
  42       INTEGER*4   NumSegPts(1)
  43       REAL*8      LineThickness, ArrowheadSize, ArrowheadAngle
  44       REAL*4      XGeomData(1), YGeomData(1), ZGeomData(1)
  45       CHARACTER*1 NULCHAR
  46       INTEGER*4   Zero
  47       POINTER     (NullPtr,Null)
  48       INTEGER*4   Null(*)
  49
  50       Debug     = 2
  51       VIsDouble = 0
  52       FileType  = 0
  53       DIsDouble = 0
  54       NULCHAR   = CHAR(0)
  55       Zero      = 0
  56       NullPtr   = 0
  57 C
  58 C Open field.plt and write the header information.
  59 C
  60       I = TECINI112('DATASET WITH 1 ORDERED ZONE, '//
  61      &              '1 QUAD ZONE OVER 2 TIME STEPS'//NULCHAR,
  62      &              'X Y P'//NULCHAR,
  63      &              'field.plt'//NULCHAR,
  64      &              '.'//NULCHAR,
  65      &               FileType,
  66      &               Debug,
  67      &               VIsDouble)
  68 C
  69 C  Open line.plt and write the header information.
  70 C
  71       VIsDouble = 1
  72       I = TECINI112('DATASET WITH ONE I-ORDERED ZONE'//NULCHAR,
  73      &              'X Y'//NULCHAR,
  74      &              'line.plt'//NULCHAR,
  75      &              '.'//NULCHAR,
  76      &               FileType,
  77      &               Debug,
  78      &               VIsDouble)
  79
  80 C
  81 C  Calculate values for the field variables.
  82 C
  83       Do 10 J = 1,5
  84       Do 10 I = 1,4
  85           X(I,J) = I
  86           Y(I,J) = J
  87           P(I,J) = I*J
  88    10 Continue
  89
  90 C
  91 C  Make sure writing to file #1.
  92 C
  93       III = 1
  94       I = TECFIL112(III)
  95
  96 C
  97 C  Write the zone header information for the ordered zone.
  98 C
  99       IMax = 4
 100       JMax = 5
 101       KMax = 1
 102       SolTime = 10.0
 103       StrandID = 1
 104       ParentZn = 0
 105         I = TECZNE112('Ordered Zone 1'//NULCHAR,
 106      &               0, ! ZONETYPE
 107      &               IMax,
 108      &               JMax,
 109      &               KMax,
 110      &               0,
 111      &               0,
 112      &               0,
 113      &               SolTime,
 114      &               StrandID,
 115      &               ParentZn,
 116      &               1,     ! ISBLOCK
 117      &               0,     ! NumFaceConnections
 118      &               0,     ! FaceNeighborMode
 119      &               0,     ! TotalNumFaceNodes
 120      &               0,     ! NumConnectedBoundaryFaces
 121      &               0,     ! TotalNumBoundaryConnections
 122      &               Null,  ! PassiveVarList
 123      &               Null,  ! ValueLocation
 124      &               Null,  ! ShareVarFromZone
 125      &               0)     ! ShareConnectivityFromZone)
 126
 127 C
 128 C  Write out the field data for the ordered zone.
 129 C
 130       III = IMax*JMax
 131       I   = TECDAT112(III,X,DIsDouble)
 132       I   = TECDAT112(III,Y,DIsDouble)
 133       I   = TECDAT112(III,P,DIsDouble)
 134
 135 C
 136 C  Calculate values for the I-ordered zone.
 137 C
 138
 139       Do 20 I = 1,50
 140          XL(I) = I
 141          YL(I) = sin(I/20.0)
 142    20 Continue
 143 C
 144 C  Switch to the 'line.plt' file (file number 2)
 145 C  and write out the line plot data.
 146 C
 147       III = 2
 148       I = TECFIL112(III)
 149 C
 150 C  Write the zone header information for the XY-data.
 151 C
 152       IMax = 50
 153       JMax = 1
 154       KMax = 1
 155       SolTime = 0.0
 156       StrandID = 0
 157       I = TECZNE112('XY Line plot'//NULCHAR,
 158      &              0,
 159      &              IMax,
 160      &              JMax,
 161      &              KMax,
 162      &              0,
 163      &              0,
 164      &              0,
 165      &              SolTime,
 166      &              StrandID,
 167      &              ParentZn,
 168      &              1,
 169      &              0,
 170      &              0,
 171      &              0,
 172      &              0,
 173      &              0,
 174      &              Null,
 175      &              Null,
 176      &              Null,
 177      &              0)
 178 C
 179 C  Write out the line plot.
 180 C
 181       DIsDouble = 1
 182       III = IMax
 183       I   = TECDAT112(III,XLDummy,DIsDouble)
 184       I   = TECDAT112(III,YLDummy,DIsDouble)
 185
 186 C
 187 C  Switch back to the field plot file and write out
 188 C  the finite-element zone.
 189 C
 190       III = 1
 191       I = TECFIL112(III)
 192 C
 193 C  Move the coordinates so this zone's not on top of the other
 194 C
 195       Do 30 J = 1,5
 196       Do 30 I = 1,4
 197           X(I,J) = I+5
 198           Y(I,J) = J
 199           P(I,J) = I*J
 200    30 Continue
 201 C
 202 C  Write the zone header information for the finite-element zone.
 203 C
 204       NPts      = 20
 205       NElm      = 12
 206       KMax      = 1
 207       SolTime   = 10.0
 208       StrandID  = 2
 209       I = TECZNE112('Finite Zone 1'//NULCHAR,
 210      &              3,  ! FEQUADRILATERAL
 211      &              NPts,
 212      &              NElm,
 213      &              KMax,
 214      &              0,
 215      &              0,
 216      &              0,
 217      &              SolTime,
 218      &              StrandID,
 219      &              ParentZn,
 220      &              1,
 221      &              0,
 222      &              0,
 223      &              0,
 224      &              0,
 225      &              0,
 226      &              Null,
 227      &              Null,
 228      &              Null,
 229      &              0)
 230 C
 231 C  Write out the field data for the finite-element zone.
 232 C
 233       IMax      = 4
 234       JMax      = 5
 235       III       = IMax*JMax
 236       DIsDouble = 0
 237       I    = TECDAT112(III,X,DIsDouble)
 238       I    = TECDAT112(III,Y,DIsDouble)
 239       I    = TECDAT112(III,P,DIsDouble)
 240
 241 C
 242 C  Calculate and then write out the connectivity list.
 243 C  Note: The NM array references cells starting with
 244 C        offset of 1.
 245 C
 246
 247       Do 40 I = 1,IMax-1
 248       Do 40 J = 1,JMax-1
 249           K = I+(J-1)*(IMax-1)
 250           L = I+(J-1)*IMax
 251           NM(1,K) = L
 252           NM(2,K) = L+1
 253           NM(3,K) = L+IMax+1
 254           NM(4,K) = L+IMax
 255    40 Continue
 256
 257       I = TECNOD112(NM)
 258 C
 259 C  Calculate vlues for the new solution variable.
 260 C
 261       Do 50 J = 1,5
 262       Do 50 I = 1,4
 263           P(I,J) = 2*I*J
 264    50 Continue
 265 C
 266 C Write the zone header information for time step 2
 267 C
 268       IMax           = 4
 269       JMax           = 5
 270       KMax           = 1
 271       SolTime        = 20.0
 272       StrandID       = 1
 273       SharingZone(1) = 1
 274       SharingZone(2) = 1
 275       SharingZone(3) = 0
 276       I = TECZNE112('Ordered Zone 2'//NULCHAR,
 277      &              0,  ! ORDERED
 278      &              IMax,
 279      &              JMax,
 280      &              KMax,
 281      &              0,
 282      &              0,
 283      &              0,
 284      &              SolTime,
 285      &              StrandID,
 286      &              ParentZn,
 287      &              1,
 288      &              0,
 289      &              0,
 290      &              0,
 291      &              0,
 292      &              0,
 293      &              Null,
 294      &              Null,
 295      &              SharingZone,
 296      &              0)
 297 C
 298 C  Write out the solution variable the grid variables are shared.
 299 C
 300       IMax      = 4
 301       JMax      = 5
 302       III       = IMax*JMax
 303       DIsDouble = 0
 304       I   = TECDAT112(III,P,DIsDouble)
 305 C
 306 C  Calculate values for the new solution variable.
 307 C
 308       Do 60 J = 1,5
 309       Do 60 I = 1,4
 310           P(I,J) = 3*I*J
 311    60 Continue
 312 C
 313 C  Write another time step for the FEZone and share from the first
 314 C
 315       SolTime = 20.0
 316       StrandID = 2
 317       KMax = 0
 318       SharingZone(1) = 2
 319       SharingZone(2) = 2
 320       SharingZone(3) = 0
 321       I = TECZNE112('Finite Zone 2'//NULCHAR,
 322      &              3,  ! FEQUADRILATERAL
 323      &              NPts,
 324      &              NElm,
 325      &              KMax,
 326      &              0,
 327      &              0,
 328      &              0,
 329      &              SolTime,
 330      &              StrandID,
 331      &              ParentZn,
 332      &              1,
 333      &              0,
 334      &              0,
 335      &              0,
 336      &              0,
 337      &              0,
 338      &              Null,
 339      &              Null,
 340      &              SharingZone,
 341      &              2)
 342 C
 343 C  Write out the solution variable the grid variables are shared.
 344 C
 345       IMax      = 4
 346       JMax      = 5
 347       III       = IMax*JMax
 348       DIsDouble = 0
 349       I   = TECDAT112(III,P,DIsDouble)
 350
 351 C
 352 C  Prepare to write out text record. Text is positioned
 353 C  at 50, 50 in frame units and has a height 5 frame units.
 354 C
 355       XP               = 50
 356       YP               = 50
 357       FH               = 5
 358       Scope            = 1
 359       Clipping         = 0
 360       PositionCoordSys = 1
 361       FontType         = 1
 362       HeightUnits      = 1
 363       AttachToZone     = 0
 364       Zone             = 0
 365       BoxType          = 0
 366       BoxMargin        = 5.0
 367       BoxLineThickness = 0.5
 368       BoxColor         = 3
 369       BoxFillColor     = 7
 370       TextAngle        = 0.0
 371       Anchor           = 0
 372       LineSpacing      = 1.5
 373       TextColor        = 0
 374
 375       III =  TECTXT112(XP,
 376      &                 YP,
 377      &                 0.0d0,
 378      &                 PositionCoordSys,
 379      &                 AttachToZone,
 380      &                 Zone,
 381      &                 FontType,
 382      &                 HeightUnits,
 383      &                 FH,
 384      &                 BoxType,
 385      &                 BoxMargin,
 386      &                 BoxLineThickness,
 387      &                 BoxColor,
 388      &                 BoxFillColor,
 389      &                 TextAngle,
 390      &                 Anchor,
 391      &                 LineSpacing,
 392      &                 TextColor,
 393      &                 Scope,
 394      &                 Clipping,
 395      &                'Hi Mom'//NULCHAR,
 396      &                NULCHAR)
 397
 398 C
 399 C  Prepare to write out geometry record (circle). Circle is
 400 C  positioned at 25, 25 in frame units and has a radius of 30.
 401 C  Circle is drawn using a dashed line pattern.
 402 C
 403
 404
 405       XP                  = 25
 406       YP                  = 25
 407       ZP                  = 0.0
 408       IsFilled            = 0
 409       Color               = 0
 410       FillColor           = 7
 411       GeomType            = 2
 412       LinePattern         = 1
 413       LineThickness       = 0.3
 414       PatternLength       = 1
 415       NumEllipsePts       = 72
 416       ArrowheadStyle      = 0
 417       ArrowheadAttachment = 0
 418       ArrowheadSize       = 0.0
 419       ArrowheadAngle      = 15.0
 420       NumSegments         = 1
 421       NumSegPts(1)        = 1
 422
 423       XGeomData(1) = 30
 424       YGeomData(1) = 0.0
 425       ZGeomData(1) = 0.0
 426
 427
 428       III =  TECGEO112(XP,
 429      &                 YP,
 430      &                 ZP,
 431      &                 PositionCoordSys,
 432      &                 AttachToZone,
 433      &                 Zone,
 434      &                 Color,
 435      &                 FillColor,
 436      &                 IsFilled,
 437      &                 GeomType,
 438      &                 LinePattern,
 439      &                 PatternLength,
 440      &                 LineThickness,
 441      &                 NumEllipsePts,
 442      &                 ArrowheadStyle,
 443      &                 ArrowheadAttachment,
 444      &                 ArrowheadSize,
 445      &                 ArrowheadAngle,
 446      &                 Scope,
 447      &                 Clipping,
 448      &                 NumSegments,
 449      &                 NumSegPts,
 450      &                 XGeomData,
 451      &                 YGeomData,
 452      &                 ZGeomData,
 453      &                 NULCHAR)
 454
 455 C
 456 C  Close out file 1.
 457 C
 458       I = TECEND112()
 459
 460 C
 461 C  Close out file 2.
 462 C
 463       III = 2
 464       I = TECFIL112(III)
 465       I = TECEND112()
 466       STOP
 467       END