class MyParticle(ptype[mode]):
sample_var = Variable('sample_var')
class MyParticle(ptype[mode]):
# JIT compilation requires a-priori knowledge of the particle
# data structure, so we define additional variables here.
p = Variable('p', dtype=np.float32, initial=0.)
p_start = Variable('p_start', dtype=np.float32, initial=fieldset.P)
class MyParticle(ptype[mode]):
sample_var = Variable('sample_var', initial=fieldset.U)
class MyParticle(ptype[mode]):
v_once = Variable('v_once',
dtype=np.float64,
initial=0.,
to_write='once')
age = Variable('age', dtype=np.float32, initial=0.)
class PumiceParticle(JITParticle):
age = Variable('age', dtype=np.float32, initial=0.)
class TempParticle(ptype[mode]):
temp = Variable('temp', dtype=np.float32)
pressure = Variable('pressure', dtype=np.float32, initial=10000)
class DensParticle(ptype[mode]):
density = Variable('density', dtype=np.float32)
class AngleParticle(ptype[mode]):
angle = Variable('angle')
class SampleParticle(ptype[mode]):
u = Variable('u', dtype=np.float32)
v = Variable('v', dtype=np.float32)
p = Variable('p', dtype=np.float32)
class TestParticle(ptype[mode]):
p = Variable('p', np.float32, initial=0.33)
p2 = Variable('p2', np.float32, initial=1, to_write=False)
p3 = Variable('p3', np.float32, to_write='once')
class PlasticParticle(JITParticle):
age = Variable('age', dtype=np.int32, initial=0)
class MyParticle(ptype[mode]):
v = Variable('v')
class SampleParticle(ptype['scipy']):
mask = Variable('mask', initial=fieldset.mask)
class TestParticle(ptype[mode]):
p = Variable('p', np.float32, initial=0.33)
n = Variable('n', np.int32, initial=2)
class TestParticle(ptype[mode]):
p = Variable('p',
dtype=np.float32 if rngfunc == 'randint' else np.float32)
class VarParticle(pclass(mode)):
a = Variable('a', dtype=np.float32, initial=fieldset.P)
class TestParticle(ptype[mode]):
age = Variable('age', dtype=np.float32)
class SampleParticle(JITParticle):
u_kernel = Variable('u_kernel', dtype=np.float32, initial=0.)
u_scipy = Variable('u_scipy', dtype=np.float32, initial=fieldset.U)
class TPressureParticle(ptype[mode]):
pressure = Variable('pressure', dtype=np.float32)
class turtle(JITParticle):
#Sampling
if param['key_alltracers']:
T = Variable('T', to_write=True, dtype=np.float32)
NPP = Variable('NPP', to_write=True, dtype=np.float32)
#All particles
active = Variable('active', to_write=True, initial=1, dtype=np.float32)
prev_lon = Variable('prev_lon',
to_write=False,
dtype=np.float64,
initial=attrgetter('lon'))
prev_lat = Variable('prev_lat',
to_write=False,
dtype=np.float64,
initial=attrgetter('lat'))
u_current = Variable('u_current', to_write=True, dtype=np.float32)
v_current = Variable('v_current', to_write=True, dtype=np.float32)
age = Variable('age', to_write=True, dtype=np.float64, initial=0.)
onland = Variable('onland',
to_write=False,
dtype=np.float32,
initial=0.) #number of beachings in a row
beached = Variable('beached',
to_write=False,
dtype=np.float32,
initial=0.) #0=ocean, 1=onland
t = Variable('t',
to_write=False,
dtype=np.float32,
initial=fieldset.tactic_factor)
SCL = Variable('SCL',
to_write=True,
dtype=np.float32,
initial=fieldset.SCL0)
if param['mode'] == 'passive':
u_swim = Variable('u_swim', to_write=False, dtype=np.float32)
v_swim = Variable('v_swim', to_write=False, dtype=np.float32)
#Mortality
if param['cold_death']:
lethargy_time = Variable('lethargy_time',
to_write=False,
dtype=np.float32,
initial=0.) #time spent under Tmin
cold_death = Variable('cold_death',
to_write=True,
dtype=np.float32,
initial=0)
if param['cold_death'] or param['mode'] == 'active':
Tmin = Variable('Tmin',
to_write=False,
dtype=np.float32,
initial=fieldset.Tmin)
#Active
if param['mode'] == 'active':
u_swim = Variable('u_swim', to_write=True, dtype=np.float32)
v_swim = Variable('v_swim', to_write=True, dtype=np.float32)
M = Variable('M', to_write=False, dtype=np.float32)
vmax = Variable('vmax', to_write=False, dtype=np.float32)
PPmax = Variable('PPmax', to_write=True, dtype=np.float32)
habT = Variable('habT', to_write=True, dtype=np.float32)
habPP = Variable('habPP', to_write=True, dtype=np.float32)
hab = Variable('hab', to_write=True, dtype=np.float32)
theta = Variable('theta', to_write=False, dtype=np.float32)
xgradh = Variable('xgradh', to_write=True, dtype=np.float32)
ygradh = Variable('ygradh', to_write=True, dtype=np.float32)
Topt = Variable('Topt',
to_write=False,
dtype=np.float32,
initial=fieldset.Topt)
if param['growth'] == 'Gompertz':
K = Variable('K',
to_write=False,
dtype=np.float32,
initial=fieldset.K0)
if param['frenzy'] or param['wave_swim']:
frenzy_speed = Variable('frenzy_speed',
to_write=False,
dtype=np.float32)
if param['wave_swim']:
theta_wave = Variable('theta_wave',
to_write=False,
dtype=np.float32)
class TestParticle(ptype[mode]):
p = Variable('p', dtype=np.float32)
class plastic_particle(JITParticle): #ScipyParticle): #
u = Variable('u', dtype=np.float32, to_write=True)
v = Variable('v', dtype=np.float32, to_write=True)
w = Variable('w', dtype=np.float32, to_write=True)
temp = Variable('temp', dtype=np.float32, to_write=False)
density = Variable('density', dtype=np.float32, to_write=True)
tpp3 = Variable('tpp3', dtype=np.float32, to_write=False)
d_phy = Variable('d_phy', dtype=np.float32, to_write=False)
nd_phy = Variable('nd_phy', dtype=np.float32, to_write=False)
a = Variable('a', dtype=np.float32, to_write=False)
kin_visc = Variable('kin_visc', dtype=np.float32, to_write=False)
sw_visc = Variable('sw_visc', dtype=np.float32, to_write=False)
vs = Variable('vs', dtype=np.float32, to_write=True)
rho_tot = Variable('rho_tot', dtype=np.float32, to_write=True)
r_tot = Variable('r_tot', dtype=np.float32, to_write=True)
delta_rho = Variable('delta_rho', dtype=np.float32, to_write=True)
vs_init = Variable('vs_init', dtype=np.float32, to_write=True)
r_pl = Variable('r_pl', dtype=np.float32, to_write='once')
rho_pl = Variable('rho_pl', dtype=np.float32, to_write='once')
class MyParticle(ptype[mode]):
sample_var = Variable('sample_var', initial=0.)
v_once = Variable('v_once',
dtype=np.float64,
initial=0.,
to_write='once')
class MyParticle(JITParticle):
psi = Variable('psi', dtype=np.float32, initial=fieldset.Psi)
class MyParticle(ptype[mode]):
temp = Variable('temp', dtype=np.float32, initial=20.)
u1 = Variable('u1', dtype=np.float32, initial=0.)
u2 = Variable('u2', dtype=np.float32, initial=0.)
v1 = Variable('v1', dtype=np.float32, initial=0.)
v2 = Variable('v2', dtype=np.float32, initial=0.)
class TestParticle(ptype[mode]):
p0 = Variable('p0', dtype=np.int32, initial=0)
p1 = Variable('p1', dtype=np.int32, initial=1)
class MyParticle(ParticleClass):
p = Variable('p', dtype=np.float32, initial=0.)
p_start = Variable('p_start', dtype=np.float32, initial=fieldset.P)
class TestParticle(ptype[mode]):
p = Variable('p', dtype=np.int32, initial=0)
class MyParticle(ptype[mode]):
temp = Variable('temp', dtype=np.float32, initial=20.)
class TestParticle(ptype[mode]):
p_float = Variable('p_float', dtype=np.float32, initial=10.)
p_double = Variable('p_double', dtype=np.float64, initial=11.)
p_int = Variable('p_int', dtype=np.int32, initial=12.)
