def set_model_constants(xx=50000.0, nx=100, va=10.0, tmax=(((60 * 360) * 24) * 3600.0), avep=(24 * 3600.0), dt=3600.0, period=(((3600 * 24) * 360) * 1), B=2.0, T0=(273.15 + 6), dT=2.0, Cs=0.001, Cp=1030.0, ra=1.5, ro=1030.0, ri=900.0, Cpo=4000.0, Cpi=2900.0, H=200.0, vo=0.2, Hb=1000.0, Li=3300000.0, Tf=(273.15 - 1.8), SW0=50.0, SW_anom=100.0, emissivity=0.99, Da=1000000.0, Do=500.0, tau_entrainment=((30 * 24) * 3600.0), **args):
    'Setup model constants. All of the constants have fixed values, but one can pass in own values or even some arbitrary values via **args.'
    C = {}
    C['xx'] = xx
    C['nx'] = nx
    C['va'] = va
    C['tmax'] = tmax
    C['dt'] = dt
    C['avep'] = avep
    C['period'] = period
    C['Cs'] = Cs
    C['Cp'] = Cp
    C['ra'] = ra
    C['ro'] = ro
    C['ri'] = ri
    C['Cpo'] = Cpo
    C['T0'] = T0
    C['dT'] = dT
    C['H'] = H
    C['vo'] = vo
    C['Hb'] = Hb
    C['Cpi'] = Cpi
    C['Li'] = Li
    C['Tf'] = Tf
    C['B'] = B
    C['emissivity'] = emissivity
    C['SW0'] = SW0
    C['SW_anom'] = SW_anom
    C['Da'] = Da
    C['Do'] = Do
    C['tau_entrainment'] = tau_entrainment
    for var in args.keys():
        C[var] = args[var]
    return C
Setup model constants. All of the constants have fixed values, but one can pass in own values or even some arbitrary values via **args.
def set_model_constants(xx=50000.0, nx=100, va=10.0, tmax=(((60 * 360) * 24) * 3600.0), avep=(24 * 3600.0), dt=3600.0, period=(((3600 * 24) * 360) * 1), B=2.0, T0=(273.15 + 6), dT=2.0, Cs=0.001, Cp=1030.0, ra=1.5, ro=1030.0, ri=900.0, Cpo=4000.0, Cpi=2900.0, H=200.0, vo=0.2, Hb=1000.0, Li=3300000.0, Tf=(273.15 - 1.8), SW0=50.0, SW_anom=100.0, emissivity=0.99, Da=1000000.0, Do=500.0, tau_entrainment=((30 * 24) * 3600.0), **args):
    
    C = {}
    C['xx'] = xx
    C['nx'] = nx
    C['va'] = va
    C['tmax'] = tmax
    C['dt'] = dt
    C['avep'] = avep
    C['period'] = period
    C['Cs'] = Cs
    C['Cp'] = Cp
    C['ra'] = ra
    C['ro'] = ro
    C['ri'] = ri
    C['Cpo'] = Cpo
    C['T0'] = T0
    C['dT'] = dT
    C['H'] = H
    C['vo'] = vo
    C['Hb'] = Hb
    C['Cpi'] = Cpi
    C['Li'] = Li
    C['Tf'] = Tf
    C['B'] = B
    C['emissivity'] = emissivity
    C['SW0'] = SW0
    C['SW_anom'] = SW_anom
    C['Da'] = Da
    C['Do'] = Do
    C['tau_entrainment'] = tau_entrainment
    for var in args.keys():
        C[var] = args[var]
    return C