Changeset 3779


Ignore:
Timestamp:
Oct 13, 2006, 3:44:12 PM (17 years ago)
Author:
sexton
Message:

refining Hobart example even further ...

Location:
anuga_work/production/hobart_2006
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/hobart_2006/project.py

    r3769 r3779  
    88from time import localtime, strftime, gmtime
    99from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area
    10 from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
    11 from anuga.coordinate_transforms.redfearn import convert_from_latlon_to_utm
     10from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_from_latlon_to_utm
    1211
    1312if sys.platform == 'win32':
     
    5857coast_name = 'coastline_points'
    5958
    60 #boundary_basename = 'puysegur' # Mw 8.7
    61 boundary_basename = 'puysegur_clip' # Mw 8.5
     59boundary_basename = 'puysegur' # Mw 8.7
     60#boundary_basename = 'puysegur_clip' # Mw 8.5
    6261
    6362#swollen/ all data output
     
    7877polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
    7978
    80 gauge_filename = gaugedir + 'hobart_gauges_final.csv'
     79#gauge_filename = gaugedir + 'hobart_gauges_final.csv'
     80gauge_filename = gaugedir + 'check_boundary.csv'
    8181
    8282codedir = getcwd()+sep                           
     
    147147plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
    148148print 'Area of bounding polygon', polygon_area(polyAll)
     149polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
     150print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
     151print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
     152
     153polyAll_refine2 = read_polygon(polygondir+'export_refine_2.csv')
     154print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
     155print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
    149156
    150157###################################################################
     
    163170northingmin25 = 5229154.555
    164171northingmax25 = 5258511.857
    165 
     172hob_clip = [[eastingmin25,northingmin25],[eastingmin25,northingmax25],[eastingmax25,northingmax25],[eastingmax25,northingmin25]]
    166173# clipping 25m data set - Bruny
    167174eastingmin25_2 = 523672.502
     
    169176northingmin25_2 = 5203120.035
    170177northingmax25_2 = 5212052.309
    171 
     178bruny_clip = [[eastingmin25_2,northingmin25_2],[eastingmin25_2,northingmax25_2],[eastingmax25_2,northingmax25_2],[eastingmax25_2,northingmin25_2]]
    172179# overall 25m clipping region - covering Hobart and Bruny
    173180eastingmin25_3 = 523672.502
     
    187194
    188195plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart4],'boundingpoly2',verbose=False)
     196print 'Area of interior regions', (polygon_area(poly_hobart1)+\
     197      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+polygon_area(poly_hobart4))/1000000.
     198
     199print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.
     200
  • anuga_work/production/hobart_2006/report/boundary_report.tex

    r3767 r3779  
    4242
    4343\begin{document}
     44\title{Comparison between ANUGA and MOST}
    4445\maketitle
    4546\begin{table} \begin{center}
Note: See TracChangeset for help on using the changeset viewer.