Ejemplo n.º 1
0
#! /usr/bin/env/python
# import necessary modulees
import polyadcirc.run_framework.subdomain as subdom
import polyadcirc.run_framework.fulldomain as fulldom
import glob
import polyadcirc.run_framework.plotADCIRC as pa

adcirc_dir = '/h1/lgraham/workspace'
fulldomain = fulldom.fulldomain(adcirc_dir+'/ADCIRC_landuse/Inlet/inputs/tides')
subdomain = subdom.subdomain(adcirc_dir+'/ADCIRC_landuse/Inlet_subdomain/inputs/tides')
subdomain.set_fulldomain(fulldomain)
num_procs = 2

fulldomain.update()
subdomain.update()

# Compare subdomain and fulldomain results
ts_names = ['fort.63', 'fort.64']
nts_names = ['maxele.63','maxvel.63', 'tinun.63','elemaxdry.63']
subdomain.update_sub2full_map()
ts_error, nts_error, time_obs, ts_data, nts_data = subdomain.compare_to_fulldomain(ts_names, 
        nts_names)
subdomain.get_Triangulation(path = subdomain.path)
pa.ts_pcolor(ts_data, time_obs, subdomain, path = subdomain.path)
pa.nts_pcolor(nts_data, subdomain, path = subdomain.path)
pa.ts_quiver(ts_data, time_obs, subdomain, path = subdomain.path)
Ejemplo n.º 2
0
#! /usr/bin/env/python
# import necessary modulees
import polyadcirc.run_framework.subdomain as subdom
import polyadcirc.run_framework.fulldomain as fulldom
import polyadcirc.run_framework.random_manningsn as rmn
import polyadcirc.pyADCIRC.volume as vol
import numpy as np
import scipy.io as sio

adcirc_dir = '/h1/lgraham/workspace'
fulldomain = fulldom.fulldomain(adcirc_dir+'/ADCIRC_landuse/Inlet/inputs/tides')
subdomain = subdom.subdomain(adcirc_dir+'/ADCIRC_landuse/Inlet_subdomain/inputs/tides')
subdomain.set_fulldomain(fulldomain)

# The following steps correspond to Table 1 in Subdomain ADICRC v.50 User Guide
# Step 1a Generate Sudomain
# subdomain.circle(0, 0, 3000)
# subdomain.setup()

# Read in fort.14, fort.15 data to subdomain and fulldomain objects
fulldomain.update()
subdomain.update()

# Read in and mapping between subdomain and fulldomain
subdomain.update_sub2full_map()

# Get list of elements in the subdomain in terms of fulldomain element numbers
elements = subdomain.sub2full_element.values()

# Get fort.63 data
fort63, time_obs = rmn.get_ts_sr(fulldomain.path, 'fort.63')
Ejemplo n.º 3
0
#! /usr/bin/env/python
# import necessary modules
import polyadcirc.run_framework.subdomain as subdom
import polyadcirc.run_framework.fulldomain as fulldom
import glob

adcirc_dir = '/work/01837/lcgraham/v50_subdomain/work'
fulldomain = fulldom.fulldomain(adcirc_dir+'/fulldomain')
subdomain = subdom.subdomain(adcirc_dir+'/subdomain')
subdomain.set_fulldomain(fulldomain)
num_procs = 2

# check to see if shape file exists, if not make it
if len(glob.glob(subdomain.path+'/shape.*14')) <= 0:
    subdomain.ellipse([40824.6, 98559.5], [98559.5, 40824,6], 60000)

# The following steps correspond to Table 1 in Subdomain ADICRC v.50 User Guide
# Step 1a Generate Sudomain
subdomain.setup()
# Step 1b Generate Full Domain Control File
subdomain.genfull()
# Step 2 Run ADCIRC on the full domain
fulldomain.update()
if subdomain.check_fulldomain():
    disp =  "Output files ``fort.06*`` exist, but running ADCIRC on fulldomain"
    print disp+"anyway."
else:
    print "Output files ``fort.06*`` do not exist, running ADCIRC on fulldomain."
fulldomain.run(num_procs, adcirc_dir)
# Step 3 Exract Subdomain Boundary Conditions
subdomain.update()