Changeset 3650 for anuga_work


Ignore:
Timestamp:
Sep 22, 2006, 12:21:09 PM (18 years ago)
Author:
sexton
Message:

updates for Hobart and want to export elevation grid for Onslow

Location:
anuga_work/production
Files:
5 edited

Legend:

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

    r3626 r3650  
    22import sys
    33
    4 from anuga.pyvolution.data_manager import sww2dem
    5 from anuga.pyvolution.ermapper_grids import read_ermapper_grid
    6 from anuga.pyvolution.util import Screen_Catcher
     4from anuga.shallow_water.data_manager import sww2dem
    75from os import sep
    86
    9 time_dir = '20060919_050831' #MSL
     7time_dir = '20060921_083846' #MSL
    108directory = project.outputdir
    119name = directory + time_dir + sep + 'source'
     
    4038sww2dem(name, basename_out = outname,
    4139            quantity = quantityname,
    42             cellsize = 20,       
     40            cellsize = 25,       
    4341            # define region for viz purposes
    4442            easting_min = project.e_min_area,
  • anuga_work/production/hobart_2006/project.py

    r3638 r3650  
    1313
    1414# data provided by Tas SES and checked by NM&I
    15 onshore_name = 'hob3_topo'
     15#onshore_name = 'hob3_topo' # original
     16onshore_name = 'hob5_topo' # 12.5m grid and clipped to 100m elevation or 3000m from coast
     17onshore_name_25 = 'hob5_topo_25' # 25m grid and clipped to 100m elevation or 3000m from coast
    1618offshore_name_tas1 = 'derwent_2m'
    1719offshore_name_tas2 = 'derwent_5m'
    18 offshore_name_tas3 = 'south_east_tas'
     20offshore_name_tas3 = 'south_east_tas' #actually this is AHO
    1921offshore_name_tas4 = 'hobart_1m'
    2022
     
    8587
    8688onshore_dem_name = datadir + onshore_name
     89onshore_dem_name_25 = datadir + onshore_name_25
     90all_onshore_dem_name = datadir + 'combined_onshore'
    8791offshore_dem_name_local1 = datadir + offshore_name_tas1
    8892offshore_dem_name_local2 = datadir + offshore_name_tas2
     
    110114outputname = outputtimedir + basename  #Used by post processing
    111115
    112 # region to export
     116# clipping 12.5m onshore data set
     117eastingmin = 520000
     118eastingmax = 536000
     119northingmin = 5245000
     120northingmax = 5260000
     121
     122# region to export for inundation map
    113123e_min_area = 500000
    114124e_max_area = 580000
     
    132142polyAll, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
    133143refzone = zone
     144NWcorner = polyAll[0]
     145SEcorner = polyAll[2]
     146
     147# region to export for Alex to make bathymetry map
     148e_min_area = 480000#NWcorner[0]
     149e_max_area = 600000#SEcorner[0]
     150n_min_area = 5156000#NWcorner[1]
     151n_max_area = 5285000#SEcorner[1]
    134152
    135153#Interior region - Hobart city area + Glenorchy, Kingston
  • anuga_work/production/hobart_2006/run_hobart.py

    r3642 r3650  
    2424                            Dirichlet_boundary, Time_boundary, File_boundary
    2525from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
     26from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files
    2627from anuga.geospatial_data.geospatial_data import *
    2728from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
     
    6162# filenames
    6263onshore_dem_name = project.onshore_dem_name
     64onshore_dem_name_25 = project.onshore_dem_name_25
    6365meshname = project.meshname+'.msh'
    6466source_dir = project.boundarydir
     
    6971convert_dem_from_ascii2netcdf(onshore_dem_name, use_cache=True, verbose=True)
    7072
    71 #creates pts file for onshore DEM
    72 dem2pts(onshore_dem_name, use_cache=True, verbose=True)
    73 
    74 print 'local offshore data sets'
     73#creates pts file for onshore DEM - 12.5
     74dem2pts(onshore_dem_name,
     75        easting_min=project.eastingmin,
     76        easting_max=project.eastingmax,
     77        northing_min=project.northingmin,
     78        northing_max= project.northingmax,
     79        use_cache=True, verbose=True)
     80
     81#creates pts file for onshore DEM - 25
     82dem2pts(onshore_dem_name_25, use_cache=True, verbose=True)
     83
     84combine_rectangular_points_files(project.onshore_dem_name + '.pts',
     85                                 project.onshore_dem_name_25 + '.pts',
     86                                 project.all_onshore_dem_name + '.pts')
     87#print 'local offshore data sets'
    7588#G1 = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya')
    7689#G2 = Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya')
     
    7891#G4 = Geospatial_data(file_name = project.offshore_dem_name_local4 + '.xya')
    7992#G11 = G1 + G2 + G3 + G4
    80 print 'aho offshore data sets'
    81 G5 = Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')
    82 G6 = Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')
    83 G7 = Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')
    84 G8 = Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')
    85 G9 = Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')
    86 G10 = Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')
    87 G11 = Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')
    88 G12 = Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')
    89 G13 = Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')
    90 G14 = Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')
    91 G15 = Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')
    92 G16 = Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')
    93 G17 = Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')
    94 G18 = Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')
    95 G19 = Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')
    96 G20 = Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')
    97 G12 = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20
    98 print 'local onshore dem'
    99 G21 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
     93#print 'aho offshore data sets'
     94#G5 = Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')
     95#G6 = Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')
     96#G7 = Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')
     97#G8 = Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')
     98#G9 = Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')
     99#G10 = Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')
     100#G11 = Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')
     101#G12 = Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')
     102#G13 = Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')
     103#G14 = Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')
     104#G15 = Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')
     105#G16 = Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')
     106#G17 = Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')
     107#G18 = Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')
     108#G19 = Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')
     109#G20 = Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')
     110#G12 = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20
     111#print 'local onshore dem'
     112#G21 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
    100113print 'coast'
    101 G22 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
    102 print 'adding data sets'
     114#G22 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
     115#print 'adding data sets'
    103116#G = G11 + G12 + G21 + G22
    104 G = G21
     117#G = G1+G2+G3+G4+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22
     118#G = G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20+G21+G22
     119#G = G21+G22+G5+G6+G7+G8+G9+G10+G11+G12+G13+G14+G15+G16+G17+G18+G19+G20
     120#G = G21+G22+G5
     121#G = Geospatial_data(file_name = project.onshore_dem_name + '.pts') #+ \
     122#    Geospatial_data(file_name = project.coast_dem_name + '.xya') #+ \
     123#    Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')
     124G = Geospatial_data(file_name = project.offshore_dem_name_local1 + '.xya')+\
     125    Geospatial_data(file_name = project.offshore_dem_name_local2 + '.xya')+\
     126    Geospatial_data(file_name = project.offshore_dem_name_local3 + '.xya')+\
     127    Geospatial_data(file_name = project.offshore_dem_name_local4 + '.xya')+\
     128    Geospatial_data(file_name = project.offshore_dem_name_aho1 + '.xya')+\
     129    Geospatial_data(file_name = project.offshore_dem_name_aho2 + '.xya')+\
     130    Geospatial_data(file_name = project.offshore_dem_name_aho3 + '.xya')+\
     131    Geospatial_data(file_name = project.offshore_dem_name_aho4 + '.xya')+\
     132    Geospatial_data(file_name = project.offshore_dem_name_aho5 + '.xya')+\
     133    Geospatial_data(file_name = project.offshore_dem_name_aho6 + '.xya')+\
     134    Geospatial_data(file_name = project.offshore_dem_name_aho7 + '.xya')+\
     135    Geospatial_data(file_name = project.offshore_dem_name_aho8 + '.xya')+\
     136    Geospatial_data(file_name = project.offshore_dem_name_aho9 + '.xya')+\
     137    Geospatial_data(file_name = project.offshore_dem_name_aho10 + '.xya')+\
     138    Geospatial_data(file_name = project.offshore_dem_name_aho11 + '.xya')+\
     139    Geospatial_data(file_name = project.offshore_dem_name_aho12 + '.xya')+\
     140    Geospatial_data(file_name = project.offshore_dem_name_aho13 + '.xya')+\
     141    Geospatial_data(file_name = project.offshore_dem_name_aho14 + '.xya')+\
     142    Geospatial_data(file_name = project.offshore_dem_name_aho15 + '.xya')+\
     143    Geospatial_data(file_name = project.offshore_dem_name_aho16 + '.xya')+\
     144    Geospatial_data(file_name = project.all_onshore_dem_name + '.pts')
    105145G.export_points_file(project.combined_dem_name + '.pts')
    106146
  • anuga_work/production/onslow_2006/export_results.py

    r3514 r3650  
    22import sys
    33
    4 from anuga.pyvolution.data_manager import sww2dem
    5 from anuga.pyvolution.ermapper_grids import read_ermapper_grid
    6 from anuga.pyvolution.util import Screen_Catcher
     4from anuga.shallow_water.data_manager import sww2dem
     5#from anuga.pyvolution.ermapper_grids import read_ermapper_grid
     6from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
    77from os import sep
    88
    99#time_dir = '20060704_063005' #HAT
    10 time_dir = '20060706_235246' #LAT
    11 #time_dir = '20060704_063234' #MSL
     10#time_dir = '20060706_235246' #LAT
     11time_dir = '20060704_063234' #MSL
    1212#time_dir = '20060515_001733' #DTED data
    1313directory = project.outputdir
     
    2525
    2626#used to catch screen output to file
    27 sys.stdout = Screen_Catcher(screen_output_name)
    28 sys.stderr = Screen_Catcher(screen_error_name)
     27#sys.stdout = Screen_Catcher(screen_output_name)
     28#sys.stderr = Screen_Catcher(screen_error_name)
    2929
    3030print 'output dir:', name
  • anuga_work/production/onslow_2006/project.py

    r3535 r3650  
    88import sys
    99
    10 from anuga.pmesh.create_mesh import convert_points_from_latlon_to_utm
     10from anuga.coordinate_transforms.redfearn import convert_points_from_latlon_to_utm
    1111
    1212from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
     
    134134n_max_area = 7610000
    135135
     136# region to export to make elevation map: JS 22/9/06
     137
     138e_min_area = 240000
     139e_max_area = 340000
     140n_min_area = 7580000
     141n_max_area = 7690000
     142
    136143#Georeferencing
    137144from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees
Note: See TracChangeset for help on using the changeset viewer.