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=['T'])

Height = ds1.T.Z
print(Height)

nx = len(ds1.T.XC)/2
print(nx)
ny = len(ds1.T.YC)/2
print(ny)

nt = -1

# Vertical Section of Temperature
plt.figure(1)
ds1['T'].where(ds1.hFacC>0)[nt,:,100,:].plot()
plt.title('Case 4 : Temperature (t=-1 ; YC = 30S)')
plt.savefig('T_Temperature_Vertical_section(x,z)_cas4'+'.png')
plt.clf()

plt.figure(2)
ds1['T'].where(ds1.hFacC>0)[nt,:,:,280].plot()
plt.title('Case 4 : Temperature (t=-1 ; XC = 0E)')
plt.savefig('T_Temperature_Vertical_section(y,z)_cas4'+'.png')
plt.clf()

@@ -501,13 +501,11 @@
 tion
-(x,z)
+_xz
 _cas
@@ -685,13 +685,11 @@
 tion
-(y,z)
+_yz
 _cas

Update C_Temperature_Vertical_sections.py