@@ -0,0 +1,1146 @@
+import numpy as np   %0Aimport matplotlib.pyplot as plt%0Aimport matplotlib.colors as colors%0Aimport cartopy.crs as ccrs%0A%0Afrom cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER%0Afrom xmitgcm import open_mdsdataset %0Aplt.ion()%0A%0Adir0 = '/homedata/bderembl/runmit/test_southatlgyre3'%0A%0Ads0 = open_mdsdataset(dir0,prefix=%5B'S'%5D)%0A%0Ant = 0%0Anz = 0%0A%0Awhile (nt %3C 150) :%0A    nt = nt+1%0A    print(nt)%0A    plt.figure(2)%0A    ax = plt.subplot(projection=ccrs.PlateCarree());%0A    ds0%5B'S'%5D%5Bnt,nz,:,:%5D.plot.pcolormesh('XC', 'YC', ax=ax,vmin=25,vmax=45,cmap='ocean');%0A    plt.title('Case 4 : Salinity ')%0A    plt.text(5,5,nt,ha='center',wrap=True)%0A    ax.coastlines()%0A    gl = ax.gridlines(draw_labels=True, alpha = 0.5, linestyle='--');%0A    gl.xlabels_top = False%0A    gl.ylabels_right = False%0A    gl.xformatter = LONGITUDE_FORMATTER%0A    gl.yformatter = LATITUDE_FORMATTER%0A    if (nt %3C 10):%0A        plt.savefig('Salinity_Cas4-'+'00'+str(nt)+'.png')%0A        plt.clf()%0A    elif (nt %3E 9) and (nt %3C 100):%0A        plt.savefig('Salinity_Cas4-'+'0'+str(nt)+'.png')%0A        plt.clf()%0A    else:%0A        plt.savefig('Salinity_Cas4-'+str(nt)+'.png')%0A        plt.clf()%0A

Create B_Salinity.py