In [1]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

import numpy as np
import netCDF4
import os
import sys 
import subprocess
import pyroms
from pyroms_toolbox import jday2date
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt 
from datetime import datetime


    

In [5]:

    

def readtext(filename):
    f = open(filename, 'r')
    # Eat first two lines
    f.readline()
    f.readline()

    ttime = []
    energy = []
    msl = []
    mass = []
    for line in f:
        a, b, c, d, e, f, g, h = re.split(',\s+', line)
        print(a, b, c, d)
        print(e, f, g, h)


    

In [6]:

    

readtext('ocean.stats.zero')


    

    




     0 0.000 0 En 4.0021589354361712E-04
CFL  0.10056 SL  4.4729E-03 Mass 7.98435E+05 Me  0.00E+00

    50 10.000 0 En 3.9369949762223531E-04
CFL  0.09720 SL  4.4731E-03 Mass 7.98436E+05 Me  1.65E-07

   100 20.000 0 En 3.8809902764580485E-04
CFL  0.09421 SL  4.3803E-03 Mass 7.98362E+05 Me -9.23E-05

   150 30.000 0 En 3.8314900914778453E-04
CFL  0.09020 SL  4.3774E-03 Mass 7.98360E+05 Me -2.94E-06

   200 40.000 0 En 3.7877958019264505E-04
CFL  0.08877 SL  4.3714E-03 Mass 7.98355E+05 Me -6.01E-06

   250 50.000 0 En 3.7462890069380199E-04
CFL  0.08700 SL  4.3659E-03 Mass 7.98350E+05 Me -5.46E-06

   300 60.000 0 En 3.7076775501069110E-04
CFL  0.08575 SL  4.3616E-03 Mass 7.98347E+05 Me -4.27E-06

   350 70.000 0 En 3.6713172170675418E-04
CFL  0.08448 SL  4.3351E-03 Mass 7.98326E+05 Me -2.64E-05

   400 80.000 0 En 3.6389093860060274E-04
CFL  0.08375 SL  4.0674E-03 Mass 7.98113E+05 Me -2.67E-04

   450 90.000 0 En 3.3805515545456814E-04
CFL  0.09962 SL  2.2727E-03 Mass 7.96687E+05 Me -1.79E-03

   500 100.000 0 En 3.1226412117449186E-04
CFL  0.14147 SL -1.0956E-03 Mass 7.94009E+05 Me -3.37E-03

   550 110.000 0 En 3.1399898595970875E-04
CFL  0.12356 SL -2.7594E-03 Mass 7.92687E+05 Me -1.67E-03

   600 120.000 0 En 2.6440728043311573E-04
CFL  0.10502 SL -3.3585E-03 Mass 7.92210E+05 Me -6.01E-04

   650 130.000 0 En 2.2203818902851661E-04
CFL  0.08741 SL -3.3730E-03 Mass 7.92199E+05 Me -1.45E-05

   700 140.000 0 En 1.8827461262810444E-04
CFL  0.06840 SL -2.8657E-03 Mass 7.92602E+05 Me  5.09E-04

   750 150.000 0 En 1.4284027464004726E-04
CFL  0.05744 SL -1.8848E-03 Mass 7.93382E+05 Me  9.83E-04

   800 160.000 0 En 1.1392134223923268E-04
CFL  0.05056 SL -1.3134E-03 Mass 7.93836E+05 Me  5.72E-04

   850 170.000 0 En 9.8748082221344259E-05
CFL  0.04077 SL -1.0410E-03 Mass 7.94053E+05 Me  2.73E-04

   900 180.000 0 En 8.7986255009010525E-05
CFL  0.03881 SL -9.3146E-04 Mass 7.94140E+05 Me  1.10E-04

   950 190.000 0 En 7.9578793244685476E-05
CFL  0.03214 SL -9.6554E-04 Mass 7.94113E+05 Me -3.41E-05

  1000 200.000 0 En 7.2443298923439316E-05
CFL  0.03055 SL -9.9291E-04 Mass 7.94091E+05 Me -2.74E-05

  1050 210.000 0 En 6.5481482746467060E-05
CFL  0.02678 SL -8.8813E-04 Mass 7.94174E+05 Me  1.05E-04

  1100 220.000 0 En 5.9345105393542501E-05
CFL  0.02419 SL -7.6736E-04 Mass 7.94270E+05 Me  1.21E-04

  1150 230.000 0 En 5.4511281640941031E-05
CFL  0.02294 SL -7.0918E-04 Mass 7.94316E+05 Me  5.82E-05

  1200 240.000 0 En 5.0074496336000540E-05
CFL  0.02441 SL -6.2274E-04 Mass 7.94385E+05 Me  8.65E-05

  1250 250.000 0 En 4.5989073227757276E-05
CFL  0.02022 SL -5.1711E-04 Mass 7.94469E+05 Me  1.06E-04

  1300 260.000 0 En 4.2533425643503937E-05
CFL  0.02017 SL -4.6668E-04 Mass 7.94509E+05 Me  5.05E-05

  1350 270.000 0 En 3.9512203548550537E-05
CFL  0.01732 SL -4.1026E-04 Mass 7.94554E+05 Me  5.64E-05

  1400 280.000 0 En 3.6812167356549476E-05
CFL  0.01789 SL -3.2087E-04 Mass 7.94625E+05 Me  8.94E-05

  1450 290.000 0 En 3.4548710121953125E-05
CFL  0.01610 SL -2.3054E-04 Mass 7.94697E+05 Me  9.04E-05

  1500 300.000 0 En 3.2956639903005890E-05
CFL  0.01501 SL -1.0576E-04 Mass 7.94796E+05 Me  1.25E-04

  1550 310.000 0 En 3.1994060809326635E-05
CFL  0.01393 SL  4.2790E-05 Mass 7.94914E+05 Me  1.49E-04

  1600 320.000 0 En 3.1549446852363636E-05
CFL  0.01349 SL  1.9492E-04 Mass 7.95035E+05 Me  1.52E-04

  1650 330.000 0 En 3.1755170602671174E-05
CFL  0.01436 SL  3.6952E-04 Mass 7.95174E+05 Me  1.75E-04

  1700 340.000 0 En 3.2093510616961276E-05
CFL  0.01589 SL  5.3071E-04 Mass 7.95302E+05 Me  1.61E-04

  1750 350.000 0 En 3.2091547813122876E-05
CFL  0.01489 SL  6.4655E-04 Mass 7.95394E+05 Me  1.16E-04

  1800 360.000 0 En 3.1606558199208056E-05
CFL  0.01404 SL  7.1507E-04 Mass 7.95448E+05 Me  6.85E-05

  1850 370.000 0 En 3.0554766806447751E-05
CFL  0.01402 SL  7.1300E-04 Mass 7.95447E+05 Me -2.08E-06

  1900 380.000 0 En 2.9187735703138164E-05
CFL  0.01350 SL  6.4836E-04 Mass 7.95395E+05 Me -6.46E-05

  1950 390.000 0 En 2.7786813751464963E-05
CFL  0.01333 SL  5.4203E-04 Mass 7.95311E+05 Me -1.06E-04

  2000 400.000 0 En 2.6573654950403144E-05
CFL  0.01294 SL  4.0288E-04 Mass 7.95200E+05 Me -1.39E-04

  2050 410.000 0 En 2.5550294648119235E-05
CFL  0.01199 SL  2.6782E-04 Mass 7.95093E+05 Me -1.35E-04

  2100 420.000 0 En 2.4724190264563750E-05
CFL  0.01282 SL  1.5926E-04 Mass 7.95007E+05 Me -1.09E-04

  2150 430.000 0 En 2.4061665167000865E-05
CFL  0.01146 SL  8.3814E-05 Mass 7.94947E+05 Me -7.54E-05

  2200 440.000 0 En 2.3523574679794141E-05
CFL  0.01008 SL  4.9214E-05 Mass 7.94919E+05 Me -3.46E-05

  2250 450.000 0 En 2.3054393746725471E-05
CFL  0.00963 SL  4.7057E-05 Mass 7.94917E+05 Me -2.16E-06

  2300 460.000 0 En 2.2583815115840778E-05
CFL  0.00923 SL  7.1449E-05 Mass 7.94937E+05 Me  2.44E-05

  2350 470.000 0 En 2.2127483153947108E-05
CFL  0.00946 SL  1.0987E-04 Mass 7.94967E+05 Me  3.84E-05

  2400 480.000 0 En 2.1703752860769137E-05
CFL  0.00957 SL  1.4276E-04 Mass 7.94993E+05 Me  3.29E-05

  2450 490.000 0 En 2.1348748177147352E-05
CFL  0.00929 SL  1.6429E-04 Mass 7.95011E+05 Me  2.15E-05

  2500 500.000 0 En 2.1054086903714983E-05
CFL  0.00902 SL  1.7046E-04 Mass 7.95015E+05 Me  6.17E-06

  2550 510.000 0 En 2.0783449114589276E-05
CFL  0.00848 SL  1.6507E-04 Mass 7.95011E+05 Me -5.39E-06

  2600 520.000 0 En 2.0502232814295275E-05
CFL  0.00798 SL  1.5905E-04 Mass 7.95006E+05 Me -6.02E-06

  2650 530.000 0 En 2.0179147438002688E-05
CFL  0.00834 SL  1.5313E-04 Mass 7.95002E+05 Me -5.91E-06

  2700 540.000 0 En 1.9808846711272145E-05
CFL  0.00817 SL  1.4794E-04 Mass 7.94998E+05 Me -5.19E-06

  2750 550.000 0 En 1.9415582852438950E-05
CFL  0.00821 SL  1.4293E-04 Mass 7.94994E+05 Me -5.00E-06

  2800 560.000 0 En 1.9023794962423388E-05
CFL  0.00804 SL  1.3281E-04 Mass 7.94986E+05 Me -1.01E-05

  2850 570.000 0 En 1.8647330512327328E-05
CFL  0.00777 SL  1.1682E-04 Mass 7.94973E+05 Me -1.60E-05

  2900 580.000 0 En 1.8294049979688880E-05
CFL  0.00732 SL  9.3697E-05 Mass 7.94954E+05 Me -2.31E-05

  2950 590.000 0 En 1.7978340948454189E-05
CFL  0.00686 SL  6.2790E-05 Mass 7.94930E+05 Me -3.09E-05

  3000 600.000 0 En 1.7716230767388042E-05
CFL  0.00712 SL  2.8566E-05 Mass 7.94903E+05 Me -3.42E-05

  3050 610.000 0 En 1.7513256324977912E-05
CFL  0.00725 SL -6.7701E-06 Mass 7.94875E+05 Me -3.53E-05

  3100 620.000 0 En 1.7362038207597260E-05
CFL  0.00727 SL -4.3220E-05 Mass 7.94846E+05 Me -3.65E-05

  3150 630.000 0 En 1.7252656508364880E-05
CFL  0.00709 SL -7.9624E-05 Mass 7.94817E+05 Me -3.64E-05

  3200 640.000 0 En 1.7188392441757075E-05
CFL  0.00689 SL -1.1478E-04 Mass 7.94789E+05 Me -3.52E-05

  3250 650.000 0 En 1.7160098879821911E-05
CFL  0.00709 SL -1.4486E-04 Mass 7.94765E+05 Me -3.01E-05

  3300 660.000 0 En 1.7140318600849792E-05
CFL  0.00724 SL -1.6627E-04 Mass 7.94748E+05 Me -2.14E-05

  3350 670.000 0 En 1.7119458424626557E-05
CFL  0.00732 SL -1.7767E-04 Mass 7.94739E+05 Me -1.14E-05

  3400 680.000 0 En 1.7104154689616464E-05
CFL  0.00730 SL -1.7817E-04 Mass 7.94738E+05 Me -5.06E-07

  3450 690.000 0 En 1.7116549130004496E-05
CFL  0.00721 SL -1.6679E-04 Mass 7.94747E+05 Me  1.14E-05

  3500 700.000 0 En 1.7178900307569013E-05
CFL  0.00707 SL -1.4167E-04 Mass 7.94767E+05 Me  2.51E-05

  3550 710.000 0 En 1.7313694157527220E-05
CFL  0.00692 SL -1.0203E-04 Mass 7.94799E+05 Me  3.96E-05

  3600 720.000 0 En 1.7526099497961747E-05
CFL  0.00680 SL -5.1057E-05 Mass 7.94839E+05 Me  5.10E-05

  3650 730.000 0 En 1.7797037665893071E-05
CFL  0.00673 SL  6.1512E-06 Mass 7.94885E+05 Me  5.72E-05

  3700 740.000 0 En 1.8089813228149357E-05
CFL  0.00672 SL  6.3360E-05 Mass 7.94930E+05 Me  5.72E-05

  3750 750.000 0 En 1.8355278710816941E-05
CFL  0.00676 SL  1.1375E-04 Mass 7.94970E+05 Me  5.04E-05

  3800 760.000 0 En 1.8568204559769947E-05
CFL  0.00682 SL  1.5234E-04 Mass 7.95001E+05 Me  3.86E-05

  3850 770.000 0 En 1.8714060447264938E-05
CFL  0.00688 SL  1.7533E-04 Mass 7.95019E+05 Me  2.30E-05

  3900 780.000 0 En 1.8799083992331018E-05
CFL  0.00691 SL  1.8229E-04 Mass 7.95025E+05 Me  6.96E-06

  3950 790.000 0 En 1.8853197772388260E-05
CFL  0.00690 SL  1.7667E-04 Mass 7.95020E+05 Me -5.61E-06

  4000 800.000 0 En 1.8907208613653862E-05
CFL  0.00688 SL  1.6237E-04 Mass 7.95009E+05 Me -1.43E-05

  4050 810.000 0 En 1.8993339782527482E-05
CFL  0.00685 SL  1.4526E-04 Mass 7.94995E+05 Me -1.71E-05

  4100 820.000 0 En 1.9134036563804566E-05
CFL  0.00683 SL  1.2980E-04 Mass 7.94983E+05 Me -1.55E-05

  4150 830.000 0 En 1.9339627104464263E-05
CFL  0.00682 SL  1.1743E-04 Mass 7.94973E+05 Me -1.24E-05

  4200 840.000 0 En 1.9623178365458754E-05
CFL  0.00684 SL  1.1312E-04 Mass 7.94970E+05 Me -4.31E-06

  4250 850.000 0 En 1.9991008109530478E-05
CFL  0.00688 SL  1.1791E-04 Mass 7.94974E+05 Me  4.79E-06

  4300 860.000 0 En 2.0431473295828010E-05
CFL  0.00693 SL  1.2360E-04 Mass 7.94978E+05 Me  5.69E-06

  4350 870.000 0 En 2.0896530686557789E-05
CFL  0.00698 SL  1.3411E-04 Mass 7.94987E+05 Me  1.05E-05

  4400 880.000 0 En 2.1339041695633583E-05
CFL  0.00703 SL  1.6213E-04 Mass 7.95009E+05 Me  2.80E-05

  4450 890.000 0 En 2.1811940208091804E-05
CFL  0.00707 SL  1.9056E-04 Mass 7.95031E+05 Me  2.84E-05

  4500 900.000 0 En 2.2310186294927576E-05
CFL  0.00729 SL  2.0353E-04 Mass 7.95042E+05 Me  1.30E-05