Changeset 7273


Ignore:
Timestamp:
Jun 29, 2009, 4:27:55 PM (14 years ago)
Author:
kristy
Message:
 
Location:
anuga_work/production/exmouth_2009
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/exmouth_2009/Arc_asc2raster_GDA94z50.py

    r7209 r7273  
    2222gp.overwriteoutput = 1
    2323
    24 scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\onslow_tsunami_scenario_2009\\"
     24scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\exmouth_tsunami_scenario_2009\\"
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090616_154908_run_final_0.0_27283_mhingee'
    28 ##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
    29 ##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
     27time_dir1 = '20090624_145617_run_final_0.0_27283_mhingee'
     28##time_dir2 = ''
     29##time_dir3 = ''
    3030
    31 ##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
     31##time_dir2 = ''
    3232
    3333time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
     
    4848   
    4949    #replication dictionary
    50     replicate = (('onslow', ''),('_', ''),('max','_M'),
     50    replicate = (('exmouth', ''),('_', ''),('max','_M'),
    5151                 ('CBD', 'CDB'),('All',''),
    5252                 ('depth','_depth'),('speed', '_speed'),
     
    5454
    5555    generate_filename = []
    56     input_ascii = glob.glob(folder + '*elevation.asc')
     56    input_ascii = glob.glob(folder + '*.asc')
    5757
    5858    for infile in input_ascii:
     
    7373        print 'Process: Define Projection'
    7474        #GDA_1994_MGA_Zone_54
    75         gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    76                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"
     75        gp.DefineProjection_management(output_DEM, "PROJCS['CM_114',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
     76                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',114.0],PARAMETER['Scale_Factor',0.9996]"
    7777                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    7878##        output_extract = output_DEM + 'E'
  • anuga_work/production/exmouth_2009/build_elevation.py

    r7209 r7273  
    4545#------------------------------------------------------------------------------
    4646
    47 print 'project.bounding_polygon', project.bounding_polygon
    48 print 'project.combined_elevation_basename', project.combined_elevation_basename
     47##print 'project.bounding_polygon', project.bounding_polygon
     48##print 'project.combined_elevation_basename', project.combined_elevation_basename
    4949
    5050# Create Geospatial data from ASCII files
  • anuga_work/production/exmouth_2009/export_results.py

    r7209 r7273  
    2222directory = project.output_folder
    2323
    24 time_dir1 = '20090616_154908_run_final_0.0_27283_mhingee'
    25 time_dir2 = '20090617_110255_run_final_0.0_27283_mhingee'
     24time_dir1 = '20090624_145617_run_final_0.0_27283_mhingee'
     25##time_dir2 = ''
    2626##time_dir3 = ''
    2727##time_dir4 = ''
     
    2929##time_dir6 = ''
    3030##
    31 time_dirs = [time_dir2]#, time_dir2]# , time_dir4, time_dir5, time_dir6]
     31time_dirs = [time_dir1]#, time_dir2]# , time_dir4, time_dir5, time_dir6]
    3232 
    3333
     
    4242######
    4343
    44 #area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
     44#area = ['Exmouth']  # strings must match keys in var_equations below
    4545area = ['All']      # 'All' means no special areas - the whole thing
    4646
     
    5656
    5757# one or more key strings from var_equations above
    58 var = ['elevation']
     58var = ['elevation','stage']
    5959
    6060######
Note: See TracChangeset for help on using the changeset viewer.