| 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'); |