Changeset 3744


Ignore:
Timestamp:
Oct 11, 2006, 12:22:40 PM (18 years ago)
Author:
ole
Message:

Made karratha production use newly resolved convert_from_latlon_to_UTM

Location:
anuga_work/production/karratha_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/karratha_2006/project.py

    r3716 r3744  
    88from anuga.coordinate_transforms.redfearn import\
    99     degminsec2decimal_degrees,\
    10      convert_points_from_latlon_to_utm
     10     convert_from_latlon_to_utm
    1111
    1212from time import localtime, strftime
    … …  
    5252outputname = outputtimedir + basename  #Used by post processing
    5353#boundaryname = outputtimedir + basename  #Used by post processing
     54
    5455
    5556
    … …  
    7576
    7677bounding_polygon, zone =\
    77                   convert_points_from_latlon_to_utm([p0, p1, p2, p3, p4,
    78                                                      p5, p6, p7, p8])
     78                  convert_from_latlon_to_utm([p0, p1, p2, p3, p4,
     79                                              p5, p6, p7, p8])
    7980refzone = zone
     81
     82
     83#Interior regions
     84
     85
     86# CIPMA point of interest
     87cipma_latitude = -20.588456
     88cipma_longitude = 116.771527
     89
     90
     91k0 = [cipma_latitude-0.02, cipma_longitude-0.02]
     92k1 = [cipma_latitude-0.02, cipma_longitude+0.02]
     93k2 = [cipma_latitude+0.02, cipma_longitude+0.02]
     94k3 = [cipma_latitude+0.02, cipma_longitude-0.02]
     95
     96cipma_polygon, zone = convert_from_latlon_to_utm([k0, k1, k2, k3])
     97assert zone == refzone
     98
    8099
    81100
    … …  
    91110k3 = [karratha_north, karratha_west]   
    92111
    93 karratha_polygon, zone = convert_points_from_latlon_to_utm([k0, k1, k2, k3])
     112karratha_polygon, zone = convert_from_latlon_to_utm([k0, k1, k2, k3])
    94113assert zone == refzone
    95114
    … …  
    106125d3 = [dampier_north, dampier_west]   
    107126
    108 dampier_polygon, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
     127dampier_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
    109128assert zone == refzone
    110129
    … …  
    121140d3 = [refinery_north, refinery_west]   
    122141
    123 refinery_polygon, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
     142refinery_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
    124143assert zone == refzone
    125144
    … …  
    138157d3 = [point_north, point_west]   
    139158
    140 point_polygon, zone = convert_points_from_latlon_to_utm([d0, d1, d2, d3])
     159point_polygon, zone = convert_from_latlon_to_utm([d0, d1, d2, d3])
    141160assert zone == refzone
    142161
    … …  
    152171                degminsec2decimal_degrees(116,46,17)]
    153172
    154 neil1_polygon, zone = convert_points_from_latlon_to_utm([neil1_point1,
    155                                                          neil1_point2,
    156                                                          neil1_point3,
    157                                                          neil1_point4])
     173neil1_polygon, zone = convert_from_latlon_to_utm([neil1_point1,
     174                                                  neil1_point2,
     175                                                  neil1_point3,
     176                                                  neil1_point4])
    158177assert zone == refzone
    159178
    … …  
    170189                degminsec2decimal_degrees(116,43,9)]
    171190
    172 neil2_polygon, zone = convert_points_from_latlon_to_utm([neil2_point1,
    173                                                          neil2_point2,
    174                                                          neil2_point3,
    175                                                          neil2_point4])
     191neil2_polygon, zone = convert_from_latlon_to_utm([neil2_point1,
     192                                                  neil2_point2,
     193                                                  neil2_point3,
     194                                                  neil2_point4])
    176195assert zone == refzone
    177196
    … …  
    190209                degminsec2decimal_degrees(116,46,17)]
    191210
    192 wb_polygon, zone = convert_points_from_latlon_to_utm([wb_point1, wb_point2,
    193                                                       wb_point3, wb_point4])
     211wb_polygon, zone = convert_from_latlon_to_utm([wb_point1, wb_point2,
     212                                               wb_point3, wb_point4])
    194213assert zone == refzone
    195214
    … …  
    208227                degminsec2decimal_degrees(116,42,50)]
    209228
    210 lwb_polygon, zone = convert_points_from_latlon_to_utm([lwb_point1, lwb_point2,
    211                                                        lwb_point3, lwb_point4])
     229lwb_polygon, zone = convert_from_latlon_to_utm([lwb_point1, lwb_point2,
     230                                                lwb_point3, lwb_point4])
    212231                                                     
    213232assert zone == refzone

  • anuga_work/production/karratha_2006/run_karratha.py

    r3716 r3744  
    7171                    #[project.refinery_polygon, 2000],
    7272                    #[project.point_polygon, 2000]]
    73                     [project.neil1_polygon, 4000],
    74                     [project.neil2_polygon, 64000]]
     73                    #[project.cipma_polygon, 20000]]
     74                    [project.cipma_polygon, 4000]]     # Caused memory error?   
     75                    #[project.neil1_polygon, 4000],
     76                    #[project.neil2_polygon, 64000]]
    7577
    7678
    … …  
    8284                                        'ocean': [1, 2, 3, 4, 5]},
    8385                         maximum_triangle_area=200000,
     86                         interior_regions=interior_regions,
    8487                         filename=meshname,
    8588                         use_cache=True,
Note: See TracChangeset for help on using the changeset viewer.