source: anuga_work/development/gong_2008/export_results.py @ 6868

Last change on this file since 6868 was 5653, checked in by sexton, 17 years ago

updates for gong study

File size: 1.5 KB
Line 
1import project, os
2import sys
3from anuga.shallow_water.data_manager import sww2dem
4from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
5from os import sep
6
7time_dir = ''
8
9cellsize = 25
10name = project.outputdir+time_dir+sep+project.basename
11
12var = [2,3,4] # depth, speed and elevation
13
14for which_var in var:
15    if which_var == 0:  # Stage
16        outname = name + '_stage'
17        quantityname = 'stage'
18
19    if which_var == 1:  # Absolute Momentum
20        outname = name + '_momentum'
21        quantityname = '(xmomentum**2 + ymomentum**2)**0.5' 
22
23    if which_var == 2:  # Depth
24        outname = name + '_depth'
25        quantityname = 'stage-elevation' 
26
27    if which_var == 3:  # Speed
28        outname = name + '_speed'
29        quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
30
31    if which_var == 4:  # Elevation
32        outname = name + '_elevation'
33        quantityname = 'elevation'  #Elevation
34
35    print 'start sww2dem'
36    sww2dem(name, basename_out = outname,
37                quantity = quantityname,
38                #timestep = timestep,
39                cellsize = cellsize,     
40                easting_min = project.eastingmin,
41                easting_max = project.eastingmax,
42                northing_min = project.northingmin,
43                northing_max = project.northingmax,       
44                reduction = max, 
45                verbose = True,
46                format = 'asc')
47
Note: See TracBrowser for help on using the repository browser.