def write2db(dbfile): ''' The function below aims at writing the resulting river inflow data and tide data into the database file ''' #Update and write the resulting data into the database file create = True #dbfile = 'GalvestonObs.db' #dire1=os.path.dirname(os.path.abspath(__file__)) dire1=os.getcwd() #dire1=os.path.dirname(dire2) folder = dire1+'/DATA/' #pdb.set_trace() ncfiles = [folder+'USGS_Rivers.nc',folder+'TCOONTide.nc'] if create: print 'Creating database: %s'%dbfile netcdfio.createObsDB(dbfile) for nc in ncfiles: print 'Inserting metadata from: %s'%nc netcdfio.netcdfObs2DB(nc,dbfile) print 'Done.' if os.path.isfile(folder+dbfile): os.remove(folder+dbfile) copy_file(dbfile, folder)
# Set the date range slightly larger just in case tstart = datetime.strptime(timestart,'%Y%m%d')-timedelta(days=30) tend = datetime.strptime(timeend,'%Y%m%d')+timedelta(days=30) latlon = [lon0-dx,lon0+dx,lat0-dx,lat0+dx] data = noaa.noaaish2nc(latlon,[tstart.year,tend.year],localdir,ncfile,shpfile) ### # Step 2: Create a database and insert the netcdf metadata into it ### print 'Creating database: %s'%dbfile if os.path.exists(dbfile): print 'Overwriting old file...' os.unlink(dbfile) netcdfio.createObsDB(dbfile) ncfiles=[ncfile] for nc in ncfiles: print 'Inserting metadata from: %s'%nc netcdfio.netcdfObs2DB(nc,dbfile) ### # Step 3: Create a SUNTANS format meteorological input file ### [coords, output, nctime] = interpWeatherStations(latlon,timestart,timeend,dt*24.0,utmzone,dbfile,maxgap=100,showplot=False) # Write to the suntans meteorological netcdf format write2NC(metncfile,coords,output,nctime)