diff --git a/modules/util/acspo_validate.py b/modules/util/acspo_validate.py index 207754f205bd359d9b9c82940780f07153e5a8e7..f58923f4e091945f2184898721f0d9d5f123849e 100644 --- a/modules/util/acspo_validate.py +++ b/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))