def populate (request):
from netCDF4 import Dataset, num2date
url = config.datasetpath
nc = Dataset(url)
lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
latc = nc.variables['latc'][:]
lonc = nc.variables['lonc'][:]
nv = nc.variables['nv'][:,:]
times = nc.variables['time']
time = num2date(times[:], units=times.units)
nodes = Node.objects.all()
nodes.delete()
cells = Cell.objects.all()
cells.delete()
times = Time.objects.all()
times.delete()
for i,d in enumerate(time):
#d = datetime.datetime(int(time[i][0]), int(time[i][1]), int(time[i][2]),\
# int(time[i][3]), int(time[i][4]), int(time[i][5]))
t = Time(date=d, index=i+1)
t.save()
for i in range(len(lat)):
n = Node(index=i+1, lat=lat[i], lon=lon[i])
n.save()
for i in range(len(latc)):
c = Cell(index=i+1, lat=latc[i], lon=lonc[i], node1=nv[0, i], node2=nv[1, i], node3=nv[2, i])
c.save()
response = HttpResponse()
response.write("done!")
return response
def Update_Geographic_Grid(self, request, queryset):
short_description = "Reload grid to populate server database."
from netCDF4 import Dataset, num2date
#import datetime
#import win32com.client
#import numpy
url = "http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3"
nc = Dataset(url)
lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
latc = nc.variables['latc'][:]
lonc = nc.variables['lonc'][:]
nv = nc.variables['nv'][:,:]
times = nc.variables['time']
time = num2date(times[:], units=times.units)
#def MakeCell(a, b, c, d, e, f):
# c = Cell(index=a, lat=b, lon=c, node1=d, node2=e, node3=f)
# c.save()
#def MakeNode(a, b, c):
# n = Node(index=a, lat=b, lon=c)
# n.save()
#def MakeTime(a, b):
# t = Time(date=a, index=b)
# t.save()
if queryset.model == Cell:
queryset.delete()
# map(MakeCell, range(len(latc))+1, latc, lonc, nv[0, :], nv[1, :], nv[2, :])
for i in range(len(latc)):
c = Cell(index=i+1, lat=latc[i], lon=lonc[i], node1=nv[0, i], node2=nv[1, i], node3=nv[2, i])
c.save()
elif queryset.model == Node:
queryset.delete()
#map(MakeNode, range(len(lat))+1, lat, lon)
for i in range(len(lat)):
n = Node(index=i+1, lat=lat[i], lon=lon[i])
n.save()
elif queryset.model == Time:
queryset.delete()
#map(MakeTime, time, range(len(time))+1)
for i,d in enumerate(time):
#d = datetime.datetime(int(time[i][0]), int(time[i][1]), int(time[i][2]),\
# int(time[i][3]), int(time[i][4]), int(time[i][5]))
t = Time(date=d, index=i+1)
t.save()
def populate(request):
from netCDF4 import Dataset, num2date
url = config.datasetpath
nc = Dataset(url)
lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
latc = nc.variables['latc'][:]
lonc = nc.variables['lonc'][:]
nv = nc.variables['nv'][:, :]
times = nc.variables['time']
time = num2date(times[:], units=times.units)
nodes = Node.objects.all()
nodes.delete()
cells = Cell.objects.all()
cells.delete()
times = Time.objects.all()
times.delete()
for i, d in enumerate(time):
#d = datetime.datetime(int(time[i][0]), int(time[i][1]), int(time[i][2]),\
# int(time[i][3]), int(time[i][4]), int(time[i][5]))
t = Time(date=d, index=i + 1)
t.save()
for i in range(len(lat)):
n = Node(index=i + 1, lat=lat[i], lon=lon[i])
n.save()
for i in range(len(latc)):
c = Cell(index=i + 1,
lat=latc[i],
lon=lonc[i],
node1=nv[0, i],
node2=nv[1, i],
node3=nv[2, i])
c.save()
response = HttpResponse()
response.write("done!")
return response
nv = nc.variables['nv'][:, :]
times = nc.variables['time']
time = num2date(times[:], units=times.units)
#matlab = win32com.client.Dispatch("matlab.application")
#matlab.Execute('''
#nc = cfdataset('http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3');
#time = nc.time('time');
#time = datevec(time);
#''')
#time = matlab.GetVariable("time", "base")
for i, d in enumerate(time):
#d = datetime.datetime(int(time[i][0]), int(time[i][1]), int(time[i][2]),\
# int(time[i][3]), int(time[i][4]), int(time[i][5]))
t = Time(date=d, index=i + 1)
t.save()
for i in range(len(lat)):
n = Node(index=i + 1, lat=lat[i], lon=lon[i])
n.save()
for i in range(len(latc)):
c = Cell(index=i + 1,
lat=latc[i],
lon=lonc[i],
node1=nv[0, i],
node2=nv[1, i],
node3=nv[2, i])
c.save()
lon = nc.variables["lon"][:]
latc = nc.variables["latc"][:]
lonc = nc.variables["lonc"][:]
nv = nc.variables["nv"][:, :]
times = nc.variables["time"]
time = num2date(times[:], units=times.units)
# matlab = win32com.client.Dispatch("matlab.application")
# matlab.Execute('''
# nc = cfdataset('http://www.smast.umassd.edu:8080/thredds/dodsC/fvcom/hindcasts/30yr_gom3');
# time = nc.time('time');
# time = datevec(time);
#''')
# time = matlab.GetVariable("time", "base")
for i, d in enumerate(time):
# d = datetime.datetime(int(time[i][0]), int(time[i][1]), int(time[i][2]),\
# int(time[i][3]), int(time[i][4]), int(time[i][5]))
t = Time(date=d, index=i + 1)
t.save()
for i in range(len(lat)):
n = Node(index=i + 1, lat=lat[i], lon=lon[i])
n.save()
for i in range(len(latc)):
c = Cell(index=i + 1, lat=latc[i], lon=lonc[i], node1=nv[0, i], node2=nv[1, i], node3=nv[2, i])
c.save()
