Changeset 4949


Ignore:
Timestamp:
Jan 17, 2008, 11:49:00 AM (16 years ago)
Author:
sexton
Message:

(a) minor updates for Cairns example in manual (b) scripts for stratification investigation for other model areas

Files:
6 added
5 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/anuga_user_manual.tex

    r4948 r4949  
    13281328\file{GetTimeseries.py} which generates the csv files for each point location. The csv files
    13291329can then be used in \code{csv2timeseries_graphs} to create the timeseries plot for each desired
    1330 quantity. \code{csv2timeseries_graphs} relies on \code{pylab} to be install which is not part
     1330quantity. \code{csv2timeseries_graphs} relies on \code{pylab} to be installed which is not part
    13311331of the standard \code{anuga} release, however it can be downloaded and installed from \code{http://matplotlib.sourceforge.net/}
    13321332
     
    13401340Figure \ref{fig:reef} shows the time series for the quantity stage for the
    13411341Elford Reef location for each scenario (the elevation at this location is negative,
    1342 therefore stage is the most appropriate to use). Note the large negative stage value when the slide was
     1342therefore stage is the more appropriate quantity to plot). Note the large negative stage value when the slide was
    13431343introduced. This is due to the double gaussian form of the initial surface
    13441344displacement of the slide. By contrast, the time series for depth is shown for the onshore location of the Cairns
  • anuga_core/documentation/user_manual/demos/cairns/GetTimeseries.py

    r4948 r4949  
    4343                            verbose=True)
    4444except ImportError:
    45     #ANUGA don't need pylab to work so the system doesn't
    46     #rely on pylab being installed
    47     print 'must have pylab install to generate plots'
     45    #ANUGA does not rely on pylab to work
     46    print 'must have pylab installed to generate plots'
    4847   
    4948
  • anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

    r4936 r4949  
    12351235                        units = 'm / sec'
    12361236                    if which_quantity == 'bearing':
    1237                         #due_east = 90.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
    1238                         #due_west = 270.0*ones(shape(model_time[0:n[j]-1,k,j],Float))
    1239                         print 'hello', bearings
    1240                         print 'east', due_east
    12411237                        plot(model_time[0:n[j]-1,k,j], bearings[0:n[j]-1,k,j], '-',
    12421238                             model_time[0:n[j]-1,k,j], due_west[0:n[j]-1], '-.',
  • anuga_work/production/pt_hedland_2006/convert2latlon.py

    r4927 r4949  
    3131
    3232timedir='20070613_044204_run_final_3.6_exmouth_nbartzis'
    33 filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
    34 #filename = project.outputdir + timedir + sep + 'stage5_interp.csv'
     33#timedir='20070613_061053_run_final_3.6_onslow_nbartzis'
     34
     35#filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
     36filename = project.outputdir + timedir + sep + 'stage5_interp.csv'
    3537#filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
    36 #filename = project.outputdir+  timedir + sep + 'stage80_interp.csv'
    37 filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
    38 #filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'
     38#filename = project.outputdir+  timedir + sep + 'stage50_interp.csv'
     39#filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
     40filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'
    3941#filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
    40 #filename_out = project.outputdir + timedir + sep + 'stage80_convert.csv'
     42#filename_out = project.outputdir + timedir + sep + 'stage50_convert.csv'
    4143alter_file(filename,filename_out)
  • anuga_work/production/pt_hedland_2006/plot_contour_data.py

    r4927 r4949  
    88directory = project.outputdir
    99timedir='20070613_044204_run_final_3.6_exmouth_nbartzis'
     10#timedir='20070613_061053_run_final_3.6_onslow_nbartzis'
    1011
    1112def get_data(filename):
     
    1920    return max(stage), min(stage)
    2021
    21 files = ['stage0_convert.csv','stage5_convert.csv','stage20_convert.csv','stage80_convert.csv']
    22 depth = [0,5,20,80]
     22files = ['stage0_convert.csv','stage5_convert.csv','stage20_convert.csv','stage50_convert.csv']
     23#files = ['stage0_convert.csv','stage20_convert.csv','stage50_convert.csv']
     24depth = [0,5,20,50]
     25#depth = [0,20,50]
    2326stage_max = []
    2427stage_min = []
     
    3336from anuga.abstract_2d_finite_volumes.util import greens_law
    3437from Numeric import arange
    35 d1 = 80.
     38d1 = 50.
    3639d2 = arange(d1,0.15,-0.1)
    3740h1 = max_stage
     
    5356#Magnitude 9.0 event (Sumba) at Pt Hedland')
    5457legend(['ANUGA','Green\'s law'])
    55 axis([-5,85,min(stage_max)*0.9,max(stage_max)*1.1])
     58axis([-5,55,min(stage_max)*0.9,max(stage_max)*1.1])
    5659grid(True)
    5760#savefig('stratification_onslow_gun')
Note: See TracChangeset for help on using the changeset viewer.