@@ -0,0 +1,1007 @@
+import numpy as np%0Aimport matplotlib.pyplot as plt %0Aimport cartopy.crs as ccrs %0A%0Afrom xmitgcm import open_mdsdataset %0Afrom cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER%0Aplt.ion()%0A%0Adir1 = '/homedata/bderembl/runmit/test_southatlgyre3'%0A%0Ads1 = open_mdsdataset(dir1,iters='all',prefix=%5B'T'%5D)%0A%0Ant = 0%0Anz = 0%0A%0A# Cartography T %0A%0Aplt.figure(1)%0Aax = plt.subplot(projection=ccrs.PlateCarree());%0Ads1%5B'T'%5D.where(ds1.hFacC%3E0)%5Bnt,nz,:,:%5D.plot.pcolormesh('XC', 'YC', ax=ax);%0Aplt.title('Cas 4 : Temperature')%0Aax.coastlines()%0Agl = ax.gridlines(draw_labels=True, alpha = 0.5, linestyle='--');%0Agl.xlabels_top = False%0Agl.ylabels_right = False%0Agl.xformatter = LONGITUDE_FORMATTER%0Agl.yformatter = LATITUDE_FORMATTER%0Aplt.savefig('T_Temperature_cas4'+'.png')%0Aplt.clf()%0A%0A# Averages%0A%0AAverage_T = ds1%5B'T'%5D.mean().values%0Aprint('Average of Temperature')%0Aprint(Average_T,'%C2%B0C')%0AAverage_T_mask = ds1%5B'T'%5D.where(ds1.hFacC%3E0).mean().values%0Aprint('Average of Temperature without continents')%0Aprint(Average_T_mask,'%C2%B0C')%0A

Create A_General_Temperature.py