Ignore:
Timestamp:
May 9, 2005, 12:33:09 PM (20 years ago)
Author:
chris
Message:
 
Location:
inundation/ga/storm_surge/analytical solutions
Files:
4 added
4 edited

Legend:

Unmodified
Added
Removed
  • inundation/ga/storm_surge/analytical solutions/Analytical solution_wave_runup.py

    r774 r1353  
    5555domain.store = True
    5656domain.format = 'sww'
    57 domain.filename = 'run-up_second_order'
     57domain.filename = 'run-up_second_order_again'
    5858
    5959print "Number of triangles = ", len(domain)
  • inundation/ga/storm_surge/analytical solutions/Analytical_solution_contracting_channel_import_mesh.py

    r774 r1353  
    66for use with the generic finite volume framework
    77
    8    Copyright 2004
    9    Christopher Zoppou, Stephen Roberts, Ole Nielsen, Duncan Gray
    10    Geoscience Australia
     8   Copyright 2005
     9   Christopher Zoppou, Stephen Roberts
     10   ANU
    1111   
    1212Specific methods pertaining to the 2D shallow water equation
     
    1818"""
    1919
    20 #-----------------------------
     20#---------------
    2121# Module imports
    2222import sys
     
    2929from pmesh2domain import pmesh_to_domain_instance
    3030
    31 
    32 #------------------------------
     31#-------
    3332# Domain
    3433filename = 'converging_channel_30846.tsh'
     
    3736print 'Number of triangles = ', len(domain)
    3837
    39 domain.default_order = 1
     38#----------------
     39# Order of scheme
     40domain.default_order = 2
    4041domain.smooth = True
    4142
     43#-------------------------------------
    4244# Provide file name for storing output
    4345domain.store = True     #Store for visualisation purposes
    4446domain.format = 'sww'   #Native netcdf visualisation format
    45 domain.filename = 'contracting_channel_first-order'
     47domain.filename = 'contracting_channel_second-order'
    4648
    47 #Decide which quantities are to be stored at each timestep
     49#----------------------------------------------------------
     50# Decide which quantities are to be stored at each timestep
    4851domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
    4952
    50 #Reduction operation for get_vertex_values             
     53#------------------------------------------
     54# Reduction operation for get_vertex_values             
    5155from util import mean
    5256domain.reduction = mean
    5357
    54 #------------------------------
    55 # Boundary Conditions
     58#------------------------
     59# Set boundary Conditions
    5660tags = {}
    5761tags['upstream'] = Dirichlet_boundary([0.2, 1.2, 0.0])
     
    6064domain.set_boundary(tags)
    6165
    62 #-----------------
    63 #Initial condition
     66#----------------------
     67# Set initial condition
    6468domain.set_quantity('elevation', 0.0)
    6569domain.set_quantity('stage', 0.2)
    6670   
    67 
    68 #-----------------
    69 #Evolution
     71#----------
     72# Evolution
    7073import time
    7174t0 = time.time()
  • inundation/ga/storm_surge/analytical solutions/MacDonald_77541.tsh

    r632 r1353  
    1221521221520 # <Hole #> <x> <y> ...Mesh Holes...
    1221531221530 # <Region #> <x> <y> <tag>...Mesh Regions...
     1221540 # <Region #> [Max Area] ...Mesh Regions...
  • inundation/ga/storm_surge/analytical solutions/Sydney.py

    r776 r1353  
    1717     Dirichlet_boundary, Transmissive_boundary
    1818from pmesh2domain import pmesh_to_domain_instance
    19 from util import Polygon_function
     19from util import Polygon_function,read_polygon
     20
    2021
    2122######################
     
    3334domain.store = True
    3435domain.format = 'sww'
    35 domain.filename = 'sydney_ubd'
     36domain.filename = 'sydney_ubd_again'
    3637
    3738#Reduction operation for get_vertex_values
Note: See TracChangeset for help on using the changeset viewer.