Ignore:
Timestamp:
Nov 30, 2006, 7:03:43 PM (18 years ago)
Author:
sexton
Message:

updates to pt hedland script (new data provided) and sydney slide scenario

Location:
anuga_work/production/pt_hedland_2006
Files:
2 edited

Legend:

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

    r3788 r4045  
    77#from anuga.utilities.polygon import read_polygon
    88import sys
    9 from anuga.coordinate_transforms.redfearn import convert_points_from_latlon_to_utm
    109from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
    1110from time import localtime, strftime
     
    2524# offshore data from AHO fairsheets
    2625offshore_name2 = 'pt_hedland_offshore_points_fairsheet_orig'
     26
     27# reissued offshore data from PMD and NM&I Nov 2006
     28offshore_name0 = 'hydrographic_charts'
     29offshore_name1 = 'xy100003817'
     30offshore_name2 = 'xy100003819'
     31offshore_name3 = 'xy100003822'
     32offshore_name4 = 'xy100003824'
     33offshore_name5 = 'xy100003825'
     34offshore_name6 = 'xy100003826'
     35offshore_name7 = 'xy100003827'
     36offshore_name8 = 'xy100003828'
     37offshore_name9 = 'xy100003829'
     38offshore_name10 = 'xy100003830'
     39offshore_name11 = 'xy100003831'
     40offshore_name12 = 'xy100003832'
     41offshore_name13 = 'xy100003837'
     42offshore_name14 = 'xy100003838'
     43offshore_name15 = 'xy100003840'
     44offshore_name16 = 'xy100003841'
     45offshore_name17 = 'xy100003843'
     46offshore_name18 = 'xy100003845'
     47offshore_name19 = 'xy100003846'
     48offshore_name20 = 'xy100003847'
     49offshore_name21 = 'xy100003945'
     50offshore_name22 = 'xy100003946'
     51offshore_name23 = 'xy100003947'
     52offshore_name24 = 'xy100003948'
     53offshore_name25 = 'xy100004071'
     54offshore_name26 = 'xy100004072'
     55offshore_name27 = 'xy100004074'
     56offshore_name28 = 'xy100004078'
     57offshore_name29 = 'xy100004079'
     58offshore_name30 = 'xy100004080'
     59offshore_name31 = 'xy100004081'
     60offshore_name32 = 'xy100004082'
     61offshore_name33 = 'xy100004083'
     62offshore_name34 = 'xy100004084'
     63offshore_name35 = 'xy100004085'
     64offshore_name36 = 'xy100004086'
     65offshore_name37 = 'xy100004087'
     66offshore_name38 = 'xy100004088'
     67offshore_name39 = 'xy100004089'
     68offshore_name40 = 'xy100004090'
     69offshore_name41 = 'xy100011753'
     70offshore_name42 = 'xy100014967'
     71offshore_name43 = 'xy100014968'
     72offshore_name44 = 'xy100014969'
     73offshore_name45 = 'xy100021030'
     74offshore_name46 = 'xy100021031'
     75offshore_name47 = 'xy100021032'
     76offshore_name48 = 'xy100021033'
     77offshore_name49 = 'xy100021035'
     78offshore_name50 = 'xy100022767'
    2779
    2880# coastline developed from aerial photography and 1.5m DLI contour
     
    86138#onshore_dem_name = datadir + onshore_name_dted
    87139onshore_dem_name = datadir + onshore_name_dli
     140offshore_dem_name0 = datadir + offshore_name0
    88141offshore_dem_name1 = datadir + offshore_name1
    89142offshore_dem_name2 = datadir + offshore_name2
     143offshore_dem_name3 = datadir + offshore_name3
     144offshore_dem_name4 = datadir + offshore_name4
     145offshore_dem_name5 = datadir + offshore_name5
     146offshore_dem_name6 = datadir + offshore_name6
     147offshore_dem_name7 = datadir + offshore_name7
     148offshore_dem_name8 = datadir + offshore_name8
     149offshore_dem_name9 = datadir + offshore_name9
     150offshore_dem_name10 = datadir + offshore_name10
     151offshore_dem_name11 = datadir + offshore_name11
     152offshore_dem_name12 = datadir + offshore_name12
     153offshore_dem_name13 = datadir + offshore_name13
     154offshore_dem_name14 = datadir + offshore_name14
     155offshore_dem_name15 = datadir + offshore_name15
     156offshore_dem_name16 = datadir + offshore_name16
     157offshore_dem_name17 = datadir + offshore_name17
     158offshore_dem_name18 = datadir + offshore_name18
     159offshore_dem_name19 = datadir + offshore_name19
     160offshore_dem_name20 = datadir + offshore_name20
     161offshore_dem_name21 = datadir + offshore_name21
     162offshore_dem_name22 = datadir + offshore_name22
     163offshore_dem_name23 = datadir + offshore_name23
     164offshore_dem_name24 = datadir + offshore_name24
     165offshore_dem_name25 = datadir + offshore_name25
     166offshore_dem_name26 = datadir + offshore_name26
     167offshore_dem_name27 = datadir + offshore_name27
     168offshore_dem_name28 = datadir + offshore_name28
     169offshore_dem_name29 = datadir + offshore_name29
     170offshore_dem_name30 = datadir + offshore_name30
     171offshore_dem_name31 = datadir + offshore_name31
     172offshore_dem_name32 = datadir + offshore_name32
     173offshore_dem_name33 = datadir + offshore_name33
     174offshore_dem_name34 = datadir + offshore_name34
     175offshore_dem_name35 = datadir + offshore_name35
     176offshore_dem_name36 = datadir + offshore_name36
     177offshore_dem_name37 = datadir + offshore_name37
     178offshore_dem_name38 = datadir + offshore_name38
     179offshore_dem_name39 = datadir + offshore_name39
     180offshore_dem_name40 = datadir + offshore_name40
     181offshore_dem_name41 = datadir + offshore_name41
     182offshore_dem_name42 = datadir + offshore_name42
     183offshore_dem_name43 = datadir + offshore_name43
     184offshore_dem_name44 = datadir + offshore_name44
     185offshore_dem_name45 = datadir + offshore_name45
     186offshore_dem_name46 = datadir + offshore_name46
     187offshore_dem_name47 = datadir + offshore_name47
     188offshore_dem_name48 = datadir + offshore_name48
     189offshore_dem_name49 = datadir + offshore_name49
     190offshore_dem_name50 = datadir + offshore_name50
     191
    90192coast_dem_name = datadir + coast_name
    91193combined_dem_name = datadir + 'pt_hedland_combined_elevation'
  • anuga_work/production/pt_hedland_2006/run_pt_hedland.py

    r3788 r4045  
    66
    77The scenario is defined by a triangular mesh created from project.polygon,
    8 the elevation data and a simulated submarine landslide.
     8the elevation data and a tsunami generated by a subduction zone earthquake.
    99
    1010Ole Nielsen and Duncan Gray, GA - 2005 and Nick Bartzis, GA - 2006
     
    2828from anuga.geospatial_data.geospatial_data import *
    2929import sys
    30 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
     30from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
    3131
    3232# Application specific imports
     
    3838#-------------------------------------------------------------------------------
    3939
    40 # creates copy of code in output dir if dir doesn't exist
    41 if access(project.outputtimedir,F_OK) == 0 :
    42     mkdir (project.outputtimedir)
    43 copy (dirname(project.__file__) +sep+ project.__name__+'.py', project.outputtimedir + project.__name__+'.py')
    44 copy (__file__, project.outputtimedir + basename(__file__))
    45 print 'project.outputtimedir',project.outputtimedir
    46 
    47 # normal screen output is stored in
    48 screen_output_name = project.outputtimedir + "screen_output.txt"
    49 screen_error_name = project.outputtimedir + "screen_error.txt"
    50 
    51 # used to catch screen output to file
    52 sys.stdout = Screen_Catcher(screen_output_name)
    53 sys.stderr = Screen_Catcher(screen_error_name)
     40copy_code_files(project.outputtimedir,__file__,dirname(project.__file__)+sep+ project.__name__+'.py' )
     41myid = 0
     42numprocs = 1
     43start_screen_catcher(project.outputtimedir, myid, numprocs)
     44
    5445print 'USER:    ', project.user
    5546
     
    7970        verbose=True)
    8071
    81 print 'create G1'
    82 G1 = Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')
    83 print 'create G2'
    84 G2 = Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')
    85 print 'create G3'
    86 G3 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
    87 print 'create G4'
    88 G4 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
    89 print 'add G1+G2+G3+G4'
    90 G = G1 + G2 + G3 + G4
     72print 'create offshore'
     73G11= Geospatial_data(file_name = project.offshore_dem_name0 + '.xya')+\
     74     Geospatial_data(file_name = project.offshore_dem_name1 + '.xya')+\
     75     Geospatial_data(file_name = project.offshore_dem_name2 + '.xya')+\
     76     Geospatial_data(file_name = project.offshore_dem_name3 + '.xya')+\
     77     Geospatial_data(file_name = project.offshore_dem_name4 + '.xya')+\
     78     Geospatial_data(file_name = project.offshore_dem_name5 + '.xya')+\
     79     Geospatial_data(file_name = project.offshore_dem_name6 + '.xya')+\
     80     Geospatial_data(file_name = project.offshore_dem_name7 + '.xya')+\
     81     Geospatial_data(file_name = project.offshore_dem_name8 + '.xya')+\
     82     Geospatial_data(file_name = project.offshore_dem_name9 + '.xya')+\
     83     Geospatial_data(file_name = project.offshore_dem_name10 + '.xya')
     84G12= Geospatial_data(file_name = project.offshore_dem_name11 + '.xya')+\
     85     Geospatial_data(file_name = project.offshore_dem_name12 + '.xya')+\
     86     Geospatial_data(file_name = project.offshore_dem_name13 + '.xya')+\
     87     Geospatial_data(file_name = project.offshore_dem_name14 + '.xya')+\
     88     Geospatial_data(file_name = project.offshore_dem_name15 + '.xya')+\
     89     Geospatial_data(file_name = project.offshore_dem_name16 + '.xya')+\
     90     Geospatial_data(file_name = project.offshore_dem_name17 + '.xya')+\
     91     Geospatial_data(file_name = project.offshore_dem_name18 + '.xya')+\
     92     Geospatial_data(file_name = project.offshore_dem_name19 + '.xya')+\
     93     Geospatial_data(file_name = project.offshore_dem_name20 + '.xya')
     94G13= Geospatial_data(file_name = project.offshore_dem_name21 + '.xya')+\
     95     Geospatial_data(file_name = project.offshore_dem_name22 + '.xya')+\
     96     Geospatial_data(file_name = project.offshore_dem_name23 + '.xya')+\
     97     Geospatial_data(file_name = project.offshore_dem_name24 + '.xya')+\
     98     Geospatial_data(file_name = project.offshore_dem_name25 + '.xya')+\
     99     Geospatial_data(file_name = project.offshore_dem_name26 + '.xya')+\
     100     Geospatial_data(file_name = project.offshore_dem_name27 + '.xya')+\
     101     Geospatial_data(file_name = project.offshore_dem_name28 + '.xya')+\
     102     Geospatial_data(file_name = project.offshore_dem_name29 + '.xya')
     103G14= Geospatial_data(file_name = project.offshore_dem_name30 + '.xya')+\
     104     Geospatial_data(file_name = project.offshore_dem_name31 + '.xya')+\
     105     Geospatial_data(file_name = project.offshore_dem_name32 + '.xya')+\
     106     Geospatial_data(file_name = project.offshore_dem_name33 + '.xya')+\
     107     Geospatial_data(file_name = project.offshore_dem_name34 + '.xya')+\
     108     Geospatial_data(file_name = project.offshore_dem_name35 + '.xya')+\
     109     Geospatial_data(file_name = project.offshore_dem_name36 + '.xya')+\
     110     Geospatial_data(file_name = project.offshore_dem_name37 + '.xya')+\
     111     Geospatial_data(file_name = project.offshore_dem_name38 + '.xya')+\
     112     Geospatial_data(file_name = project.offshore_dem_name39 + '.xya')
     113G15= Geospatial_data(file_name = project.offshore_dem_name40 + '.xya')+\
     114     Geospatial_data(file_name = project.offshore_dem_name41 + '.xya')+\
     115     Geospatial_data(file_name = project.offshore_dem_name42 + '.xya')+\
     116     Geospatial_data(file_name = project.offshore_dem_name43 + '.xya')+\
     117     Geospatial_data(file_name = project.offshore_dem_name44 + '.xya')+\
     118     Geospatial_data(file_name = project.offshore_dem_name45 + '.xya')+\
     119     Geospatial_data(file_name = project.offshore_dem_name46 + '.xya')+\
     120     Geospatial_data(file_name = project.offshore_dem_name47 + '.xya')+\
     121     Geospatial_data(file_name = project.offshore_dem_name48 + '.xya')+\
     122     Geospatial_data(file_name = project.offshore_dem_name49 + '.xya')+\
     123     Geospatial_data(file_name = project.offshore_dem_name50 + '.xya')
     124   
     125print 'create onshore'
     126G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
     127print 'create coast'
     128G3 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
     129print 'add'
     130G = G11 + G12 + G13 + G14 + G15 + G2 + G3
    91131print 'export G'
    92132G.export_points_file(project.combined_dem_name + '.pts')
Note: See TracChangeset for help on using the changeset viewer.