示例#1
0
home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter7/'
out = '/Dropbox/Thesis/IMG/chapter7/'
NCIhome = '/short/fr3/ndh401'
NCIfolder = '/chapter7/'

#==========================================
# NCI general case - initialisation
#==========================================

try:
    arg1 = sys.argv[1]
except IndexError:
    print "Provide arguments <runnum> <numofjobs> <scenario>"

para = Para()
para.central_case(N=100)
para.randomize()
para.set_property_rights(scenario='CS')

run_no = int(arg1)

print '============================================================'
print 'Initialisation for run no: ' + str(run_no)
print '============================================================'

mod = Model(para, ch7=True, turn_off_env=True)
E_lambda = mod.chapter7_initialise()
print '============================================================'
print 'E_lambda: ' + str(E_lambda)
print '============================================================'
示例#2
0
# Decentralised water storage model: master file - Neal Hughes

# Model testing

from __future__ import division
import numpy as np
from para import Para
import model
import time
from econlearn.samplegrid import test
import pylab

if __name__ == '__main__':

    para = Para(rebuild=True)
    para.set_property_rights(scenario='RS')

    para.ch7['inflow_share'] = 0
    para.ch7['capacity_share'] = 0

    mod = model.Model(para, ch7=True, turn_off_env=True)
    #mod.plannerQV_ch7(T=125000, stage2=True, d=0.2, simulate=True, envoff=True)
    #print mod.sim.ITEROLD
    #stats_envoff = mod.sim.stats
    #series_envoff = mod.sim.series

    #del mod
    #mod = model.Model(para, ch7=True, turn_off_env=False)
    #mod.plannerQV_ch7(T=125000, stage2=True, d=0.2, simulate=True, envoff=False)
    #print mod.sim.ITEROLD
    #stats = mod.sim.stats
示例#3
0
# Decentralised water storage model: master file - for control of multiple model runs.   Neal Hughes

from __future__ import division
import numpy as np
import model
from para import Para
from results import chapter8
from results.chartbuilder import *

para = Para(rebuild=True, charts=False)
para.central_case(N=100)
para.solve_para()
para.set_property_rights(scenario='CS')

mod = model.Model(para)

home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter8/'
out = '/Dropbox/Thesis/IMG/chapter8/'

SW = []
SWb = []
S = []
TIME = []

m = 5
sp = [False for i in range(m)]
t1 = [5000, 10000, 20000, 50000, 80000]
t2 = [5000, 10000, 20000, 50000, 80000]
d = [0.2] * m
"""