old_contents: import numpy as np
  import cython_hazel
  import matplotlib.pyplot as p
  
  synModeInput = 5
  nSlabsInput = 1
  B1Input = np.asarray([3.0,80.0,41.0])
  B2Input = np.asarray([0.0,0.0,0.0])
  hInput = 3.e0
  tau1Input = 1.e0
  tau2Input = 0.e0
  boundaryInput  = np.asarray([0.0,0.0,0.0,0.0])
  transInput = 1
  atomicPolInput = 1
  anglesInput = np.asarray([90.0,0.0,90.0])
  lambdaAxisInput = np.asarray([-1.5e0,2.5e0])
  nLambdaInput = 150
  dopplerWidthInput = 6.e0
  dopplerWidth2Input = 0.e0
  dampingInput = 0.e0
  dopplerVelocityInput = 0.e0
  dopplerVelocity2Input = 0.e0
  ffInput = 0.e0
  nbarInput = np.asarray([0.0,0.0,0.0,0.0])
  omegaInput = np.asarray([0.0,0.0,0.0,0.0])
  
  cython_hazel.cyHazelInit()
  
  # Compute the Stokes parameters using many default parameters, using Allen's data
  [l, stokes, etaOutput, epsOutput] = cython_hazel.cyHazelSynth(synModeInput, nSlabsInput, B1Input, B2Input, hInput, tau1Input, 
  	tau2Input, boundaryInput, transInput, atomicPolInput, anglesInput, lambdaAxisInput, nLambdaInput, dopplerWidthInput, 
  	dopplerWidth2Input, dampingInput, dopplerVelocityInput, dopplerVelocity2Input, ffInput, nbarInput, omegaInput)
  
  # Now plot the Stokes parameters
  labels = ['I/Imax','Q/Imax','U/Imax','V/Imax']
  
  for i in range(4):
  	p.subplot(2,2,i+1)
  	p.plot(l - 10829.0911, stokes[i,:])
  	p.xlabel('Wavelength [A]')
  	p.ylabel(labels[i])
  	
  p.tight_layout()
  
  p.show()
  
new_contents: import numpy as np
  import pyhazel
  import matplotlib.pyplot as pl
  
  synModeInput = 5
  nSlabsInput = 1
  B1Input = np.asarray([3.0, 80.0, 41.0])
  B2Input = np.asarray([0.0, 0.0, 0.0])
  hInput = 3.e0
  tau1Input = 1.e0
  tau2Input = 0.e0
  boundaryInput = np.asarray([0.0, 0.0, 0.0, 0.0])
  transInput = 1
  atomicPolInput = 1
  anglesInput = np.asarray([90.0, 0.0, 90.0])
  lambdaAxisInput = np.asarray([-1.5e0, 2.5e0])
  nLambdaInput = 150
  dopplerWidthInput = 6.e0
  dopplerWidth2Input = 0.e0
  dampingInput = 0.e0
  dopplerVelocityInput = 0.e0
  dopplerVelocity2Input = 0.e0
  ffInput = 0.e0
  beta = 1
  nbarInput = np.asarray([0.1, 0.1, 0.1, 0.1])
  omegaInput = np.asarray([0.21, 0.1, 0.1, 0.1])
  
  pyhazel.init()
  
  # Compute the Stokes parameters using many default parameters, using
  # Allen's data
  [l, stokes, etaOutput, epsOutput] = pyhazel.synth(synModeInput, nSlabsInput, B1Input, B2Input, hInput,
                                                    tau1Input, tau2Input, boundaryInput, transInput, atomicPolInput, anglesInput,
                                                    lambdaAxisInput, nLambdaInput, dopplerWidthInput, dopplerWidth2Input, dampingInput,
                                                    dopplerVelocityInput, dopplerVelocity2Input, beta, ffInput, nbarInput, omegaInput)
  
  # Now plot the Stokes parameters
  labels = ['I/Imax', 'Q/Imax', 'U/Imax', 'V/Imax']
  
  for i in range(4):
      pl.subplot(2, 2, i + 1)
      pl.plot(l - 10829.0911, stokes[i, :])
      pl.xlabel(r'Wavelength [$\AA$]')
      pl.ylabel(labels[i])
  
  pl.tight_layout()
  
  pl.show()
  
subject: Update test.py