Ejemplos de Para en Python

Ejemplo n.º 1

# 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

Ejemplo n.º 2

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 '============================================================'

Ejemplo n.º 3

Archivo: chapter5.py Proyecto: nealbob/regrivermod

# Decentralised water storage model: master file - Neal Hughes

# Chapter 5 model runs

from __future__ import division
import numpy as np
from para import Para
import model
import pickle

para = Para()
para.central_case(N=100, printp=False)
para.set_property_rights(scenario='OA')
para.solve_para()

home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter5/'

scenarios = ['CS', 'SWA', 'OA', 'NS', 'CS-SL', 'SWA-SL', 'CS-SWA']
results = {scen: 0 for scen in scenarios}
policies = {scen: 0 for scen in scenarios}

for i in range(1):
    #try:

    para.central_case(N = 100)
    para.t_cost = 100000000000
    para.aproximate_shares(nonoise=True)
    #if i > 0:
    #    para.randomize(N = 100)
    #    para.aproximate_shares()

Ejemplo n.º 4

from __future__ import division
import numpy as np
from para import Para
from model import Model
from results import chapter6
from results.chartbuilder import *
import multiprocessing
from multiprocessing.queues import Queue
import sys

home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter6/'
NCI = '/short/fr3/ndh401/chapter6/'
out = '/Dropbox/Thesis/IMG/chapter6/'

para = Para()

scenarios = ['CS'] #['RS-HL-O', 'RS-HL', 'RS-O', 'RS', 'CS', 'CS-O', 'CS-HL', 'CS-HL-O', 'CS-U']
results = {scen: 0 for scen in scenarios}
Lambda = {scen: 0 for scen in scenarios}
LambdaK = {scen: 0 for scen in scenarios}

#==========================================
# Central case (with trade)
#==========================================

for scen in scenarios:
    para.set_property_rights(scenario=scen)
    para.aproximate_shares(nonoise=True)

    mod = Model(para)

Ejemplo n.º 5

# 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
"""

Ejemplo n.º 6

Archivo: chapter7_init.py Proyecto: nealbob/regrivermod

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 '============================================================'

Ejemplo n.º 7

Archivo: chapter8.py Proyecto: nealbob/regrivermod

# 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
"""

Ejemplo n.º 8

Archivo: chapter3.py Proyecto: fagan2888/regrivermod

from __future__ import division
import numpy as np
from para import Para
from model import Model
from results import chapter3
import pickle

# Store results here

home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter3/'
NCI = '/short/fr3/ndh401/chapter3/'

# Initialise parameters

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

# Create model instance

mod = Model(para)

# Model runs

result = {'paras' : [], 'stats' : [], 'series' : []}
#temp = [0,0]

for i in range(1000):

Ejemplo n.º 9

Archivo: test.py Proyecto: nealbob/regrivermod

# 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

Ejemplo n.º 10

Archivo: chapter6.py Proyecto: nealbob/regrivermod

from __future__ import division
import numpy as np
from para import Para
from model import Model
from results import chapter6
from results.chartbuilder import *
import multiprocessing
from multiprocessing.queues import Queue
import sys

home = "/home/nealbob"
folder = "/Dropbox/Model/results/chapter6/"
NCI = "/short/fr3/ndh401/chapter6/"
out = "/Dropbox/Thesis/IMG/chapter6/"

para = Para()

scenarios = ["CS"]  # ['RS-HL-O', 'RS-HL', 'RS-O', 'RS', 'CS', 'CS-O', 'CS-HL', 'CS-HL-O', 'CS-U']
results = {scen: 0 for scen in scenarios}
Lambda = {scen: 0 for scen in scenarios}
LambdaK = {scen: 0 for scen in scenarios}

# ==========================================
# Central case (with trade)
# ==========================================

for scen in scenarios:
    para.set_property_rights(scenario=scen)
    para.aproximate_shares(nonoise=True)

    mod = Model(para)