Changeset 6041


Ignore:
Timestamp:
Dec 1, 2008, 11:48:14 AM (15 years ago)
Author:
kristy
Message:

change of name

File:
1 moved

Legend:

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

    r6027 r6041  
    2626
    2727
    28 time_dir1 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
    29 ##time_dir2 = ''
    30 ##time_dir3 = ''
    31 ##time_dir4 = ''
    32 ##time_dir5 = ''
    33 ##time_dir6 = ''
     28time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
     29time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
     30time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
     31time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
     32time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
     33time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
    3434
    35 time_dirs = [time_dir1] #, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    36  
     35time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     36
    3737for time_dir in time_dirs:
    3838
     
    4343
    4444    # Process: Create File GDB...
    45     gp.CreateFileGDB_management(folder, "raster")
     45    #gp.CreateFileGDB_management(folder, "raster")
    4646
    4747    gp.Workspace = raster_gbd
     
    5656
    5757    generate_filename = []
    58     input_ascii = glob.glob(folder + '*.asc')
     58    input_ascii = glob.glob(folder + '*max.asc')
    5959
    6060    for infile in input_ascii:
     
    7070        print 'Output DEM ',output_DEM
    7171        output_extract = output_DEM + 'E'
    72 ##        print 'Output Extract ',output_extract
     72        print 'Output Extract ',output_extract
    7373       
    7474        print 'Process: ASCII to Raster'
     
    8080                                                   ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
    8181
    82 ##        print 'Process: Extract by Mask'
    83 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
     82        print 'Process: Extract by Mask'
     83        gp.ExtractByMask_sa(output_DEM, land, output_extract)
    8484
    8585
Note: See TracChangeset for help on using the changeset viewer.