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