Ignore:
Timestamp:
Mar 10, 2008, 4:00:52 PM (16 years ago)
Author:
nick
Message:

update geraldton

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/geraldton/run_geraldton.py

    r5007 r5150  
    168168    print 'domain id', id(domain)
    169169    #print 'Reading Boundary file',project.boundaries_dir_namea + '.sww'
     170    print'set_boundary'
     171
     172    Br = Reflective_boundary(domain)
     173    Bd = Dirichlet_boundary([kwargs['tide'],0,0])
     174    Bo = Dirichlet_boundary([kwargs['tide']+5.0,0,0])
    170175
    171176    if project.source != 'other':
     
    173178                    domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'],
    174179                    use_cache=True, verbose=True)
    175                    
    176     kwargs['input_start_time']=domain.starttime
    177 
    178     print 'finished reading boundary file'
    179 
    180     Br = Reflective_boundary(domain)
    181     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
    182     Bo = Dirichlet_boundary([kwargs['tide']+5.0,0,0])
    183 
    184 
    185     print'set_boundary'
    186 
    187     domain.set_boundary({'back': Br,
    188                          'side': Bd,
    189                          'ocean': Bf})
     180        print 'finished reading boundary file'
     181        domain.set_boundary({'back': Br,
     182                             'side': Bd,
     183                             'ocean': Bf})
     184    else:
     185        print 'set ocean'               
     186        domain.set_boundary({'back': Br,
     187                             'side': Bd,
     188                             'ocean': Bd})
     189
     190#    kwargs['input_start_time']=domain.starttime
     191
     192
     193
    190194    print'finish set boundary'
    191195
     
    195199    t0 = time.time()
    196200
    197     for t in domain.evolve(yieldstep = 240, finaltime = kwargs['starttime']):
     201    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']):
    198202        domain.write_time()
    199203        domain.write_boundary_statistics(tags = 'ocean')     
    200204
    201         if allclose(t, 120):
     205        if allclose(t, 240):
    202206            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
    203207
    204         if allclose(t, 720):
     208        if allclose(t, 1440):
    205209            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
    206210
     
    274278   
    275279    kwargs['output_dir']=project.output_run_time_dir
    276     kwargs['bathy_file']=project.combined_dir_name
    277 #    kwargs['bathy_file']=project.combined_small_dir_name + '.pts'
     280#    kwargs['bathy_file']=project.combined_dir_name+'.txt'
     281    kwargs['bathy_file']=project.combined_dir_name_small + '.txt'
    278282    kwargs['boundary_file']=project.boundaries_in_dir_name + '.sww'
    279283    kwargs['file_name']=project.home+'detail.csv'
Note: See TracChangeset for help on using the changeset viewer.