Ignore:
Timestamp:
Aug 20, 2009, 9:20:31 AM (16 years ago)
Author:
kristy
Message:
 
File:
1 edited

Legend:

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

    r7355 r7385  
    2121
    2222directory = project.output_folder
    23 time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt'
     23time_dir1 = '20090504_124158_run_final_0_58068_small_kvanputt'
    2424time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt'
    2525time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
    2626time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
    2727
    28 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4
     28time_dirs = [time_dir1
    2929
    3030# sww filename extensions ie. hobart_time_17640_0.sww, input into list 17640
    3131# make sure numbers are in sequential order
    32 times = [16980, 33960, 50940]
     32times = [22680, 45360, 68040]
    3333           
    3434#Modify the cellsize value to set the size of the raster you require
    3535#Take into account mesh size when aplying this paramater
    36 cellsize = 20 #250
     36cellsize = 250 #250
    3737
    3838#Now set the timestep at which you want the raster generated.
     
    120120            print 'which var=%s' % which_var
    121121            #print 'asc_name ', str(asc_name)
    122             if times == []:     
    123                 print 'Only one sww file, no Max file required'
     122            if len(names) >1:
     123                maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
     124
     125                print 'max asc outname ', maxasc_outname
     126                print 'asc_name ', str(asc_name)
     127               
     128                MaxAsc(maxasc_outname, asc_name)
    124129            else:
    125                 print 'Multiple sww files'
    126                 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    127                 print 'maxasc_outname = ', maxasc_outname
    128                 MaxAsc(maxasc_outname, asc_name)
     130                print 'Only one sww file, no Max file required'
     131
Note: See TracChangeset for help on using the changeset viewer.