forked from zhixin-cheng/TrackMPD-v.2
-
Notifications
You must be signed in to change notification settings - Fork 0
/
input_conf.m
194 lines (142 loc) · 6.94 KB
/
input_conf.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
function conf = input_conf(varargin)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Inputs data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
conf.Data.BaseDir = [pwd '\InputDataPOM']; %Folder containing inputs
conf.Data.Domain = [conf.Data.BaseDir '\coast_clean.dat'] ; % domain file (full path)
conf.Data.ParticlesFile = [conf.Data.BaseDir '\particles_JB_scenario1.csv']; % particle release file (full path)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Ocean General Ciculation model (OGCM)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
conf.OGCM.Model_name = 'POM'; %Select POM or FVCOM
conf.OGCM.DomainName = 'JervisBay';
conf.OGCM.BaseDir = [pwd '\InputDataPOM']; %Folder containing OGCM inputs
conf.OGCM.TimeStep = 1/24; %in days
conf.OGCM.VerticalLayer = 'sigma2depthVar'; % Only for 3D mode
% Options: 'Rectangular'
% 'Sigma2depthCte': cte over time, independient of Elevation
% 'Sigma2depthVar': variable over time,dependient of Elevation
if strcmpi(conf.OGCM.Model_name,'POM') % When using POM outputs
%%%Particular inputs for POM (SARCCM version)
conf.OGCM.POM_Prefix = 'u_h';
conf.OGCM.POM_Suffix = '.nc';
conf.OGCM.t0 = datenum(1998,06,26,00,0,0); %POM reference time variable
%(date of the first time step:
%datenum(y,m,d,H,M,S)
conf.OGCM.Hmin = 1; %Depth for land
%%%
elseif strcmpi(conf.OGCM.Model_name,'FVCOM') % When using FVCOM outputs
%%%Particular inputs for FVCOM
conf.OGCM.FVCOMFile = [conf.OGCM.BaseDir '\YRE_0001.nc'];
conf.OGCM.FVCOMGrid = [conf.OGCM.BaseDir '\YRE_gdr.nc'];
conf.OGCM.NumLonGrid = 200; %number of points in the new rectangular grid (lon dimension)
conf.OGCM.NumLatGrid = 200; %number of points in the new rectangular grid (lat dimension)
conf.OGCM.NameTime = 't'; %name of variable time
conf.OGCM.NameW = 'ww'; %name of variable for vertical velocity (3D mode)
%%%
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Trajectory
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Mode
conf.Traj.Mode = '3D'; %'2D' or '3D' tracking
% Time parameters
conf.Traj.ReleaseTime = datenum(1998,6,26,0,0,0); %datenum(y,m,d,H,M,S)
conf.Traj.TrajectoryDuration = 1; % in days
conf.Traj.TimeStep = 60/24/60; % in days
conf.Traj.Direction = 'forward'; %forward or backward
% Output
conf.Traj.ScenarioName = [conf.OGCM.DomainName '_ReleasePoint1']; %Free choice
conf.Traj.BaseDir = [pwd '\Outputs']; %Folder to save outputs
% Chunk Method: to avoid memory problems, the time domain can be divided in
% different partitions (chunks) with a duration of duration of "chunklen"
% (e.g. 5 days) to calculate advection. Necessary for long simulations
% and/or high resolution grids.
conf.Traj.chunklen = 2; % duration of partitions (in days)
% Dispersion
conf.Traj.Kh = 1; %m2/s
conf.Traj.Kv = 0.00001; %0.00001; %m2/s For 3D mode
% Beaching
conf.Traj.Beaching = 'yes'; %'yes' or 'no'
% Refloating
conf.Traj.Refloating = 'no'; %'yes' or 'no'
if strcmpi(conf.Traj.Refloating,'yes')
conf.Traj.Tw = 2; % half-live of particle to remain on the beach before washing off again (in days)
conf.Traj.xcoord_elev = 38; %grid cell (x-coord) to calculate the time series of elevation
conf.Traj.ycoord_elev = 73; %grid cell (y-coord) to calculate the time seriss of elevation
end
% % Windage
% conf.Traj.Windage = 'no'; %'yes' or 'no'
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Behaviour (3D mode)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%BehaviourType:
% 1: macro or microplastic, plastic density<1 g cm^-3, no biofouling, no degradation
% 2: macroplastic (or microplastic with a known ws), plastic density>1 g cm^-3, no biofouling, no degradation
% 3: microplastic (shape,size,density), plastic density>1 g cm^-3, no biofouling, no degradation
% 4: microplastic (shape,size,density), Biofouling (stationary), no degradation
% 5: microplastic (shape,size,density), Biofouling (non-stationary), no degradation
% 6: microplastic (shape,size,density), No Biofouling, degradation
BehaviourType = 1;
if BehaviourType==1
conf.beh=behaviour1;
%%%Inputs (only useful if windage)
conf.beh.PolymerDensity = 0.9; %g cm^-3 (default value: 0.9)
%%%
elseif BehaviourType==2
conf.beh=behaviour2;
%%%Inputs
conf.beh.Ws = 0.0005;
%%%
elseif BehaviourType==3
conf.beh=behaviour3;
%%%Inputs
conf.beh.PolymerType = 'acrylonitrile-butadiene-styrene'; %Not mandatory
conf.beh.PolymerAcronym = 'ABS'; %Not mandatory
conf.beh.PolymerDensity = 1.665; %g cm^-3
conf.beh.Category = 'LongCylinder'; % (Sphere, LongCylinder,ShortCylinder)
conf.beh.Shape = Cylinder(0.15*10^(-3),0.3*10^(-3)); %Sphere(R) or Cylinder(R,L) (m)
conf.beh.WaterDensity = 1.025; %g cm^-3)
conf.beh.WaterViscosity = 10^(-6); % m^2/s
%%%
elseif BehaviourType==4
conf.beh=behaviour4;
%%% Inputs
conf.beh.PolymerType = 'acrylonitrile-butadiene-styrene'; %Not mandatory
conf.beh.PolymerAcronym = 'ABS'; %Not mandatory
conf.beh.PolymerDensity = 1.026; %g cm^-3
conf.beh.Category = 'Sphere'; % (Sphere, LongCylinder,ShortCylinder)
conf.beh.Shape = Sphere(0.15*10^(-3)); %Sphere(R) or Cylinder(R,L) (m)
conf.beh.BiofoulingThickness = 0.01*10^(-3); %meters
conf.beh.FilmDensity = 1.05; %g cm^-3
conf.beh.WaterDensity = 1.025; %g cm^-3
conf.beh.WaterViscosity = 10^(-6); % m^2/s
%%%
elseif BehaviourType==5
conf.beh=behaviour5;
%%% Inputs
conf.beh.PolymerType = 'acrylonitrile-butadiene-styrene'; %Not mandatory
conf.beh.PolymerAcronym = 'ABS'; %Not mandatory
conf.beh.PolymerDensity = 1.026; %g cm^-3
conf.beh.Category = 'LongCylinder'; % (Sphere, LongCylinder,ShortCylinder)
conf.beh.Shape = Cylinder(0.15*10^(-3),0.3*10^(-3)); %Sphere(R) or Cylinder(R,L) (m)
conf.beh.BiofoulingThickness0 = 0*10^(-3); %meters
conf.beh.FilmDensity = 1.05; %g cm^-3
conf.beh.BiofoulingRate = 0.005*10^(-3); %meters/day;
conf.beh.WaterDensity = 1.025; %g cm^-3
conf.beh.WaterViscosity = 10^(-6); % m^2/s
%%%
elseif BehaviourType==6
conf.beh=behaviour6;
%%%Inputs
conf.beh.PolymerType = 'acrylonitrile-butadiene-styrene'; %Not mandatory
conf.beh.PolymerAcronym = 'ABS'; %Not mandatory
conf.beh.PolymerDensity = 1.026; %g cm^-3
conf.beh.Category = 'Sphere'; % (Sphere, LongCylinder,ShortCylinder)
conf.beh.Shape0 = Sphere(1*10^(-3)); %Sphere(R) or Cylinder(R,L) (m)
conf.beh.DegradationRate = 30; %Percentage of size decrease per day
conf.beh.WaterDensity = 1.025; %g cm^-3)
conf.beh.WaterViscosity = 10^(-6); % m^2/s
%%%
end
end