Notebook Conversion between .dat and .csv file.

In [1]:

import matplotlib.pyplot as plt
import numpy as np
%matplotlib inline

from itertools import izip
import csv

In [2]:

def radians(degrees):
    rad = np.pi/180.*degrees
    return rad

In [3]:

l = range (98)
lnew=l[0::3]

nstn = 32
# pathname for data - all of the tide runs are stored in this directory
fp = open('/ocean/imachuca/MEOPAR/NEMO-forcing/fromMike/vigf7_pres_ele_harm_obs2.dat')

In [4]:

name=['Neah Bay','Port Renfrew','Port Angeles', 'Victoria Harbour','Port Townsend', 'Bangor','Seattle',
       'Tacoma', 'Cherry Point', 'Friday Harbor', 'Hanbury Point', 'Sidney', 'Fulford Harbour',
       'Patos Island', 'Tsawwassen', 'Point Atkinson', 'Winchelsea Islands', 'Little River',
       'Twin Islets', 'Campbell River', 'Seymour Narrows', 'Owen Bay', 'Big Bay', 'Chatham Point',
       'Yorke Island', 'Alert Bay', 'Port Hardy', 'Montagu Point', 'Siwash Bay', 'Winter Harbour',
       'Bella Bella', 'Tofino']
lons=np.zeros(nstn); lats=np.zeros(nstn)
K1amp = np.zeros(nstn); Q1amp = np.zeros(nstn); O1amp = np.zeros(nstn); P1amp = np.zeros(nstn)
K1phase = np.zeros(nstn); Q1phase = np.zeros(nstn); O1phase = np.zeros(nstn); P1phase = np.zeros(nstn)
M2amp = np.zeros(nstn); N2amp = np.zeros(nstn); S2amp = np.zeros(nstn); K2amp = np.zeros(nstn)
M2phase = np.zeros(nstn); N2phase = np.zeros(nstn); S2phase = np.zeros(nstn); K2phase = np.zeros(nstn)
for i in range(nstn):
      
    #Details
    line = fp.readline()
    
    lons[i] = round(float(line[1:10])-360,3)
    lats[i] = round(float(line[11:18]),3)
    
    #Diurnals
    line = fp.readline()
    values = line.split(' ')
    diurnals = [round(float(e),2) for e in values if (e != '' and e != '\n')]
 
    K1amp[i] = diurnals[6]
    Q1amp[i] = diurnals[0]
    O1amp[i] = diurnals[2]
    P1amp[i] = diurnals[4]
    K1phase[i] = diurnals[7]
    Q1phase[i] = diurnals[1]
    O1phase[i] = diurnals[3]
    P1phase[i] = diurnals[5]   
    #Semi-diurnals
    line = fp.readline()
    values = line.split(' ')
    semis = [float(e) for e in values if (e != '' and e != '\n')]
    M2amp[i] = semis[2]
    N2amp[i] = semis[0]
    S2amp[i] = semis[4]
    K2amp[i] = semis[6]
    M2phase[i] = semis[3]
    N2phase[i] = semis[1]
    S2phase[i] = semis[5]
    K2phase[i] = semis[7]

In [5]:

print K1amp

[ 49.7   45.33  66.7   63.6   76.4   83.5   83.1   83.8   81.7   75.8   75.3
  76.6   75.3   79.    83.4   86.2   87.5   90.19  90.37  84.58  69.19
  67.76  83.27  65.41  55.75  51.63  49.9   52.5   52.6   43.1   46.5
  38.81]

In [6]:

outfile ='other_constituents.csv'

with open(outfile, 'wb') as csvfile:
    writer = csv.writer(csvfile, delimiter=',')
    writer.writerow([
            'Site', 'Lat', 'Lon',
            'O1 amp', 'O1 phase (deg UT)', 'P1 amp', 'P1 phase (deg UT)', 'Q1 amp', 'Q1 phase (deg UT)', 
            'S2 amp', 'S2 phase (deg UT)', 'N2 amp', 'N2 phase (deg UT)', 'K2 amp', 'K2 phase (deg UT)' 
        ])
    for stn in range(nstn):
        location = name[stn]
        writer.writerow([location, lats[stn], lons[stn],
                        O1amp[stn], O1phase[stn], P1amp[stn], P1phase[stn],Q1amp[stn], Q1phase[stn],
                        S2amp[stn], S2phase[stn], N2amp[stn], N2phase[stn],K2amp[stn], K2phase[stn]])

In [6]: