1 |
edhill |
1.1 |
clear |
2 |
|
|
clear path |
3 |
|
|
|
4 |
|
|
global Nx Ny Nz |
5 |
|
|
global lat long dz dm mdep |
6 |
|
|
global delt_su su_its t_su delt |
7 |
|
|
global descriptor this_path |
8 |
|
|
global f deltaf Q beta r_expt r_heat H |
9 |
|
|
global time rots it |
10 |
|
|
global g Cp rho_bar alpha |
11 |
|
|
global u v t w |
12 |
|
|
global iterations |
13 |
|
|
|
14 |
|
|
|
15 |
|
|
param_file_name = ... |
16 |
|
|
input(' Please enter the name of the m-file with the parameters for this run : ','s') ; |
17 |
|
|
feval(param_file_name) ; |
18 |
|
|
|
19 |
|
|
iterations |
20 |
|
|
|
21 |
|
|
it = input(' Please enter iteration : ','s') |
22 |
|
|
|
23 |
|
|
path = this_path |
24 |
|
|
cmdstr=['cd ' path ]; |
25 |
|
|
eval(cmdstr); |
26 |
|
|
path=pwd |
27 |
|
|
|
28 |
|
|
wfilename=(['wVeltave.' it ]); |
29 |
|
|
w=rdmeta(wfilename,'b'); |