Changeset 4145


Ignore:
Timestamp:
Jan 8, 2007, 3:36:39 PM (17 years ago)
Author:
sexton
Message:

report updates and fixes to check_list

Files:
39 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/planning/Inundation Scenarios - 06 07 Plan.Gnt

    r4071 r4145  
    11<?xml version="1.0" encoding="utf-8"?>
    22<!--****WARNING: DO NOT MODIFY THIS FILE DIRECTLY******-->
    3 <Gantt Version="0.902" Checksum="+R8DXMfKYgs=">
     3<Gantt Version="0.902" Checksum="FkpD2Gk1zVI=">
    44  <Display Width="0" Height="0" Split="178" />
    5   <Miscellaneous Saved="2006-12-11 18:35:09Z" Created="Monday, 11 December 2006 7:35:09 " />
     5  <Miscellaneous Saved="2006-12-19 10:48:35Z" Created="Monday, 18 December 2006 11:48:35 " />
    66  <Globalization>
    77    <Culture LCID="3081" ISO="ENA" DN="English (Australia)" Cal="System.Globalization.GregorianCalendar" Parent="English">English (Australia)</Culture>
     
    1010  <FirstDay Date="2006-08-28" />
    1111  <LastDay Date="2007-06-30" />
    12   <Tasks Size="38">
     12  <Tasks Size="40">
    1313    <!--S=Start D=Duration P=Predecessor H=Height L=Lag VA=Valign R=Remarks F=Flag I=Importance-->
    1414    <!--B=BarType BC=BarColor BTx=BarText-->
     
    1717    <Task ID="3" Type="1" S="2006-09-20" SH="0" P="2" In="0" B="2" BC="-16776961" D="17" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Hobart scenario</Task>
    1818    <Task ID="4" Type="1" S="2006-09-29" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Karratha/Dampier data to NM&amp;I</Task>
    19     <Task ID="5" Type="1" S="2006-10-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of K/D data</Task>
     19    <Task ID="5" Type="1" S="2006-10-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of K/D data for CIPMA scenario</Task>
    2020    <Task ID="6" Type="1" S="2006-09-15" SH="12" P="0" In="0" B="2" BC="-16776961" D="50.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D scenario</Task>
    2121    <Task ID="7" Type="1" S="2006-10-06" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome offshore data to NM&amp;I</Task>
    22     <Task ID="8" Type="1" S="2006-10-26" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome data</Task>
    23     <Task ID="9" Type="1" S="2006-10-27" SH="12" P="8" In="0" B="2" BC="-16776961" D="15.0" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome setup</Task>
     22    <Task ID="8" Type="1" S="2006-12-15" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome data</Task>
     23    <Task ID="9" Type="1" S="2007-01-08" SH="0" P="8" In="0" B="2" BC="-16776961" D="3.0" H="35" L="15.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome setup</Task>
    2424    <Task ID="10" Type="1" S="2006-11-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Perth offshore data to NM&amp;I</Task>
    2525    <Task ID="11" Type="1" S="2006-12-07" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Perth data</Task>
    26     <Task ID="12" Type="1" S="2006-12-08" SH="0" P="11" In="0" B="2" BC="-16776961" D="3" H="35" L="0.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth setup</Task>
     26    <Task ID="12" Type="1" S="2006-12-08" SH="0" P="11" In="0" B="2" BC="-16776961" D="7" H="35" L="0.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth setup</Task>
    2727    <Task ID="13" Type="1" S="2006-11-24" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton offshore data to NM&amp;I</Task>
    28     <Task ID="14" Type="1" S="2006-12-14" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton data</Task>
    29     <Task ID="15" Type="1" S="2006-12-18" SH="0" P="14" In="0" B="2" BC="-16776961" D="2.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton setup</Task>
     28    <Task ID="14" Type="1" S="2006-12-18" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton data</Task>
     29    <Task ID="15" Type="1" S="2006-12-19" SH="12" P="14" In="0" B="2" BC="-16776961" D="3.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton setup</Task>
    3030    <Task ID="16" Type="1" S="2006-11-27" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth offshore data to NM&amp;I</Task>
    31     <Task ID="17" Type="1" S="2006-12-11" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth data</Task>
    32     <Task ID="18" Type="1" S="2006-12-12" SH="12" P="17" In="0" B="2" BC="-16776961" D="3.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth setup</Task>
    33     <Task ID="19" Type="1" S="2007-01-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay offshore data to NM&amp;I</Task>
    34     <Task ID="20" Type="1" S="2007-02-02" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay data</Task>
    35     <Task ID="21" Type="1" S="2007-01-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of BC for Shark Bay</Task>
    36     <Task ID="22" Type="1" S="2007-02-05" SH="0" P="20" In="0" B="2" BC="-16776961" D="15" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Shark Bay validation exercise</Task>
    37     <Task ID="23" Type="1" S="2007-01-17" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Decision of number and type of events</Task>
    38     <Task ID="24" Type="1" S="2007-01-26" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of events</Task>
    39     <Task ID="25" Type="1" S="2007-01-29" SH="12" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D run</Task>
    40     <Task ID="26" Type="1" S="2007-02-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome run</Task>
    41     <Task ID="27" Type="1" S="2007-02-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth run</Task>
    42     <Task ID="28" Type="1" S="2007-02-19" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton run</Task>
    43     <Task ID="29" Type="1" S="2007-02-26" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth run</Task>
    44     <Task ID="30" Type="1" S="2007-03-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Onslow run</Task>
    45     <Task ID="31" Type="1" S="2007-03-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Pt Hedland run</Task>
    46     <Task ID="32" Type="1" S="2006-12-05" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of offshore data for NSW slump scenarios</Task>
    47     <Task ID="33" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Wollongong slump scenario</Task>
    48     <Task ID="34" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Sydney slump scenario</Task>
    49     <Task ID="35" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Newcastle slump scenario</Task>
    50     <Task ID="36" Type="1" S="2006-12-15" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">ANUGA release</Task>
    51     <Task ID="37" Type="1" S="2007-01-19" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Beowulf operational</Task>
    52     <Task ID="38" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="2" BC="-16776961" D="11" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Christmas Break</Task>
     31    <Task ID="17" Type="1" S="2006-12-21" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth data</Task>
     32    <Task ID="18" Type="1" S="2007-01-11" SH="0" P="17" In="0" B="2" BC="-16776961" D="3" H="35" L="15" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth setup</Task>
     33    <Task ID="19" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Deliver of updated K/D data</Task>
     34    <Task ID="20" Type="1" S="2007-01-03" SH="12" P="19" In="0" B="2" BC="-16776961" D="2.5" H="35" L="8.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D resetup</Task>
     35    <Task ID="21" Type="1" S="2007-01-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay offshore data to NM&amp;I</Task>
     36    <Task ID="22" Type="1" S="2007-02-02" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay data</Task>
     37    <Task ID="23" Type="1" S="2007-01-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of BC for Shark Bay</Task>
     38    <Task ID="24" Type="1" S="2007-02-05" SH="0" P="22" In="0" B="2" BC="-16776961" D="15" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Shark Bay validation exercise</Task>
     39    <Task ID="25" Type="1" S="2007-01-17" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Decision of number and type of events</Task>
     40    <Task ID="26" Type="1" S="2007-01-26" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of events</Task>
     41    <Task ID="27" Type="1" S="2007-01-29" SH="12" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D run</Task>
     42    <Task ID="28" Type="1" S="2007-02-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome run</Task>
     43    <Task ID="29" Type="1" S="2007-02-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth run</Task>
     44    <Task ID="30" Type="1" S="2007-02-19" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton run</Task>
     45    <Task ID="31" Type="1" S="2007-02-26" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth run</Task>
     46    <Task ID="32" Type="1" S="2007-03-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Onslow run</Task>
     47    <Task ID="33" Type="1" S="2007-03-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Pt Hedland run</Task>
     48    <Task ID="34" Type="1" S="2006-12-05" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of offshore data for NSW slump scenarios</Task>
     49    <Task ID="35" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Wollongong slump scenario</Task>
     50    <Task ID="36" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Sydney slump scenario</Task>
     51    <Task ID="37" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Newcastle slump scenario</Task>
     52    <Task ID="38" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">ANUGA release</Task>
     53    <Task ID="39" Type="1" S="2007-01-19" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Beowulf operational</Task>
     54    <Task ID="40" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="2" BC="-16776961" D="11" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Christmas Break</Task>
    5355  </Tasks>
    5456  <TextStyles Size="8">
  • anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

    r4131 r4145  
    833833                   
    834834    import sys
    835     if not sys.version.startswith('2.4'):
     835    #if not sys.version.startswith('2.4'):
    836836        # Backwards compatibility
    837         from sets import Set as set
    838        
    839    
    840     if sys.platform == 'win32':
    841         for i,j in enumerate(quantity):
    842             quantity[i] = quantity[i].lower()
    843         p = list(set(quantity).difference(set(all_quantity)))
    844         if len(p) <> 0:
    845             msg = 'Quantities %s do not exist - please try again' %p
    846             raise Exception, msg
     837    #   from sets import Set as set
     838
     839    from sets import Set as set
     840           
     841    for i,j in enumerate(quantity):
     842        quantity[i] = quantity[i].lower()
     843    p = list(set(quantity).difference(set(all_quantity)))
     844    if len(p) <> 0:
     845        msg = 'Quantities %s do not exist - please try again' %p
     846        raise Exception, msg
    847847       
    848848    return
     
    928928    max_momentums = []
    929929    max_speeds = []
    930     c = 0
     930    max_depths = []
    931931    model_time_plot3d = zeros((n0,m), Float)
    932932    stages_plot3d = zeros((n0,m), Float)
     
    944944            max_stage = 0
    945945            max_momentum = 0
    946             max_speed = 0   
     946            max_speed = 0
     947            max_depth = 0
    947948            gaugeloc = locations[k]
    948949            thisfile = file_loc[j]+sep+'gauges_time_series'+'_'+gaugeloc+'.csv'
     
    983984                    if m > max_momentum: max_momentum = m
    984985                    if vel > max_speed: max_speed = vel
     986                    if z > 0 and depth > max_depth: max_depth = depth
    985987                   
    986988                   
     
    990992            min_stages.append(min_stage)
    991993            max_momentums.append(max_momentum)
    992             max_speeds.append(max_speed)     
     994            max_speeds.append(max_speed)
     995            max_depths.append(max_depth)
    993996            #### finished generating quantities for each swwfile #####
    994997       
     
    10261029        profilefig = 'solution_xprofile'
    10271030        savefig('profilefig')
    1028                
    1029     #stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1])
    1030     #stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0])   
    1031     #vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1])
    1032     #mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 
     1031
     1032    depth_axis = axis([time_min/60.0, time_max/60.0, 0, max(max_depths)*1.1])
     1033    stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1])
     1034    stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0])   
     1035    vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1])
     1036    mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 
    10331037   
    10341038    cstr = ['g', 'r', 'b', 'c', 'm', 'y', 'k']
     
    10371041    elev_output = []
    10381042    if len(label_id) > 1: graphname_report = []
     1043    pp = 1
     1044    div = 11.
    10391045    for k in gauge_index:
    10401046        g = gauges[k]
    10411047        count1 = 0
    10421048        if report == True and len(label_id) > 1:
    1043             s = '\\begin{figure}[hbt] \n \\centering \n \\begin{tabular}{cc} \n'
     1049            s = '\\begin{figure}[ht] \n \\centering \n \\begin{tabular}{cc} \n'
    10441050            fid.write(s)
    10451051        if len(label_id) > 1: graphname_report = []
     
    10631069                    plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j])
    10641070                    units = 'm'
     1071                    axis(depth_axis)
    10651072                if which_quantity == 'stage':
    1066                     if elevations[0:n[j]-1,k,j] < 0:
     1073                    if elevations[0,k,j] < 0:
    10671074                        plot(model_time[0:n[j]-1,k,j], stages[0:n[j]-1,k,j], '-', c = cstr[j])
     1075                        axis(stage_axis)
    10681076                    else:
    10691077                        plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j])
    1070                     #axis(stage_axis)
     1078                        axis(depth_axis)                 
    10711079                    units = 'm'
    10721080                if which_quantity == 'momentum':
    10731081                    plot(model_time[0:n[j]-1,k,j], momenta[0:n[j]-1,k,j], '-', c = cstr[j])
    1074                     #axis(mom_axis)
     1082                    axis(mom_axis)
    10751083                    units = 'm^2 / sec'
    10761084                if which_quantity == 'xmomentum':
    10771085                    plot(model_time[0:n[j]-1,k,j], xmom[0:n[j]-1,k,j], '-', c = cstr[j])
    1078                     #axis(mom_axis)
     1086                    axis(mom_axis)
    10791087                    units = 'm^2 / sec'
    10801088                if which_quantity == 'ymomentum':
    10811089                    plot(model_time[0:n[j]-1,k,j], ymom[0:n[j]-1,k,j], '-', c = cstr[j])
    1082                     #axis(mom_axis)
     1090                    axis(mom_axis)
    10831091                    units = 'm^2 / sec'
    10841092                if which_quantity == 'speed':
    10851093                    plot(model_time[0:n[j]-1,k,j], speed[0:n[j]-1,k,j], '-', c = cstr[j])
    1086                     #axis(vel_axis)
     1094                    axis(vel_axis)
    10871095                    units = 'm / sec'
    10881096                if which_quantity == 'bearing':
     
    11701178                s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label)
    11711179                fid.write(s)
    1172                 c += 1
    1173                 if c % 6 == 0: fid.write('\\clearpage \n')
     1180                cc += 1
     1181                if cc % 6 == 0: fid.write('\\clearpage \n')
    11741182                savefig(graphname_latex)               
    11751183               
     
    11861194                    else:
    11871195                        word_quantity += plot_quantity[i]
    1188                 print 'hello', elevations[0,k,j]
    11891196                where1 = 0
    11901197                count1 += 1
     
    12131220            s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label)
    12141221            fid.write(s)
    1215             c += 1
    1216             if c % 6 == 0: fid.write('\\clearpage \n')         
     1222            if float((k+1)/div - pp) == 0.:
     1223                fid.write('\\clearpage \n')
     1224                pp += 1
    12171225           
    12181226            #### finished generating figures ###
  • anuga_validation/solitary_waves/solitary_wave_runup.py

    r4025 r4145  
    5454
    5555# Structured mesh
    56 dx = 30           # Resolution: Length of subdivisions on x axis (length)
    57 dy = 30           # Resolution: Length of subdivisions on y axis (width)
     56dx = 5           # Resolution: Length of subdivisions on x axis (length)
     57dy = 5           # Resolution: Length of subdivisions on y axis (width)
    5858
    5959length = east-west
  • anuga_work/production/onslow_2006/make_report.py

    r4134 r4145  
    187187  \input{execsum}
    188188
     189\clearpage
     190
    189191  \\tableofcontents
    190192 
     
    210212"""
    211213fid.write(s)
    212 
    213 # Generate latex output for location points
    214 s = '\\begin{table} \\begin{center} \n'
    215 fid.write(s)
    216 s = '\caption{Defined point locations for %s study area.}' %scenario_name
    217 fid.write(s)
    218 s = """
    219 \label{table:locations}
    220 \\begin{tabular}{|l|l|l|l|}\hline
    221 \\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
    222 """
    223 fid.write(s)
    224 
    225 gauges, locations, elevation = get_gauges_from_file(project.gauge_filename)
    226 
    227 for name, gauges, elev in zip(locations, gauges, elevation):
    228     east = gauges[0]
    229     north = gauges[1]
    230     s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev)
    231     fid.write(s)
    232 
    233 s = '\\end{tabular} \n  \end{center} \n \end{table} \n \n'
    234 fid.write(s)
    235 
    236 s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map
    237 fid.write(s)
    238 
    239 s  = """
    240 \caption{Point locations used for Onslow study.} 
    241 \label{fig:points}
    242 \end{figure}
    243 """
    244 fid.write(s)
    245214   
    246215s = '\input{interpretation} \n'
     
    257226    fid.write(s)
    258227
     228# Generate latex output for location points
     229s = '\\begin{table} \\begin{center} \n'
     230fid.write(s)
     231s = '\caption{Defined point locations for %s study area.}' %scenario_name
     232fid.write(s)
     233s = """
     234\label{table:locations}
     235\\begin{tabular}{|l|l|l|l|}\hline
     236\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
     237"""
     238fid.write(s)
     239
     240gauges, locations, elevation = get_gauges_from_file(project.gauge_filename)
     241
     242for name, gauges, elev in zip(locations, gauges, elevation):
     243    east = gauges[0]
     244    north = gauges[1]
     245    s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev)
     246    fid.write(s)
     247
     248s = '\\end{tabular} \n  \end{center} \n \end{table} \n \n'
     249fid.write(s)
     250
     251#s = '\\begin{figure}[h] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map
     252s = '\\begin{figure}[h] \n \centerline{ \includegraphics[scale=0.7]{../report_figures/%s}}' %gauge_map
     253fid.write(s)
     254
     255s  = """
     256\caption{Point locations used for Onslow study.} 
     257\label{fig:points}
     258\end{figure}
     259
     260\clearpage
     261"""
     262fid.write(s)
     263
    259264# Closing
    260265
     
    265270"""
    266271fid.write(s)
    267 
    268 #for i, name in enumerate(production_dirs.keys()):
    269 
    270 #    s = '\input{%s} \n \clearpage \n \n' %damage_maps[production_dirs[name]]
    271 #    fid.write(s)
    272272
    273273s = """
     
    285285    \input{references}
    286286
     287\clearpage
     288
    287289    \\appendix
    288290
     
    302304     \label{sec:timeseries}
    303305     \input{timeseriesdiscussion}
     306     \clearpage
    304307"""
    305308fid.write(s)
     
    309312
    310313s="""
    311 
    312 \pagebreak
    313314
    314315   \section{Damage modelling inputs}
  • anuga_work/production/onslow_2006/report/anuga_setup.tex

    r4134 r4145  
    5555\begin{figure}[hbt]
    5656
    57   \centerline{ \includegraphics[scale=0.15]{../report_figures/onslow_resolution_zones.jpg}}
     57  \centerline{ \includegraphics[scale=0.125]{../report_figures/onslow_resolution_zones.jpg}}
    5858
    5959  \caption{Study area for Onslow scenario highlighting four regions of increased refinement.
     
    8787\caption{Parameters used in ANUGA for the Onslow scenario.}
    8888\begin{tabular}{|l|l|l|}\hline
    89 Mesh & & \hline
    90 & resolution in Region 1 & 500 m$^2$ \hline
    91 & resolution in Region 2 & 2500 m$^2$ \hline
    92 & resolution in Region 3 & 20000 m$^2$ \hline
    93 & remaining resolution & 100 000 m$^2$ \hline
    94 Model parameters & & \hline
    95 & friction & 0.017
    96 & minimum stored height & 0.1 m \hline
     89Mesh & & \\ \hline
     90& resolution in Region 1 & 500 m$^2$ \\ \hline
     91& resolution in Region 2 & 2500 m$^2$ \\ \hline
     92& resolution in Region 3 & 20000 m$^2$ \\ \hline
     93& remaining resolution & 100 000 m$^2$ \\ \hline
     94Model parameters & & \\ \hline
     95& friction & 0.017 \\ \hline
     96& minimum stored height & 0.1 m \\ \hline
    9797\end{tabular}
    9898\end{center}
  • anuga_work/production/onslow_2006/report/damage_inputs.tex

    r4134 r4145  
    1 \begin{table}[p]
     1\begin{table}[h]
    22\begin{center}
    33\caption{Framed residential building collapse probability. $h$ is the
     
    1616\end{table}
    1717
    18 \begin{table}
     18\begin{table}[h]
    1919\begin{center}
    2020\caption{Mortality and injury state probability}
     
    3434\end{table}
    3535
    36 \begin{table}
     36\begin{table}[h]
    3737\begin{center}
    3838\caption{Injury level classificationse. Floor height is assumed to be 30cm}
  • anuga_work/production/onslow_2006/report/introduction.tex

    r3407 r4145  
    2222detailing critical infrastructure as well as damage modelling estimates.
    2323
    24 This report is the first in a series of tsunami assessments
    25 of the North West Shelf. The scenario used for this study has
    26 an unknown return period, but is considered a plausible event (see
    27 Section \ref{sec:methodology}).
    28 Subsequent assessments will use refined hazard models with
    29 associate return periods. A suite of assessments will be
    30 made for Onslow and other localities, as advised by FESA.
     24This report details the impact assessments for a range of tsunami events.
     25These events are based on the probabilistic hazard assessment conducted
     26for the Western Australian coastline. A number of events are selected for
     27return periods of 500, 1000 and 2000 years, see Section \ref{sec:tsunamiscenario}.
     28%This report is the first in a series of tsunami assessments
     29%of the North West Shelf. The scenario used for this study has
     30%an unknown return period, but is considered a plausible event (see
     31%Section \ref{sec:methodology}).
     32%Subsequent assessments will use refined hazard models with
     33%associate return periods. A suite of assessments will be
     34%made for Onslow and other localities, as advised by FESA.
    3135
    3236Onslow has a population of around 800 and
     
    3842
    3943The modelling technique to simulate the
    40 impact ashore is discussed in Section \ref{sec:methodology} and data inputs
    41 discussed in Section \ref{sec:data}.
     44impact ashore is discussed in Section \ref{sec:methodology} with
     45event and data inputs
     46discussed in Sections \ref{sec:tsunamiscenario} and \ref{sec:data} respectively.
    4247The inundation results are presented and discussed in Section \ref{sec:results}
    4348and the impact modelling results outlined in Section \ref{sec:impact}.
  • anuga_work/production/onslow_2006/report/latexoutput.tex

    r4134 r4145  
    1 \begin{figure}[hbt]
     1\begin{figure}[ht]
     2 \centering
     3 \begin{tabular}{cc}
     4\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthstage.png}&
     5\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthspeed.png}\\
     6\end{tabular}
     7 \caption{Time series for stage and speed at Beadon Point Loading Berth location (elevation -8.70m)}
     8 \label{fig:gaugeBeadonPointLoadingBerth}
     9 \end{figure}
     10 
     11\begin{figure}[ht]
     12 \centering
     13 \begin{tabular}{cc}
     14\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalstage.png}&
     15\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalspeed.png}\\
     16\end{tabular}
     17 \caption{Time series for depth and speed at Hospital location (elevation 6.02m)}
     18 \label{fig:gaugeHospital}
     19 \end{figure}
     20 
     21\begin{figure}[ht]
    222 \centering
    323 \begin{tabular}{cc}
     
    929 \end{figure}
    1030 
     31\begin{figure}[ht]
     32 \centering
     33 \begin{tabular}{cc}
     34\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationstage.png}&
     35\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationspeed.png}\\
     36\end{tabular}
     37 \caption{Time series for depth and speed at Power Station location (elevation 5.81m)}
     38 \label{fig:gaugePowerStation}
     39 \end{figure}
     40 
     41\begin{figure}[ht]
     42 \centering
     43 \begin{tabular}{cc}
     44\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwaystage.png}&
     45\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwayspeed.png}\\
     46\end{tabular}
     47 \caption{Time series for depth and speed at Airport Runway location (elevation 5.49m)}
     48 \label{fig:gaugeAirportRunway}
     49 \end{figure}
     50 
     51\begin{figure}[ht]
     52 \centering
     53 \begin{tabular}{cc}
     54\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksstage.png}&
     55\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksspeed.png}\\
     56\end{tabular}
     57 \caption{Time series for depth and speed at Beadon Creek Docks location (elevation 1.65m)}
     58 \label{fig:gaugeBeadonCreekDocks}
     59 \end{figure}
     60 
     61\begin{figure}[ht]
     62 \centering
     63 \begin{tabular}{cc}
     64\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynestage.png}&
     65\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynespeed.png}\\
     66\end{tabular}
     67 \caption{Time series for stage and speed at West of Groyne location (elevation -2.11m)}
     68 \label{fig:gaugeWestofGroyne}
     69 \end{figure}
     70 
     71\begin{figure}[ht]
     72 \centering
     73 \begin{tabular}{cc}
     74\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthstage.png}&
     75\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthspeed.png}\\
     76\end{tabular}
     77 \caption{Time series for stage and speed at Beadon Creek mouth location (elevation -2.90m)}
     78 \label{fig:gaugeBeadonCreekmouth}
     79 \end{figure}
     80 
     81\begin{figure}[ht]
     82 \centering
     83 \begin{tabular}{cc}
     84\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockstage.png}&
     85\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockspeed.png}\\
     86\end{tabular}
     87 \caption{Time series for stage and speed at Beadon Creek south of dock location (elevation -1.81m)}
     88 \label{fig:gaugeBeadonCreeksouthofdock}
     89 \end{figure}
     90 
     91\begin{figure}[ht]
     92 \centering
     93 \begin{tabular}{cc}
     94\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallstage.png}&
     95\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallspeed.png}\\
     96\end{tabular}
     97 \caption{Time series for depth and speed at Centre dam wall location (elevation 2.09m)}
     98 \label{fig:gaugeCentredamwall}
     99 \end{figure}
     100 
     101\begin{figure}[ht]
     102 \centering
     103 \begin{tabular}{cc}
     104\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowstage.png}&
     105\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowspeed.png}\\
     106\end{tabular}
     107 \caption{Time series for depth and speed at Dam overflow location (elevation 0.88m)}
     108 \label{fig:gaugeDamoverflow}
     109 \end{figure}
     110 
     111\clearpage
     112\begin{figure}[ht]
     113 \centering
     114 \begin{tabular}{cc}
     115\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerstage.png}&
     116\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerspeed.png}\\
     117\end{tabular}
     118 \caption{Time series for depth and speed at Light Tower location (elevation 3.86m)}
     119 \label{fig:gaugeLightTower}
     120 \end{figure}
     121 
     122\begin{figure}[ht]
     123 \centering
     124 \begin{tabular}{cc}
     125\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayweststage.png}&
     126\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBaywestspeed.png}\\
     127\end{tabular}
     128 \caption{Time series for stage and speed at Beadon Bay west location (elevation -4.62m)}
     129 \label{fig:gaugeBeadonBaywest}
     130 \end{figure}
     131 
     132\begin{figure}[ht]
     133 \centering
     134 \begin{tabular}{cc}
     135\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeaststage.png}&
     136\includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeastspeed.png}\\
     137\end{tabular}
     138 \caption{Time series for stage and speed at Beadon Bay east location (elevation -3.56m)}
     139 \label{fig:gaugeBeadonBayeast}
     140 \end{figure}
     141 
  • anuga_work/production/onslow_2006/report/modelling_methodology.tex

    r4134 r4145  
    2929evenly spaced along the Sunda Arc subduction zone.
    3030The assessment used the Method of Splitting Tsunamis (MOST)
    31 \cite{VT:MOST} model. A detailed probabilistic hazard map has now been completed
    32 for the WA coastline \cite{prob:fesa} which is based on the paradigm used by the URS corporation's
    33 Probabilistic Tsunami Hazard Analysis \cite{somerville:urs}.
     31\cite{VT:MOST} model. This preliminary hazard map had no probability attached to the event which is
     32required to conduct a tsunami risk assessment. Using
     33the Probabilistic Tsunami Hazard Analysis \cite{somerville:urs} paradigm
     34used by the URS corporation, a detailed probabilistic hazard map has now been completed
     35for the WA coastline \cite{prob:fesa}.
    3436
    3537While MOST and URS are suitable for generating and propagating the tsunami wave from its source,
     
    5759Further details
    5860regarding the inundation modelling requirements for this study can be found in
    59 Appendix \ref{anugasetup}.
     61Appendix \ref{sec:anugasetup}.
    6062
    6163The risk of a given tsunami scenario can only be determined when the likelihood
     
    7779for a range of probabilities (or return periods). As Figure \ref{fig:probonslow}
    7880shows, for a given probability, a number of events are possible. The resulting
    79 impact to Onslow would then vary depending on the source of the event. 
     81impact to Onslow would then vary depending on the source of the event.  Further detail
     82on the tsunami scenarios are outlined in Section \ref{sec:tsunamiscenario}.
    8083
    8184% used for the 2005 report when looking at one event
  • anuga_work/production/onslow_2006/report/onslow_2006_report.tex

    r4134 r4145  
    5353  \input{execsum}
    5454
     55\clearpage
     56
    5557  \tableofcontents
    5658 
     
    7476     \label{sec:results}
    7577         
     78\input{interpretation}
     79\input{HAT_map}
     80 \clearpage
     81 
     82\input{MSL_map}
     83 \clearpage
     84 
     85\input{LAT_map}
     86 \clearpage
     87 
    7688\begin{table} \begin{center}
    7789\caption{Defined point locations for onslow study area.}
     
    97109 \end{table}
    98110 
    99 \begin{figure}[hbt]
    100  \centerline{ \includegraphics[width=\paperwidth]{../report_figures/onslow_dli_gauge.jpg}}
     111\begin{figure}[h]
     112 \centerline{ \includegraphics[scale=0.7]{../report_figures/onslow_dli_gauge.jpg}}
    101113\caption{Point locations used for Onslow study.} 
    102114\label{fig:points}
    103115\end{figure}
    104 \input{interpretation}
    105 \input{HAT_map}
    106  \clearpage
    107  
    108 \input{MSL_map}
    109  \clearpage
    110  
    111 \input{LAT_map}
    112  \clearpage
    113  
     116
     117\clearpage
    114118
    115119   \section{Impact modelling}
     
    130134    \input{references}
    131135
     136\clearpage
     137
    132138    \appendix
    133139
     
    147153     \label{sec:timeseries}
    148154     \input{timeseriesdiscussion}
     155     \clearpage
    149156\input{latexoutput}
    150157 \clearpage
    151158 
    152 
    153 
    154 \pagebreak
    155159
    156160   \section{Damage modelling inputs}
Note: See TracChangeset for help on using the changeset viewer.