source: anuga_validation/convergence_study/animatesww2d_alt.py @ 4695

def animatesww2d_alt(sww_filename, movie_filename, range):
    """Plot cross section of model output
    """

    # Read SWW file     
    from Scientific.IO.NetCDF import NetCDFFile
    fid = NetCDFFile(sww_filename, 'r')
    
    x = fid.variables['x']
    y = fid.variables['y']
    volumes = fid.variables['volumes']
    
    elevation = fid.variables['elevation']
    time = fid.variables['time']
    stage = fid.variables['stage']
    xmomentum = fid.variables['xmomentum']
    ymomentum = fid.variables['ymomentum']
    
    
    # Set up plot

    spacing=10 # grid interval in metres
    hmax=5     # ??

    
#figure(1);

#if nargin > 1
#    figure(1);set(gcf,'Position',[200 400 800 300]);
#    mov=avifile(movie_filename,'FPS',5)
#end

#i=find(V_y==0);

#[xi,yi]=meshgrid(min(V_x):spacing:max(V_x),min(V_y):spacing:max(V_y));

#elev=griddata(V_x,V_y,V_elevation,xi,yi);

#if nargin < 3
#    range=[min(V_x) max(V_x) min(min(V_stage))-1 max(max(V_stage))+1];
#end

#m=V_stage(1,i);

#for t=1:length(V_time)
    
#    m=max(m,V_stage(t,i));
#    plot(V_x(i),V_stage(t,i),V_x(i),V_elevation(i),V_x(i),m,'.');

#    axis(range)
    

#    title(num2str(V_time(t)));drawnow;

#    if nargin>1
#        F=getframe(gcf);
#        mov=addframe(mov,F);   
#    end
   
#end

#if nargin > 1
#    mov=close(mov);
#end


if __name__ == '__main__':

    import sys
    animatesww2d_alt(sys.argv[1], 'cross_section_plot', 5)