# -*- coding: utf-8 -*- """ Created on Wed Dec 4 12:07:23 2013 @author: team """ import oilspill_wrapper oilspill_wrapper.mul_GNOME_inputs(12, 2013, 12, 21, 0, 48) oilspill_wrapper.run_mul_GNOME(12, 666393, 3076780, 2013, 12, 21, 48, 900) oilspill_wrapper.GNOME_GM_visualization(12) oilspill_wrapper.GNOME_GE_animation(12, 13, 2013, 12, 21) # oilspill_wrapper.rk4(12,900,900,48,666393,3076780)
# -*- coding: utf-8 -*- """ Created on Wed Dec 4 12:07:23 2013 @author: team """ import oilspill_wrapper oilspill_wrapper.mul_GNOME_inputs(12,2013,12,21,0,48) oilspill_wrapper.run_mul_GNOME(12,666393,3076780,2013,12,21,48,900) oilspill_wrapper.GNOME_GM_visualization(12) oilspill_wrapper.GNOME_GE_animation(12,13,2013,12,21) # oilspill_wrapper.rk4(12,900,900,48,666393,3076780)
def HyosPy(starttime, endtime, period, mode=1):
"""
input:
starttime='2014-02-21'
endtime='2014-02-22'
period=24 oil spill simulation duration;
mode = 1
Mode 1:
download ROMS output from the server and run GNOME, have a quicker look (takes a few seconds)
Mode 2:
download river, wind, tidal, ROMS initial data, run SUNTANS and GNOME, have an accurate look
(takes a few minutes)
Mode 3:
a combination of mode 1 and 2, use the output of both ROMS and SUNTANS, duplicate the particles
that enters the overlap region
Some suggestions:
1) if the initial location is out of SUNTANS domain, run both mode 1 and mode 3
2) if the initial location is in SUNTANS domain, simply run mode 2
"""
#starttime='2014-02-21'
#endtime='2014-02-22'
dire=os.getcwd()
if mode==1:
downloadROMS(starttime,endtime,ROMSsrc='hindcast')
ROMS_file=dire+'/DATA/'+'txla_subset_HIS.nc'
ROMS_out=dire+'/GNOME/'+'hiroms_ss_rho.nc'
oilspill_wrapper.init_model(opt='ROMS')
oilspill_wrapper.HIROMS(ROMS_file,ROMS_out)
oilspill_wrapper.run_mul_GNOME(311584.1,3113650.2,starttime,endtime,period,900, opt='ROMS')
oilspill_wrapper.GNOME_GM_visualization(opt='ROMS')
oilspill_wrapper.GNOME_GE_animation(13,starttime,opt='ROMS')
elif mode==2:
infile=dire+'/CoarseTri/rundata/'+'GalvCoarse_0000.nc'
outfile=dire+'/GNOME/'+'txsuntans.nc'
hydro_wrapper.runSUNTANS(starttime, endtime)
oilspill_wrapper.init_model()
oilspill_wrapper.Tx_SUNTANS(infile,outfile) ##input: SUNTANS input; output: GNOME input
oilspill_wrapper.run_mul_GNOME(321947.94,3256260.05,starttime,endtime,period,900)
oilspill_wrapper.GNOME_GM_visualization()
oilspill_wrapper.GNOME_GE_animation(13,starttime)
print 'end'
elif mode==3:
print "under developing!!\n"
####use the ROMS current data to run GNOME####
downloadROMS(starttime,endtime, ROMSsrc='hindcast')
ROMS_file=dire+'/DATA/'+'txla_subset_HIS.nc'
ROMS_out=dire+'/GNOME/'+'hiroms_ss_rho.nc'
oilspill_wrapper.init_model(opt='ROMS') #The option is 'ROMS', since will run GNOME using ROMS output first
oilspill_wrapper.HIROMS(ROMS_file,ROMS_out)
#oilspill_wrapper.run_mul_GNOME(311584.1,3113650.2,starttime,endtime,period,900, opt='ROMS')
#(utm_x,utm_y)=utm.from_latlon(28.600572, -94.728385)[0:2]
(utm_x,utm_y)=utm.from_latlon(28.353786, -95.315109)[0:2]
oilspill_wrapper.run_mul_GNOME(utm_x,utm_y,starttime,endtime,period,900, opt='ROMS')
oilspill_wrapper.GNOME_GM_visualization(opt='ROMS')
####duplicate the particles####
oilspill_wrapper.duplicate(starttime,endtime,900)
oilspill_wrapper.visualization2()
else:
print "There is no such mode, check HyosPy input mode!!!\n"
