File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -53,7 +53,6 @@ def init_data(my_data, rp):
5353 # set the density to be stratified in the y-direction
5454 myg = my_data .grid
5555
56- xcenter = 0.5 * (myg .xmin + myg .xmax )
5756 ycenter = 0.5 * (myg .ymin + myg .ymax )
5857
5958 p = myg .scratch_array ()
@@ -73,10 +72,10 @@ def init_data(my_data, rp):
7372 idx_l = myg .x2d < (myg .xmax - myg .xmin )/ 3.0
7473 idx_r = myg .x2d >= (myg .xmax - myg .xmin )/ 3.0
7574
76- ymom [idx_l ] = amp * np .sin (4.0 * np .pi * f_l * myg .x2d [idx_l ]/ ( myg . xmax - myg . xmin )) * np . exp ( - ( myg . y2d [ idx_l ] - ycenter ) ** 2 / sigma ** 2 )
77- ymom [ idx_r ] = amp * np . sin ( 4.0 * np . pi * f_r * myg .x2d [ idx_r ] / ( myg . xmax - myg .xmin ))* np .exp (- (myg .y2d [idx_r ]- ycenter )** 2 / sigma ** 2 )
78-
79- print ( ymom [:, 100 ] )
75+ ymom [idx_l ] = amp * np .sin (4.0 * np .pi * f_l * myg .x2d [idx_l ] /
76+ ( myg .xmax - myg .xmin ))* np .exp (- (myg .y2d [idx_l ]- ycenter )** 2 / sigma ** 2 )
77+ ymom [ idx_r ] = amp * np . sin ( 4.0 * np . pi * f_r * myg . x2d [ idx_r ] /
78+ ( myg . xmax - myg . xmin )) * np . exp ( - ( myg . y2d [ idx_r ] - ycenter ) ** 2 / sigma ** 2 )
8079
8180 ymom *= dens
8281
You can’t perform that action at this time.
0 commit comments