Changeset 5139 for anuga_work/production/tonga/project.py

Timestamp:

Mar 7, 2008, 4:37:36 PM (17 years ago)

Author:

herve

Message:

deleting data files

File:

• anuga_work/production/tonga/project.py (modified) (3 diffs)

anuga_work/production/tonga/project.py

@@ -2 +2 @@
 """Common filenames and locations for topographic data, meshes and outputs.
 """
-print 'hello'
+
 from os import sep, environ, getenv, getcwd
 from os.path import expanduser
 import sys
 from time import localtime, strftime, gmtime
-from anuga.utilities.polygon import read_polygon, plot_polygons, \
-                                    polygon_area, is_inside_polygon
+from anuga.utilities.polygon import read_polygon, plot_polygons, is_inside_polygon, number_mesh_triangles
+#from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
+from anuga.utilities.system_tools import get_user_name, get_host_name
+
+# file and system info
+#---------------------------------
+#codename = 'project.py'
+
+home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent diruser = get_user_name()
+user = get_user_name()
+host = get_host_name()
+
+# INUNDATIONHOME is the inundation directory, not the data directory.
+
+#time stuff
+time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
+gtime = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
+build_time = time+'_build'
+run_time = time+'_run'
+print 'gtime: ', gtime
+
+#Making assumptions about the location of scenario data
+state = 'sw_pacific'
+scenario_name = 'Tonga'
+scenario = 'fixed_wave'
+
+tide = 0
+
+alpha = 0.1
+friction=0.01
+starttime=10000
+midtime=21600
+finaltime=10000
+export_cellsize=50
+setup='trial'
+
+if setup =='trial':
+    print'trial'
+    res_factor=10
+    time_thinning=48
+    yieldstep=240
+if setup =='basic': 
+    print'basic'
+    res_factor=4
+    time_thinning=12
+    yieldstep=120
+if setup =='final': 
+    print'final'
+    res_factor=1
+    time_thinning=4
+    yieldstep=60
+
+dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(user)
+
+
+topo_gridfile ='tongatapu_10mgrid'
+
+anuga_dir = home+state+sep+scenario+sep+'anuga'+sep
+
+topographies_in_dir = home+sep+state+sep+scenario+sep+'elevation_final'+sep+'points'+sep
+topographies_dir = anuga_dir+'topographies'+sep
+#topographies_time_dir = topographies_dir+build_time+sep
+
+meshes_dir = anuga_dir+'meshes'+sep
+meshes_dir_name = meshes_dir + scenario_name
+
+polygons_dir = anuga_dir+'polygons'+sep
+tide_dir = anuga_dir+'tide_data'+sep
+
+boundaries_in_dir_name = boundaries_in_dir + scenario_name
+boundaries_dir = anuga_dir+'boundaries'+sep
+boundaries_dir_name = boundaries_dir + scenario_name
+
+#output locations
+output_dir = anuga_dir+'outputs'+sep
+output_build_time_dir = anuga_dir+'outputs'+sep+build_time+dir_comment+sep
+output_run_time_dir = anuga_dir+'outputs'+sep+run_time+dir_comment+sep
+output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post pr
 
 ###############################
 # Domain definitions
 ###############################
-
-# bounding polygon for study area
-bounding_polygon = read_polygon('extent.txt')
-
-print 'Area of bounding polygon', polygon_area(bounding_polygon)/1000000.0
+from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
+
+poly_all = read_polygon(polygons_dir+'extent.txt')
+res_poly_all = 100000*res_factor
+
 
 ###############################
@@ -39 +115 @@
 poly_island13= read_polygon('islands13.txt')
 
-#plot_polygons([bounding_polygon,poly_cairns,poly_island0,poly_island1,\
-#               poly_island2,poly_island3,poly_shallow],\
-#              'boundingpoly',verbose=False)
+
+islands_res = 5000
+cairns_res = 5000 
+interior_regions = [[project.poly_tongatapu, tongatapu_res],
+                    [project.poly_island1, islands_res],
+                    [project.poly_island2, islands_res],
+                    [project.poly_island3, islands_res],
+                    [project.poly_island4, islands_res],
+                    [project.poly_island5, islands_res],
+                    [project.poly_island6, islands_res],
+                    [project.poly_island7, islands_res],
+                    [project.poly_island8, islands_res],
+                    [project.poly_island9, islands_res],
+                    [project.poly_island10, islands_res],
+                    [project.poly_island11, islands_res],
+                    [project.poly_island12, islands_res],
+                    [project.poly_island13, islands_res]]
+
+boundary_tags={'ocean_east': [0],'ocean_north': [1],'ocean_west': [2],'land': [3,4,5,6,7],'ocean_southeast':[8,9]}
+
+trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+
+print 'min number triangles', trigs_min
 
 ###################################################################
@@ -57 +153 @@
 slide_depth = 207.
 
-gauge_filename = 'gauges.csv'
+
+
+
+
+
+
+
+
+
+
+# topo Filenames
+dem_name = 'tongatapu_10mgrid' 
+meshname = 'tongatapu.msh'
+
+# Create DEM from asc data
+convert_dem_from_ascii2netcdf(dem_name, use_cache=True, verbose=True)
+
+# Create pts file for onshore DEM
+dem2pts(dem_name, use_cache=True, verbose=True)
+
+
+
+
+
+
+# -*- coding: cp1252 -*-
+"""Common filenames and locations for topographic data, meshes and outputs.
+"""
+
+from os import sep, environ, getenv, getcwd
+from os.path import expanduser
+import sys
+from time import localtime, strftime, gmtime
+from anuga.utilities.polygon import read_polygon, plot_polygons, \
+                                    polygon_area, is_inside_polygon
+
+###############################
+# Domain definitions
+###############################
+
+# bounding polygon for study area
+bounding_polygon = read_polygon('Y:\data\sw_pacific\tonga\anuga\boundaries\extent.txt')
+
+print 'Area of bounding polygon', polygon_area(bounding_polygon)/1000000.0
+
+###############################
+# Interior region definitions
+###############################
+
+# interior polygons
+poly_tongatapu = read_polygon('tongatapu_finemesharea.txt')
+poly_island1 = read_polygon('island1.txt')
+poly_island2 = read_polygon('island2.txt')
+poly_island3 = read_polygon('islands3.txt')
+poly_island4 = read_polygon('islands4.txt')
+poly_island5= read_polygon('islands5.txt')
+poly_island6= read_polygon('islands6.txt')
+poly_island7= read_polygon('islands7.txt')
+poly_island8= read_polygon('islands8.txt')
+poly_island9= read_polygon('islands9.txt')
+poly_island10= read_polygon('islands10.txt')
+poly_island11= read_polygon('islands11.txt')
+poly_island12= read_polygon('islands12.txt')
+poly_island13= read_polygon('islands13.txt')
+
+#plot_polygons([bounding_polygon,poly_cairns,poly_island0,poly_island1,\
+#               poly_island2,poly_island3,poly_shallow],\
+#              'boundingpoly',verbose=False)
+
+###################################################################
+# Clipping regions for export to asc and regions for clipping data
+###################################################################
+
+# exporting asc grid
+eastingmin = 670500
+eastingmax = 712750
+northingmin = 7646000
+northingmax = 7677000
+
+
+slide_origin = [701290, 7665750] # move onto the continental shelf, depth = 500
+slide_depth = 207.
+
+#gauge_filename = 'gauges.csv'