Change dimension ordering in writing out to netCDF

app/regrid_domain_ceto/file/grid_interp.py

@@ -88,8 +88,8 @@ if rank == 0:
 
     all_attributes = {'lon':{'standard_name':'longitude', 'units':'degrees_east', 'long_name':'longitude'},
             'lat':{'standard_name':'latitude', 'units':'degrees_north', 'long_name':'latitude'},
-            'depth':{'standard_name':'depth', 'units':'m', 'long_name':'depth, measured downwards from free surface', 'axis':'Z', 'positive':'down'},
-            'h':{'standard_name':'h', 'units':'m', 'long_name':'model bathymetry depth, measured downwards from geoid', 'axis':'Z', 'positive':'down', 'missing_value':-32768},
+            'depth':{'standard_name':'depth', 'units':'m', 'long_name':'depth, measured downwards from free surface', 'axis':'depth', 'positive':'down'},
+            'h':{'standard_name':'h', 'units':'m', 'long_name':'model bathymetry depth, measured downwards from geoid', 'axis':'depth', 'positive':'down', 'missing_value':-32768},
             'temp':{'standard_name':'sea_water_potential_temperature','units':'C','missing_value':-32768, 'long_name':'Sea Water Potential Temperature'} ,
             'salinity':{'standard_name':'sea_water_salinity','units':'psu', 'missing_value':-32768, 'long_name':'Sea Water Salinity'},
             'u':{'standard_name':'eastward_sea_water_velocity','units':'m s-1', 'missing_value':-32768, 'long_name':'Eastward Current Velocity'},
@@ -126,10 +126,9 @@ if rank == 0:
 
         for this_var, this_mode in varlist.items():
             if this_mode == 'surface':
-                outfile.add_variable(this_var, np.swapaxes(collected_interp_data[this_var],1,2), ['time', 'lon', 'lat'],
-                                         attributes=all_attributes[this_var], ncopts=ncopts)
+                outfile.add_variable(this_var, np.swapaxes(collected_interp_data[this_var],1,2), ['time', 'lat', 'lon'],  attributes=all_attributes[this_var], ncopts=ncopts)
             else:
-                outfile.add_variable(this_var, np.swapaxes(collected_interp_data[this_var],1,3), ['time', 'lon', 'lat', 'depth'],
+                outfile.add_variable(this_var, np.swapaxes(collected_interp_data[this_var],1,3), ['time', 'depth', 'lat','lon'],
                                          attributes=all_attributes[this_var], ncopts=ncopts)
             
 else: