Commit 6f3772f2 authored by Julio's avatar Julio
Browse files

started changing code to consider 3100 ICs but did not yet succeeded.

parent a12c454b
No preview for this file type
......@@ -5,6 +5,8 @@ import numpy as np
import h5py
import scipy.io
from a03_FIELD_BICKLY import velocity_field_Bickley
"""
Created on Thu Mar 24 15:47:47 2022
@author: Julio Urizarna Carasa
......@@ -47,7 +49,7 @@ save_plot_to = '/home/julio/Documentos/PhD/02_projects/GitLab/maxey-riley/scr
################################
Case_v = np.array(["Bickley", "Vortex", "Faraday", "Test"])
Case_elem = 1
Case_elem = 0
......@@ -55,15 +57,15 @@ Case_elem = 1
# Define time grid #
####################
tini = 0.0
tini = 10.0
if Case_v[Case_elem] == Case_v[2]:
# Time domain for Faraday field is restricted according to available data.
tend = 3.44
nt = 200
else:
tend = 5.00
nt = 1000
tend = 30.0
nt = 3001
taxis = np.linspace(tini, tend, nt)
......@@ -107,19 +109,27 @@ x_fd_v = -c * np.log(1.0 - xi_fd_v)
# Decide whether to apply parallel computing #
##############################################
parallel_flag = False #True
parallel_flag = True
########################################
# Define particle's initial conditions #
########################################
'''
mat = scipy.io.loadmat('IniCond.mat')
x0 = mat['X'][0]
y0 = mat['Y'][0]
u0 = mat['vx'][0]
v0 = mat['vy'][0]
#u0 = mat['vx'][0]
#v0 = mat['vy'][0]
if Case_v[Case_elem] == Case_v[0]:
'''Bickley Jet velocity field'''
vel = velocity_field_Bickley()
u0 = vel.get_velocity(x0, y0, tini)
v0 = vel.get_velocity(x0, y0, tini)
'''
......@@ -129,7 +139,7 @@ if Case_v[Case_elem] == Case_v[2]:
else:
x0, y0 = 1.0, 0.0
u0, v0 = 0.15, 0.15
'''
############################################
......
Markdown is supported
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