Changeset 6498


Ignore:
Timestamp:
Mar 12, 2009, 10:57:08 AM (15 years ago)
Author:
jgriffin
Message:

changed event number; added Arc_asc2raster_GDA94z56.py

Location:
anuga_work/production/australia_ph2/sydney
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/sydney/Arc_asc2raster_GDA94z56.py

    r6288 r6498  
    3131##time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
    3232##time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
    33 time_dir3 = '20090130_160955_run_final_0_7875_250m_vlarge_jgriffin'
     33time_dir3 = '20090211_135851_run_final_0_7875_jgriffin'
    3434#time_dir4 = '20090123_154945_run_final_0_7875_250m_small_jgriffin'
    3535##time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
     
    6262
    6363    generate_filename = []
    64     input_ascii = glob.glob(folder + '*stage_max.asc')
     64    input_ascii = glob.glob(folder + '*max.asc')
    6565    print time_dir
    6666
  • anuga_work/production/australia_ph2/sydney/comparisons/export_results_max.py

    r6298 r6498  
    2020directory = project_250m.output_dir
    2121
    22 time_dir1 = '20090206_152643_run_final_0_250m_large_jgriffin'
     22time_dir1 = '20090210_162935_run_final_0_250m_vlarge_jgriffin'
    2323#time_dir2 = '20090123_155045_run_final_0_7875_250m_large_jgriffin'
    2424#time_dir3 = '20090123_154945_run_final_0_7875_250m_small_jgriffin'
     
    8080
    8181            name1 = directory+time_dir+sep+project_250m.scenario_name
    82  
    83             names = [name1]
     82            name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_16680_0'
     83            name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_33360_0'
     84            name4 = directory+time_dir+sep+project_250m.scenario_name+'_time_50040_0'
     85            name5 = directory+time_dir+sep+project_250m.scenario_name+'_time_66720_0'
     86
     87            names = [name1, name2 , name3, name4, name5]
    8488     
    8589            asc_name = []   
  • anuga_work/production/australia_ph2/sydney/comparisons/project_250m.py

    r6295 r6498  
    4040scenario_name = 'sydney'
    4141scenario = 'sydney'
    42 area = 'vlarge' #'small_b10'
     42area = 'large' #'small_b10'
    4343
    4444# Model specific parameters. One or all can be changed each time the
     
    104104# Used in run_model.py
    105105# Format for points easting,northing (no header)                   
    106 interior_regions_data = [['coast_3km_buffer_' + area + '.csv', 10000]]
     106interior_regions_data = [['coast_3km_buffer_' + area + '.csv', 50000]]
    107107
    108108#------------------------------------------------------------------------------
  • anuga_work/production/australia_ph2/sydney/comparisons/run_sydney_250m.py

    r6295 r6498  
    221221
    222222    domain.set_boundary({'back': Br,
    223                          'side1': Bt,
    224                          'side2':Bt,
     223                         'side1': Bd,
     224                         'side2':Bd,
    225225                         'ocean':Rb})
    226226
  • anuga_work/production/australia_ph2/sydney/export_results_max.py

    r6298 r6498  
    22Generates ascii grids of nominated areas -
    33Input: sww file from run_pt_hedland.py
    4        boundaries for grids from project_250m.py
     4       boundaries for grids from project.py
    55Outputs: ascii grids of specified variables
    66Stored in the 'outputs_dir' folder for respective .sww file
     
    1111"""
    1212
    13 import project_250m, os
     13import project, os
    1414import sys
    1515from anuga.lib.maxasc.maxasc import MaxAsc
     
    1717from anuga.shallow_water.data_manager import sww2dem
    1818from os import sep
     19from os.path import join
    1920
    20 directory = project_250m.output_dir
     21directory = project.output_folder
     22print directory
    2123
    22 time_dir1 = '20090206_150250_run_final_0_250m_large_jgriffin'
    23 #time_dir2 = '20090123_155045_run_final_0_7875_250m_large_jgriffin'
     24time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
     25time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
    2426#time_dir3 = '20090123_154945_run_final_0_7875_250m_small_jgriffin'
    25 time_dirs = [time_dir1]#,time_dir2,time_dir3]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     27time_dirs = [time_dir1,time_dir2]#,time_dir2,time_dir3]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    2628
    2729#cellsize = 20
     
    6971        else:
    7072            try:
    71                 easting_min = eval('project_250m.xmin%s' % which_area)
    72                 easting_max = eval('project_250m.xmax%s' % which_area)
    73                 northing_min = eval('project_250m.ymin%s' % which_area)
    74                 northing_max = eval('project_250m.ymax%s' % which_area)
     73                easting_min = eval('project.xmin%s' % which_area)
     74                easting_max = eval('project.xmax%s' % which_area)
     75                northing_min = eval('project.ymin%s' % which_area)
     76                northing_max = eval('project.ymax%s' % which_area)
    7577            except AttributeError:
    7678                print 'Unrecognized area name: %s' % which_area
     
    7981        for time_dir in time_dirs:
    8082
    81             name1 = directory+time_dir+sep+project_250m.scenario_name
    82             #name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_52620_0'
    83            # name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_78360_0'
     83            name1 = join(directory, time_dir, project.scenario_name)
     84            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
     85            #name3 = join(directory, time_dir, project.scenario_name)+'_time_33360_0'
     86            #name4 = join(directory, time_dir, project.scenario_name)+'_time_50040_0'
     87            #name5 = join(directory, time_dir, project.scenario_name)+'_time_66720_0'
    8488
    85             names = [name1]#, name2] #, name3]
     89            names = [name1, name2]# , name3, name4, name5]
    8690     
    8791            asc_name = []   
     
    108112                asc_name.append(outname + '.asc')
    109113   
    110             maxasc_outname = directory+time_dir+sep+project_250m.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
     114            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
    111115
    112116            print 'max asc outname ', maxasc_outname
  • anuga_work/production/australia_ph2/sydney/project.py

    r6466 r6498  
    2929tide = 0                # difference between MSL and HAT
    3030zone = 50               # specify zone of model
    31 event_number = 6368    # the event number or the mux file name
     31event_number = 58326    # the event number or the mux file name
    3232alpha = 0.1             # smoothing parameter for mesh
    3333friction=0.01           # manning's friction coefficient
     
    3636
    3737# index is only used when wave = Tb
    38 index = 2971            # index from the PTHA
    39 wave = 'Tb'             # Bf (sts wave) Tb (index wave)
     38index = 2956            # index from the PTHA
     39wave = 'Bf'             # Bf (sts wave) Tb (index wave)
    4040
    4141setup = 'final'         # This can be one of three values
Note: See TracChangeset for help on using the changeset viewer.