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
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)