forc = forcing.get(dt, days) # Load forcing
#------ Change the defaults ---------------
par['chi'] = 0.0001 # No exchange of water in tank
par['X'] = 1 # Basin length
par['Y'] = 1 # Basin width
par['Z'] = 1 # Basin depth
par['V'] = par['X'] * par['Y'] * par['Z']
# ----------------------------------------
Mussel_levels = [0, 1, 2, 4, 8, 16, 30,
60] # Units: number of mussels in tank (or bay)
Tank1_multioutput = {}
for level in Mussel_levels:
InitCond['n_muss'] = float(level)
Tank1_multioutput[str(level)] = model.run(days, dt, InitCond, par, forc)
pickle.dump(Tank1_multioutput, open('Tank1_multioutput.p', 'wb'))
print 'Experiment is DONE!'
fig1 = pexp.densityVSproduction(Tank1_multioutput, optimum=5)
fig1.savefig('Exp2_Tank_tinyFlow_densVsprod.png')
fig2 = pexp.densityVSoxygen(Tank1_multioutput, optimum=5)
fig2.savefig('Exp2_Tank_tinyFlow_densVsoxy.png')
fig3 = pexp.densityVSammonia(Tank1_multioutput, optimum=5)
fig3.savefig('Exp2_Tank_tinyFlow_densVsammonia.png')
import model_NPZD2_SHELLE_openBoundary_withForcing as model
import forcing
import plot_experiment as pexp
experiment_name = 'Exp1_Tank_noFlow'
days, dt, par, InitCond = model.load_defaults()
days = 365 * 1.5
forc = forcing.get(dt, days) # Load forcing
#------ Change the defaults ---------------
par['chi'] = 0.2 # No exchange of water in tank
par['X'] = 1. # Basin length
par['Y'] = 1. # Basin width
par['Z'] = 1. # Basin depth
par['V'] = par['X'] * par['Y'] * par['Z'] # <<<<<<<<<<<THIS IS NEW
# ----------------------------------------
Mussel_levels = [0, 5, 10, 20, 30, 50, 100,
150] # Units: number of mussels in tank (or bay)
multioutput = {}
for level in Mussel_levels:
InitCond['n_muss'] = float(level)
multioutput[str(level)] = model.run(days, dt, InitCond, par, forc)
print 'Experiment is DONE!'
pexp.densityVSproduction(multioutput, optimum=0.)
#------ Change the defaults ---------------
par['chi'] = 0.001 # No exchange of water in tank
par['X'] = 2000 # Basin length
par['Y'] = 100 # Basin width
par['Z'] = 10 # Basin depth
par['V'] = par['X'] * par['Y'] * par['Z']
# ----------------------------------------
Mussel_levels = [
10000, 100000, 1000000, 5000000, 10000000, 50000000, 100000000, 100000000
] # Units: number of mussels in tank (or bay)
Tank1_multioutput = {}
for level in Mussel_levels:
InitCond['n_muss'] = float(level)
Tank1_multioutput[str(level)] = model.run(days, dt, InitCond, par, forc)
pickle.dump(Tank1_multioutput, open('Tank1_multioutput.p', 'wb'))
print 'Experiment is DONE!'
pexp.densityVSproduction(Tank1_multioutput)
#fig1.savefig('Exp4_Embayment_HigherFlow_densVsprod.png')
pexp.densityVSoxygen(Tank1_multioutput)
#fig2.savefig('Exp4_Embayment_HigherFlow_densVsoxy.png')
pexp.densityVSammonia(Tank1_multioutput)
#fig3.savefig('Exp4_Embayment_HigherFlow_densVsammonia.png')