Help reading nefcdf output from an ocean model

Mon Dec 18 12:12:21 EST 2006

Good Day,

I work with Regional Ocean Model (ROMS) whose outputs are Netcdf files. I

attached the ncdump-c output file for your convienience, as I was unable to open the files with gradsnc (sdfopen - xdfopen).

I currently run the GrADS version 1.9b4 on Linux cluster. Is there a way to create a descriptor file to read this Netcdf file?

Thank you!

Deboarh

> > ----------------------------------------------------------------------------
> > -----------
> > bash-2.05b$ ncdump -c ocean_avg.nc
> > netcdf ocean_avg {
> > dimensions:
> > xi_rho = 43 ;
> > xi_u = 42 ;
> > xi_v = 43 ;
> > xi_psi = 42 ;
> > eta_rho = 82 ;
> > eta_u = 82 ;
> > eta_v = 81 ;
> > eta_psi = 81 ;
> > s_rho = 16 ;
> > s_w = 17 ;
> > tracer = 13 ;
> > boundary = 4 ;
> > ocean_time = UNLIMITED ; // (20 currently)
> > variables:
> > int ntimes ;
> > ntimes:long_name = "number of long time-steps" ;
> > int ndtfast ;
> > ndtfast:long_name = "number of short time-steps" ;
> > double dt ;
> > dt:long_name = "size of long time-steps" ;
> > dt:units = "second" ;
> > double dtfast ;
> > dtfast:long_name = "size of short time-steps" ;
> > dtfast:units = "second" ;
> > double dstart ;
> > dstart:long_name = "time stamp assigned to model
> > initilization" ;
> > dstart:units = "days since 0000-01-01 00:00:00" ;
> > dstart:calendar = "365.25_day" ;
> > int nHIS ;
> > nHIS:long_name = "number of time-steps between history
> > records" ;
> > int nRST ;
> > nRST:long_name = "number of time-steps between restart
> > records" ;
> > nRST:cycle = "only latest two records are maintained" ;
> > int ntsAVG ;
> > ntsAVG:long_name = "starting time-step for accumulation of
> > time- averaged fields" ;
> > int nAVG ;
> > nAVG:long_name = "number of time-steps between
> > time-averaged rec ords" ;
> > int ntsDIA ;
> > ntsDIA:long_name = "starting time-step for accumulation of
> > diagn ostic fields" ;
> > int nDIA ;
> > nDIA:long_name = "number of time-steps between diagnostic
> > record s" ;
> > double Falpha ;
> > Falpha:long_name = "Power-law shape barotropic filter
> > parameter"  ;
> > double Fbeta ;
> > Fbeta:long_name = "Power-law shape barotropic filter
> > parameter" ;
> > double Fgamma ;
> > Fgamma:long_name = "Power-law shape barotropic filter
> > parameter"  ;
> > double tnu2(tracer) ;
> > tnu2:long_name = "Laplacian mixing coefficient for tracers"
> > ;
> > tnu2:units = "meter2 second-1" ;
> > double visc2 ;
> > visc2:long_name = "Laplacian mixing coefficient for
> > momentum" ;
> > visc2:units = "meter2 second-1" ;
> > double Akt_bak(tracer) ;
> > Akt_bak:long_name = "background vertical mixing
> > coefficient for tracers" ;
> > Akt_bak:units = "meter2 second-1" ;
> > double Akv_bak ;
> > Akv_bak:long_name = "background vertical mixing
> > coefficient for momentum" ;
> > Akv_bak:units = "meter2 second-1" ;
> > double rdrg ;
> > rdrg:long_name = "linear drag coefficient" ;
> > rdrg:units = "meter second-1" ;
> > double rdrg2 ;
> > rdrg2:long_name = "quadratic drag coefficient" ;
> > double Zob ;
> > Zob:long_name = "bottom roughness" ;
> > Zob:units = "meter" ;
> > double Zos ;
> > Zos:long_name = "surface roughness" ;
> > Zos:units = "meter" ;
> > double Znudg ;
> > Znudg:long_name = "free-surface nudging/relaxation inverse
> > time scale" ;
> > Znudg:units = "day-1" ;
> > double M2nudg ;
> > M2nudg:long_name = "2D momentum nudging/relaxation inverse
> > time scale" ;
> > M2nudg:units = "day-1" ;
> > double M3nudg ;
> > M3nudg:long_name = "3D momentum nudging/relaxation inverse
> > time scale" ;
> > M3nudg:units = "day-1" ;
> > double Tnudg(tracer) ;
> > Tnudg:long_name = "Tracers nudging/relaxation inverse time
> > scale " ;
> > Tnudg:units = "day-1" ;
> > double rho0 ;
> > rho0:long_name = "mean density used in Boussinesq
> > approximation"  ;
> > rho0:units = "kilogram meter-3" ;
> > double R0 ;
> > R0:long_name = "background density used in linear equation
> > of st ate" ;
> > R0:units = "kilogram meter-3" ;
> > double Tcoef ;
> > Tcoef:long_name = "thermal expansion coefficient" ;
> > Tcoef:units = "Celsius-1" ;
> > double Scoef ;
> > Scoef:long_name = "Saline contraction coefficient" ;
> > Scoef:units = "PSU-1" ;
> > double gamma2 ;
> > gamma2:long_name = "Slipperiness parameter" ;
> > int BioIter ;
> > BioIter:long_name = "number of iterations to achieve
> > convergence " ;
> > double AttSW ;
> > AttSW:long_name = "light attenuation by seawater" ;
> > AttSW:units = "meter-1" ;
> > double AttChl ;
> > AttChl:long_name = "light attenuation by chlorophyll" ;
> > AttChl:units = "meter-2 milligram_Chl-1" ;
> > double PARfrac ;
> > PARfrac:long_name = "photosynthetically available
> > radiation frac tion" ;
> > double Vp0 ;
> > Vp0:long_name = "Eppley temperature-limited growth
> > parameter" ;
> > double I_thNH4 ;
> > I_thNH4:long_name = "radiation threshold for nitrification"
> > ;
> > I_thNH4:units = "watt meter-2" ;
> > double D_p5NH4 ;
> > D_p5NH4:long_name = "half-saturation radiation for
> > nitrification " ;
> > D_p5NH4:units = "watt meter-2" ;
> > double NitriR ;
> > NitriR:long_name = "nitrification rate" ;
> > NitriR:units = "day-1" ;
> > double K_NO3 ;
> > K_NO3:long_name = "inverse half-saturation for
> > phytoplankton NO3  uptake" ;
> > K_NO3:units = "meter3 millimole_N-1" ;
> > double K_NH4 ;
> > K_NH4:long_name = "inverse half-saturation for
> > phytoplankton NH4  uptake" ;
> > K_NH4:units = "meter3 millimole_N-1" ;
> > double K_Phy ;
> > K_Phy:long_name = "zooplankton half-saturation constant
> > for inge stion" ;
> > K_Phy:units = "day-1" ;
> > double Chl2C_m ;
> > Chl2C_m:long_name = "maximum chlorophyll to carbon ratio" ;
> > Chl2C_m:units = "milligram_chl milligram_carbon-1" ;
> > double ChlMin ;
> > ChlMin:long_name = "minimum chlorophyll threshold" ;
> > ChlMin:units = "milligram_chl meter-3" ;
> > double PhyCN ;
> > PhyCN:long_name = "phytoplankton Carbon:Nitrogen ratio" ;
> > PhyCN:units = "mole_C mole_N-1" ;
> > double PhyIP ;
> > PhyIP:long_name = "phytoplankton, NH4 inhibition parameter"
> > ;
> > PhyIP:units = "millimole_N-1" ;
> > double PhyIS ;
> > PhyIS:long_name = "phytoplankton, initial slope of P-I
> > curve" ;
> > PhyIS:units = "milligram_C milligram_Chl-1 watt-1 meter2
> > day-1" ;
> > double PhyMin ;
> > PhyMin:long_name = "minimum phytoplankton threshold" ;
> > PhyMin:units = "millimole_N meter-3" ;
> > double PhyMR ;
> > PhyMR:long_name = "phytoplankton mortality rate" ;
> > PhyMR:units = "day-1" ;
> > double ZooAE_N ;
> > ZooAE_N:long_name = "zooplankton mitrogen assimilation
> > efficienc y" ;
> > double ZooBM ;
> > ZooBM:long_name = "rate for zooplankton basal metabolism" ;
> > ZooBM:units = "day-1" ;
> > double ZooCN ;
> > ZooCN:long_name = "zooplankton Carbon:Nitrogen ratio" ;
> > ZooCN:units = "mole_C mole_N-1" ;
> > double ZooER ;
> > ZooER:long_name = "zooplankton specific excretion rate" ;
> > ZooER:units = "day-1" ;
> > double ZooGR ;
> > ZooGR:long_name = "zooplankton maximum growth rate" ;
> > ZooGR:units = "day-1" ;
> > double ZooMin ;
> > ZooMin:long_name = "minimum zooplankton threshold" ;
> > ZooMin:units = "millimole_N meter-3" ;
> > double ZooMR ;
> > ZooMR:long_name = "zooplankton mortality rate" ;
> > ZooMR:units = "day-1" ;
> > double LDeRRN ;
> > LDeRRN:long_name = "rate of large detritus nitrogen
> > re-mineraliz ation" ;
> > LDeRRN:units = "day-1" ;
> > double LDeRRC ;
> > LDeRRC:long_name = "rate of large detritus carbon
> > re-mineralizat ion" ;
> > LDeRRC:units = "day-1" ;
> > double CoagR ;
> > CoagR:long_name = "coagulation rate" ;
> > CoagR:units = "day-1" ;
> > double SDeRRN ;
> > SDeRRN:long_name = "remineralization rate for small
> > detritus nit rogen" ;
> > SDeRRN:units = "day-1" ;
> > double SDeRRC ;
> > SDeRRC:long_name = "remineralization rate for small
> > detritus car bon" ;
> > SDeRRC:units = "day-1" ;
> > double wPhy ;
> > wPhy:long_name = "vertical sinking velocity for
> > phytoplankton" ;
> > wPhy:units = "meter day-1" ;
> > double wLDet ;
> > wLDet:long_name = "vertical sinking velocity for large
> > detritus"  ;
> > wLDet:units = "meter day-1" ;
> > double wSDet ;
> > wSDet:long_name = "vertical sinking velocity for small
> > detritus"  ;
> > wSDet:units = "meter day-1" ;
> > double pCO2air ;
> > pCO2air:long_name = "partial pressure of CO2 in the air" ;
> > pCO2air:units = "parts per million by volume" ;
> > char spherical ;
> > spherical:long_name = "grid type logical switch" ;
> > spherical:option_T = "spherical" ;
> > spherical:option_F = "Cartesian" ;
> > double xl ;
> > xl:long_name = "domain length in the XI-direction" ;
> > xl:units = "meter" ;
> > double el ;
> > el:long_name = "domain length in the ETA-direction" ;
> > el:units = "meter" ;
> > double theta_s ;
> > theta_s:long_name = "S-coordinate surface control
> > parameter" ;
> > double theta_b ;
> > theta_b:long_name = "S-coordinate bottom control parameter"
> > ;
> > double Tcline ;
> > Tcline:long_name = "S-coordinate surface/bottom layer
> > width" ;
> > Tcline:units = "meter" ;
> > double hc ;
> > hc:long_name = "S-coordinate parameter, critical depth" ;
> > hc:units = "meter" ;
> > double s_rho(s_rho) ;
> > s_rho:long_name = "S-coordinate at RHO-points" ;
> > s_rho:valid_min = -1. ;
> > s_rho:valid_max = 0. ;
> > s_rho:standard_name = "ocean_s_coordinate" ;
> > s_rho:formula_terms = "s: s_rho eta: zeta depth: h a:
> > theta_s b:  theta_b depth_c: hc" ;
> > s_rho:field = "s_rho, scalar" ;
> > double s_w(s_w) ;
> > s_w:long_name = "S-coordinate at W-points" ;
> > s_w:valid_min = -1. ;
> > s_w:valid_max = 0. ;
> > s_w:standard_name = "ocean_s_coordinate" ;
> > s_w:formula_terms = "s: s_w eta: zeta depth: h a: theta_s
> > b: the ta_b depth_c: hc" ;
> > s_w:field = "s_w, scalar" ;
> > double Cs_r(s_rho) ;
> > Cs_r:long_name = "S-coordinate stretching curves at
> > RHO-points" ;
> > Cs_r:valid_min = -1. ;
> > Cs_r:valid_max = 0. ;
> > Cs_r:field = "Cs_r, scalar" ;
> > double Cs_w(s_w) ;
> > Cs_w:long_name = "S-coordinate stretching curves at
> > W-points" ;
> > Cs_w:valid_min = -1. ;
> > Cs_w:valid_max = 0. ;
> > Cs_w:field = "Cs_w, scalar" ;
> > double h(eta_rho, xi_rho) ;
> > h:long_name = "bathymetry at RHO-points" ;
> > h:units = "meter" ;
> > h:coordinates = "x_rho y_rho" ;
> > h:field = "bath, scalar" ;
> > double f(eta_rho, xi_rho) ;
> > f:long_name = "Coriolis parameter at RHO-points" ;
> > f:units = "second-1" ;
> > f:coordinates = "x_rho y_rho" ;
> > f:field = "coriolis, scalar" ;
> > double pm(eta_rho, xi_rho) ;
> > pm:long_name = "curvilinear coordinate metric in XI" ;
> > pm:units = "meter-1" ;
> > pm:coordinates = "x_rho y_rho" ;
> > pm:field = "pm, scalar" ;
> > double pn(eta_rho, xi_rho) ;
> > pn:long_name = "curvilinear coordinate metric in ETA" ;
> > pn:units = "meter-1" ;
> > pn:coordinates = "x_rho y_rho" ;
> > pn:field = "pn, scalar" ;
> > double x_rho(eta_rho, xi_rho) ;
> > x_rho:long_name = "x-locations of RHO-points" ;
> > x_rho:units = "meter" ;
> > x_rho:field = "x_rho, scalar" ;
> > double y_rho(eta_rho, xi_rho) ;
> > y_rho:long_name = "y-locations of RHO-points" ;
> > y_rho:units = "meter" ;
> > y_rho:field = "y_rho, scalar" ;
> > double x_u(eta_u, xi_u) ;
> > x_u:long_name = "x-locations of U-points" ;
> > x_u:units = "meter" ;
> > x_u:field = "x_u, scalar" ;
> > double y_u(eta_u, xi_u) ;
> > y_u:long_name = "y-locations of U-points" ;
> > y_u:units = "meter" ;
> > y_u:field = "y_u, scalar" ;
> > double x_v(eta_v, xi_v) ;
> > x_v:long_name = "x-locations of V-points" ;
> > x_v:units = "meter" ;
> > x_v:field = "x_v, scalar" ;
> > double y_v(eta_v, xi_v) ;
> > y_v:long_name = "y-locations of V-points" ;
> > y_v:units = "meter" ;
> > y_v:field = "y_v, scalar" ;
> > double x_psi(eta_psi, xi_psi) ;
> > x_psi:long_name = "x-locations of PSI-points" ;
> > x_psi:units = "meter" ;
> > x_psi:field = "x_psi, scalar" ;
> > double y_psi(eta_psi, xi_psi) ;
> > y_psi:long_name = "y-locations of PSI-points" ;
> > y_psi:units = "meter" ;
> > y_psi:field = "y_psi, scalar" ;
> > double ocean_time(ocean_time) ;
> > ocean_time:long_name = "averaged time since initialization"
> > ;
> > ocean_time:units = "seconds since 0000-01-01 00:00:00" ;
> > ocean_time:calendar = "365.25_day" ;
> > ocean_time:field = "time, scalar, series" ;
> > float zeta(ocean_time, eta_rho, xi_rho) ;
> > zeta:long_name = "averaged free-surface" ;
> > zeta:units = "meter" ;
> > zeta:time = "ocean_time" ;
> > zeta:coordinates = "x_rho y_rho" ;
> > zeta:field = "free-surface, scalar, series" ;
> > float ubar(ocean_time, eta_u, xi_u) ;
> > ubar:long_name = "averaged vertically integrated
> > u-momentum comp onent" ;
> > ubar:units = "meter second-1" ;
> > ubar:time = "ocean_time" ;
> > ubar:coordinates = "x_u y_u" ;
> > ubar:field = "ubar-velocity, scalar, series" ;
> > float vbar(ocean_time, eta_v, xi_v) ;
> > vbar:long_name = "averaged vertically integrated
> > v-momentum comp onent" ;
> > vbar:units = "meter second-1" ;
> > vbar:time = "ocean_time" ;
> > vbar:coordinates = "x_v y_v" ;
> > vbar:field = "vbar-velocity, scalar, series" ;
> > float u(ocean_time, s_rho, eta_u, xi_u) ;
> > u:long_name = "averaged u-momentum component" ;
> > u:units = "meter second-1" ;
> > u:time = "ocean_time" ;
> > u:coordinates = "x_u y_u" ;
> > u:field = "u-velocity, scalar, series" ;
> > float v(ocean_time, s_rho, eta_v, xi_v) ;
> > v:long_name = "averaged v-momentum component" ;
> > v:units = "meter second-1" ;
> > v:time = "ocean_time" ;
> > v:coordinates = "x_v y_v" ;
> > v:field = "v-velocity, scalar, series" ;
> > float omega(ocean_time, s_w, eta_rho, xi_rho) ;
> > omega:long_name = "averaged S-coordinate vertical momentum
> > compo nent" ;
> > omega:units = "meter3 second-1" ;
> > omega:time = "ocean_time" ;
> > omega:coordinates = "x_rho y_rho" ;
> > omega:field = "omega, scalar, series" ;
> > float temp(ocean_time, s_rho, eta_rho, xi_rho) ;
> > temp:long_name = "averaged potential temperature" ;
> > temp:units = "Celsius" ;
> > temp:time = "ocean_time" ;
> > temp:coordinates = "x_rho y_rho" ;
> > temp:field = "temperature, scalar, series" ;
> > float salt(ocean_time, s_rho, eta_rho, xi_rho) ;
> > salt:long_name = "averaged salinity" ;
> > salt:units = "PSU" ;
> > salt:time = "ocean_time" ;
> > salt:coordinates = "x_rho y_rho" ;
> > salt:field = "salinity, scalar, series" ;
> > float NO3(ocean_time, s_rho, eta_rho, xi_rho) ;
> > NO3:long_name = "averaged nitrate concentration" ;
> > NO3:units = "millimole_N03 meter-3" ;
> > NO3:time = "ocean_time" ;
> > NO3:coordinates = "x_rho y_rho" ;
> > NO3:field = "NO3, scalar, series" ;
> > float NH4(ocean_time, s_rho, eta_rho, xi_rho) ;
> > NH4:long_name = "averaged ammonium concentration" ;
> > NH4:units = "millimole_NH4 meter-3" ;
> > NH4:time = "ocean_time" ;
> > NH4:coordinates = "x_rho y_rho" ;
> > NH4:field = "NH4, scalar, series" ;
> > float chlorophyll(ocean_time, s_rho, eta_rho, xi_rho) ;
> > chlorophyll:long_name = "averaged chlorophyll
> > concentration" ;
> > chlorophyll:units = "milligrams_chlorophyll meter-3" ;
> > chlorophyll:time = "ocean_time" ;
> > chlorophyll:coordinates = "x_rho y_rho" ;
> > chlorophyll:field = "chlorophyll, scalar, series" ;
> > float phytoplankton(ocean_time, s_rho, eta_rho, xi_rho) ;
> > phytoplankton:long_name = "averaged phytoplankton
> > concentration"  ;
> > phytoplankton:units = "millimole_nitrogen meter-3" ;
> > phytoplankton:time = "ocean_time" ;
> > phytoplankton:coordinates = "x_rho y_rho" ;
> > phytoplankton:field = "phytoplankton, scalar, series" ;
> > float zooplankton(ocean_time, s_rho, eta_rho, xi_rho) ;
> > zooplankton:long_name = "averaged zooplankton
> > concentration" ;
> > zooplankton:units = "millimole_nitrogen meter-3" ;
> > zooplankton:time = "ocean_time" ;
> > zooplankton:coordinates = "x_rho y_rho" ;
> > zooplankton:field = "zooplankton, scalar, series" ;
> > float LdetritusN(ocean_time, s_rho, eta_rho, xi_rho) ;
> > LdetritusN:long_name = "averaged large fraction nitrogen
> > detritu s concentration" ;
> > LdetritusN:units = "millimole_nitrogen meter-3" ;
> > LdetritusN:time = "ocean_time" ;
> > LdetritusN:coordinates = "x_rho y_rho" ;
> > LdetritusN:field = "LdetritusN, scalar, series" ;
> > float SdetritusN(ocean_time, s_rho, eta_rho, xi_rho) ;
> > SdetritusN:long_name = "averaged small fraction nitrogen
> > detritu s concentration" ;
> > SdetritusN:units = "millimole_nitrogen meter-3" ;
> > SdetritusN:time = "ocean_time" ;
> > SdetritusN:coordinates = "x_rho y_rho" ;
> > SdetritusN:field = "SdetritusN, scalar, series" ;
> > float LdetritusC(ocean_time, s_rho, eta_rho, xi_rho) ;
> > LdetritusC:long_name = "averaged large fraction carbon
> > detritus concentration" ;
> > LdetritusC:units = "millimole_carbon meter-3" ;
> > LdetritusC:time = "ocean_time" ;
> > LdetritusC:coordinates = "x_rho y_rho" ;
> > LdetritusC:field = "LdetritusC, scalar, series" ;
> > float SdetritusC(ocean_time, s_rho, eta_rho, xi_rho) ;
> > SdetritusC:long_name = "averaged small fraction carbon
> > detritus concentration" ;
> > SdetritusC:units = "millimole_carbon meter-3" ;
> > SdetritusC:time = "ocean_time" ;
> > SdetritusC:coordinates = "x_rho y_rho" ;
> > SdetritusC:field = "SdetritusC, scalar, series" ;
> > float TIC(ocean_time, s_rho, eta_rho, xi_rho) ;
> > TIC:long_name = "averaged total inorganic carbon" ;
> > TIC:units = "millimole_carbon meter-3" ;
> > TIC:time = "ocean_time" ;
> > TIC:coordinates = "x_rho y_rho" ;
> > TIC:field = "TIC, scalar, series" ;
> > float alkalinity(ocean_time, s_rho, eta_rho, xi_rho) ;
> > alkalinity:long_name = "averaged total alkalinity" ;
> > alkalinity:units = "milliequivalents meter-3" ;
> > alkalinity:time = "ocean_time" ;
> > alkalinity:coordinates = "x_rho y_rho" ;
> > alkalinity:field = "alkalinity, scalar, series" ;
> > global attributes:
> > :type = "ROMS/TOMS averages file" ;
> > :Conventions = "CF-1.0" ;
> > :title = "ROMS/TOMS 2.2 - Wind-Driven
> > Upwelling/Downwelling over  a Periodic Channel" ;
> > :var_info = "External/varinfo.dat" ;
> > :rst_file = "ocean_rst.nc" ;
> > :avg_file = "ocean_avg.nc" ;
> > :dia_file = "ocean_dia.nc" ;
> > :bpar_file = "External/bioFasham.in" ;
> > :os = "Linux" ;
> > :cpu = "x86_64" ;
> > :compiler_system = "pgi" ;
> > :compiler_command = "pgf90" ;
> > :compiler_flags = " -u -Bstatic -fastsse -Mipa=fast" ;
> > :history = "ROMS/TOMS, Version 2.2, Tuesday - November 14,
> > 2006 -  2:07:36 PM" ;
> > :tiling = "001x001" ;
> > :CPP_options = " ANA_BIOLOGY, ANA_BPFLUX, ANA_BSFLUX,
> > ANA_BTFLUX , ANA_GRID, ANA_INITIAL, ANA_SMFLUX,
> > ANA_SPFLUX, ANA_SRFLUX, ANA_SSFLUX, ANA_STF LUX, ANA_VMIX,
> > ASSUMED_SHAPE, AVERAGES, BIO_FASHAM, BIO_SEDIMENT, CARBON,
> > DENITR IFICATION, DIAGNOSTICS_BIO, DIAGNOSTICS_TS,
> > DIAGNOSTICS_UV, DJ_GRADPS, DOUBLE_PR ECISION, EW_PERIODIC,
> > MIX_S_TS, MIX_S_UV, NONLINEAR, POWER_LAW, PROFILE, !RST_SI
> > NGLE, SALINITY, SOLVE3D, SPLINES, !TALK_PROGNOSTIC,
> > TS_U3HADVECTION, TS_C4VADVEC TION, TS_DIF2, UPWELLING,
> > UV_ADV, UV_COR, UV_U3HADVECTION, UV_LDRAG, UV_VIS2, VA
> > R_RHO_2D," ;
> > data:
> > s_rho = -0.96875, -0.90625, -0.84375, -0.78125, -0.71875, -0.65625,
> > -0.59375, -0.53125, -0.46875, -0.40625, -0.34375, -0.28125, -0.21875,
> > -0.15625, -0.09375, -0.03125 ;
> > s_w = -1, -0.9375, -0.875, -0.8125, -0.75, -0.6875, -0.625, -0.5625, -0.5,
> > -0.4375, -0.375, -0.3125, -0.25, -0.1875, -0.125, -0.0625, 0 ;
> > ocean_time = 10800, 32400, 54000, 75600, 97200, 118800, 140400, 162000,
> > 183600, 205200, 226800, 248400, 270000, 291600, 313200, 334800, 356400,
> > 378000, 399600, 421200 ;
> > }
>

Summary of resource utilization
-------------------------------
 CPU time:        0.020 sec
 Overhead CPU:    0.000 sec
 CPU model:         Ultra-4 (1536M)
 Job origin:      herceg at MARINE.RUTGERS.EDU