> expname = "EKF400_v1.3_assim_ghcn_d_better_precip_R30_timeloc"

> fsyr = 1602

> feyr = 2004

> syr_cru = 1901

> eyr_cru = 2004

> syr_recon = 1750

> eyr_recon = 1900

> syr_ncep = 1948

> eyr_ncep = 2009

> sixmonstatevector = T

> if (sixmonstatevector) {
+     s <- 2
+ } else {
+     s <- 12
+ }

> season = c(3, 9)

> nmem = 30

> generate_ECHAM = F

> generate_ECHAM_1901_70 = F

> generate_ECHAM_103 = F

> generate_ECHAM_covar = F

> generate_ECHAM_anom = F

> generate_NCEP = F

> generate_20CR = F

> generate_CRUALLVAR = F

> generate_HadCRU4 = F

> generate_LUTPAULKUT = F

> generate_t_yuri = F

> generate_slp_yuri = F

> generate_GHCN = F

> generate_GHCN_precip = F

> generate_DOCUM = F

> generate_PROXIES = F

> generate_PROXIESnew = F

> if (generate_PROXIESnew == T) {
+     regression_months = c("t.second", "t.third", "t.fourth", 
+         "t.fifth", "p.first", "p.second", "p.third")
+     TRW = T
+     MXD = T
+     SCHWEINGR = T
+     PAGES = F
+     pages_lm_fit = "CRU"
+     type = c("tree")
+     NTREND = F
+ }

> generate_PAGES = F

> generate_NTREND = F

> old_statvec = F

> new_statvec = T

> yuri_temp = F

> yuri_slp = F

> ghcn_temp = F

> isti_instead_ghcn = F

> ghcn_prec = T

> precip_ratio = F

> gauss_ana = F

> check_norm = F

> ghcn_wday = F

> import_luca = F

> trw_only = F

> mxd_only = F

> schweingr_only = F

> if (generate_PROXIESnew) {
+     if ((generate_PAGES & PAGES) | (generate_NTREND & NTREND) | 
+         (generate_PROXIES & generate_PROXIESnew) | (trw_only) | 
+         (mxd_only) | (schweingr_only)) {
+         stop("WARNING! These switches should not be set to TRUE simultaneously: \n\n       generate_PROXIES & generate_PROXIESnew \n       generate_PAGES   & PAGES\n       generate_NTRED   & NTREND\n       trw_only         & generate_PROXIESnew\n       mxd_only         & generate_PROXIESnew\n       schweingr_only   & generate_PROXIESnew\n         ")
+     }
+ }

> no_forc_big_ens = F

> covarclim = 0

> cov_inflate = F

> inflate_fac = 1.02

> state = "changing"

> n_covar = 0

> PHclim_loc = F

> PHclim_lvec_factor = 1

> mixed_loc = F

> update_PHclim = F

> save_ananomallts = F

> calc_decorr_dist = F

> region = "global"

> cor_length_period = "annual"

> loc = T

> if (loc) {
+     l_dist_temp2 = 1000 * 1.5
+     l_dist_slp = 1800 * 1.5
+     l_dist_precip = 300 * 1.5
+     l_dist_gph500 = 1800 * 1.5
+     l_dist_gph100 = 2500 * 1.5
+     l_dist_u850 = 1200 * 1.5
+     l_dist_u200 = 1200 * 1.5
+     l_dist_v850 = 1000 * 1.5
+     l_dist_v200 = 1000 * 1.5
+     l_dist_omega500 = 300 * 1.5
+     l_dist_t500 = 1000 * 1.5
+     l_dist_ind = 999999
+     l_dist_wdays = 300 * 1.5
+     l_dist_blocks = 1800 * 1.5
+     l_dist_cycfreq = 1800 * 1.5
+ } else {
+     l_dist_temp2 = 999999
+     l_dist_slp = 999999
+     l_dist_precip = 999999
+     l_dist_gph500 = 999999
+     l_dist_gph100 = 999999
+     l_dist_u850 = 999999
+     l_dist_u200 = 999999
+     l_dist_v850 = 999999
+     l_dist_v200 = 999999
+     l_dist_omega500 = 999999
+     l_dist_t850 = 999999
+     l_dist_ind = 999999
+ }

> shape_wgt = "circle"

> landcorr = F

> first_prox_per_grid = F

> firstproxres = 10

> avg_prox_per_grid = T

> ins_tim_loc = T

> instmaskprox = F

> reduced_proxies = F

> every2grid = T

> land_only = F

> fasttest = F

> tps_only = T

> tpsw_only = F

> no_stream = F

> loo = F

> if (loo) {
+     tps_only = T
+     no_stream = F
+ }

> check_assimdata = T

> if (no_stream & tps_only) {
+     tps_only = F
+     print("ACHTUNG: tps_only was set to FALSE")
+ }

> scaleprox = T

> anomaly_assim = T

> nseas <- 12

> check_dist = F

> ana.enssize = F

> NCEP_SOCOL = F

> twcr_vali = F

> ncep_vali = F

> monthly_out = F

> calc_prepplot = F

> write_coor = F

> write_netcdf = F

> version = "v1.3"

> if (!monthly_out & write_netcdf) {
+     write_netcdf = F
+     print("ACHTUNG: write_netcdf set to FALSE because monthly_out=F")
+ }

> load_prepplot = F

> statyr = 1955

> load_image = T

> calc_vali_stat = T

> CRPS = TRUE

> vali_plots = F

> ind_ECHAM = F

> ind_recon = F

> validation_set = "cru_vali"

> monthly = F

> pseudoproxy = F

> plot_dweights = F

> write_nc = F

> recalc <- F

> reload <- F

> plstat <- NULL

> countseries <- F