old_contents: import numpy as np
  import matplotlib.pyplot as plt 
  
  from xmitgcm import open_mdsdataset 
  plt.ion()
  
  dir1 = '/homedata/bderembl/runmit/test_southatlgyre3'
  
  ds1 = open_mdsdataset(dir1,iters='all',prefix=['S'])
  
  Height = ds1.S.Z
  print(Height)
  
  nx = len(ds1.S.XC)/2
  print(nx)
  ny = len(ds1.S.YC)/2
  print(ny)
  
  nt = -1
  
  # Vertical Section of Salinity
  plt.figure(1)
  ds1['S'].where(ds1.hFacC>0)[nt,:,100,:].plot()
  plt.title('Case 4 : Salinity (t=-1 ; YC = 30S)')
  plt.savefig('S_Salinity_Vertical_section(x,z)_cas4'+'.png')
  plt.clf()
  
  plt.figure(2)
  ds1['S'].where(ds1.hFacC>0)[nt,:,:,280].plot()
  plt.title('Case 4 : Salinity (t=-1 ; XC = 0E)')
  plt.savefig('S_Salinity_Vertical_section(y,z)_cas4'+'.png')
  plt.clf()
  
new_contents: import numpy as np
  import matplotlib.pyplot as plt 
  
  from xmitgcm import open_mdsdataset 
  plt.ion()
  
  dir1 = '/homedata/bderembl/runmit/test_southatlgyre3'
  
  ds1 = open_mdsdataset(dir1,iters='all',prefix=['S'])
  
  Height = ds1.S.Z
  print(Height)
  
  nx = len(ds1.S.XC)/2
  print(nx)
  ny = len(ds1.S.YC)/2
  print(ny)
  
  nt = -1
  
  # Vertical Section of Salinity
  plt.figure(1)
  ds1['S'].where(ds1.hFacC>0)[nt,:,100,:].plot()
  plt.title('Case 4 : Salinity (t=-1 ; YC = 30S)')
  plt.savefig('S_Salinity_Vertical_section_xz_cas4'+'.png')
  plt.clf()
  
  plt.figure(2)
  ds1['S'].where(ds1.hFacC>0)[nt,:,:,280].plot()
  plt.title('Case 4 : Salinity (t=-1 ; XC = 0E)')
  plt.savefig('S_Salinity_Vertical_section_yz_cas4'+'.png')
  plt.clf()
  
subject: Update C_Salinity_Vertical_sections.py