Changeset 6516


Ignore:
Timestamp:
Mar 16, 2009, 10:27:09 AM (15 years ago)
Author:
nick2009
Message:

Change boundary condition, so not to double-dip on the tide.

Location:
anuga_work/production/geraldton
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/geraldton/NewScripts/export_results_max.py

    r6420 r6516  
    2929##time_dirs = [time_dir1, time_dir2]#, time_dir3]
    3030
    31 time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
     31time_dir1 = '20090225_171439_run_final_0_wa04_val_kvanputt'
    3232##time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
    3333##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
     
    3939 
    4040
    41 cellsize = 250
     41cellsize = 20
    4242##cellsize = 5
    4343
     
    6363
    6464# one or more key strings from var_equations above
    65 var = ['elevation']
     65var = ['elevation','stage']
    6666
    6767######
  • anuga_work/production/geraldton/NewScripts/project.py

    r6421 r6516  
    2626# Model specific parameters.
    2727# One or all can be changed each time the run_model script is executed
    28 tide = 0                # difference between MSL and HAT
     28tide = 0.0               # difference between MSL and HAT
    2929zone = 50               # specify zone of model
    3030event_number = 'wa04_val'   # the event number or the mux file name
     
    3232friction=0.01           # manning's friction coefficient
    3333starttime=0             # start time for simulation
    34 finaltime=1000          # final time for simulation
     34finaltime=80000          # final time for simulation
    3535
    3636setup = 'final'         # This can be one of three values
     
    9292# Used in get_timeseries.py. 
    9393# Format easting,northing,name,elevation (with header)
    94 gauges_filename = 'gauges.csv'
     94gauges_filename = 'nick_geraldton_gauges.csv'
    9595
    9696# BUILDINGS EXPOSURE - for identifying inundated houses
  • anuga_work/production/geraldton/NewScripts/run_model.py

    r6420 r6516  
    141141                    domain, mean_stage=project.tide,
    142142                    time_thinning=1,
    143                     default_boundary=Bd,
     143                    default_boundary=Dirichlet_boundary([0, 0, 0]),
    144144                    boundary_polygon=bounding_polygon_sts,                   
    145145                    use_cache=True,
  • anuga_work/production/geraldton/export_results_max.py

    r6065 r6516  
    2020
    2121directory = project.output_dir
    22 ##time_dir1 = '20081209_155330_run_final_0_27255_250m_none_lfountai'
     22time_dir1 = '20090225_171439_run_final_0_wa04_val_kvanputt'
    2323##time_dir2 = '20081209_161229_run_final_0_68693_250m_none_lfountai'
    2424##time_dirs = [time_dir1, time_dir2]
     
    2626##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
    2727##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
    28 time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
     28##time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
    2929##time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
    30 time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
     30##time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
    3131##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
    3232##
    33 time_dirs = [time_dir3, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
     33time_dirs = [time_dir1] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3434
    3535 
     
    3838cellsize = 250
    3939
    40 timestep = None    # None means no timestep!
    41 #timestep = 0
     40##timestep = None    # None means no timestep!
     41timestep = 0
    4242
    4343######
     
    5959
    6060# one or more key strings from var_equations above
    61 var = ['stage']
     61var = ['stage,elevation']
    6262
    6363######
     
    8989
    9090            name1 = directory+time_dir+sep+project.scenario_name
    91             name2 = directory+time_dir+sep+project.scenario_name+'_time_41820_0'
     91#            name2 = directory+time_dir+sep+project.scenario_name+'_time_41820_0'
    9292           
    93             names = [name1, name2]
     93            names = [name1]
    9494     
    9595            asc_name = []   
Note: See TracChangeset for help on using the changeset viewer.