Ignore:
Timestamp:
Apr 17, 2009, 4:06:53 PM (15 years ago)
Author:
kristy
Message:

exported all events

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/lockhart_river/export_results_max.py

    r6746 r6835  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090306_131720_run_final_0_64448_1889_Tb_mhingee'
    25 time_dir2 = '20090306_165354_run_final_0_58331_1889_Tb_mhingee'
    26 time_dir3 = '20090306_165528_run_final_0_27347_1889_Tb_mhingee'
    27 ##time_dir4 = ''
    28 ##time_dir5 = ''
    29 ##time_dir6 = ''
    30 ##
    31 time_dirs = [time_dir1, time_dir2, time_dir3] #time_dir2]# , time_dir4, time_dir5, time_dir6]
    32  
     24time_dir1 = '20090408_142152_run_final_0_51316_1994_Tb__mhingee'
     25time_dir2 = '20090408_152731_run_final_0_51316_1994_Tb__mhingee'
     26time_dir3 = '20090408_200614_run_final_0_62786_1994_Tb__mhingee'
     27
     28time_dirs = [time_dir1, time_dir2, time_dir3]
    3329
    3430cellsize = 250
    3531##cellsize = 5
    3632
    37 #timestep = None    # None means no timestep!
    38 timestep = 0
     33timestep = None    # None means no timestep!
     34#timestep = 0
    3935
    4036######
     
    5652
    5753# one or more key strings from var_equations above
    58 var = ['elevation','stage']
     54var = ['stage']
    5955
    6056######
     
    8581        for time_dir in time_dirs:
    8682
    87             name1 = join(directory, time_dir, project.scenario_name)
    88             #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0'
    89             #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0'
     83            name = join(directory, time_dir, project.scenario_name)
    9084
    91             names = [name1] #, name2]#, name3]
    92      
    93            # asc_name = []   
    94                
    95             for name in names:
    96                          
    97                 outname = name + '_' + which_area + '_' + which_var
    98                 quantityname = var_equations[which_var]
     85            outname = name + '_' + which_area + '_' + which_var
     86            quantityname = var_equations[which_var]
    9987
    100                 print 'start sww2dem: time_dir=%s' % time_dir
    101                
    102                 sww2dem(name, basename_out = outname,
    103                             quantity = quantityname,
    104                             timestep = timestep,
    105                             cellsize = cellsize,     
    106                             easting_min = easting_min,
    107                             easting_max = easting_max,
    108                             northing_min = northing_min,
    109                             northing_max = northing_max,       
    110                             reduction = max,
    111                             verbose = True,
    112                             format = 'asc')
     88            print 'start sww2dem: time_dir=%s' % time_dir
     89           
     90            sww2dem(name, basename_out = outname,
     91                        quantity = quantityname,
     92                        timestep = timestep,
     93                        cellsize = cellsize,     
     94                        easting_min = easting_min,
     95                        easting_max = easting_max,
     96                        northing_min = northing_min,
     97                        northing_max = northing_max,       
     98                        reduction = max,
     99                        verbose = True,
     100                        format = 'asc')
    113101
    114                 asc_name.append(outname + '.asc')
    115    
    116             maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
    117 
    118             print 'max asc outname ', maxasc_outname
    119             print 'asc_name ', str(asc_name)
    120            
    121             MaxAsc(maxasc_outname, asc_name)
     102         
Note: See TracChangeset for help on using the changeset viewer.