GridTool 3.3, VGRID 3.6, POSTGRID 3.6 and 3.7

****************************************************************************

1)  project.mapbc (created by GridTool): Patch/flow-boundary-condition file

                    ----------------------------
       parameter(mpatch=***)
c
       integer bcpch(mpatch)
c
       open(13,file='project.mapbc',form='formatted')
c
c      4 header lines
       read(13,*)
       read(13,*)
       read(13,*)
       read(13,*)
       do 1 i=1,mpatch
         read(13,*,end=100)ipatch,bcpch(ipatch)
         npatch=ipatch
 1     continue
 100   continue
                   ----------------------------

  where:
     npatch = number of surface patches defining the geometry
      bcpch = flow boundary condition assigned to each surface patch

****************************************************************************

2)  project.bc (created by VGRID/POSTGRID): patch/surface-triangle file

                   ----------------------------
       parameter(mbf=***)
c
       integer patch(mbf),node(3,mbf)
c
       open(unit=12,file='project.bc',form='formatted')
c
       read(12,91)nbf,nb1,npatch,igrid
       read(12,*)
       do 92 if=1,nbf
         read(12,93)if,patch(if),(node(i,if),i=1,3)
 92    continue
 91    format(4i8)
 93    format(5i8)
                   ----------------------------

  where:
         nbf = number of boundary triangular faces
         nb1 = number of surface grid points along patch boundaries
      npatch = number of surface patches
       igrid = 1 for inviscid grids; 2 for viscous grids
   patch(if) = surface patch number
  node(3,if) = nodes of face if

****************************************************************************

3)  project.cogsg (created by VGRID/POSTGRID): grid coordinates and
                                               connectivity file

                   ----------------------------
       parameter(mp=***, me=***)
c
       integer int(me,4)
       real*8  crd(mp,3),t

c
       open(10,file='project.cogsg',form='unformatted')
c
       read(10)inew,ne,np,nb,npv,nev,t,
     &            ((int(ie,in),ie=1,ne),in=1,4)
       read(10)((crd(ip,i),ip=1,np),i=1,3)
                   ----------------------------

  where:
     int = connectivity array
     crd = coordinates array
      ne = total number of cells (tetrahedra)
      np = total number of grid points (including nb)
    inew = a dummy variable (= 0)
      nb = number of grid points on the boundaries
     npv = number of grid points in the viscous layers
           (for Euler grids = 0)
     nev = number of cells in the viscous layers
           (for Euler grids = 0)
       t = a dummy variable

****************************************************************************

4)  project.poin1: field nodes in the B.L. linked to the surface mesh nodes

                     --------------------------
       parameter(mp=***)
c
       integer lpoin(mp)
c
       open(14,file='project.poin1',form='formatted')
c
       read(14,*)nbn,npv
       do ip=1,npv
         read(14,*)jp,lpoin(ip)
       enddo
                        ---------------------

 where

         nbn = number of grid nodes on the boundaries
         npv = number of grid nodes inside the boundary layer (including nbn)
          jp = node index
   lpoin(ip) = surface node linked to the field node "ip"

  Note 1: "viscous" field nodes are generated along rays emanating from
          the surface mesh nodes. lpoin gives the corresponding surface
          node index for each viscous field node along the ray.

  Note 2: lpoin(ip) = ip, for ip=1,nbn (surface mesh nodes).

****************************************************************************