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