Changeset 1868


Ignore:
Timestamp:
Oct 5, 2005, 5:40:34 PM (18 years ago)
Author:
ole
Message:

Ongoing work with karratha

Location:
production/karratha_2005
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • production/karratha_2005/create_mesh.py

    r1855 r1868  
    9292    """
    9393
    94     from pmesh.mesh import *
     94    from pmesh.mesh import Mesh
    9595    from pyvolution.coordinate_transforms.redfearn import redfearn
    9696    from pyvolution.util import populate_polygon
     
    107107
    108108
     109
     110   
    109111    print "***********************"
    110112    print "geo ref", geo
  • production/karratha_2005/project.py

    r1857 r1868  
    6363
    6464#Interior regions
    65 karratha_south = degminsec2decimal_degrees(-20,45,0)
    66 karratha_north = degminsec2decimal_degrees(-20,37,0)
    67 karratha_west = degminsec2decimal_degrees(116,45,0)
    68 karratha_east = degminsec2decimal_degrees(116,55,0)
     65karratha_south = degminsec2decimal_degrees(-20,44,0)
     66karratha_north = degminsec2decimal_degrees(-20,42,0)
     67karratha_west = degminsec2decimal_degrees(116,48,0)
     68karratha_east = degminsec2decimal_degrees(116,53,30)
    6969
    7070k0 = [karratha_south, karratha_west]
     
    7575karratha_polygon = [k0, k1, k2, k3]
    7676
     77
     78#Interior regions
     79dampier_south = degminsec2decimal_degrees(-20,40,0)
     80dampier_north = degminsec2decimal_degrees(-20,38,10)
     81dampier_west = degminsec2decimal_degrees(116,43,0)
     82dampier_east = degminsec2decimal_degrees(116,45,0)
     83
     84d0 = [dampier_south, dampier_west]
     85d1 = [dampier_south, dampier_east]
     86d2 = [dampier_north, dampier_east]
     87d3 = [dampier_north, dampier_west]   
     88
     89dampier_polygon = [d0, d1, d2, d3]
     90
     91
     92#Interior regions
     93refinery_south = degminsec2decimal_degrees(-20,38,0)
     94refinery_north = degminsec2decimal_degrees(-20,36,0)
     95refinery_west = degminsec2decimal_degrees(116,44,0)
     96refinery_east = degminsec2decimal_degrees(116,46,10)
     97
     98d0 = [refinery_south, refinery_west]
     99d1 = [refinery_south, refinery_east]
     100d2 = [refinery_north, refinery_east]
     101d3 = [refinery_north, refinery_west]   
     102
     103refinery_polygon = [d0, d1, d2, d3]
  • production/karratha_2005/run_karratha.py

    r1855 r1868  
    1111"""
    1212
    13 tide = 0.75   #HMWS estimate by Colin French, GA
     13#tide = 0.75   #HMWS estimate by Colin French, GA
     14tide = 0       #HMW
    1415
    1516
     
    8384from create_mesh import create_mesh
    8485
    85 interior_regions = [[project.karratha_polygon, 10000]]
     86interior_regions = [[project.karratha_polygon, 25000],
     87                    [project.dampier_polygon, 8000],
     88                    [project.refinery_polygon, 8000]]
     89
    8690m = cache(create_mesh,
    8791          project.polygon,
    8892          {'boundary_tags': {'back': [7, 8], 'side': [0, 6], 'ocean': [1, 2, 3, 4, 5]},
    89           'resolution': 80000,
     93          'resolution': 100000,
    9094          'filename': project.meshname + '.msh',
    9195          'interior_regions': interior_regions},
     
    142146
    143147
    144 #Run
    145 #for t in domain.evolve(yieldstep = 600, finaltime = 15000):
    146 #    domain.write_time()
    147 #    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
    148 #
    149 #for t in domain.evolve(yieldstep = 10, finaltime = 35000):
    150 #    domain.write_time()
    151 #    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
     148#Evolve
     149import time
     150t0 = time.time()
    152151
    153 for t in domain.evolve(yieldstep = 60, finaltime = 40000):
     152for t in domain.evolve(yieldstep = 60, finaltime = 15000):
     153    domain.write_time()
     154    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
     155
     156for t in domain.evolve(yieldstep = 20, finaltime = 35000,
     157                       skip_initial_step = True):
     158    domain.write_time()
     159    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
     160
     161for t in domain.evolve(yieldstep = 60, finaltime = 40000,
     162                       skip_initial_step = True):
    154163    domain.write_time()
    155164    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')       
     165
     166print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracChangeset for help on using the changeset viewer.