Changeset 5377


Ignore:
Timestamp:
May 29, 2008, 8:55:23 AM (16 years ago)
Author:
kristy
Message:

updated

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/busselton/run_boundary_points.py

    r5000 r5377  
    2121from anuga.shallow_water.data_manager import URS_points_needed_to_file, \
    2222     urs_ungridded2sww
    23 import project_new                 # Definition of file names and polygons
     23import project                 # Definition of file names and polygons
     24
     25#-----------------------------------------------------------------------------
     26# run_boundary_points
     27#-----------------------------------------------------------------------------
     28 
     29LL_LAT = -50.0          #ll_lat - lower left latitude, in decimal degrees
     30LL_LONG = 80.0          #ll-long - lower left longitude, in decimal degrees
     31GRID_SPACING = 2.0/60.0 #grid_spacing - in deciamal degrees
     32LAT_AMOUNT = 4800       #lat_amount - number of latitudes
     33LONG_AMOUNT = 3600      #long_amount- number of longs
     34ZONE = 50
    2435
    2536#------------------------------------------------------------------------------
     
    2738# output to file
    2839#------------------------------------------------------------------------------
    29 poly = project_new.poly_all
     40poly = project.poly_all
    3041print "poly", poly
    31 print "project.boundaries_in_dir_name",project_new.boundaries_in_dir_name
    32 URS_points_needed_to_file(project_new.boundaries_in_dir_name,poly, export_csv=True, zone=50)
     42print "project.boundaries_in_dir_name",project.boundaries_in_dir_name
     43URS_points_needed_to_file(project.boundaries_in_dir_name,poly, ZONE,
     44                          LL_LAT, LL_LONG,GRID_SPACING,LAT_AMOUNT, LONG_AMOUNT,export_csv=True)
     45#URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True, zone=50)
     46#URS_points_needed_to_file(project.boundaries_dir_name5,poly, export_csv=True)
    3347#urs_ungridded2sww(project.boundaries_in_dir_name2, project.boundaries_dir_name5, mean_stage= project.tide,  verbose=True)
Note: See TracChangeset for help on using the changeset viewer.