Changeset 6733


Ignore:
Timestamp:
Apr 7, 2009, 8:28:20 AM (15 years ago)
Author:
Leharne
Message:

updates to Gold Coast tsunami scenario

Location:
anuga_work/production/gold_coast_2009
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/gold_coast_2009/export_results_max.py

    r6657 r6733  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
     24#time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
    2525time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
    2626time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
     
    2828time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
    2929time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
    30 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
     30time_dirs = [time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
    3131
    3232##cellsize = 250
    33 cellsize = 5
     33cellsize = 20
    3434
    3535#timestep = None    # None means no timestep!
     
    5454
    5555# one or more key strings from var_equations above
    56 var = ['depth', 'elevation']
     56var = ['depth']
    5757
    5858######
     
    8484
    8585            name1 = join(directory, time_dir, project.scenario_name)
    86             #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
    87             #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
     86            name2 = join(directory, time_dir, project.scenario_name)+'_time_21960_0'
     87            name3 = join(directory, time_dir, project.scenario_name)+'_time_43920_0'
     88            #name4 = join(directory, time_dir, project.scenario_name)+'_time_65880_0'
    8889
    89             names = [name1] #, name2]#, name3]
     90            names = [name1, name2, name3]
    9091     
    91            # asc_name = []   
     92            asc_name = []   
    9293               
    9394            for name in names:
     
    110111                            format = 'asc')
    111112
    112 ##                asc_name.append(outname + '.asc')
    113 ##   
    114 ##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    115 ##
    116 ##            print 'max asc outname ', maxasc_outname
    117 ##            print 'asc_name ', str(asc_name)
    118 ##           
    119 ##            MaxAsc(maxasc_outname, asc_name)
     113                asc_name.append(outname + '.asc')
     114   
     115            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     116
     117            print 'max asc outname ', maxasc_outname
     118            print 'asc_name ', str(asc_name)
     119           
     120            MaxAsc(maxasc_outname, asc_name)
  • anuga_work/production/gold_coast_2009/get_timeseries.py

    r6657 r6733  
    1313
    1414from os import sep, rename, listdir, system
     15from os.path import join
    1516from anuga.abstract_2d_finite_volumes.util import sww2csv_gauges,csv2timeseries_graphs
    1617import project
    1718
    18 directory = project.output_dir
     19directory = project.output_folder
    1920
    20 time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
     21#time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
    2122time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
    2223time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
     
    2425time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
    2526time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
    26 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     27time_dirs = [time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    2728
    2829for time_dir in time_dirs:
    29     name = directory+time_dir+sep+project.scenario_name
    30     gauge = project.gauges_folder+sep+project.gauges_filename
     30    name = join(directory,time_dir,project.scenario_name)
     31    gauge = project.gauges
    3132    sww2csv_gauges(name+'.sww',gauge,
    3233                   quantities = ['stage','speed','depth','elevation'],
    3334                   verbose=True)
    3435
    35    
    36 ##    name1 = directory+time_dir+sep+project.scenario_name
    37 ####    name2 = directory+time_dir+sep+project.scenario_name+'_time_43380_0'
    38 ####    name3 = directory+time_dir+sep+project.scenario_name+'_time_78360_0'
    39 ##
    40 ##    names = [name1]#, name2, name3]
    41 ##
    42 ##    for i, name in enumerate(names):
    43 ##        gauge = project.gauges_dir_name
    44 ####        out_name = 'gauge_%d_' % i
    45 ##        print name
    46 ##        sww2csv_gauges(name+'.sww',gauge,#out_name = out_name,
    47 ##                        quantities = ['stage','speed','depth','elevation'],
    48 ##                        verbose=True)
    49 ##
Note: See TracChangeset for help on using the changeset viewer.