1 subroutine initphcommonvars() bind(C, name='initPhastaCommonVars') 2 use iso_c_binding 3 include "common.h" 4 5 character(8), dimension(13) :: names 6 & = (/ 'Input ', 'PrProces', 'Rezoning', 'Elm_Form', 7 & 'Solver ', 'Bnd_Flux', 'Output ', 'Mapping ', 8 & 'Gather ', 'Scatter ', 'Begin ', 'End ', 9 & 'Back ' /) 10 ccode = reshape(names,shape(names)) 11 12 intmax = 3 13 master = 0 14 icd = 11 15 16 indsym = reshape((/ 1, 2, 4, 7, 11, 17 & 2, 3, 5, 8, 12, 18 & 4, 5, 6, 9, 13, 19 & 7, 8, 9, 10, 14, 20 & 11, 12, 13, 14, 15 /), 21 & shape(indsym)) 22 23 resfrt = 0.00000000000000d+0 24 25 mbeg = 1 26 mend = 100000 27 mprec = 2 28 29 fin = 'input.dat' 30 fgeom = 'geombc.dat' 31 fpar = 'partition.dat' 32 fbndc = 'bc.dat' 33 fmat = 'material.dat' 34 fecho = 'echo.dat' 35 frstin = 'restart' 36 frstou = 'restart' 37 fhist = 'histor.dat' 38 ferror = 'error.dat' 39 ftable = 'table.dat' 40 fforce = 'forces.dat' 41 fgraph = 'graph.dat' 42 ftime = 'time.out' 43 44 iin = 10 45 igeom = 11 46 ipar = 12 47 ibndc = 13 48 imat = 14 49 iecho = 15 50 iout = 16 51 ichmou = 17 52 irstin = 18 53 irstou = 19 54 ihist = 20 55 iflux = 21 56 ierror = 22 57 itable = 23 58 iforce = 24 59 igraph = 25 60 itime = 26 61 62 mcsyst = 4 63 melCat = 8 64 nenCat = reshape((/ 2, 2, 2, 2, 3, 3, 3, 3, ! 1D 65 & 4, 3, 3, 4, 9, 6, 6, 9, ! 2D 66 & 8, 4, 6, 6, 27, 10, 18, 18 /), 67 & shape(nenCat)) ! 3D 68 nfaCat = reshape((/ 2, 2, 2, 2, 2, 2, 2, 2, ! 1D 69 & 4, 3, 3, 4, 4, 3, 3, 4, ! 2D 70 & 6, 4, 5, 5, 6, 4, 5, 5 /), ! 3D 71 & shape(nfaCat)) 72 73 74 mnodeb = reshape((/ 1, 0, 0, 0, 0, 0, 0, 0, 0, 75 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 76 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 77 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 78 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 79 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 80 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 81 & 1, 0, 0, 0, 0, 0, 0, 0, 0, ! 1D 82 & 1, 2, 0, 0, 0, 0, 0, 0, 0, 83 & 1, 2, 0, 0, 0, 0, 0, 0, 0, 84 & 1, 2, 0, 0, 0, 0, 0, 0, 0, 85 & 1, 2, 0, 0, 0, 0, 0, 0, 0, 86 & 1, 2, 5, 0, 0, 0, 0, 0, 0, 87 & 1, 2, 4, 0, 0, 0, 0, 0, 0, 88 & 1, 2, 4, 0, 0, 0, 0, 0, 0, 89 & 1, 2, 4, 0, 0, 0, 0, 0, 0, ! 2D 90 & 1, 2, 3, 4, 0, 0, 0, 0, 0, 91 & 1, 2, 3, 0, 0, 0, 0, 0, 0, 92 & 1, 2, 3, 0, 0, 0, 0, 0, 0, 93 & 1, 2, 5, 4, 0, 0, 0, 0, 0, 94 & 1, 2, 3, 4, 9, 10, 11, 12, 21, 95 & 1, 2, 3, 5, 6, 9, 0, 0, 0, 96 & 1, 2, 3, 7, 9, 8, 0, 0, 0, 97 & 1, 2, 5, 4, 7, 10, 13, 14, 16 /), 98 & shape(mnodeb)) ! 3D 99 100 101 Planck = 6.62617600000000d-34 102 Stefan = 5.66970000000000d-08 103 Nh = 6.02204500000000d+23 104 Rh = 8.31441000000000d+0 105 gamma = 1.40000000000000d+0 106 gamma1 = 0.40000000000000d+0 107 xN2 = 0.79000000000000d+0 108 xO2 = 0.21000000000000d+0 109 Msh = reshape((/ 2.80000000000000d-2, 3.20000000000000d-2, 110 & 3.00000000000000d-2, 1.40000000000000d-2, 111 & 1.60000000000000d-2 /), 112 & shape(Msh)) 113 h0sh = reshape((/ 0.00000000000000d+0, 0.00000000000000d+0, 114 & 8.97750000000000d+4, 4.70820000000000d+5, 115 & 2.46790000000000d+5 /), 116 & shape(h0sh)) 117 Trot = reshape((/ 2.87000000000000d+0, 2.08000000000000d+0, 118 & 2.45000000000000d+0, 0.00000000000000d+0, 119 & 0.00000000000000d+0 /), 120 & shape(Trot)) 121 sigs = reshape((/ 2.00000000000000d+0, 2.00000000000000d+0, 122 & 1.00000000000000d+0, 0.00000000000000d+0, 123 & 0.00000000000000d+0 /), 124 & shape(sigs)) 125 Tvib = reshape((/ 3.39350000000000d+3, 2.27356000000000d+3, 126 & 2.73887000000000d+3, 0.00000000000000d+0, 127 & 0.00000000000000d+0 /), 128 & shape(Tvib)) 129 g0s = reshape((/ 1.00000000000000d+0, 3.00000000000000d+0, 130 & 4.00000000000000d+0, 4.00000000000000d+0, 131 & 9.00000000000000d+0 /), 132 & shape(g0s)) 133 dofs = reshape((/ 5.00000000000000d+0, 5.00000000000000d+0, 134 & 5.00000000000000d+0, 3.00000000000000d+0, 135 & 3.00000000000000d+0 /), 136 & shape(dofs)) 137 end subroutine initphcommonvars 138