Trying to convert the analysis to netcdf.

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}"

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
 

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()