Changeset 5393


Ignore:
Timestamp:
Jun 4, 2008, 3:20:10 PM (17 years ago)
Author:
jakeman
Message:

updating scripts for boxingday validation

Location:
anuga_work/development/boxingday08
Files:
5 added
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/boxingday08/combine_good_data.py

    r5247 r5393  
    8383if not exists(dir+sep+filename+'_clipped'+'.csv'):
    8484    clip_csv(dir+sep+filename,project.bounding_polygon)
     85else:
     86    print dir+sep+filename+'_clipped'+'.csv already exists ensure the .csv file has not changed'
    8587
    8688print 'Creating geospatial data'
     
    8991else:
    9092    G0 = Geospatial_data(file_name=dir+sep+filenames[0]+'.csv')
     93    print 'read bay'
    9194    G1=Geospatial_data(file_name=dir+sep+filenames[1]+'_clipped.csv')
     95    print 'read andaman'
    9296    G1.clip(project.bounding_polygon)
     97    print 'clipped andaman'
    9398    G1.clip_outside(G0)
    9499
  • anuga_work/development/boxingday08/export_results.py

    r5247 r5393  
    55from os import sep
    66
     7#name='good_simulation/good_simulation_polylinefriction20080529_115703'
     8#name='good_simulation/good_simulation_ursfriction20080529_160110'
     9name='good_simulation/good_simulationsourcefriction20080529_160144'
    710
    8 #time_dir = '20070322_035135' # birubi
    9 
    10 #directory = project_slide.outputdir
    11 #name = directory + time_dir + sep + project_slide.basename
    12 #name = 'test500000.035000'
    13 #name = 'test'
    14 name = 'good_simulationsourcefriction20080404_140016'
    15 is_parallel = False
    16 if is_parallel == True: nodes = 4
    1711print 'output dir:', name
    1812
     
    3832    quantityname = 'elevation'  #Elevation
    3933
    40 if is_parallel == True:
    41     for i in range(0,nodes):
    42         namei = name + '_P%d_%d' %(i,nodes)
    43         outnamei = outname + '_P%d_%d' %(i,nodes)
    44         print 'start sww2dem for sww file %d' %(i)
    45         sww2dem(namei, basename_out = outnamei,
    46                     quantity = quantityname,
    47                     timestep = 1,
    48                     cellsize = 100,     
    49                     #easting_min = project.e_min_area,
    50                     #easting_max = project.e_max_area,
    51                     #northing_min = project.n_min_area,
    52                     #northing_max = project.n_max_area,       
    53                     reduction = max,
    54                     verbose = True,
    55                     format = 'asc')
    56 else:
    57     print 'start sww2dem'
    58     sww2dem(name, basename_out = outname,
    59                 quantity = quantityname,
    60                 #timestep = 0,
    61                 cellsize = 50,     
    62                 #easting_min = 416920,
    63                 #easting_max = 425310,
    64                 #northing_min = 877864,
    65                 #northing_max = 869570,       
    66                 reduction = max,
    67                 verbose = True,
    68                 format = 'asc')
     34print 'start sww2dem'
     35sww2dem(name, basename_out = outname,
     36        quantity = quantityname,
     37        cellsize = 20,
     38        easting_min = 417348,
     39        easting_max = 425656,
     40        northing_min = 870565,
     41        northing_max = 877008,   
     42        reduction = max,
     43        verbose = True,
     44        format = 'asc')
  • anuga_work/development/boxingday08/project.py

    r5247 r5393  
    4444
    4545dir='mesh_polygons'
    46 #extent = 'extent.csv'
    47 extent = 'boundary_new_pts.csv'
     46extent = 'extent.csv'
     47#extent = 'boundary_new_pts.csv'#old boundary
    4848# bounding polygon for study area
    4949bounding_polygon = read_polygon(dir+sep+extent)
     
    7575#south = -5.0
    7676
    77 #max_extent of ol simulation is why does this work and not below
     77#max_extent of old simulation is why does this work and not below
    7878#7.40911081272 8.71484358635 99.1683687224 97.513856322
    7979
    8080#max extent of clipped MOST file
    81 east  = 100
    82 west  = 96.0
     81
     82east  = 97.7
     83west  = 96.7
    8384north = 9.0
    84 south = 6.0
    85 
    86 
     85south = 5.9
     86     
     87     
Note: See TracChangeset for help on using the changeset viewer.