def runSimulation():
conf = GppConfig()
gm = GliderModel(conf.myDataDir+'RiskMap.shelf',conf.myDataDir+'/roms5/')
u,v,time,depth,lat,lon = gm.GetRomsData(yy,mm,dd,numDays)
plt.figure(),plt.imshow(gm.riskMapArray,origin='lower')
FullSimulation,HoldValsOffMap=False,True
for i in range(0,200):
'''
xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
gm.SimulateDive(gm.lat_pts[1]+random.gauss(0,0.1),gm.lon_pts[6]+random.gauss(0,0.1),gm.lat_pts[5],gm.lon_pts[1],80,u,v,lat,lon,depth,i,False)
'''
gm.InitSim(gm.lat_pts[1]+random.gauss(0,0.1),gm.lon_pts[6]+random.gauss(0,0.1),gm.lat_pts[5],gm.lon_pts[1],80,i,FullSimulation,HoldValsOffMap)
xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
gm.SimulateDive_R(20)
print 'Flag for Simulation Completion = ',gm.doneSimulating
#xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
# gm.SimulateDive_R(True,10)
#gm.SimulateDive(gm.lat_pts[1]+random.gauss(0,0.1),gm.lon_pts[6]+random.gauss(0,0.1),gm.lat_pts[5],gm.lon_pts[1],80,u,v,lat,lon,depth,i,False)
if possibleCollision == False:
tempX,tempY = gm.GetPxPyFromLatLon(np.array(latArray),np.array(lonArray))
x_sims,y_sims = tempX[-1:],tempY[-1:] # TODO: This might be wrong!
plt.plot(tempX,tempY)
else:
tempX,tempY = gm.GetPxPyFromLatLon(np.array(latArray),np.array(lonArray))
plt.plot(tempX,tempY,'r.-')
x_sims,y_sims = 0,0
def MultiProcessedSimulator(queue):
conf = GppConfig()
gm = GliderModel(conf.myDataDir + 'RiskMap.shelf',
conf.myDataDir + '/roms/')
u, v, time, depth, lat, lon = gm.GetRomsData(yy, mm, dd, numDays)
simRuns = []
FullSimulation, HoldValsOffMap = False, True
for i in range(0, 50):
'''
xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
gm.SimulateDive(gm.lat_pts[1]+random.gauss(0,0.1),gm.lon_pts[6]+random.gauss(0,0.1),gm.lat_pts[5],gm.lon_pts[1],80,u,v,lat,lon,depth,i,False)
'''
gm.InitSim(gm.lat_pts[1] + random.gauss(0, 0.1),
gm.lon_pts[6] + random.gauss(0, 0.1), gm.lat_pts[5],
gm.lon_pts[1], 80, i, FullSimulation, HoldValsOffMap)
xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
gm.SimulateDive_R(20)
print 'Flag for Simulation Completion = ', gm.doneSimulating
tempX, tempY = None, None
#xFin,yFin,latFin,lonFin,latArray,lonArray,depthArray,tArray,possibleCollision,CollisionReason,totalDist=\
# gm.SimulateDive_R(True,10)
#gm.SimulateDive(gm.lat_pts[1]+random.gauss(0,0.1),gm.lon_pts[6]+random.gauss(0,0.1),gm.lat_pts[5],gm.lon_pts[1],80,u,v,lat,lon,depth,i,False)
if possibleCollision == False:
tempX, tempY = gm.GetPxPyFromLatLon(np.array(latArray),
np.array(lonArray))
x_sims, y_sims = tempX[-1:], tempY[
-1:] # TODO: This might be wrong!
#plt.plot(tempX,tempY)
else:
tempX, tempY = gm.GetPxPyFromLatLon(np.array(latArray),
np.array(lonArray))
#plt.plot(tempX,tempY,'r.-')
x_sims, y_sims = 0, 0
simRuns.append((tempX, tempY, possibleCollision))
queue.put(simRuns)