Changeset 2986


Ignore:
Timestamp:
May 26, 2006, 3:53:42 PM (18 years ago)
Author:
ole
Message:

More stochastic plotting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • development/stochastic_study/plot_spread.py

    r2979 r2986  
    5050# Simple plot of timeseries for different realisations
    5151for j in range(number_of_realisations):
    52     #print j, data[:,j]
     52    print j, data[:,j]
    5353    plot(time, data[:,j], 'k-')
    5454   
    5555    xlabel('time(s)')
    5656    ylabel('stage (m)')
    57     title('Realisation %d of %d' %(j, number_of_realisations-1))   
     57    title('Realisation %d of %d' %(j, number_of_realisations-1))    #
     58#raw_input('Next')
    5859
     60
     61hold(False)
     62# Plot spread of stage values for each timestep
     63for i in range(200,project.number_of_timesteps):
     64    # Plot histogram
     65
     66    w = data[i,:]
     67    plot(w, 'k.')
     68    xlabel('realisations')
     69    ylabel('stage (m)')
     70    title('Timestep %d of %d (t=%.2f)'\
     71          %(i, project.number_of_timesteps-1, time[i]))       
    5972    raw_input('Next')
    6073
    6174
    62 
     75hold(False)
    6376# Plot histogram of stage values for each timestep
    64 for i in range(project.number_of_timesteps):
     77for i in range(200,project.number_of_timesteps):
    6578    # Plot histogram
    6679
     
    7891    title('Timestep %d of %d (t=%.2f)'\
    7992          %(i, project.number_of_timesteps-1, time[i]))       
    80     raw_input('Next')
     93    #raw_input('Next')
    8194
    8295
Note: See TracChangeset for help on using the changeset viewer.