import sys
sys.path.insert(1, "../")
import os
import numpy as np
import pickle
import shutil
import matplotlib.pyplot as plt
try:
Lorenz96
Coupled_Lorenz
except (NameError):
import dates_setup
dates_setup.initialize_dates()
from lorenz_models import Lorenz96, Coupled_Lorenz
# Define environment variables and update Python search path;
# this is a necessary call that must be inserted in the beginning of any driver.
#
import dates_utility as utility # import DATeS utility module(s)
def get_model_info(timespan,
ensemble_size=25,
model=None,
load_ensemble=False,
ignore_existing_repo=False,
repo_file=None):
"""
* > python qg1p5_mcclhmc_test_driver.py *
* *
********************************************************************************************
"""
import sys
sys.path.insert(1, "../")
import numpy as np # this is just to create
from matplotlib import pyplot as plt
import matplotlib.animation as animation
# Define environment variables and update Python search path;
# this is a necessary call that must be inserted in the beginning of any driver.
import dates_setup
dates_setup.initialize_dates(random_seed=1)
#
import dates_utility as utility # import DATeS utility module(s)
# Create a model object
# ---------------------
from qg_1p5_model import QG1p5
model = QG1p5(model_configs=dict(
MREFIN=5,
model_name='QG-1.5',
model_grid_type='cartesian',
observation_operator_type='wind-magnitude',
observation_vector_size=100,
observation_error_variances=4.0,
observation_errors_covariance_method='diagonal',
background_error_variances=5.0,
