def testEspVICwithPerturb(self):
"""Test ESP forecast VIC simulation, with initialization
from perturbed model simulations."""
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
self.options['vic']['initialize'] = "perturb"
forecast.execute(self.dbname, self.options)
def testNmmeVIC(self):
"""Test forecast VIC simulations using IRI forecast data."""
self.options['forecast']['startdate'] = "2017-7-1"
self.options['forecast']['enddate'] = "2017-7-31"
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "nmme"
forecast.execute(self.dbname, self.options)
def testEspDSSAT(self):
self.options['forecast']['model'] = 'vic, dssat'
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
self.options['forecast']['startdate'] = "2011-2-1"
self.options['forecast']['enddate'] = "2011-4-30"
forecast.execute(self.dbname, self.options)
def testIriVIC(self):
"""Test forecast VIC simulations using IRI forecast data."""
self.options['forecast']['startdate'] = "2012-2-1"
self.options['forecast']['enddate'] = "2012-4-30"
self.options['forecast']['ensemble size'] = 1
self.options['forecast']['method'] = "iri"
forecast.execute(self.dbname, self.options)
def run():
"""Main RHEAS routine."""
config_filename, dbname, db_update, verbose, logfile = parseArgs()
if verbose:
log_level = logging.DEBUG
else:
log_level = logging.INFO
if logfile is None:
logging.basicConfig(level=log_level, format='%(levelname)s: %(message)s')
else:
logging.basicConfig(filename=logfile, level=log_level, format='%(levelname)s: %(message)s')
log = logging.getLogger(__name__)
if dbname is None:
dbname = "rheas"
dbio.connect(dbname)
# check if database update is requested
if db_update:
log.info("Updating database!")
update(dbname, config_filename)
else:
options = config.loadFromFile(config_filename)
# check what simulations have been requested
if "nowcast" in options:
nowcast.execute(dbname, options)
if "forecast" in options:
forecast.execute(dbname, options)
def testEspDSSAT(self):
self.options['forecast']['model'] = 'vic, dssat'
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
self.options['forecast']['startdate'] = "2011-2-1"
self.options['forecast']['enddate'] = "2011-4-30"
forecast.execute(self.dbname, self.options)
def testEspVICwithAssimilation(self):
self.options['forecast']['startdate'] = "2011-1-1"
self.options['forecast']['enddate'] = "2011-1-31"
self.options['vic']['ensemble size'] = 3
self.options['forecast']['ensemble size'] = 3
self.options['vic']['observations'] = "smos"
forecast.execute(self.dbname, self.options)
def testEspVICwithAssimilation(self):
self.options['forecast']['startdate'] = "2010-2-1"
self.options['forecast']['enddate'] = "2010-4-30"
self.options['vic']['ensemble size'] = 3
self.options['forecast']['ensemble size'] = 3
self.options['vic']['observations'] = "smos"
forecast.execute(self.dbname, self.options)
def run():
"""Main RHEAS routine."""
config_filename, dbname, db_update, verbose, logfile = parseArgs()
if verbose:
log_level = logging.DEBUG
else:
log_level = logging.INFO
if logfile is None:
logging.basicConfig(level=log_level,
format='%(levelname)s: %(message)s')
else:
logging.basicConfig(filename=logfile,
level=log_level,
format='%(levelname)s: %(message)s')
log = logging.getLogger(__name__)
if dbname is None:
dbname = "rheas"
dbio.connect(dbname)
# check if database update is requested
if db_update:
log.info("Updating database!")
update(dbname, config_filename)
else:
options = config.loadFromFile(config_filename)
# check what simulations have been requested
if "nowcast" in options:
nowcast.execute(dbname, options)
if "forecast" in options:
forecast.execute(dbname, options)
def testNmmeVIC(self):
"""Test forecast VIC simulations using IRI forecast data."""
self.options['forecast']['startdate'] = "2017-7-1"
self.options['forecast']['enddate'] = "2017-7-31"
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "nmme"
forecast.execute(self.dbname, self.options)
def testEspVICwithAssimilation(self):
"""Test ESP forecast VIC simulation with data assimilation."""
self.options['forecast']['startdate'] = "2011-4-1"
self.options['forecast']['enddate'] = "2011-4-30"
self.options['vic']['ensemble size'] = 3
self.options['forecast']['ensemble size'] = 3
self.options['vic']['observations'] = "smos"
forecast.execute(self.dbname, self.options)
def testEspVICwithAssimilation(self):
"""Test ESP forecast VIC simulation with data assimilation."""
self.options['forecast']['startdate'] = "2011-1-1"
self.options['forecast']['enddate'] = "2011-1-31"
self.options['vic']['ensemble size'] = 3
self.options['forecast']['ensemble size'] = 3
self.options['vic']['observations'] = "smos"
forecast.execute(self.dbname, self.options)
def run():
"""Main RHEAS routine."""
config_filename, dbname, db_update = parseArgs()
if dbname is None:
dbname = "rheas"
dbio.connect(dbname)
# check if database update is requested
if db_update:
print "Updating database!"
update(dbname, config_filename)
else:
options = config.loadFromFile(config_filename)
# check what simulations have been requested
if "nowcast" in options:
nowcast.execute(dbname, options)
if "forecast" in options:
forecast.execute(dbname, options)
def run():
"""Main RHEAS routine."""
config_filename, dbname, db_update = parseArgs()
if dbname is None:
dbname = "rheas"
dbio.connect(dbname)
# check if database update is requested
if db_update:
print "Updating database!"
update(dbname, config_filename)
else:
options = config.loadFromFile(config_filename)
# check what simulations have been requested
if "nowcast" in options:
nowcast.execute(dbname, options)
if "forecast" in options:
forecast.execute(dbname, options)
def run():
"""Main RHEAS routine."""
config_filename, dbname, db_update = parseArgs()
if dbname is None:
dbname = "rheas"
try:
pg.connect("dbname={0}".format(dbname))
except:
print("Cannot connect to database {0}. Please restart it by running \n {1}/pg_ctl -D {2}/postgres restart".format(
dbname, rpath.bins, rpath.data))
sys.exit()
# check if database update is requested
if db_update:
print "Updating database!"
update(dbname, config_filename)
else:
options = config.loadFromFile(config_filename)
# check what simulations have been requested
if "nowcast" in options:
nowcast.execute(dbname, options)
if "forecast" in options:
forecast.execute(dbname, options)
def testEspVIC(self):
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
forecast.execute(self.dbname, self.options)
def testEspDSSAT(self):
self.options['forecast']['model'] = 'vic, dssat'
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
forecast.execute(self.dbname, self.options)
def testIriVIC(self):
self.options['forecast']['startdate'] = "2012-2-1"
self.options['forecast']['enddate'] = "2012-4-30"
self.options['forecast']['ensemble size'] = 1
self.options['forecast']['method'] = "iri"
forecast.execute(self.dbname, self.options)
def testEspVIC(self):
"""Test ESP forecast VIC simulation, with random initialization."""
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
forecast.execute(self.dbname, self.options)
def testEspVICwithPerturb(self):
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
self.options['vic']['initialize'] = "perturb"
forecast.execute(self.dbname, self.options)
def testEspVICwithPerturb(self):
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
self.options['vic']['initialize'] = "perturb"
forecast.execute(self.dbname, self.options)
def testIriVIC(self):
self.options['forecast']['startdate'] = "2001-4-1"
self.options['forecast']['enddate'] = "2001-6-30"
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "iri"
forecast.execute(self.dbname, self.options)
def testEspVIC(self):
self.options['forecast']['ensemble size'] = 2
self.options['forecast']['method'] = "esp"
forecast.execute(self.dbname, self.options)
