Commit ec8e16d2 authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

fix parameter tweaking

parent e94fa37b
Pipeline #22425 passed with stage
in 6 minutes and 37 seconds
...@@ -694,15 +694,17 @@ class DNS(_code): ...@@ -694,15 +694,17 @@ class DNS(_code):
if type(opt.dealias_type) == type(None): if type(opt.dealias_type) == type(None):
opt.dealias_type = self.parameters['dealias_type'] opt.dealias_type = self.parameters['dealias_type']
self.parameters['nu'] = (opt.kMeta * 2 / opt.n)**(4./3) self.parameters['nu'] = (opt.kMeta * 2 / opt.n)**(4./3)
# check value of kMax
kM = opt.n * 0.5
if opt.dealias_type == 1:
kM *= 0.8
# tweak forcing/viscosity based on forcint type
if opt.forcing_type == 'linear': if opt.forcing_type == 'linear':
# custom famplitude for 288 and 576 # custom famplitude for 288 and 576
if opt.n == 288: if opt.n == 288:
self.parameters['famplitude'] = 0.45 self.parameters['famplitude'] = 0.45
elif opt.n == 576: elif opt.n == 576:
self.parameters['famplitude'] = 0.47 self.parameters['famplitude'] = 0.47
kM = opt.n * 0.5
if opt.dealias_type == 1:
kM *= 0.8
elif opt.forcing_type == 'fixed_energy_injection_rate': elif opt.forcing_type == 'fixed_energy_injection_rate':
# use the fact that mean dissipation rate is equal to injection rate # use the fact that mean dissipation rate is equal to injection rate
self.parameters['nu'] = ( self.parameters['nu'] = (
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment