Changeset 4021


Ignore:
Timestamp:
Nov 21, 2006, 3:02:03 PM (17 years ago)
Author:
sexton
Message:

(i) updates for Dampier report (ii) update cairns demo

Files:
1 added
1 deleted
12 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/demos/cairns/runcairns.py

    r4003 r4021  
    3535# Define scenario as either slump or fixed_wave.
    3636#-------------------------------------------------------------------------------
    37 scenario = 'slump' # OR 'fixed_wave'
     37#scenario = 'slump' # OR
     38scenario = 'fixed_wave'
    3839if access(scenario,F_OK) == 0:
    3940    mkdir (scenario)
     
    193194        domain.write_boundary_statistics(tags = 'ocean_east')     
    194195
    195         # save every 30 secs as wave starts inundating ashore
    196         for t in domain.evolve(yieldstep = 30, finaltime = 10000,
    197                                skip_initial_step = True):
    198             domain.write_time()
    199             domain.write_boundary_statistics(tags = 'ocean_east')
     196    # save every 30 secs as wave starts inundating ashore
     197    for t in domain.evolve(yieldstep = 30, finaltime = 10000,
     198                           skip_initial_step = True):
     199        domain.write_time()
     200        domain.write_boundary_statistics(tags = 'ocean_east')
    200201           
    201202print 'That took %.2f seconds' %(time.time()-t0)
  • anuga_core/source/anuga/shallow_water/smf.py

    r4009 r4021  
    9292    a3D = a2D / (1 + (15.5*sqrt(depth/(length*sint))))
    9393
     94    if verbose is True:
     95       
    9496    #a few temporary print statements
    9597    if verbose is True:
     
    111113        print '\t2D amplitude: ', a2D
    112114        print '\t3D amplitude: ', a3D
     115        print '\t HELLO', a0, ut, s0, t0, w, a2D, a3D
    113116
    114117    #keep an eye on some of the assumptions built into the maths
     
    117120        if verbose is True:
    118121            print 'WARNING: slope out of range (5 - 30 degrees) ', slope
     122            print 'hello'
    119123    if ((depth/length < 0.06) or (depth/length > 1.5)):
    120124        if verbose is True:
  • anuga_work/production/dampier_2006/export_results.py

    r4009 r4021  
    1919
    2020if which_var == 1:  # Absolute Momentum
    21     outname = name + '_momentum'
     21    outname = name + '_momentum_i0'
    2222    quantityname = '(xmomentum**2 + ymomentum**2)**0.5' 
    2323
    2424if which_var == 2:  # Depth
    25     outname = name + '_depth'
     25    outname = name + '_depth_test100'
    2626    quantityname = 'stage-elevation' 
    2727
    2828if which_var == 3:  # Speed
    29     outname = name + '_speed'
     29    outname = name + '_speed_i8'
    3030    #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
    3131    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
     
    4242        sww2dem(namei, basename_out = outnamei,
    4343                    quantity = quantityname,
    44                     cellsize = 20,     
     44                    #timestep = 0,
     45                    cellsize = 100,     
    4546                    easting_min = project.e_min_area,
    4647                    easting_max = project.e_max_area,
  • anuga_work/production/dampier_2006/make_report_cipma.py

    r4009 r4021  
    8888gauge_map1 = 'dampier_gauges_1.jpg'
    8989gauge_map2 = 'dampier_gauges_2.jpg'
    90 gauge_map3 = 'dampier_gauges_3.jpg'
    9190
    9291# Create sections and graphs for each designated production directory
     
    111110                                                  report = True,
    112111                                                  reportname = reportname,
    113                                                   plot_quantity = ['stage', 'speed'],
     112                                                  plot_quantity = ['stage', 'momentum'],
    114113                                                  surface = False,
    115114                                                  time_min = None,
     
    133132                                          report = True,
    134133                                          reportname = report_name,
    135                                           plot_quantity = ['stage', 'speed'],
     134                                          plot_quantity = ['stage', 'momentum'],
    136135                                          surface = False,
    137136                                          time_min = None,
     
    257256s  = """
    258257\caption{Point locations used for Dampier study; surrounding the facility.
    259 See Appendix \\ref{sec:timeseries} for stage and speed results.
     258See Appendix \\ref{sec:timeseries} for stage and momentum results.
    260259
    261260\label{fig:points2}
     
    263262"""
    264263fid.write(s)
    265 s = '\\begin{figure}[hbt] \centerline{\includegraphics[scale=0.5]{../report_figures/%s}}' %gauge_map1
    266 fid.write(s)
    267 s  = """
    268 \caption{Point location used for Dampier study; surrounding the middle section of the pipeline.
    269 See Appendix \\ref{sec:timeseries} for stage and speed results.}
    270 \label{fig:points1}
    271 \end{figure}
    272 """
    273 fid.write(s)
    274264s = '\\begin{figure}[hbt] \centerline{\includegraphics[scale=0.5]{../report_figures/%s}}' %gauge_map3
    275265fid.write(s)
     
    277267\caption{Point locations used for Dampier study; section
    278268of the pipeline towards the boundary of the region and other boundary gauges.
    279 See Appendix \\ref{sec:timeseries} for stage and speed results.} 
     269See Appendix \\ref{sec:timeseries} for stage and momentum results.} 
    280270\label{fig:points3}
    281271\end{figure}
  • anuga_work/production/dampier_2006/report/anuga_setup.tex

    r4009 r4021  
    1919
    2020To set up a model for the tsunami scenario, a study area is first
    21 determined. Here, we define the study area to cover around km of coastline
    22 and extending offshore to around the 70m depth contour. Onshore data covers the
     21determined. Here, we define the study area to
     22extend offshore to around the 70m depth contour with the onshore data covering the
    2323area of interest to CIPMA. The study area is around 6250 km$^2$.
    2424
     
    3030area will be the maximum cell area within the defined region and that each
    3131cell in the region does not necessarily have the same area.}.
    32 The cell areas should not be too small as this will cause unrealisticly long computational time,
    33 and not too great as this may inadequately capture important behaviour.
    34 %There are no gains in choosing the area to be less than the supporting data.
    35 Figure \ref{fig:dampier_area} shows the study area with regions of difference cell areas. The total number
    36 of cells is 177 000.
     32\Figure \ref{fig:dampier_area} shows the study area with regions of difference cell areas. The total number
     33of cells is 177 191.
    3734Lateral accuracy refers to the distance at which we are confident in stating a region is inundated.
    3835Figure \ref{fig:dampier_area} shows the maximum triangular cell area and lateral accuracy for each region.
     
    4643
    4744  \caption{Study area for Dampier scenario highlighting four regions of increased refinement.
    48 Region 1: Surrounds onshore gas infrastructure at Dampier with a cell area of 500 m$^2$ (lateral accuracy 30 m).
     45Region 1: Surrounds onshore gas infrastructure at Dampier with a cell area of 500 m$^2$ (lateral accuracy 30 m) - not shown here.
    4946Region 2: Surrounds the gas pipelines with a cell area of 2000 m$^2$ (lateral accuracy 100 m).
    5047Region 3: Surrounds the coastal regions with a cell area of 10000 m$^2$ (lateral accuracy  m).
  • anuga_work/production/dampier_2006/report/data.tex

    r3988 r4021  
    1313as the vertical datum. Mean Sea Level (MSL) is approximately equal to
    14140m AHD with the Highest Astronomical Tide (HAT)
    15 and Lowest Astronomical Tide (LAT) defined as 2.4m AHD
    16 and -m AHD respectively for Dampier \cite{antt:06}.
     15defined as 2.4m AHD for Dampier \cite{antt:06}.
    1716
    1817Data for this study have been sourced from a number of agencies. With
  • anuga_work/production/dampier_2006/report/execsum.tex

    r3957 r4021  
    33Risk Assessment Methods Project (RAMP). The purpose of this
    44study is to determine the potential tsunami inundation extent and
    5 particle velocity as the tsunami
     5particle momentum as the tsunami
    66impacts the major gas facilities and infrastructure in Dampier.
    77
     
    1010and results (Section \ref{sec:results})
    1111for a specific tsunami-genic event as it impacts the Dampier gas facilities and infrastructure.
    12 In particular, maximum inundation and velocity maps are shown
     12In particular, maximum inundation and depth and velocity time series are shown
    1313for the event occurring at highest astronomical tide.
    1414Impact assessments will be conducted by CIPMA and are not discussed here.
  • anuga_work/production/dampier_2006/report/metadata.tex

    r3955 r4021  
    1 \includepdf[pages={1-5}]{MetadataforATWSDampierScenario}
     1\includepdf[pages={1-6}]{MetadataforATWSDampierScenario}
Note: See TracChangeset for help on using the changeset viewer.