File tree Expand file tree Collapse file tree 4 files changed +13
-13
lines changed Expand file tree Collapse file tree 4 files changed +13
-13
lines changed Original file line number Diff line number Diff line change 4949 "mhd" : "T" ,
5050 # Patch 1 - Analytical for v and B
5151 # gamma = 5/3
52- # rho = 25/(36π )
53- # p = 5/(12π )
54- # v = (-sin(2π y), sin(2π x), 0)
55- # B = (-sin(2π y)/sqrt(4π ), sin(4π x)/sqrt(4π ), 0)
52+ # rho = 25/(36*pi )
53+ # p = 5/(12*pi )
54+ # v = (-sin(2*pi* y), sin(2*pi* x), 0)
55+ # B = (-sin(2*pi* y)/sqrt(4*pi ), sin(4*pi* x)/sqrt(4*pi ), 0)
5656 "patch_icpp(1)%hcid" : 250 ,
5757 "patch_icpp(1)%geometry" : 7 ,
5858 "patch_icpp(1)%x_centroid" : 0.5 ,
Original file line number Diff line number Diff line change 5151 "fd_order" : 2 ,
5252 # Patch 1 - Analytical for v and B
5353 # gamma = 5/3
54- # rho = 25/(36π )
55- # p = 5/(12π )
56- # v = (-sin(2π y), sin(2π x), 0)
57- # B = (-sin(2π y)/sqrt(4π ), sin(4π x)/sqrt(4π ), 0)
54+ # rho = 25/(36*pi )
55+ # p = 5/(12*pi )
56+ # v = (-sin(2*pi* y), sin(2*pi* x), 0)
57+ # B = (-sin(2*pi* y)/sqrt(4*pi ), sin(4*pi* x)/sqrt(4*pi ), 0)
5858 "patch_icpp(1)%hcid" : 250 ,
5959 "patch_icpp(1)%geometry" : 7 ,
6060 "patch_icpp(1)%x_centroid" : 0.5 ,
Original file line number Diff line number Diff line change 132132
133133 case (250 ) ! MHD Orszag- Tang vortex
134134 ! gamma = 5 / 3
135- ! rho = 25 / (36 π )
136- ! p = 5 / (12 π )
137- ! v = (- sin (2 π y), sin (2 π x), 0 )
138- ! B = (- sin (2 π y)/ sqrt (4 π ), sin (4 π x)/ sqrt (4 π ), 0 )
135+ ! rho = 25 / (36 * pi )
136+ ! p = 5 / (12 * pi )
137+ ! v = (- sin (2 * pi * y), sin (2 * pi * x), 0 )
138+ ! B = (- sin (2 * pi * y)/ sqrt (4 * pi ), sin (4 * pi * x)/ sqrt (4 * pi ), 0 )
139139
140140 q_prim_vf(momxb)%sf(i, j, 0 ) = - sin (2._wp * pi* y_cc(j))
141141 q_prim_vf(momxb + 1 )%sf(i, j, 0 ) = sin (2._wp * pi* x_cc(i))
Original file line number Diff line number Diff line change @@ -112,7 +112,7 @@ contains
112112 ! 2 : vx - > - (divB) * Bx
113113 ! 3 : vy - > - (divB) * By
114114 ! 4 : vz - > - (divB) * Bz
115- ! 5 : E - > - (divB) * (v· B)
115+ ! 5 : E - > - (divB) * (v dot B)
116116 ! 6 : Bx - > - (divB) * vx
117117 ! 7 : By - > - (divB) * vy
118118 ! 8 : Bz - > - (divB) * vz
You can’t perform that action at this time.
0 commit comments