Changeset 8452


Ignore:
Timestamp:
Jun 21, 2012, 11:30:04 AM (12 years ago)
Author:
martins
Message:
 
Location:
trunk/anuga_work/development/Busselton
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_work/development/Busselton/Analysis/export_results_max.py

    r8449 r8452  
    2424rootDir = r"/short/w85/Shane/Output"
    2525scenario_name = "Busselton"
    26 event = "B4"
    27 #directory = join(rootDir,event)
     26event = "B2Full"
    2827
    29 time_dir1 = "final_"+event+"_sfm547_np1"
    30 time_dir2 = "final_"+event+"_sfm547_np1_From21600"
    31 
    32 time_dirs = [time_dir1,time_dir2] 
     28time_dirs = [#"final_"+event+"_sfm547_np1",\
     29#             "final_"+event+"_sfm547_np1_From28800",\
     30#             "final_"+event+"_sfm547_np1_From57600",\
     31#             "final_"+event+"_sfm547_np1_From86400",\
     32#             "final_"+event+"_sfm547_np1_From108000",\
     33             "final_"+event+"_sfm547_np1_From115200"]
    3334
    3435#Modify the cellsize value to set the size of the raster you require
     
    6667######
    6768
    68 if myid == 0:
    69     time_dir = time_dirs[0]
    70     var = vars[0]
     69proc = 0
     70Needed = True
     71
     72for dirID in range(len(time_dirs)):   
    7173   
    72 elif myid == 1:
    73     time_dir = time_dirs[1]
    74     var = vars[1]
     74    if myid == proc:
     75        time_dir = time_dirs[dirID]
     76        var = vars[0]
     77       
     78    elif myid == proc + 1:
     79        time_dir = time_dirs[dirID]
     80        var = vars[1]
     81       
     82    proc+=2 
     83
     84if myid >= len(time_dirs)*2:
     85    print "I am not needed: Proc ", myid
     86    Needed = False
     87
     88if Needed:
     89       
     90    if var not in var_equations:
     91        print 'Unrecognized variable name: %s' % which_var
     92        sys.exit()
    7593   
    76 elif myid == 2:
    77     time_dir = time_dirs[0]
    78     var = vars[1]
     94    for which_area in area:
     95        if which_area == 'All':
     96            easting_min = None
     97            easting_max = None
     98            northing_min = None
     99            northing_max = None
     100        else:
     101            try:
     102                easting_min = eval('project.xmin%s' % which_area)
     103                easting_max = eval('project.xmax%s' % which_area)
     104                northing_min = eval('project.ymin%s' % which_area)
     105                northing_max = eval('project.ymax%s' % which_area)
     106            except AttributeError:
     107                print 'Unrecognized area name: %s' % which_area
     108                break
    79109   
    80 elif myid == 3:
    81     time_dir = time_dirs[1]
    82     var = vars[0]
     110        names = []
     111        Name = glob.glob(join(rootDir,event,time_dir,scenario_name) +'*.sww')
     112        names.append(Name)
    83113   
    84 if var not in var_equations:
    85     print 'Unrecognized variable name: %s' % which_var
    86     sys.exit()
    87 
    88 for which_area in area:
    89     if which_area == 'All':
    90         easting_min = None
    91         easting_max = None
    92         northing_min = None
    93         northing_max = None
    94     else:
    95         try:
    96             easting_min = eval('project.xmin%s' % which_area)
    97             easting_max = eval('project.xmax%s' % which_area)
    98             northing_min = eval('project.ymin%s' % which_area)
    99             northing_max = eval('project.ymax%s' % which_area)
    100         except AttributeError:
    101             print 'Unrecognized area name: %s' % which_area
    102             break
    103 
    104     names = []
    105     Name = glob.glob(join(rootDir,event,time_dir,scenario_name) +'*.sww')
    106     names.append(Name)
    107 
    108     for names in names:
    109        
    110         asc_name = []   
    111         print 'start sww2dem: time_dir = %s' % time_dir               
    112         for name in names:
    113 
    114             test = name
    115             name = test[:-4]
    116             outname = name + '_' + which_area + '_' + var + str(cellsize) + "m"
    117             quantityname = var_equations[var]
    118 
    119             print 'start sww2dem: name = %s' % os.path.basename(name)
    120             print outname
    121             anuga.sww2dem(name+'.sww', outname+'.asc',
    122                         quantity = quantityname,
    123                         cellsize = cellsize,     
    124                         easting_min = easting_min,
    125                         easting_max = easting_max,
    126                         northing_min = northing_min,
    127                         northing_max = northing_max,       
    128                         reduction = timestep,
    129                         verbose = True)
    130 
    131             asc_name.append(outname + '.asc')
     114        for names in names:
    132115           
     116            asc_name = []   
     117            print 'start sww2dem: time_dir = %s' % time_dir               
     118            for name in names:
     119   
     120                test = name
     121                name = test[:-4]
     122                outname = name + '_' + which_area + '_' + var + str(cellsize) + "m"
     123                quantityname = var_equations[var]
     124   
     125                print 'start sww2dem: name = %s' % os.path.basename(name)
     126                print outname
     127                anuga.sww2dem(name+'.sww', outname+'.asc',
     128                            quantity = quantityname,
     129                            cellsize = cellsize,     
     130                            easting_min = easting_min,
     131                            easting_max = easting_max,
     132                            northing_min = northing_min,
     133                            northing_max = northing_max,       
     134                            reduction = timestep,
     135                            verbose = True)
     136   
     137                asc_name.append(outname + '.asc')
     138               
    133139barrier()
    134140
  • trunk/anuga_work/development/Busselton/Analysis/master.py

    r8449 r8452  
    2727wall[32] = "07:00:00"
    2828
    29 np = 4
     29np = 16
    3030
    3131ScriptName = "/home/547/sfm547/Busselton/Analysis/export_results_max.py"
  • trunk/anuga_work/development/Busselton/B1Full/anuga.log

    r8449 r8452  
    1 2012-06-16 18:13:52,079 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
    2 2012-06-16 18:13:52,210 INFO                       polygon:959 |--------------------------------------------------------------------------------
    3 2012-06-16 18:13:52,210 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
    4 2012-06-16 18:13:52,211 INFO                       polygon:962 |--------------------------------------------------------------------------------
    5 2012-06-16 18:13:52,211 INFO                       polygon:976 |Interior 100                      56.48              564762
    6 2012-06-16 18:13:52,212 INFO                       polygon:984 |Bounding 10000                    352.99             35299
    7 2012-06-16 18:13:52,212 INFO                       polygon:989 |Estimated total number of triangles: 857230
    8 2012-06-16 18:13:52,212 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
     12012-06-18 16:19:21,783 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
     22012-06-18 16:19:21,816 INFO                       polygon:959 |--------------------------------------------------------------------------------
     32012-06-18 16:19:21,816 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
     42012-06-18 16:19:21,816 INFO                       polygon:962 |--------------------------------------------------------------------------------
     52012-06-18 16:19:21,817 INFO                       polygon:976 |Interior 100                      56.48              564762
     62012-06-18 16:19:21,817 INFO                       polygon:984 |Bounding 10000                    352.99             35299
     72012-06-18 16:19:21,818 INFO                       polygon:989 |Estimated total number of triangles: 857230
     82012-06-18 16:19:21,818 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
  • trunk/anuga_work/development/Busselton/B2Full/anuga.log

    r8449 r8452  
    1 2012-06-15 15:13:55,701 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
    2 2012-06-15 15:13:55,798 INFO                       polygon:959 |--------------------------------------------------------------------------------
    3 2012-06-15 15:13:55,798 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
    4 2012-06-15 15:13:55,798 INFO                       polygon:962 |--------------------------------------------------------------------------------
    5 2012-06-15 15:13:55,802 INFO                       polygon:976 |Interior 100                      56.48              564762
    6 2012-06-15 15:13:55,802 INFO                       polygon:984 |Bounding 10000                    352.99             35299
    7 2012-06-15 15:13:55,802 INFO                       polygon:989 |Estimated total number of triangles: 857230
    8 2012-06-15 15:13:55,802 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
     12012-06-18 16:31:17,239 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
     22012-06-18 16:31:17,265 INFO                       polygon:959 |--------------------------------------------------------------------------------
     32012-06-18 16:31:17,265 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
     42012-06-18 16:31:17,265 INFO                       polygon:962 |--------------------------------------------------------------------------------
     52012-06-18 16:31:17,266 INFO                       polygon:976 |Interior 100                      56.48              564762
     62012-06-18 16:31:17,266 INFO                       polygon:984 |Bounding 10000                    352.99             35299
     72012-06-18 16:31:17,266 INFO                       polygon:989 |Estimated total number of triangles: 857230
     82012-06-18 16:31:17,266 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
  • trunk/anuga_work/development/Busselton/B3Full/anuga.log

    r8449 r8452  
    1 2012-06-13 09:57:27,185 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
    2 2012-06-13 09:57:27,296 INFO                       polygon:959 |--------------------------------------------------------------------------------
    3 2012-06-13 09:57:27,297 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
    4 2012-06-13 09:57:27,297 INFO                       polygon:962 |--------------------------------------------------------------------------------
    5 2012-06-13 09:57:27,298 INFO                       polygon:976 |Interior 100                      56.48              564762
    6 2012-06-13 09:57:27,298 INFO                       polygon:984 |Bounding 10000                    352.99             35299
    7 2012-06-13 09:57:27,298 INFO                       polygon:989 |Estimated total number of triangles: 857230
    8 2012-06-13 09:57:27,299 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
     12012-06-18 16:31:17,284 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
     22012-06-18 16:31:17,306 INFO                       polygon:959 |--------------------------------------------------------------------------------
     32012-06-18 16:31:17,306 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
     42012-06-18 16:31:17,306 INFO                       polygon:962 |--------------------------------------------------------------------------------
     52012-06-18 16:31:17,307 INFO                       polygon:976 |Interior 100                      56.48              564762
     62012-06-18 16:31:17,307 INFO                       polygon:984 |Bounding 10000                    352.99             35299
     72012-06-18 16:31:17,307 INFO                       polygon:989 |Estimated total number of triangles: 857230
     82012-06-18 16:31:17,307 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
  • trunk/anuga_work/development/Busselton/B4Full/anuga.log

    r8449 r8452  
    1 2012-06-16 18:14:42,834 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
    2 2012-06-16 18:14:42,868 INFO                       polygon:959 |--------------------------------------------------------------------------------
    3 2012-06-16 18:14:42,868 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
    4 2012-06-16 18:14:42,868 INFO                       polygon:962 |--------------------------------------------------------------------------------
    5 2012-06-16 18:14:42,869 INFO                       polygon:976 |Interior 100                      56.48              564762
    6 2012-06-16 18:14:42,870 INFO                       polygon:984 |Bounding 10000                    352.99             35299
    7 2012-06-16 18:14:42,870 INFO                       polygon:989 |Estimated total number of triangles: 857230
    8 2012-06-16 18:14:42,870 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
     12012-06-20 22:53:16,886 INFO           anuga.utilities.log:0   |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL
     22012-06-20 22:53:16,952 INFO                       polygon:959 |--------------------------------------------------------------------------------
     32012-06-20 22:53:16,953 INFO                       polygon:960 |Polygon  Max triangle area (m^2)  Total area (km^2)  Estimated #triangles
     42012-06-20 22:53:16,953 INFO                       polygon:962 |--------------------------------------------------------------------------------
     52012-06-20 22:53:16,953 INFO                       polygon:976 |Interior 100                      56.48              564762
     62012-06-20 22:53:16,953 INFO                       polygon:984 |Bounding 10000                    352.99             35299
     72012-06-20 22:53:16,954 INFO                       polygon:989 |Estimated total number of triangles: 857230
     82012-06-20 22:53:16,954 INFO                       polygon:990 |Note: This is generally about 20% less than the final amount
Note: See TracChangeset for help on using the changeset viewer.