1 | include(dom.inc)
   2 |       PROGRAM avbptest
   3 | 
   4 | !       ================================================================!
   5 | !                                                                       !
   6 | !       avbptest.F: This program tests the data structure used for the  !
   7 | !                   geometrical pre-processing in Prissma. An AVBP      !
   8 | !                   complete geometry is tested.                        !
   9 | !                                                                       !
  10 | !       nota      : For any implementation of the DATAStructure (DATAS) !
  11 | !                   follow the steps used in this example.              !
  12 | !                                                                       !
  13 | !                                                                       !
  14 | !       author    : J. AMAYA (avril 2007)                               !
  15 | !                                                                       !
  16 | !       ================================================================!
  17 | 
  18 |         USE datas
  19 | 
  20 |         IMPLICIT NONE
  21 | 
  22 |         DOM_INT          :: ndirs
  23 |         DOM_REAL :: s(4,3)
  24 |         CHARACTER*80     :: file1, file2, path
  25 |         REAL             :: t1, t2
  26 | 
  27 |         WRITE(*,*) "--- Testing a real 3D geometry: ---"
  28 |         CALL CPU_TIME(t1)
  29 |         WRITE(*,*) "    Starting time: ", t1
  30 | 
  31 | !       ------------------------------------!
  32 | !       Creation of vertex, faces and cells !
  33 | !       ------------------------------------!
  34 |         WRITE(*,*) "--- Adding faces and cells: ---"
  35 |         file1 = "MESH/datastest.coor"
  36 |         file2 = "MESH/datastest.conn"
  37 |         CALL avbp2dom(file1, file2)
  38 | 
  39 | !       -------------------------!
  40 | !       Print faces' information !
  41 | !       -------------------------!
  42 | 
  43 |         CALL printfacesinfo
  44 | 
  45 | !       ---------------!
  46 | !       Set directions !
  47 | !       ---------------!
  48 | 
  49 |         ndirs = 4
  50 |         s(1,:) = (/ 1, 0, 0 /)
  51 |         s(2,:) = (/ 0, 1, 0 /)
  52 |         s(3,:) = (/ 0, 0, 1 /)
  53 |         s(4,:) = (/ 1, 1, 0 /)
  54 | 
  55 | !       ------------------------!
  56 | !       Link neighbouring cells !
  57 | !       ------------------------!
  58 |         WRITE(*,*) "--- Creating all links: ---"
  59 | 
  60 |         CALL create_cell_link(ndirs, s)
  61 | 
  62 | !       ----------------------------!
  63 | !       Printing cells' information !
  64 | !       ----------------------------!
  65 | 
  66 |         CALL printcellsinfo(ndirs)
  67 | 
  68 |         CALL CPU_TIME(t2)
  69 |         WRITE(*,*) "Test finished in ", t2 - t1, " seconds."
  70 | 
  71 | !       ---------------------------------!
  72 | !       Write geometry in Prissma format !
  73 | !       ---------------------------------!
  74 |         path="./"
  75 |         CALL writeinfiles(path)
  76 | 
  77 |       END PROGRAM avbptest