Changeset 6459


Ignore:
Timestamp:
Mar 5, 2009, 9:45:32 AM (16 years ago)
Author:
kristy
Message:

Changed wave to Tb but put if statements in so that you can change back to Bf

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/dampier/run_model.py

    r6445 r6459  
    141141# Prepare time boundary
    142142TB.prepare_timeboundary(project.boundary_csv)
    143 f = file_function(project.boundary_csv[:-4] + '.tms',domain)
     143f = file_function(project.boundary_csv[:-4] + '.tms')
    144144
    145145Br = Reflective_boundary(domain)
    146146Bt = Transmissive_stage_zero_momentum_boundary(domain)
    147147Bd = Dirichlet_boundary([project.tide, 0, 0])
    148 Bf = Field_boundary(project.event_sts+'.sts',
    149                     domain, mean_stage=project.tide,
    150                     time_thinning=1,
    151                     default_boundary=Bd,
    152                     boundary_polygon=bounding_polygon_sts,                   
    153                     use_cache=True,
    154                     verbose=True)
    155 Tb = Time_boundary(domain,f)
    156148
    157 domain.set_boundary({'back': Br,
     149if project.wave == 'Bf':
     150    Bf = Field_boundary(project.event_sts+'.sts',
     151                        domain, mean_stage=project.tide,
     152                        time_thinning=1,
     153                        default_boundary=Bd,
     154                        boundary_polygon=bounding_polygon_sts,                   
     155                        use_cache=True,
     156                        verbose=True)
     157    domain.set_boundary({'back': Br,
    158158                     'side': Bd,
    159                      'ocean': Tb})
     159                     'ocean': Bf})
     160   
     161elif project.wave == 'Tb':
     162    Tb = Time_boundary(domain,f,default_boundary=Bd )
     163
     164    domain.set_boundary({'back': Br,
     165                         'side': Bd,
     166                         'ocean': Tb})
     167else:
     168    print 'No wave specified in project script (Bf or Tb)'
     169   
    160170
    161171#-------------------------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.