Changeset 7702


Ignore:
Timestamp:
Apr 29, 2010, 3:50:10 PM (14 years ago)
Author:
steve
Message:

Using small file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

    r7676 r7702  
    1616#------------------------------------------------------------------------------
    1717
    18 import unittest
    1918import os
    2019import sys
    2120import time
    2221import pypar
    23 
    2422import numpy as num
    2523
    26 
    27 
    28 
     24#------------------------
     25# ANUGA Modules
     26#------------------------
    2927from anuga.utilities.numerical_tools import ensure_numeric
    3028from anuga.utilities.util_ext        import double_precision
     
    4846#--------------------------------------------------------------------------
    4947
    50 #mesh_filename = "merimbula_10785_1.tsh"
    51 mesh_filename = "merimbula_43200.tsh"
    52 #mesh_filename = "test-100.tsh"
    53 yieldstep = 50
     48#mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.0
     49#mesh_filename = "merimbula_43200.tsh"   ; x0 = 756000.0 ; x1 = 756500.0
     50mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5
     51yieldstep = 5
    5452finaltime = 200
    55 quantity = 'stage'
    5653verbose = True
    5754
     
    6057#--------------------------------------------------------------------------
    6158class Set_Stage:
    62     """Set an initial condition with constant water height, for x<x0
     59    """Set an initial condition with constant water height, for x0<x<x1
    6360    """
    6461
     
    7673if myid == 0:
    7774    domain = create_domain_from_file(mesh_filename)
    78     domain.set_quantity('stage', Set_Stage(756000.0, 756500.0, 2.0))
     75    domain.set_quantity('stage', Set_Stage(x0, x1, 2.0))
    7976else:
    8077    domain = None
     
    8683if myid == 0 and verbose: print 'DISTRIBUTING DOMAIN'
    8784domain = distribute(domain)
     85
     86domain.smooth = False
    8887
    8988#------------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.