function plotfields( iter )
clear 
clear path 

global Nx Ny Nz
global lat long dz dm mdep
global delt_su su_its t_su delt
global descriptor this_path
global f deltaf Q beta r_expt r_heat H
global time rots it
global g Cp rho_bar alpha
global u v t w
global iterations


for j=1:5:Ny
j
figure(1)
subplot(2,1,1)
pcolor(flipud(transpose(squeeze(t(:,j,:)))));shading flat; colorbar;
title(['t timestep ' num2str(eval(it)) '    y section ' int2str(j) ]);
drawnow
subplot(2,1,2)
pcolor(flipud(transpose(squeeze(w(:,j,:)))));shading flat; colorbar; 
title(['w timestep ' num2str(eval(it)) '    y section ' int2str(j) ]);
drawnow

pause

end

return