Ignore:
Timestamp:
May 22, 2009, 4:40:11 PM (11 years ago)
Author:
rwilson
Message:

Fiddling with layout of user guide.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/demos/cairns/ExportResults.py

    r5408 r7064  
    1 import project, os
     1import os
    22import sys
     3import project
    34
    45from anuga.shallow_water.data_manager import sww2dem
     
    78
    89scenario = 'slide'
    9 #scenario = 'fixed_wave'
    1010
    1111name = scenario + sep + scenario + 'source'
     
    1313print 'output dir:', name
    1414which_var = 4
    15 if which_var == 0:  # Stage
     15
     16if which_var == 0:    # Stage
    1617    outname = name + '_stage'
    1718    quantityname = 'stage'
    1819
    19 if which_var == 1:  # Absolute Momentum
     20if which_var == 1:    # Absolute Momentum
    2021    outname = name + '_momentum'
    21     quantityname = '(xmomentum**2 + ymomentum**2)**0.5'  #Absolute momentum
     22    quantityname = '(xmomentum**2 + ymomentum**2)**0.5'    #Absolute momentum
    2223
    23 if which_var == 2:  # Depth
     24if which_var == 2:    # Depth
    2425    outname = name + '_depth'
    2526    quantityname = 'stage-elevation'  #Depth
    2627
    27 if which_var == 3:  # Speed
     28if which_var == 3:    # Speed
    2829    outname = name + '_speed'
    2930    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
    3031
    31 if which_var == 4:  # Elevation
     32if which_var == 4:    # Elevation
    3233    outname = name + '_elevation'
    3334    quantityname = 'elevation'  #Elevation
     
    3536print 'start sww2dem'
    3637
    37 sww2dem(name, basename_out = outname,
    38             quantity = quantityname,
    39             cellsize = 100,     
    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 = 'ers')
    47 
     38sww2dem(name,
     39        basename_out=outname,
     40        quantity=quantityname,
     41        cellsize=100,     
     42        easting_min=project.eastingmin,
     43        easting_max=project.eastingmax,
     44        northing_min=project.northingmin,
     45        northing_max=project.northingmax,       
     46        reduction=max,
     47        verbose=True,
     48        format='ers')
Note: See TracChangeset for help on using the changeset viewer.