Changeset 5396

Timestamp:

Jun 10, 2008, 4:05:44 PM (15 years ago)

Author:

Leharne

Message:

Updating output directory in project and run scripts to automatically include run time, user and mesh resolution. Still in progress.

Location:

anuga_work/development/convergence_okushiri_2008

Files:

• project_truescale.py (modified) (1 diff)
• run_okushiri_truescale.py (modified) (2 diffs)

anuga_work/development/convergence_okushiri_2008/project_truescale.py

@@ -3 +3 @@
 
 """
+from os import sep, environ, getenv, getcwd, umask
+from os.path import expanduser, basename, join
+from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon, number_mesh_triangles
+import sys
+from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
+from time import localtime, strftime, gmtime
+from anuga.utilities.system_tools import get_user_name, get_host_name
+
+codename = 'project_grad.py' # FIXME can be obtained automatically as __file__
+
+home = join(getenv('INUNDATIONHOME'), 'data', 'anuga_validation',
+            'convergence_okushiri_2008')# Location of Data   
+user = get_user_name()
+host = get_host_name()
+#needed when running using mpirun, mpirun doesn't inherit umask from .bashrc
+umask(002)
+
+#time stuff
+time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
+run_time = time+'_run'
+
+#Set anuga directories
+anuga_dir = join(home, 'anuga')
+
+dir_comment='_'+setup+'_'+str(tide)+'_'+\
+             str(user)+'_'+boundary_event+'_'+which_data
+
+mesh_dir = join(anuga_dir, 'meshes')
+mesh_name = join(mesh_dir, 'okushiri_truescale')
+
+polygons_dir = join(anuga_dir, 'polygons') # Created with ArcGIS (csv files)
+tide_dir = join(anuga_dir, 'tide_data')
+
+#output locations
+output_dir = join(anuga_dir, 'outputs')+sep
+output_run_time_dir = output_dir +run_time+dir_comment+sep
+output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
+
+#gauges
+gauges_dir = join(anuga_dir,'gauges')
+gauge_name = 'gauge_location_onslow.csv'
+gauges_dir_name = gauges_dir + gauge_name
 
 # Given boundary wave

anuga_work/development/convergence_okushiri_2008/run_okushiri_truescale.py

@@ -40 +40 @@
 print domain.statistics()
 
+z = domain.get_vertex_coordinates()
 
+# Write vertex coordinates to file
+filename='okushiri_mesh_vertex_coordinates_truescale.txt'
+fid=open(filename,'w')
+fid.write('x (m), y (m)\n')
+for i in range(len(z)):
+    pt=z[i]
+    x=pt[0]
+    y=pt[1]
+    fid.write('%.6f, %.6f\n' %(x, y))
+        
 #-------------------------
 # Initial Conditions
@@ -57 +68 @@
 #-------------------------
 domain.set_name(project_truescale.output_filename)  # Name of output sww file
+domain.set_datadir(project_truescale.output_dir)
 domain.set_default_order(2)               # Apply second order scheme 
 domain.set_all_limiters(0.9)              # Max second order scheme (old lim)