snapshot...

modules/util/acspo_validate.py

@@ -95,15 +95,15 @@ def acspo_validate(oper_file, cspp_file, rel_tol=0.001, outfile_nc=None):
     valid_lat_oper = lat_oper[both_valid]
 
     xarray_data = xr.Dataset({
-        'sst_cspp': xr.DataArray([valid_sst_cspp], coords=None, dims=None, name='sst_cspp'),
-        'sst_oper': xr.DataArray([valid_sst_oper], coords=None, dims=None, name='sst_oper'),
-        'cspp_lat': xr.DataArray([valid_lat_cspp], coords=None, dims=None, name='cspp_lat'),
-        'cspp_lon': xr.DataArray([valid_lon_cspp], coords=None, dims=None, name='cspp_lon'),
-        'oper_lat': xr.DataArray([valid_lat_oper], coords=None, dims=None, name='oper_lat'),
-        'oper_lon': xr.DataArray([valid_lon_oper], coords=None, dims=None, name='oper_lon'),
+        'sst_cspp': xr.DataArray(valid_sst_cspp, coords=None, dims=None, name='sst_cspp'),
+        'sst_oper': xr.DataArray(valid_sst_oper, coords=None, dims=None, name='sst_oper'),
+        'cspp_lat': xr.DataArray(valid_lat_cspp, coords=None, dims=None, name='cspp_lat'),
+        'cspp_lon': xr.DataArray(valid_lon_cspp, coords=None, dims=None, name='cspp_lon'),
+        'oper_lat': xr.DataArray(valid_lat_oper, coords=None, dims=None, name='oper_lat'),
+        'oper_lon': xr.DataArray(valid_lon_oper, coords=None, dims=None, name='oper_lon'),
     })
     if outfile_nc is not None:
-        xarray_data.to_netcdf('/Users/tomrink/outfile_nc')
+        xarray_data.to_netcdf(outfile_nc)
 
     print('fraction approx equal: ',
           np.sum(np.isclose(valid_sst_cspp, valid_sst_oper, rtol=rel_tol))/np.sum(both_valid))