[cig-commits] r14115 - seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Sun Feb 22 06:30:43 PST 2009


Author: dkomati1
Date: 2009-02-22 06:30:42 -0800 (Sun, 22 Feb 2009)
New Revision: 14115

Added:
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_elements_OpenDX_to_CUBIT.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_points_OpenDX_to_CUBIT.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_CUBIT.inp
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_OpenDX.dx
Log:
added doubling brick in CUBIT format, and OpenDX to CUBIT converters written in Fortran


Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_elements_OpenDX_to_CUBIT.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_elements_OpenDX_to_CUBIT.f90	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_elements_OpenDX_to_CUBIT.f90	2009-02-22 14:30:42 UTC (rev 14115)
@@ -0,0 +1,19 @@
+
+program convert_elements
+
+implicit none
+
+integer, parameter :: NELEM = 32
+
+integer i,i1,i2,i3,i4,i5,i6,i7,i8
+
+  do i = 1,NELEM
+! point order is inverted in OpenDX and numbers start at 0, not 1
+    read(*,*) i4,i1,i8,i5,i3,i2,i7,i6
+    write(*,230) i,i1+1,i2+1,i3+1,i4+1,i5+1,i6+1,i7+1,i8+1
+  enddo
+
+ 230 format(i4,',',i4,',',i4,',',i4,',',i4,',',i4,',',i4,',',i4,',',i4)
+
+end program convert_elements
+

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_points_OpenDX_to_CUBIT.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_points_OpenDX_to_CUBIT.f90	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/convert_points_OpenDX_to_CUBIT.f90	2009-02-22 14:30:42 UTC (rev 14115)
@@ -0,0 +1,17 @@
+
+program convert_points
+
+implicit none
+
+integer, parameter :: NPOIN = 67
+
+integer i
+real a,b,c
+
+  do i = 1,NPOIN
+    read(*,*) a,b,c
+    write(*,*) i,',',a,',',b,',',c
+  enddo
+
+end program convert_points
+

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_CUBIT.inp
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_CUBIT.inp	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_CUBIT.inp	2009-02-22 14:30:42 UTC (rev 14115)
@@ -0,0 +1,103 @@
+*HEADING
+cubit(/home/user/pbasini/CUBIT/prova.inp): 02/18/2009: 14:36:21
+*NODE
+           1 ,   1.500000     ,   1.000000     ,   2.000000    
+           2 ,   1.500000     ,   1.000000     ,   1.500000    
+           3 ,   1.500000     ,   1.500000     ,   1.500000    
+           4 ,   1.500000     ,   1.500000     ,   2.000000    
+           5 ,   2.000000     ,   1.000000     ,   2.000000    
+           6 ,   2.000000     ,   1.000000     ,   1.000000    
+           7 ,   2.000000     ,   1.500000     ,   1.000000    
+           8 ,   2.000000     ,   1.500000     ,   2.000000    
+           9 ,   1.500000     ,   2.000000     ,   1.000000    
+          10 ,   1.500000     ,   2.000000     ,   2.000000    
+          11 ,   2.000000     ,   2.000000     ,  0.5000000    
+          12 ,   2.000000     ,   2.000000     ,   2.000000    
+          13 ,   1.000000     ,   1.000000     ,   1.000000    
+          14 ,   1.000000     ,   1.000000     ,  0.5000000    
+          15 ,   1.000000     ,   2.000000     ,  0.5000000    
+          16 ,   1.000000     ,   2.000000     ,   1.000000    
+          17 ,   1.500000     ,   1.000000     ,   1.000000    
+          18 ,   2.000000     ,   1.000000     ,  0.5000000    
+          19 ,   1.000000     ,   1.000000     ,   1.500000    
+          20 ,   1.000000     ,   1.000000     ,   2.000000    
+          21 ,   1.000000     ,   1.500000     ,   1.500000    
+          22 ,   1.000000     ,   1.500000     ,   2.000000    
+          23 ,   1.000000     ,   2.000000     ,   2.000000    
+          24 ,   1.000000     ,   1.000000     ,  0.0000000E+00
+          25 ,   2.000000     ,   1.000000     ,  0.0000000E+00
+          26 ,   2.000000     ,   2.000000     ,  0.0000000E+00
+          27 ,   1.000000     ,   2.000000     ,  0.0000000E+00
+          28 ,   1.500000     ,  0.5000000     ,   1.500000    
+          29 ,   1.500000     ,  0.5000000     ,   2.000000    
+          30 ,   2.000000     ,  0.5000000     ,   1.000000    
+          31 ,   2.000000     ,  0.5000000     ,   2.000000    
+          32 ,   1.500000     ,  0.0000000E+00 ,   1.000000    
+          33 ,   1.500000     ,  0.0000000E+00 ,   2.000000    
+          34 ,   2.000000     ,  0.0000000E+00 ,  0.5000000    
+          35 ,   2.000000     ,  0.0000000E+00 ,   2.000000    
+          36 ,   1.000000     ,  0.0000000E+00 ,  0.5000000    
+          37 ,   1.000000     ,  0.0000000E+00 ,   1.000000    
+          38 ,   1.000000     ,  0.5000000     ,   1.500000    
+          39 ,   1.000000     ,  0.5000000     ,   2.000000    
+          40 ,   1.000000     ,  0.0000000E+00 ,   2.000000    
+          41 ,   2.000000     ,  0.0000000E+00 ,  0.0000000E+00
+          42 ,   1.000000     ,  0.0000000E+00 ,  0.0000000E+00
+          43 ,  0.5000000     ,   1.000000     ,   2.000000    
+          44 ,  0.5000000     ,   1.000000     ,   1.500000    
+          45 ,  0.5000000     ,   1.500000     ,   1.500000    
+          46 ,  0.5000000     ,   1.500000     ,   2.000000    
+          47 ,  0.0000000E+00 ,   1.000000     ,   2.000000    
+          48 ,  0.0000000E+00 ,   1.000000     ,   1.000000    
+          49 ,  0.0000000E+00 ,   1.500000     ,   1.000000    
+          50 ,  0.0000000E+00 ,   1.500000     ,   2.000000    
+          51 ,  0.5000000     ,   2.000000     ,   1.000000    
+          52 ,  0.5000000     ,   2.000000     ,   2.000000    
+          53 ,  0.0000000E+00 ,   2.000000     ,  0.5000000    
+          54 ,  0.0000000E+00 ,   2.000000     ,   2.000000    
+          55 ,  0.5000000     ,   1.000000     ,   1.000000    
+          56 ,  0.0000000E+00 ,   1.000000     ,  0.5000000    
+          57 ,  0.0000000E+00 ,   1.000000     ,  0.0000000E+00
+          58 ,  0.0000000E+00 ,   2.000000     ,  0.0000000E+00
+          59 ,  0.5000000     ,  0.5000000     ,   1.500000    
+          60 ,  0.5000000     ,  0.5000000     ,   2.000000    
+          61 ,  0.0000000E+00 ,  0.5000000     ,   1.000000    
+          62 ,  0.0000000E+00 ,  0.5000000     ,   2.000000    
+          63 ,  0.5000000     ,  0.0000000E+00 ,   1.000000    
+          64 ,  0.5000000     ,  0.0000000E+00 ,   2.000000    
+          65 ,  0.0000000E+00 ,  0.0000000E+00 ,  0.5000000    
+          66 ,  0.0000000E+00 ,  0.0000000E+00 ,   2.000000    
+          67 ,  0.0000000E+00 ,  0.0000000E+00 ,  0.0000000E+00
+*ELEMENT, TYPE=C3D8R, ELSET=EB1
+   1,   1,   2,   3,   4,   5,   6,   7,   8
+   2,   4,   3,   9,  10,   8,   7,  11,  12
+   3,  13,  14,  15,  16,  17,  18,  11,   9
+   4,  20,  19,  21,  22,   1,   2,   3,   4
+   5,   2,  17,   9,   3,   6,  18,  11,   7
+   6,  22,  21,  16,  23,   4,   3,   9,  10
+   7,  19,  13,  16,  21,   2,  17,   9,   3
+   8,  24,  25,  26,  27,  14,  18,  11,  15
+   9,   1,   2,  28,  29,   5,   6,  30,  31
+  10,  29,  28,  32,  33,  31,  30,  34,  35
+  11,  13,  14,  36,  37,  17,  18,  34,  32
+  12,  20,  19,  38,  39,   1,   2,  28,  29
+  13,   2,  17,  32,  28,   6,  18,  34,  30
+  14,  39,  38,  37,  40,  29,  28,  32,  33
+  15,  19,  13,  37,  38,   2,  17,  32,  28
+  16,  24,  25,  41,  42,  14,  18,  34,  36
+  17,  43,  44,  45,  46,  47,  48,  49,  50
+  18,  46,  45,  51,  52,  50,  49,  53,  54
+  19,  13,  14,  15,  16,  55,  56,  53,  51
+  20,  20,  19,  21,  22,  43,  44,  45,  46
+  21,  44,  55,  51,  45,  48,  56,  53,  49
+  22,  22,  21,  16,  23,  46,  45,  51,  52
+  23,  19,  13,  16,  21,  44,  55,  51,  45
+  24,  24,  57,  58,  27,  14,  56,  53,  15
+  25,  43,  44,  59,  60,  47,  48,  61,  62
+  26,  60,  59,  63,  64,  62,  61,  65,  66
+  27,  13,  14,  36,  37,  55,  56,  65,  63
+  28,  20,  19,  38,  39,  43,  44,  59,  60
+  29,  44,  55,  63,  59,  48,  56,  65,  61
+  30,  39,  38,  37,  40,  60,  59,  63,  64
+  31,  19,  13,  37,  38,  44,  55,  63,  59
+  32,  24,  57,  67,  42,  14,  56,  65,  36

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_OpenDX.dx
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_OpenDX.dx	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/doubling_brick/dimitri_superbrick_without_duplicates_OpenDX.dx	2009-02-22 14:30:42 UTC (rev 14115)
@@ -0,0 +1,142 @@
+ object 1 class array type float rank 1 shape 3 items           67  data follows
+   1.500000       1.000000       2.000000    
+   1.500000       1.000000       1.500000    
+   1.500000       1.500000       1.500000    
+   1.500000       1.500000       2.000000    
+   2.000000       1.000000       2.000000    
+   2.000000       1.000000       1.000000    
+   2.000000       1.500000       1.000000    
+   2.000000       1.500000       2.000000    
+   1.500000       2.000000       1.000000    
+   1.500000       2.000000       2.000000    
+   2.000000       2.000000      0.5000000    
+   2.000000       2.000000       2.000000    
+   1.000000       1.000000       1.000000    
+   1.000000       1.000000      0.5000000    
+   1.000000       2.000000      0.5000000    
+   1.000000       2.000000       1.000000    
+   1.500000       1.000000       1.000000    
+   2.000000       1.000000      0.5000000    
+   1.000000       1.000000       1.500000    
+   1.000000       1.000000       2.000000    
+   1.000000       1.500000       1.500000    
+   1.000000       1.500000       2.000000    
+   1.000000       2.000000       2.000000    
+   1.000000       1.000000       0.000000    
+   2.000000       1.000000       0.000000    
+   2.000000       2.000000       0.000000    
+   1.000000       2.000000       0.000000    
+   1.500000      0.5000000       1.500000    
+   1.500000      0.5000000       2.000000    
+   2.000000      0.5000000       1.000000    
+   2.000000      0.5000000       2.000000    
+   1.500000       0.000000       1.000000    
+   1.500000       0.000000       2.000000    
+   2.000000       0.000000      0.5000000    
+   2.000000       0.000000       2.000000    
+   1.000000       0.000000      0.5000000    
+   1.000000       0.000000       1.000000    
+   1.000000      0.5000000       1.500000    
+   1.000000      0.5000000       2.000000    
+   1.000000       0.000000       2.000000    
+   2.000000       0.000000       0.000000    
+   1.000000       0.000000       0.000000    
+  0.5000000       1.000000       2.000000    
+  0.5000000       1.000000       1.500000    
+  0.5000000       1.500000       1.500000    
+  0.5000000       1.500000       2.000000    
+   0.000000       1.000000       2.000000    
+   0.000000       1.000000       1.000000    
+   0.000000       1.500000       1.000000    
+   0.000000       1.500000       2.000000    
+  0.5000000       2.000000       1.000000    
+  0.5000000       2.000000       2.000000    
+   0.000000       2.000000      0.5000000    
+   0.000000       2.000000       2.000000    
+  0.5000000       1.000000       1.000000    
+   0.000000       1.000000      0.5000000    
+   0.000000       1.000000       0.000000    
+   0.000000       2.000000       0.000000    
+  0.5000000      0.5000000       1.500000    
+  0.5000000      0.5000000       2.000000    
+   0.000000      0.5000000       1.000000    
+   0.000000      0.5000000       2.000000    
+  0.5000000       0.000000       1.000000    
+  0.5000000       0.000000       2.000000    
+   0.000000       0.000000      0.5000000    
+   0.000000       0.000000       2.000000    
+   0.000000       0.000000       0.000000    
+ object 2 class array type int rank 1 shape 8 items           32  data follows
+     3      0      7      4      2      1      6      5
+     9      3     11      7      8      2     10      6
+    15     12      8     16     14     13     10     17
+    21     19      3      0     20     18      2      1
+     2      1      6      5      8     16     10     17
+    22     21      9      3     15     20      8      2
+    20     18      2      1     15     12      8     16
+    26     23     14     13     25     24     10     17
+    28      0     30      4     27      1     29      5
+    32     28     34     30     31     27     33     29
+    36     12     31     16     35     13     33     17
+    38     19     28      0     37     18     27      1
+    27      1     29      5     31     16     33     17
+    39     38     32     28     36     37     31     27
+    37     18     27      1     36     12     31     16
+    41     23     35     13     40     24     33     17
+    45     42     49     46     44     43     48     47
+    51     45     53     49     50     44     52     48
+    15     12     50     54     14     13     52     55
+    21     19     45     42     20     18     44     43
+    44     43     48     47     50     54     52     55
+    22     21     51     45     15     20     50     44
+    20     18     44     43     15     12     50     54
+    26     23     14     13     57     56     52     55
+    59     42     61     46     58     43     60     47
+    63     59     65     61     62     58     64     60
+    36     12     62     54     35     13     64     55
+    38     19     59     42     37     18     58     43
+    58     43     60     47     62     54     64     55
+    39     38     63     59     36     37     62     58
+    37     18     58     43     36     12     62     54
+    41     23     35     13     66     56     64     55
+ attribute "element type" string "cubes"
+ attribute "ref" string "positions"
+ object 3 class array type float rank 0 items           32  data follows
+           1
+           2
+           3
+           4
+           5
+           6
+           7
+           8
+           9
+          10
+          11
+          12
+          13
+          14
+          15
+          16
+          17
+          18
+          19
+          20
+          21
+          22
+          23
+          24
+          25
+          26
+          27
+          28
+          29
+          30
+          31
+          32
+ attribute "dep" string "connections"
+ object "irregular positions irregular connections" class field
+ component "positions" value 1
+ component "connections" value 2
+ component "data" value 3
+ end



More information about the CIG-COMMITS mailing list