v2.17.0, v2.18.0dev

enkf/CHANGELOG

@@ -1,3 +1,9 @@
+v.2.17.0
+ PS 09022023
+ -- Added UPDATE option "--calculate-spread-only" to calculate forecast and
+ ensemble spread without updating the ensemble. Note that it will still
+ do the full update without saving the analysis (rather than reading
+ the analysis when available).
 v.2.15.0
  PS 05012023
  -- Implemented detection and elimination of bad observation batches

enkf/calc/obsstats.c

@@ -664,7 +664,7 @@ hashtable* das_processbatches(dasystem* das, hashtable* batches)
  }
 
  /*
- * detect bad batches and write thier stats to FNAME_BADBATCHES
+ * detect bad batches and write their stats to FNAME_BADBATCHES
  */
  if (das->nbadbatchspecs > 0) {
  FILE* f = NULL;

enkf/common/version.c

@@ -13,4 +13,4 @@
  *
  *****************************************************************************/
 
-char* ENKF_VERSION = "2.15.0";
+char* ENKF_VERSION = "2.17.0";

enkf/doc/enkf-userguide.tex

@@ -39,7 +39,7 @@
 
 \begin{document}
 
-\title{EnKF-C user guide\\{\normalsize version 2.15.0}}
+\title{EnKF-C user guide\\{\normalsize version 2.17.0}}
 
 \author{Pavel Sakov}
 \date{June 19, 2014 -- \today}
@@ -1699,6 +1699,8 @@ \section{UPDATE}
  confirm that static ensemble is conditioned for using log space
  --calculate-spread
  calculate forecast and analysis ensemble spread and write to spread.nc
+ --calculate-spread-only
+ calculate forecast and analysis ensemble spread only and write to spread.nc
  --calculate-forecast-spread
  calculate forecast ensemble spread and write to spread.nc
  --calculate-forecast-spread-only

enkf/update/diags.c

@@ -43,9 +43,6 @@ void das_allocatespread(dasystem* das, char fname[])
  if (rank != 0)
  return;
 
- if (file_exists(fname))
- return;
-
  ncw_create(fname, NC_CLOBBER | das->ncformat, &ncid);
  for (vid = 0; vid < nvar; ++vid) {
  char* varname_src = model_getvarname(m, vid);

enkf/update/enkf_update.c

@@ -33,6 +33,8 @@ static void usage()
  enkf_printf(" confirm that static ensemble is conditioned for using log space\n");
  enkf_printf(" --calculate-spread\n");
  enkf_printf(" calculate forecast and analysis ensemble spread and write to %s\n", FNAME_SPREAD);
+ enkf_printf(" --calculate-spread-only\n");
+ enkf_printf(" calculate forecast and analysis ensemble spread only and write to %s\n", FNAME_SPREAD);
  enkf_printf(" --calculate-forecast-spread\n");
  enkf_printf(" calculate forecast ensemble spread and write to %s\n", FNAME_SPREAD);
  enkf_printf(" --calculate-forecast-spread-only\n");
@@ -69,6 +71,7 @@ static void parse_commandline(int argc, char* argv[], char** fname, int* updates
  int no_update = 0;
  int vcorrs_only = 0;
  int fspread_only = 0;
+ int spread_only = 0;
  int i;
 
  if (argc < 2)
@@ -91,12 +94,15 @@ static void parse_commandline(int argc, char* argv[], char** fname, int* updates
  *updatespec |= (UPDATE_DOFORECASTSPREAD | UPDATE_DOANALYSISSPREAD);
  i++;
  continue;
+ } else if (strcmp(argv[i], "--calculate-spread-only") == 0) {
+ spread_only = 1;
+ i++;
+ continue;
  } else if (strcmp(argv[i], "--calculate-forecast-spread") == 0) {
  *updatespec |= UPDATE_DOFORECASTSPREAD;
  i++;
  continue;
  } else if (strcmp(argv[i], "--calculate-forecast-spread-only") == 0) {
- *updatespec |= UPDATE_DOFORECASTSPREAD;
  fspread_only = 1;
  i++;
  continue;
@@ -160,12 +166,17 @@ static void parse_commandline(int argc, char* argv[], char** fname, int* updates
  if (*fname == NULL)
  enkf_quit("command line: parameter file not specified");
 
+ if (vcorrs_only + fspread_only + spread_only > 1)
+ enkf_quit("more than one argument with \"only\"");
+
  if (no_update)
  *updatespec &= ~UPDATE_NEEDAN;
  if (vcorrs_only)
  *updatespec = UPDATE_DOVERTCORRS;
  if (fspread_only)
  *updatespec = UPDATE_DOFORECASTSPREAD;
+ if (spread_only)
+ *updatespec = (UPDATE_DOFORECASTSPREAD | UPDATE_DOANALYSISSPREAD);
 }
 
 /**