Changes between Version 2 and Version 3 of PolarizationMaps

Timestamp:

05/01/18 16:17:08 (7 years ago)

Author:

Jonathan

Comment:

—

PolarizationMaps

@@ -39 +39 @@
 
 {{{
-
 #run it at command line with: 
 #visit -assume_format Chombo -cli -s /Users/mhuartee/Box\ Sync/py/polarisation_xmdv_3.py
@@ -71 +70 @@
 dy=(ylim[1]-ylim[0])/float(pix[1])
 dx=(xlim[1]-xlim[0])/float(pix[0])
-where='/scratch/bliu17/Laurence/run_dir_3d_mhd/out/'
+dz=0
 ###############################
 
@@ -83 +82 @@
 DefineScalarExpression("fact", "(px*px+py*py+pz*pz)/rho/2.")#Ekinetic
 
-if mode == MODE3D #Create expressions in advance
+if mode == MODE3D: #Create expressions in advance
      DefineScalarExpression("Bproj", "By*cos(angjet)+Bz*sin(angjet)")
      DefineScalarExpression("I", "fact*(Bx^2+Bproj^2)")
      DefineScalarExpression("Q", "fact*(Bx^2-Bproj^2)")
      DefineScalarExpression("U", "2.*fact*Bproj*Bx")
-else if mode == MODE25D
+elif mode == MODE25D:
      DefineScalarExpression("xp", "1.0")
      DefineScalarExpression("Bzz", "Bz*xp/x")
-     DefineScalarExpression("Bxx", Bz*sqrt(x^2-xp^2)/x")
+     DefineScalarExpression("Bxx", "Bz*sqrt(x^2-xp^2)/x")
      DefineScalarExpression("Bproj", "By*cos(angjet)+Bzz*sin(angjet)")
      DefineScalarExpression("dzdx", "x/sqrt(x^2-xp^2)")   
      DefineScalarExpression("I","2*fact*( Bxx^2+Bproj^2) * dzdx")
-     DefineScalarExpression("Q","2*fact*( (Bxx^2-Bproj^2)*dxdz")
+     DefineScalarExpression("Q","2*fact*( Bxx^2-Bproj^2)*dzdx")
      DefineScalarExpression("U","2*2.*fact*Bxx*Bproj")    
 
@@ -130 +129 @@
 #initial IO
 lines=pix[0]*pix[1]
-filename='/scratch/bliu17/Laurence/run_dir_3d_mhd/polarMaps/'+hoy+'pol.okc'
+filename='/scratch/jcarrol5/'+hoy+'pol.okc'
 f = open(filename, 'w+')
 print f.write("%g %g %g\n" % (6,lines,6))
@@ -143 +142 @@
 #compute (integral in z)
 yshift=.5*tt*(zlim[1]-zlim[0])
-
+z=0;
 x=xlim[0]
 while x<=xlim[1]:
@@ -157 +156 @@
             LineoutAtts.point1 = (x,y,0)
             LineoutAtts.point2 = (xlim[1],y+dz,0)
-      elsif mode == MODE3D:
+      elif mode == MODE3D:
            LineoutAtts.point1 = (x,y-yshift,zlim[0])
            LineoutAtts.point2 = (x,y+yshift,zlim[1])
@@ -181 +180 @@
             math.pow(Q0,2.))/I0 
             )
+      else:
+         degree=0
 
       if Q0 != 0.:
          psi=0.5*(math.atan(U0/Q0)+0.*math.pi)
+      else:
+         psi=0
 
       f.write("%g %g %g %g %g %g\n" % (x,y,z,degree*math.sin(psi),