diff --git a/templates/event-generator/adapt_member.sh b/templates/event-generator/adapt_member.sh index 5e21c3064f1461e6a730eb9a8eebd12d1a467d66..48bad02db8c9b3c13cfcf345a34de84a2fd42d87 100644 --- a/templates/event-generator/adapt_member.sh +++ b/templates/event-generator/adapt_member.sh @@ -53,8 +53,7 @@ interpolate_SST "${DESTINATION_GRID}" "${SST_FORCING}" "${INTERPOLATED_SST}" # Integrate sst to extpar integrate_sst_to_extpar "${INTERPOLATED_SST}" "${EXTERNAL_PARAMETERS_FILE}" -#cp ${EXTERNAL_PARAMETERS_FILE} extpar.nc # Integrate sst to analysis -# integrate_sst_to_analysis "${INTERPOLATED_SST}" "${ANALYSIS_FILE}" +integrate_sst_to_analysis "${INTERPOLATED_SST}" "${ANALYSIS_FILE}" diff --git a/templates/event-generator/adapt_member_utils.sh b/templates/event-generator/adapt_member_utils.sh index b70f7cee006669a061bf0dc730cda3da689e3914..cbbe6254a9a3859a0a7de49a7327e2b44debaaa8 100644 --- a/templates/event-generator/adapt_member_utils.sh +++ b/templates/event-generator/adapt_member_utils.sh @@ -123,10 +123,10 @@ import numpy as np os.environ['ECCODES_GRIB_NO_INDEX'] = '1' # Open files 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. - ds["T_SEA"]= sst_ds["T_SEA"].isel(time=day_of_year+1) +# 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. +# ds["T_SEA"]= sst_ds["T_SEA"].isel(time=day_of_year+1) ds.to_netcdf("analysis.nc") END diff --git a/templates/event-generator/prepare_namelist.py b/templates/event-generator/prepare_namelist.py index ffeb64297e5230e1ed71e8e4a9ecd09cf3453a07..0223de68d7eddb2f3a21e4c7ed6cddeb26ec2e0a 100644 --- a/templates/event-generator/prepare_namelist.py +++ b/templates/event-generator/prepare_namelist.py @@ -19,10 +19,10 @@ ATMOSPHERE_NAMELIST_PATH = Path("%simulation.namelist_paths.atmosphere%") MASTER_NAMELIST_PATH = Path("%simulation.namelist_paths.master%") # Get first-guess and analysis names: -analysis_filename = (RUNDIR / "igaf*.m*.grb") -analysis_filename = analysis_filename.parent.glob(analysis_filename.name).__next__() -analysis_filename = analysis_filename.name -# analysis_filename = "analysis.nc" +# analysis_filename = (RUNDIR / "igaf*.m*.grb") +# analysis_filename = analysis_filename.parent.glob(analysis_filename.name).__next__() +# analysis_filename = analysis_filename.name +analysis_filename = "analysis.nc" first_guess_filename = (RUNDIR / "igfff00030000.m*.grb").resolve()