Changeset 4521


Ignore:
Timestamp:
May 30, 2007, 2:35:04 PM (17 years ago)
Author:
sexton
Message:

update to script for looking at "cannot find triangle for point"

Location:
anuga_work/production/onslow_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/onslow_2006/export_results.py

    r4504 r4521  
    55from os import sep
    66
    7 #time_dir = '20070507_215717_run_final_1.5_nbartzis' # HAT 1 in Dampier 10000 yr
     7time_dir = '20070507_215717_run_final_1.5_nbartzis' # HAT 1 in Dampier 10000 yr
    88#time_dir = '20070517_061654_run_final_0.0_nbartzis' # MSL 1 in Dampier 10000 yr
    99#time_dir = '' # HAT 1 in Broome 10000 yr
     
    1515#time_dir = '' # HAT 1 in Exmouth 10000 yr
    1616#time_dir = '' # MSL 1 in Exmouth 10000 yr
    17 time_dir = '20060704_063005' # 2006 simulation HAT Mw9
    18 #cellsize = 30
    19 cellsize = 250
     17#time_dir = '20060704_063005' # 2006 simulation HAT Mw9
     18cellsize = 25
     19#cellsize = 150
    2020timestep = None
    2121directory = project.outputdir
    2222#name = directory + time_dir + sep + project.scenario_name
    23 #name = directory+time_dir+sep+project.scenario_name
    24 name = directory+time_dir+sep+'source'
     23name = directory+time_dir+sep+project.scenario_name
     24#name = directory+time_dir+sep+'source'
    2525
    26 is_parallel = False
    27 #is_parallel = True
     26#is_parallel = False
     27is_parallel = True
    2828if is_parallel == True: nodes = 4
    2929print 'output dir:', name
    3030
    31 #var = [2,3,4] # depth and speed
    32 var = [3] # depth and speed
     31var = [2,3,4] # depth and speed
     32#var = [3] # depth and speed
    3333
    3434for which_var in var:
     
    4646
    4747    if which_var == 3:  # Speed
    48         outname = name + '_speed'
     48        outname = name + '_speed_150'
    4949        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
    5050
  • anuga_work/production/onslow_2006/test_export_results.py

    r4518 r4521  
    3939            timestep = timestep,
    4040            cellsize = cellsize,     
    41             easting_min = project.e_min_area,
    42             easting_max = project.e_max_area,
    43             northing_min = project.n_min_area,
    44             northing_max = project.n_max_area,       
     41            #easting_min = project.e_min_area,
     42            #easting_max = project.e_max_area,
     43            #northing_min = project.n_min_area,
     44            #northing_max = project.n_max_area,       
    4545            reduction = max,
    4646            verbose = True,
Note: See TracChangeset for help on using the changeset viewer.