old_contents: import numpy as np
  import matplotlib.pyplot as plt 
  import cartopy.crs as ccrs 
  
  from xmitgcm import open_mdsdataset 
  from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
  plt.ion()
  
  dir1 = '/homedata/bderembl/runmit/test_southatlgyre6'
  
  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 = 0
  while (nt < 1000) :
    nt=nt+1
    print(nt)
    plt.figure(1)
    ax = plt.subplot(projection=ccrs.PlateCarree());
    ds1['S'].where(ds1.hFacC>0)[nt,:,:,280].plot()
    plt.title('Vertical Section (yz) of Salinity (XC = 0E)')
    plt.text(5,5,nt,ha='center',wrap=True)
    ax.coastlines()
    gl = ax.gridlines(draw_labels=True, alpha = 0.5, linestyle='--');
    gl.xlabels_top = False
    gl.ylabels_right = False
    gl.xformatter = LONGITUDE_FORMATTER
    gl.yformatter = LATITUDE_FORMATTER
    if (nt < 10):
        plt.savefig('Salinity_Vertical_section_yz_Cas6-'+'00'+str(nt)+'.png')
        plt.clf()
     elif (nt > 9) and (nt < 100):
        plt.savefig('Salinity_Vertical_section_yz_Cas6'+'0'+str(nt)+'.png')
        plt.clf()
     else:
        plt.savefig('Salinity_Vertical_section_yz_Cas6'+str(nt)+'.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_southatlgyre6'
  
  ds1 = open_mdsdataset(dir1,iters='all',prefix=['S'])
  
  nt = 0
  nx = 280
  
  while (nt < 1000) :
    nt=nt+1
    print(nt)
    plt.figure(1)
    ds1['S'].where(ds1.hFacC>0)[nt,:,:,nx].plot()
    #plt.pcolormesh('XC','YC',,vmin=30,vmax=39)
    plt.text(-60,-4000,nt)
    plt.title('Vertical Section (xz) of Salinity (YC = 30S)')
    if (nt < 10):
        plt.savefig('Salinity_Vertical_section_xz_Cas6-'+'00'+str(nt)+'.png')
        plt.clf()
    elif (nt > 9) and (nt < 100):
        plt.savefig('Salinity_Vertical_section_xz_Cas6-'+'0'+str(nt)+'.png')
        plt.clf()
    else:
        plt.savefig('Salinity_Vertical_section_xz_Cas6-'+str(nt)+'.png')
        plt.clf()
  
subject: Update A_Salinity_section_yz_movie.py