Ignore:
Timestamp:
Feb 3, 2009, 11:51:59 AM (10 years ago)
Author:
ole
Message:

Looped elevation data conversions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/standardised_version/run_busselton.py

    r6259 r6261  
    2525
    2626# Standard modules
    27 from os import sep
    2827import os
    29 from os.path import dirname, basename
    30 from os import mkdir, access, F_OK
    31 from shutil import copy
     28from os.path import dirname
    3229import time
    33 import sys
    3430
    3531# Related major packages
     
    4844import project  # Definition of file names and polygons
    4945
    50    
     46
    5147#-----------------------------------------------------------------------
    5248# Copy scripts to time stamped output directory and capture screen
     
    5652#copy script must be before screen_catcher
    5753#copy_code_files(project.output_run_time_dir, __file__,
    58 #         dirname(project.__file__)+sep+ project.__name__+'.py' )
     54#         dirname(project.__file__)+os.sep+ project.__name__+'.py' )
    5955#start_screen_catcher(project.output_run_time_dir, myid, numprocs)
    6056
     
    9591print domain.statistics()
    9692
    97 
    9893domain.set_name(project.scenario_name)
    9994domain.set_datadir(project.output_run_time_dir)
     
    111106                      geo_reference=domain.geo_reference)
    112107domain.set_quantity('stage', IC, use_cache=True, verbose=True)
    113 
    114108domain.set_quantity('friction', project.friction)
    115 
    116109domain.set_quantity('elevation',
    117110                    filename=project.combined_dir_name+'.pts',
     
    124117# Setup boundary conditions
    125118#-------------------------------------------------------------------------
    126 
    127119print 'Set boundary - available tags:', domain.get_boundary_tags()
    128 
    129 boundary_urs_out=project.boundaries_dir_event + sep + project.scenario_name
    130120
    131121Br = Reflective_boundary(domain)
    132122Bd = Dirichlet_boundary([project.tide,0,0])
    133 
    134 
    135 Bf = Field_boundary(boundary_urs_out+'.sts',
     123Bf = Field_boundary(project.urs_boundary_name+'.sts',
    136124                    domain, mean_stage=project.tide,
    137125                    time_thinning=1,
     
    140128                    use_cache=True,
    141129                    verbose=True)
    142 
    143130
    144131domain.set_boundary({'back': Br,
Note: See TracChangeset for help on using the changeset viewer.