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