Ignore:
Timestamp:
Dec 6, 2006, 5:35:00 PM (18 years ago)
Author:
sexton
Message:

updates for nsw slide modelling and cairns demo

File:
1 edited

Legend:

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

    r3514 r4063  
    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
    9 time_dir = '20060706_235036' #test with coarse grid
    10 #time_dir = '20060707_001859' #MSL DLI data
     9#time_dir = '20060706_235036' #test with coarse grid
     10time_dir = '20060707_001859' #MSL DLI data
    1111#time_dir = '20060707_003301' #HAT DLI data
    1212#time_dir = '20060707_003424' #LAT DLI data
     
    4949
    5050if which_var == 4:  # Elevation
    51     outname = name + '_elevation'
     51    outname = name + '_elevation_50'
    5252    quantityname = 'elevation'  #Elevation
    5353
     
    5656sww2dem(name, basename_out = outname,
    5757            quantity = quantityname,
    58             cellsize = 20,       # Trevor would like this at 25
     58            cellsize = 50,       # Trevor would like this at 25
    5959            # define region for viz purposes
    60             easting_min = project.e_min_area,
    61             easting_max = project.e_max_area,
    62             northing_min = project.n_min_area,
    63             northing_max = project.n_max_area,       
     60            #easting_min = project.e_min_area,
     61            #easting_max = project.e_max_area,
     62            #northing_min = project.n_min_area,
     63            #northing_max = project.n_max_area,       
    6464            reduction = max, #this is because we want max quantityname
    6565            verbose = True,
Note: See TracChangeset for help on using the changeset viewer.