Changeset 6949


Ignore:
Timestamp:
May 5, 2009, 9:35:20 AM (11 years ago)
Author:
myall
Message:

re running carnarvon, and exporting results

Location:
anuga_work/production/australia_ph2
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py

    r6741 r6949  
    2525output_dir="anuga\\outputs\\"
    2626
    27 time_dir1 = '20090324_150124_run_final_0_27285_109_Bf_mhingee'
    28 time_dir2 = '20090324_161452_run_final_0_27285_109_Tb_mhingee'
     27time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt'
     28time_dir2 = ''
    2929
    30 time_dirs = [time_dir2] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
     30time_dirs = [time_dir1] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
    3131 
    3232for time_dir in time_dirs:
     
    5252
    5353    generate_filename = []
    54     input_ascii = glob.glob(folder + '*.asc')
     54    input_ascii = glob.glob(folder + '*elevation.asc')
    5555
    5656    for infile in input_ascii:
     
    7070
    7171        print 'Process: Define Projection'
    72         gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_49',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
    73                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',111.0],PARAMETER['Scale_Factor',0.9996]"
     72        gp.DefineProjection_management(output_DEM, "PROJCS['CM_113.5',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
     73                                                   ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',113.5],PARAMETER['Scale_Factor',0.9996]"
    7474                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    7575##        output_extract = output_DEM + '_E'
  • anuga_work/production/australia_ph2/carnarvon/project.py

    r6856 r6949  
    4040friction=0.01           # manning's friction coefficient
    4141starttime=0             # start time for simulation
    42 finaltime=1000 # final time for simulation
    43 setup = 'trial'         # This can be one of three values
     42finaltime=60000 # final time for simulation
     43setup = 'final'         # This can be one of three values
    4444                        #    trial - coarsest mesh, fast
    4545                        #    basic - coarse mesh
     
    128128# Thinned ordering file from Hazard Map (geographic)
    129129# Format is index,latitude,longitude (with header)
    130 urs_order_filename = 'urs_order_simple.csv'
     130urs_order_filename = 'urs_order.csv'
    131131
    132132# Landward bounding points
  • anuga_work/production/australia_ph2/gold_coast/Arc_asc2raster_GDA94z50.py

    r6842 r6949  
    2929time_dir3 = '20090409_175154_run_final_0_51453_3727_Tb__kvanputt'
    3030
    31 time_dirs = [time_dir1, time_dir2, time_dir3]
     31##time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee'
     32
     33time_dirs = [time_dir1]#, time_dir2, time_dir3]
    3234
    3335for time_dir in time_dirs:
     
    3941##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean"
    4042
    41     print 'Process: Create File GDB'
    42     gp.CreateFileGDB_management(folder, "raster")
     43##    print 'Process: Create File GDB'
     44##    gp.CreateFileGDB_management(folder, "raster")
    4345
    4446    gp.Workspace = raster_gbd
     
    5355
    5456    generate_filename = []
    55     input_ascii = glob.glob(folder + '*max.asc')
     57    input_ascii = glob.glob(folder + '*elevation.asc')
    5658
    5759    for infile in input_ascii:
Note: See TracChangeset for help on using the changeset viewer.