function experiment

global Nx Ny Nz
global dx dy lat long 
global dz depth
global delt iteration it
global descriptor 
global this_path
global g Cp rho_bar alpha
global f0 B0 Q0 beta r_heat r_expt H 

% prescribed parameters
Nx = 124; Ny = 124; Nz = 30;
dx = 0.01; dy = 0.01; dz = 0.005;  
delt = 0.1; 
iteration = ['0000250209'; ...
'0000251210'; ...
'0000253609'; ...
'0000254610'; ...
'0000255610'; ...
'0000256610'; ...
'0000257009'; ...
'0000258010'; ...
'0000259010'; ...
'0000260010'; ...
'0000260409'; ...
'0000261410'; ...
'0000262410'; ...
'0000263410'; ...
'0000264411'; ...
'0000265411'; ...
'0000266411'; ...
'0000266810'; ...
'0000267811'; ...
'0000268811'; ...
'0000269811'; ...
'0000270210'; ...
'0000274611'; ...
'0000275611'; ...
'0000276611'; ...
'0000277010'; ...
'0000278011'; ...
'0000279011'; ...
'0000280011'; ...
'0000280410'; ...
'0000281411'; ...
'0000282411'; ...
'0000283411'; ...
'0000283910'; ...
'0000284911'; ...
'0000285911'; ...
'0000286911'; ...
'0000287410'; ...
'0000288411'; ...
'0000289411'; ...
'0000290411'; ...
'0000290910'; ...
'0000291911'; ...
'0000292911'; ...
'0000293911'; ...
'0000294410'; ...
'0000295411'; ...
'0000296411'; ...
'0000297411'; ...
'0000297910'; ...
'0000298911'; ...
'0000299911'; ...
'0000300911'; ...
'0000301410'; ...
'0000302411'; ...
'0000303411'; ...
'0000304411'; ...
'0000304810'; ...
'0000305811'; ...
'0000306811'; ...
'0000307811'; ...
'0000308310'; ...
'0000310311'; ...
'0000311810'; ... 
'0000312811'; ...
'0000313811'; ...
'0000314811'; ...
'0000316211'; ...
'0000317211'; ...
'0000318211'; ...
'0000318510'; ...
'0000319511'; ...
'0000320511'; ...
'0000321511'; ...
'0000321810'; ...
'0000322811'; ...
'0000323811'; ...
'0000324811'; ...
'0000325010'; ...
'0000326011'; ...
'0000327011'; ...
'0000328011'; ...
'0000329511'; ...
'0000330511'; ...
'0000331511'; ...
'0000331910'; ...
'0000332911'; ...
'0000334911'; ...
'0000335310'; ...
'0000336311'; ...
'0000337311'; ...
'0000338610'; ...
'0000339611'];
descriptor = 'reference donut';
this_path  = ['/d04/helen/tank/R15donutprep03/dumps'];
g = 9.81; Cp = 4000; rho_bar = 999.8; alpha = 2e-4;
f0 = 0.75; Q0 = 500; beta = 0; r_heat = 0.2; r_expt = 0.48; 

% derived  parameters
lat = [1:1:Nx]; 
long = [1:1:Ny];
% depth = -cumsum(dz)+dz(1)/2;
%b0 = g*alpha*Q0/rho_bar*Cp;
% H = depth(Nz)-dz(Nz)/2;