Skip to content
Snippets Groups Projects
Commit c2c2b6d4 authored by Oriol Tintó's avatar Oriol Tintó
Browse files

Fix bug

parent f98fad56
No related branches found
No related tags found
No related merge requests found
Pipeline #19181 failed
......@@ -169,7 +169,7 @@ os.environ['ECCODES_GRIB_NO_INDEX'] = '1'
with read("${ANALYSIS_FILE}") as ds, read("${INTERPOLATED_SST}") as sst_ds:
date=ds.time.values
day_of_year = int((date - date.astype('datetime64[Y]')) / np.timedelta64(1, 'D'))
Replace analysis T_SEA with info from the sst_clim file.
# Replace analysis T_SEA with info from the sst_clim file.
ds["T_SEA"]= sst_ds["T_SEA"].isel(time=day_of_year+1)
ds.to_netcdf("analysis.nc")
END
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment