writeMesh.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	writemesh (volid)

Function/Subroutine Documentation

writemesh()

subroutine writemesh

(

integer, intent(in)

volid

)

Definition at line 29 of file writeMesh.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE fvmbag_meshcontrol_mod
      USE inoutfile_mod
      use element_mod , only : nixs
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: VOLID
C-----------------------------------------------
C     Local variable
C-----------------------------------------------
      CHARACTER(LEN=1024) :: FILENAME
      CHARACTER(LEN = 100) :: RAD_INTER, RAD_HEADER, RAD_BEGIN, RAD_END
      INTEGER :: II, N1, N2, N3, N4
 
      INTEGER :: LEN_TMP_NAME
      CHARACTER(len=4096) :: TMP_NAME
 
 
      rad_header = "#RADIOSS STARTER"
      rad_inter = "#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|"
      rad_begin = "/begin"
      RAD_END = "/END"
 
      WRITE(FILENAME, "(A6,I4.4,A4)") "MONVOL",VOLID,".inc"
      
      TMP_NAME=OUTFILE_NAME(1:OUTFILE_NAME_LEN)//FILENAME(1:len_trim(FILENAME))    
      LEN_TMP_NAME = OUTFILE_NAME_LEN+len_trim(FILENAME)
      OPEN(UNIT=10,FILE=TMP_NAME(1:LEN_TMP_NAME),FORM='formatted')
 
      WRITE(10, "(A100)") RAD_INTER
      WRITE(10, "(A5)") "/NODE"
      WRITE(10, "(A100)") RAD_INTER
      DO II = NB_TOTAL_NODE_OLD + 1, NB_TOTAL_NODE
         WRITE(10, "(I10, 1PG20.13, 1PG20.13, 1PG20.13)") ITAB_TMP(II), 
     .        NODE_COORD(1, II), NODE_COORD(2, II), NODE_COORD(3, II)
      ENDDO
C     VOID MATERIAL
      WRITE(10, "(A100)") RAD_INTER
      WRITE(10, "(A11)") "/MAT/VOID/1"
      WRITE(10, "(A5)") "tetra"
      WRITE(10, "(A12)") "/PROP/VOID/1"
      WRITE(10, "(A5)") "tetra"
      WRITE(10, "(A7)") "/PART/1"
      WRITE(10, "(A5)") "tetra"
      WRITE(10, "(I10, I10)") 1, 1
      WRITE(10, "(A14)") "/GRBRIC/PART/1"
      WRITE(10, "(A5)") "tetra"
      WRITE(10, "(I10)") 1
      WRITE(10, "(A100)") RAD_INTER
C     WRITE TETRAS
      WRITE(10, "(A100)") RAD_INTER
      WRITE(10, "(A9)") "/TETRA4/1"
      WRITE(10, "(A100)") RAD_INTER
      DO II = NB_TOTAL_SOLIDS_OLD + 1, NB_TOTAL_SOLIDS
         N1 = IXS_TEMP(NIXS * (II - 1) + 2)
         N2 = IXS_TEMP(NIXS * (II - 1) + 4)
         N3 = IXS_TEMP(NIXS * (II - 1) + 7)
         N4 = IXS_TEMP(NIXS * (II - 1) + 9)
         WRITE(10, "(5I10)") II, ITAB_TMP(N1), ITAB_TMP(N2), ITAB_TMP(N3), ITAB_TMP(N4)
      ENDDO
C     END OF FILE
      WRITE(10, "(A8)") "#enddata"
      CLOSE(10)