Added 'coordinates' attribute to basin/strait vars.

Author: jkrasting

tools/analysis/refineDiag_ocean_month.py

@@ -119,6 +119,7 @@ def main(args):
     hfbasin = np.ma.array(hfbasin,fill_value=nc_misval)
     hfbasin.long_name = 'Northward Ocean Heat Transport'
     hfbasin.units = 'W'
+    hfbasin.coordinates = 'region'
     hfbasin.cell_methods = 'yq:sum time:mean'
     hfbasin.comment = 'Indo-Pacific heat transport begins at 34 S'
     hfbasin.time_avg_info = 'average_T1,average_T2,average_DT'
@@ -199,6 +200,8 @@ def main(args):
     for k in hfbasin.__dict__.keys():
       if k[0] != '_': hfbasin_out.setncattr(k,hfbasin.__dict__[k])
 
+    region_out.setncattr('standard_name','region')
+
     average_T1_out.setncatts(average_T1.__dict__)
     average_T2_out.setncatts(average_T2.__dict__)
     average_DT_out.setncatts(average_DT.__dict__)

tools/analysis/refineDiag_ocean_month_rho2.py

@@ -64,6 +64,7 @@ def main(args):
     msftyrho = np.ma.array(msftyrho,fill_value=1.e20)
     msftyrho.long_name = 'Ocean Y Overturning Mass Streamfunction'
     msftyrho.units = 'kg s-1'
+    msftyrho.coordinates = 'region'
     msftyrho.cell_methods = 'rho2_i:sum yq:sum time:mean'
     msftyrho.time_avg_info = 'average_T1,average_T2,average_DT'
     msftyrho.standard_name = 'ocean_y_overturning_mass_streamfunction'
@@ -78,6 +79,7 @@ def main(args):
     msftyrhompa = np.ma.array(msftyrhompa,fill_value=1.e20)
     msftyrhompa.long_name = 'ocean Y overturning mass streamfunction due to parameterized mesoscale advection'
     msftyrhompa.units = 'kg s-1'
+    msftyrhompa.coordinates = 'region'
     msftyrhompa.cell_methods = 'rho2_i:sum yq:sum time:mean'
     msftyrhompa.time_avg_info = 'average_T1,average_T2,average_DT'
     msftyrhompa.standard_name = 'ocean_y_overturning_mass_streamfunction_due_to_parameterized_'+\
@@ -138,6 +140,8 @@ def main(args):
 
     msftyrhompa_out = f_out.createVariable('msftyrhompa', np.float32, ('time', 'basin', 'rho2_i', 'yq'), fill_value=1.e20)
     msftyrhompa_out.missing_value = 1.e20
+
+    region_out.setncattr('standard_name','region')
 
     average_T1_out = f_out.createVariable('average_T1', np.float64, ('time'))
     average_T2_out = f_out.createVariable('average_T2', np.float64, ('time'))

tools/analysis/refineDiag_ocean_month_z.py

@@ -110,6 +110,7 @@ def main(args):
     msftyyz = np.ma.array(msftyyz,fill_value=nc_misval)
     msftyyz.long_name = 'Ocean Y Overturning Mass Streamfunction'
     msftyyz.units = 'kg s-1'
+    msftyyz.coordinates = 'region'
     msftyyz.cell_methods = 'z_i:sum yq:sum time:mean'
     msftyyz.time_avg_info = 'average_T1,average_T2,average_DT'
     msftyyz.standard_name = 'ocean_y_overturning_mass_streamfunction'
@@ -124,6 +125,7 @@ def main(args):
     msftyzmpa = np.ma.array(msftyzmpa,fill_value=nc_misval)
     msftyzmpa.long_name = 'ocean Y overturning mass streamfunction due to parameterized mesoscale advection'
     msftyzmpa.units = 'kg s-1'
+    msftyzmpa.coordinates = 'region'
     msftyzmpa.cell_methods = 'z_i:sum yq:sum time:mean'
     msftyzmpa.time_avg_info = 'average_T1,average_T2,average_DT'
     msftyzmpa.standard_name = 'ocean_y_overturning_mass_streamfunction_due_to_parameterized_'+\
@@ -139,6 +141,7 @@ def main(args):
     msftyzsmpa = np.ma.array(msftyzsmpa,fill_value=nc_misval)
     msftyzsmpa.long_name = 'ocean Y overturning mass streamfunction due to parameterized submesoscale advection'
     msftyzsmpa.units = 'kg s-1'
+    msftyzsmpa.coordinates = 'region'
     msftyzsmpa.cell_methods = 'z_i:sum yq:sum time:mean'
     msftyzsmpa.time_avg_info = 'average_T1,average_T2,average_DT'
     msftyzsmpa.standard_name = 'ocean_meridional_overturning_mass_streamfunction_due_to_parameterized_'+\
@@ -162,6 +165,7 @@ def main(args):
     strait_names = np.array( [strait.cmor_name for strait in straits] )
     mfo.long_name = 'Sea Water Transport'
     mfo.units = 'kg s-1'
+    mfo.coordinates = 'strait'
     mfo.cell_methods = 'time:mean'
     mfo.time_avg_info = 'average_T1,average_T2,average_DT'
     mfo.standard_name = 'sea_water_transport_across_line'
@@ -264,6 +268,9 @@ def main(args):
     for k in wmo.__dict__.keys():
       if k[0] != '_': wmo_out.setncattr(k,wmo.__dict__[k])
 
+    region_out.setncattr('standard_name','region')
+    strait_out.setncattr('standard_name','region')
+
     average_T1_out.setncatts(average_T1.__dict__)
     average_T2_out.setncatts(average_T2.__dict__)
     average_DT_out.setncatts(average_DT.__dict__)