Ignore:
Timestamp:
Oct 17, 2006, 12:21:30 PM (18 years ago)
Author:
sexton
Message:

introduce refined grids around specified sites for Hobart

File:
1 edited

Legend:

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

    r3805 r3809  
    124124#outputname = outputtimedir + basename  #Used by post processing
    125125
    126 
    127126###############################
    128127# Domain definitions
     
    146145polyAll = read_polygon(polygondir+'new_extent_2.csv')
    147146plot_polygons([polyAll, polyAll2],'boundingpoly',verbose=False)
    148 print 'Area of bounding polygon', polygon_area(polyAll)
     147#print 'Area of bounding polygon', polygon_area(polyAll)
    149148polyAll_refine = read_polygon(polygondir+'extent_refine.csv')
    150 print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
    151 print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
     149#print 'Area of refined bounding polygon', polygon_area(polyAll_refine)
     150#print 'change', polygon_area(polyAll_refine)/polygon_area(polyAll)
    152151
    153152polyAll_refine2 = read_polygon(polygondir+'export_refine_2.csv')
    154 print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
    155 print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
     153#print 'Area of refined bounding polygon 2', polygon_area(polyAll_refine2)
     154#print 'change 2', polygon_area(polyAll_refine2)/polygon_area(polyAll)
    156155
    157156###################################################################
     
    164163northingmin = 5245000
    165164northingmax = 5260000
    166 
     165# for original model setup
    167166# clipping 25m data set - Hobart
    168167eastingmin25 = 524208.387
     
    183182northingmax25_3 = 5258511.857
    184183
     184# export for refined grids at paleo sites
     185eastingmin_site13 = 532835.94
     186eastingmax_site13 = 533848.23
     187northingmin_site13 = 5237384.18
     188northingmax_site13 = 5238195.03
     189
     190eastingmin_kingston = 525905.49
     191eastingmax_kingston = 527773.42
     192northingmin_kingston = 5240256.71
     193northingmax_kingston = 5242148.58
     194
     195eastingmin_bruny = 523638.9
     196eastingmax_bruny = 533000.88
     197northingmin_bruny = 5204967.41
     198northingmax_bruny = 5212177.69
     199
    185200###############################
    186201# Interior region definitions
     
    194209poly_hobart5 = read_polygon(polygondir+'new_interior_region.csv')
    195210
    196 plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'boundingpoly2',verbose=False)
    197 print 'Area of interior regions', (polygon_area(poly_hobart1)+\
    198       polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+\
    199       polygon_area(poly_hobart4)+polygon_area(poly_hobart5))/1000000.
    200 
    201 print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.
    202 print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.
     211plot_polygons([polyAll, poly_hobart1,poly_hobart2,poly_hobart3,poly_hobart5],'orig_model',verbose=False)
     212
     213# to refine around paleo sites
     214poly_kingston = read_polygon(polygondir+'Kingston.csv')
     215poly_bruny = read_polygon(polygondir+'bruny_refine.csv')
     216poly_site13 = read_polygon(polygondir+'Site13.csv')
     217
     218plot_polygons([polyAll, poly_kingston,poly_bruny,poly_site13,poly_hobart5],'orig_model',verbose=False)
     219#print 'Area of interior regions', (polygon_area(poly_hobart1)+\
     220#      polygon_area(poly_hobart2)+polygon_area(poly_hobart3)+\
     221#      polygon_area(poly_hobart4)+polygon_area(poly_hobart5))/1000000.
     222
     223#print 'Area of original 25m data clipping', (polygon_area(hob_clip)+polygon_area(bruny_clip))/1000000.
     224#print 'Area of new interior region', polygon_area(poly_hobart5)/1000000.
Note: See TracChangeset for help on using the changeset viewer.