Changeset 5896


Ignore:
Timestamp:
Nov 6, 2008, 12:08:50 PM (15 years ago)
Author:
kristy
Message:

Updated script to combine the output asc grids into one maximum asc grid for the corresponding which_area + which_var, due to the problem of an .sww file only 2GB big (hence more than one for a scenario).

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/export_results_all.py

    r5850 r5896  
    2121
    2222
    23 time_dir1 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
    24 time_dir2 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
    25 time_dir3 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
     23##time_dir1 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
     24##time_dir2 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
     25##time_dir3 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
    2626
    27 time_dirs = [time_dir1, time_dir2, time_dir3]
     27time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
     28time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
     29time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
     30time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
     31time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
     32time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
    2833
    29 #cellsize = 20
    30 cellsize = 250
     34
     35time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     36
     37cellsize = 20
     38#cellsize = 250
    3139
    3240timestep = None    # None means no timestep!
     
    5058
    5159# one or more key strings from var_equations above
    52 var = ['stage', 'elevation']
     60var = ['depth', 'speed']
    5361
    5462######
     
    5765for time_dir in time_dirs:
    5866    name1 = directory+time_dir+sep+project.scenario_name
    59 #    name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything
    60     names = [name1]
    61    
     67    name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything
     68    name3 = directory+time_dir+sep++project.scenario_name+'_time_79800_0' #need to get assistance on how to make this into anything
     69
     70    names = [name1, name2, name3]
     71
     72    asc_name = []
     73     
    6274    for name in names:
    6375        for which_area in area:
     
    8294                    break
    8395
    84                 outname = name + which_area + '_' + which_var
     96                outname = name + '_' + which_area + '_' + which_var
    8597                quantityname = var_equations[which_var]
    8698
     
    98110                            verbose = True,
    99111                            format = 'asc')
     112
     113                asc_name.append(name)
     114
     115    maxasc_outname = directory+time_dir + which_area+'_'+which_var+'_max.asc'
     116   
     117    MaxAsc(maxasc_outname, asc_name)
     118
     119
Note: See TracChangeset for help on using the changeset viewer.