Changeset 4927


Ignore:
Timestamp:
Jan 10, 2008, 6:04:35 PM (16 years ago)
Author:
sexton
Message:

scripts for stratification investigation

Location:
anuga_work/production
Files:
2 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/onslow_2006/convert2latlon.py

    r4887 r4927  
    3030from os import sep
    3131
    32 timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
    33 #filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
     32#timedir = '20070608_060316_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
     33timedir = '20070613_034754_run_final_1.5_dampier_nbartzis'
     34#timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
     35filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
    3436#filename = project.outputdir + timedir + sep + 'stage5_interp.csv'
    35 filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
     37#filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
    3638#filename = project.outputdir+  timedir + sep + 'stage80_interp.csv'
    37 #filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
     39filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
    3840#filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'
    39 filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
     41#filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
    4042#filename_out = project.outputdir + timedir + sep + 'stage80_convert.csv'
    4143alter_file(filename,filename_out)
  • anuga_work/production/onslow_2006/plot_contour_data.py

    r4890 r4927  
    33
    44import project
    5 from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend
     5from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend, grid
    66from os import sep
    77 
    88directory = project.outputdir
    99timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
     10#timedir = '20070608_060316_run_final_1.5_onslow_nbartzis'
     11#timedir = '20070613_034754_run_final_1.5_dampier_nbartzis'
    1012
    1113def get_data(filename):
     
    3234
    3335from anuga.abstract_2d_finite_volumes.util import greens_law
     36from Numeric import arange
    3437d1 = 80.
    35 d2 = range(d1,5,-1)
     38d2 = arange(d1,0.15,-0.1)
    3639h1 = max_stage
     40#h1 = 4.0
    3741green = []
    3842for d in d2:
     
    4650title('ANUGA modelled maximum stage versus Green\'s approximation \n \
    4751Magnitude 9.3 event (Java) at Onslow')
     52#title('ANUGA modelled maximum stage versus Green\'s approximation \n \
     53#Magnitude 9.1 event (Sumba) at Onslow')
     54#title('ANUGA modelled maximum stage versus Green\'s approximation \n \
     55#Magnitude 9.0 event (Sumba) at Onslow')
    4856legend(['ANUGA','Green\'s law'])
    4957axis([-5,85,min(stage_max)*0.9,max(stage_max)*1.1])
    50 savefig('stratification')
     58grid(True)
     59#savefig('stratification_onslow_gun')
     60savefig('stratification_exmouth_gun')
     61#savefig('stratification_dampier_gun')
    5162
    5263close('all')
Note: See TracChangeset for help on using the changeset viewer.