Changeset 4927

Timestamp:

Jan 10, 2008, 6:04:35 PM (17 years ago)

Author:

sexton

Message:

scripts for stratification investigation

Location:

anuga_work/production

Files:

• onslow_2006/convert2latlon.py (modified) (1 diff)
• onslow_2006/plot_contour_data.py (modified) (3 diffs)
• pt_hedland_2006/convert2latlon.py (added)
• pt_hedland_2006/plot_contour_data.py (added)

anuga_work/production/onslow_2006/convert2latlon.py

@@ -30 +30 @@
 from os import sep
 
-timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
-#filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
+#timedir = '20070608_060316_run_final_1.5_onslow_nbartzis' # HAT 1 in Onslow 10000 yr
+timedir = '20070613_034754_run_final_1.5_dampier_nbartzis'
+#timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
+filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
 #filename = project.outputdir + timedir + sep + 'stage5_interp.csv'
-filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
+#filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
 #filename = project.outputdir+  timedir + sep + 'stage80_interp.csv'
-#filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
+filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
 #filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'
-filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
+#filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
 #filename_out = project.outputdir + timedir + sep + 'stage80_convert.csv'
 alter_file(filename,filename_out)

anuga_work/production/onslow_2006/plot_contour_data.py

@@ -3 +3 @@
 
 import project
-from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend
+from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend, grid
 from os import sep
  
 directory = project.outputdir
 timedir = '20070906_061418_run_final_1.5_jsexton_exmouth_original'
+#timedir = '20070608_060316_run_final_1.5_onslow_nbartzis'
+#timedir = '20070613_034754_run_final_1.5_dampier_nbartzis'
 
 def get_data(filename):
@@ -32 +34 @@
 
 from anuga.abstract_2d_finite_volumes.util import greens_law
+from Numeric import arange
 d1 = 80.
-d2 = range(d1,5,-1)
+d2 = arange(d1,0.15,-0.1)
 h1 = max_stage
+#h1 = 4.0
 green = []
 for d in d2:
@@ -46 +50 @@
 title('ANUGA modelled maximum stage versus Green\'s approximation \n \
 Magnitude 9.3 event (Java) at Onslow')
+#title('ANUGA modelled maximum stage versus Green\'s approximation \n \
+#Magnitude 9.1 event (Sumba) at Onslow')
+#title('ANUGA modelled maximum stage versus Green\'s approximation \n \
+#Magnitude 9.0 event (Sumba) at Onslow')
 legend(['ANUGA','Green\'s law'])
 axis([-5,85,min(stage_max)*0.9,max(stage_max)*1.1])
-savefig('stratification')
+grid(True)
+#savefig('stratification_onslow_gun')
+savefig('stratification_exmouth_gun')
+#savefig('stratification_dampier_gun')
 
 close('all')