Changeset 3507


Ignore:
Timestamp:
Aug 17, 2006, 6:20:53 PM (18 years ago)
Author:
ole
Message:

Work post suresh kumar visit

File:
1 edited

Legend:

Unmodified
Added
Removed
  • development/stochastic_study/plot_spread.py

    r3035 r3507  
    4040#hold(False)
    4141
    42 
    43 for j in range(1):
     42N=250
     43for j in range(N):
    4444    #print j, filenames[j]
    4545    plot(time, data[:,j], 'k-')
     
    5454    #raw_input('Next')
    5555
    56 
    57 raw_input('Stats')   
     56title('%d realisations' %N)
     57savefig('realisations', dpi = 300)
     58#raw_input('Stats')   
    5859 
    5960
     
    6162hold(False)
    6263# Plot spread of stage values for each timestep
    63 for i in range(0,project.number_of_timesteps):
     64for i in range(0, project.number_of_timesteps, 10):
    6465    # Plot histogram#
    6566
     
    6869    xlabel('realisations')
    6970    ylabel('stage (m)')
    70     title('Study %s: spread at timestep %d of %d (t=%.2f)'\
    71           %(study, i, project.number_of_timesteps-1, time[i]))       
    72 raw_input('Next')
     71    #title('Study %s: spread at timestep %d of %d (t=%.2f)'\
     72    #      %(study, i, project.number_of_timesteps-1, time[i]))
     73    title('Spread at timestep %d of %d (t=%.2f)'\
     74          %(i, project.number_of_timesteps-1, time[i]))
     75
     76    savefig('spread_%d' %i, dpi = 300)   
     77    #raw_input('Next')
    7378
    7479show()
Note: See TracChangeset for help on using the changeset viewer.