netcdf ECB_ROMS/ECB_ROMS_best.ncd {
  dimensions:
    xi_rho = 100;
    xi_u = 99;
    xi_v = 100;
    xi_psi = 99;
    eta_rho = 150;
    eta_u = 150;
    eta_v = 149;
    eta_psi = 149;
    N = 20;
    s_rho = 20;
    s_w = 21;
    tracer = 21;
    boundary = 4;
    ocean_time = UNLIMITED;   // (13 currently)
    time = 1981;
    time1 = 3002;
    time2 = 3369;
    time3 = 6362;
  variables:
    int ntimes;
      :long_name = "number of long time-steps";

    int ndtfast;
      :long_name = "number of short time-steps";

    double dt;
      :units = "second";
      :long_name = "size of long time-steps";

    double dtfast;
      :units = "second";
      :long_name = "size of short time-steps";

    double dstart;
      :units = "days since 2009-01-01 00:00:00";
      :long_name = "time stamp assigned to model initilization";
      :calendar = "gregorian_proleptic";

    int nHIS;
      :long_name = "number of time-steps between history records";

    int ndefHIS;
      :long_name = "number of time-steps between the creation of history files";

    int nRST;
      :long_name = "number of time-steps between restart records";

    int ntsAVG;
      :long_name = "starting time-step for accumulation of time-averaged fields";

    int nAVG;
      :long_name = "number of time-steps between time-averaged records";

    int ndefAVG;
      :long_name = "number of time-steps between the creation of average files";

    int nSTA;
      :long_name = "number of time-steps between stations records";

    double Falpha;
      :long_name = "Power-law shape barotropic filter parameter";

    double Fbeta;
      :long_name = "Power-law shape barotropic filter parameter";

    double Fgamma;
      :long_name = "Power-law shape barotropic filter parameter";

    double Akt_bak(tracer=21);
      :units = "meter2 second-1";
      :long_name = "background vertical mixing coefficient for tracers";

    double Akv_bak;
      :units = "meter2 second-1";
      :long_name = "background vertical mixing coefficient for momentum";

    double Akk_bak;
      :units = "meter2 second-1";
      :long_name = "background vertical mixing coefficient for turbulent energy";

    double Akp_bak;
      :units = "meter2 second-1";
      :long_name = "background vertical mixing coefficient for length scale";

    double rdrg;
      :units = "meter second-1";
      :long_name = "linear drag coefficient";

    double rdrg2;
      :long_name = "quadratic drag coefficient";

    double Zob;
      :units = "meter";
      :long_name = "bottom roughness";

    double Zos;
      :units = "meter";
      :long_name = "surface roughness";

    double gls_p;
      :long_name = "stability exponent";

    double gls_m;
      :long_name = "turbulent kinetic energy exponent";

    double gls_n;
      :long_name = "turbulent length scale exponent";

    double gls_cmu0;
      :long_name = "stability coefficient";

    double gls_c1;
      :long_name = "shear production coefficient";

    double gls_c2;
      :long_name = "dissipation coefficient";

    double gls_c3m;
      :long_name = "buoyancy production coefficient (minus)";

    double gls_c3p;
      :long_name = "buoyancy production coefficient (plus)";

    double gls_sigk;
      :long_name = "constant Schmidt number for TKE";

    double gls_sigp;
      :long_name = "constant Schmidt number for PSI";

    double gls_Kmin;
      :long_name = "minimum value of specific turbulent kinetic energy";

    double gls_Pmin;
      :long_name = "minimum Value of dissipation";

    double Charnok_alpha;
      :long_name = "Charnok factor for surface roughness";

    double Zos_hsig_alpha;
      :long_name = "wave amplitude factor for surface roughness";

    double sz_alpha;
      :long_name = "surface flux from wave dissipation";

    double CrgBan_cw;
      :long_name = "surface flux due to Craig and Banner wave breaking";

    double Znudg;
      :units = "day-1";
      :long_name = "free-surface nudging/relaxation inverse time scale";

    double M2nudg;
      :units = "day-1";
      :long_name = "2D momentum nudging/relaxation inverse time scale";

    double M3nudg;
      :units = "day-1";
      :long_name = "3D momentum nudging/relaxation inverse time scale";

    double Tnudg(tracer=21);
      :units = "day-1";
      :long_name = "Tracers nudging/relaxation inverse time scale";

    double FSobc_in(boundary=4);
      :units = "second-1";
      :long_name = "free-surface inflow, nudging inverse time scale";

    double FSobc_out(boundary=4);
      :units = "second-1";
      :long_name = "free-surface outflow, nudging inverse time scale";

    double M2obc_in(boundary=4);
      :units = "second-1";
      :long_name = "2D momentum inflow, nudging inverse time scale";

    double M2obc_out(boundary=4);
      :units = "second-1";
      :long_name = "2D momentum outflow, nudging inverse time scale";

    double Tobc_in(boundary=4, tracer=21);
      :units = "second-1";
      :long_name = "tracers inflow, nudging inverse time scale";

    double Tobc_out(boundary=4, tracer=21);
      :units = "second-1";
      :long_name = "tracers outflow, nudging inverse time scale";

    double M3obc_in(boundary=4);
      :units = "second-1";
      :long_name = "3D momentum inflow, nudging inverse time scale";

    double M3obc_out(boundary=4);
      :units = "second-1";
      :long_name = "3D momentum outflow, nudging inverse time scale";

    double rho0;
      :units = "kilogram meter-3";
      :long_name = "mean density used in Boussinesq approximation";

    double gamma2;
      :long_name = "slipperiness parameter";

    int LuvSrc;
      :long_name = "momentum point sources and sink activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LwSrc;
      :long_name = "mass point sources and sink activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LtracerSrc(tracer=21);
      :long_name = "tracer point sources and sink activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LsshCLM;
      :long_name = "sea surface height climatology processing switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int Lm2CLM;
      :long_name = "2D momentum climatology processing switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int Lm3CLM;
      :long_name = "3D momentum climatology processing switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LtracerCLM(tracer=21);
      :long_name = "tracer climatology processing switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LnudgeM2CLM;
      :long_name = "2D momentum climatology nudging activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LnudgeM3CLM;
      :long_name = "3D momentum climatology nudging activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int LnudgeTCLM(tracer=21);
      :long_name = "tracer climatology nudging activation switch";
      :flag_values = 0, 1; // int
      :flag_meanings = ".FALSE. .TRUE.";

    int BioIter;
      :long_name = "number of iterations to achieve convergence";

    double AttSW;
      :units = "meter-1";
      :long_name = "light attenuation by seawater";

    double AttChl;
      :units = "meter-2 milligram_Chl-1";
      :long_name = "light attenuation by chlorophyll";

    double PARfrac;
      :long_name = "photosynthetically available radiation fraction";

    double Vp0;
      :long_name = "Phytoplankton growth rate for temperature < 20C";

    double I_thNH4;
      :units = "watt meter-2";
      :long_name = "radiation threshold for nitrification";

    double D_p5NH4;
      :units = "watt meter-2";
      :long_name = "half-saturation radiation for nitrification";

    double NitriR;
      :units = "day-1";
      :long_name = "nitrification rate";

    double K_NO3;
      :units = "meter3 millimole_N-1";
      :long_name = "inverse half-saturation for phytoplankton NO3 uptake";

    double K_NH4;
      :units = "meter3 millimole_N-1";
      :long_name = "inverse half-saturation for phytoplankton NH4 uptake";

    double K_Phy;
      :units = "day-1";
      :long_name = "zooplankton half-saturation constant for ingestion";

    double Chl2C_m;
      :units = "milligram_chl milligram_carbon-1";
      :long_name = "maximum chlorophyll to carbon ratio";

    double ChlMin;
      :units = "milligram_chl meter-3";
      :long_name = "minimum chlorophyll threshold";

    double PhyCN;
      :units = "mole_C mole_N-1";
      :long_name = "phytoplankton Carbon:Nitrogen ratio";

    double PhyIP;
      :units = "millimole_N-1";
      :long_name = "phytoplankton, NH4 inhibition parameter";

    double PhyIS;
      :units = "milligram_C milligram_Chl-1 watt-1 meter2 day-1";
      :long_name = "phytoplankton, initial slope of P-I curve";

    double PhyMin;
      :units = "millimole_N meter-3";
      :long_name = "minimum phytoplankton threshold";

    double PhyMR;
      :units = "day-1";
      :long_name = "phytoplankton mortality rate";

    double ZooAE_N;
      :long_name = "zooplankton mitrogen assimilation efficiency";

    double ZooBM;
      :units = "day-1";
      :long_name = "rate for zooplankton basal metabolism";

    double ZooCN;
      :units = "mole_C mole_N-1";
      :long_name = "zooplankton Carbon:Nitrogen ratio";

    double ZooER;
      :units = "day-1";
      :long_name = "zooplankton specific excretion rate";

    double ZooGR;
      :units = "day-1";
      :long_name = "zooplankton maximum growth rate at 0C";

    double ZooMin;
      :units = "millimole_N meter-3";
      :long_name = "minimum zooplankton threshold";

    double ZooMR;
      :units = "day-1";
      :long_name = "zooplankton mortality rate";

    double LDeRRN;
      :units = "day-1";
      :long_name = "rate of large detritus nitrogen re-mineralization";

    double LDeRRC;
      :units = "day-1";
      :long_name = "rate of large detritus carbon re-mineralization";

    double CoagR;
      :units = "day-1";
      :long_name = "coagulation rate";

    double SDeRRN;
      :units = "day-1";
      :long_name = "remineralization rate for small detritus nitrogen";

    double SDeRRC;
      :units = "day-1";
      :long_name = "remineralization rate for small detritus carbon";

    double wPhy;
      :units = "meter day-1";
      :long_name = "vertical sinking velocity for phytoplankton";

    double wLDet;
      :units = "meter day-1";
      :long_name = "vertical sinking velocity for large detritus";

    double wSDet;
      :units = "meter day-1";
      :long_name = "vertical sinking velocity for small detritus";

    double wISS1;
      :units = "meter day-1";
      :long_name = "vertical sinking velocity for Inorg. Susp. Solids";

    double wISS2;
      :units = "meter day-1";
      :long_name = "vertical sinking velocity for Inorg. Susp. Solids";

    double tau_c;
      :units = "Pa";
      :long_name = "critical stress for Inorg. Susp. Solids";

    double Merosion;
      :units = "g m-2 day-1 Pa-1";
      :long_name = "vertical erosion rate for Inorg. Susp. Solids";

    double rkd1;
      :units = "meter-1";
      :long_name = "light attenuation by low salinity seawater in Chespeake Bay";

    double rkd2;
      :units = "meter-1";
      :long_name = "light attenuation by high salinity seawater in Chespeake Bay";

    double rkdChl1;
      :units = "meter^2 milligram_Chl-1";
      :long_name = "light attenuation by chlorophyll at low salinity seawater in Chespeake Bay";

    double rkdChl2;
      :units = "meter^2 milligram_Chl-1";
      :long_name = "light attenuation by chlorophyll at high salinity seawater in Chespeake Bay";

    double rkdTSS1;
      :units = "meter^2 g-1";
      :long_name = "light attenuation by TSS at low salinity seawater in Chespeake Bay";

    double rkdTSS2;
      :units = "meter^2 g-1";
      :long_name = "light attenuation by TSS at high salinity seawater in Chespeake Bay";

    double rkdS1;
      :units = "meter-1 PSU-1";
      :long_name = "light attenuation by CDOM (use salinity instead) at low salinity seawater in Chespeake Bay";

    double rkdS2;
      :units = "meter-1 PSU-1";
      :long_name = "light attenuation by CDOM (use salinity instead) at high salinity seawater in Chespeake Bay";

    double pCO2air;
      :units = "parts per million by volume";
      :long_name = "partial pressure of CO2 in the air";

    double ElDON;
      :units = "day-1";
      :long_name = "phytoplankton exudation rate of labile DON";

    double EsDON;
      :units = "day-1";
      :long_name = "phytoplankton exudation rate of semilabile DON";

    double deltN;
      :long_name = "ratio of semilabile DON to total DON in cell";

    double deltC;
      :long_name = "ratio of semilabile DOC to total DOC in cell";

    double gammaC;
      :long_name = "parameter of carbon excess-based DOC exudation";

    double slCexc;
      :long_name = "fraction of semilabile DOC produced by C excess uptake";

    double a0N;
      :units = "day-1";
      :long_name = "Remineralization rate of DON at 0 deg.C [1/day]";

    double a0C;
      :units = "day-1";
      :long_name = "Remineralization rate of DOC at 0 deg.C [1/day]";

    double LDeNSR;
      :units = "day-1";
      :long_name = "Large detritus solubilization rate of N at 0C [1/day]";

    double LDeCSR;
      :units = "day-1";
      :long_name = "Large detritus solubilization rate of C at 0C [1/day]";

    double SDeNSR;
      :units = "day-1";
      :long_name = "Small detritus solubilization rate of N at 0C [1/day]";

    double SDeCSR;
      :units = "day-1";
      :long_name = "Small detritus solubilization rate of C at 0C [1/day]";

    double AttDOC;
      :units = "meter2 milligram_Chl-1";
      :long_name = "light attenuation by semilabile and refractory DOC";

    double DOCatt_m;
      :units = "millimole_C meter-3";
      :long_name = "offset to light attenuation by DOC";

    int spherical;
      :long_name = "grid type logical switch";
      :flag_values = 0, 1; // int
      :flag_meanings = "Cartesian spherical";

    double xl;
      :units = "meter";
      :long_name = "domain length in the XI-direction";

    double el;
      :units = "meter";
      :long_name = "domain length in the ETA-direction";

    int Vtransform;
      :long_name = "vertical terrain-following transformation equation";

    int Vstretching;
      :long_name = "vertical terrain-following stretching function";

    double theta_s;
      :long_name = "S-coordinate surface control parameter";

    double theta_b;
      :long_name = "S-coordinate bottom control parameter";

    double Tcline;
      :units = "meter";
      :long_name = "S-coordinate surface/bottom layer width";

    double hc;
      :units = "meter";
      :long_name = "S-coordinate parameter, critical depth";

    int grid;
      :cf_role = "grid_topology";
      :topology_dimension = 2; // int
      :node_dimensions = "xi_psi eta_psi";
      :face_dimensions = "xi_rho: xi_psi (padding: both) eta_rho: eta_psi (padding: both)";
      :edge1_dimensions = "xi_u: xi_psi eta_u: eta_psi (padding: both)";
      :edge2_dimensions = "xi_v: xi_psi (padding: both) eta_v: eta_psi";
      :node_coordinates = "lon_psi lat_psi";
      :face_coordinates = "lon_rho lat_rho";
      :edge1_coordinates = "lon_u lat_u";
      :edge2_coordinates = "lon_v lat_v";
      :vertical_dimensions = "s_rho: s_w (padding: none)";

    double Cs_r(s_rho=20);
      :long_name = "S-coordinate stretching curves at RHO-points";
      :valid_min = -1.0; // double
      :valid_max = 0.0; // double
      :field = "Cs_r, scalar";

    double Cs_w(s_w=21);
      :long_name = "S-coordinate stretching curves at W-points";
      :valid_min = -1.0; // double
      :valid_max = 0.0; // double
      :field = "Cs_w, scalar";

    double h(eta_rho=150, xi_rho=100);
      :units = "meter";
      :long_name = "bathymetry at RHO-points";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";
      :field = "bath, scalar";

    double f(eta_rho=150, xi_rho=100);
      :units = "second-1";
      :long_name = "Coriolis parameter at RHO-points";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";
      :field = "coriolis, scalar";

    double pm(eta_rho=150, xi_rho=100);
      :units = "meter-1";
      :long_name = "curvilinear coordinate metric in XI";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";
      :field = "pm, scalar";

    double pn(eta_rho=150, xi_rho=100);
      :units = "meter-1";
      :long_name = "curvilinear coordinate metric in ETA";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";
      :field = "pn, scalar";

    double angle(eta_rho=150, xi_rho=100);
      :units = "radians";
      :long_name = "angle between XI-axis and EAST";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";
      :field = "angle, scalar";

    double mask_rho(eta_rho=150, xi_rho=100);
      :long_name = "mask on RHO-points";
      :flag_values = 0.0, 1.0; // double
      :flag_meanings = "land water";
      :grid = "grid";
      :location = "face";
      :coordinates = "lat_rho lon_rho ";

    double mask_u(eta_u=150, xi_u=99);
      :long_name = "mask on U-points";
      :flag_values = 0.0, 1.0; // double
      :flag_meanings = "land water";
      :grid = "grid";
      :location = "edge1";
      :coordinates = "lat_u lon_u ";

    double mask_v(eta_v=149, xi_v=100);
      :long_name = "mask on V-points";
      :flag_values = 0.0, 1.0; // double
      :flag_meanings = "land water";
      :grid = "grid";
      :location = "edge2";
      :coordinates = "lat_v lon_v ";

    double mask_psi(eta_psi=149, xi_psi=99);
      :long_name = "mask on psi-points";
      :flag_values = 0.0, 1.0; // double
      :flag_meanings = "land water";
      :grid = "grid";
      :location = "node";
      :coordinates = "lat_psi lon_psi ";

    float zeta(time3=6362, eta_rho=150, xi_rho=100);
      :units = "meter";
      :long_name = "free-surface";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "free-surface, scalar, series";
      :_FillValue = 1.0E37f; // float

    float ubar(time3=6362, eta_u=150, xi_u=99);
      :units = "meter second-1";
      :long_name = "vertically integrated u-momentum component";
      :time = "ocean_time";
      :grid = "grid";
      :location = "edge1";
      :coordinates = "time3_run time3 lat_u lon_u ";
      :field = "ubar-velocity, scalar, series";
      :_FillValue = 1.0E37f; // float

    float vbar(time3=6362, eta_v=149, xi_v=100);
      :units = "meter second-1";
      :long_name = "vertically integrated v-momentum component";
      :time = "ocean_time";
      :grid = "grid";
      :location = "edge2";
      :coordinates = "time3_run time3 lat_v lon_v ";
      :field = "vbar-velocity, scalar, series";
      :_FillValue = 1.0E37f; // float

    float temp(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "Celsius";
      :long_name = "potential temperature";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "temperature, scalar, series";
      :_FillValue = 1.0E37f; // float

    float salt(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :long_name = "salinity";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "salinity, scalar, series";
      :_FillValue = 1.0E37f; // float

    float NO3(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_N03 meter-3";
      :long_name = "nitrate concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "NO3, scalar, series";
      :_FillValue = 1.0E37f; // float

    float NH4(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_NH4 meter-3";
      :long_name = "ammonium concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "NH4, scalar, series";
      :_FillValue = 1.0E37f; // float

    float chlorophyll(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "milligrams_chlorophyll meter-3";
      :long_name = "chlorophyll concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "chlorophyll, scalar, series";
      :_FillValue = 1.0E37f; // float

    float phytoplankton(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "phytoplankton concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "phytoplankton, scalar, series";
      :_FillValue = 1.0E37f; // float

    float zooplankton(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "zooplankton concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "zooplankton, scalar, series";
      :_FillValue = 1.0E37f; // float

    float LdetritusN(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "large fraction nitrogen detritus concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "LdetritusN, scalar, series";
      :_FillValue = 1.0E37f; // float

    float SdetritusN(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "small fraction nitrogen detritus concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "SdetritusN, scalar, series";
      :_FillValue = 1.0E37f; // float

    float LdetritusC(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_carbon meter-3";
      :long_name = "large fraction carbon detritus concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "LdetritusC, scalar, series";
      :_FillValue = 1.0E37f; // float

    float SdetritusC(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_carbon meter-3";
      :long_name = "small fraction carbon detritus concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "SdetritusC, scalar, series";
      :_FillValue = 1.0E37f; // float

    float TIC(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_carbon meter-3";
      :long_name = "total inorganic carbon";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "TIC, scalar, series";
      :_FillValue = 1.0E37f; // float

    float alkalinity(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "milliequivalents meter-3";
      :long_name = "total alkalinity";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "alkalinity, scalar, series";
      :_FillValue = 1.0E37f; // float

    float oxygen(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_oxygen meter-3";
      :long_name = "dissolved oxygen concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "Oxygen, scalar, series";
      :_FillValue = 1.0E37f; // float

    float semilabileDON(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "semilabile dissolved organic nitrogen concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "semilabileDON, scalar, series";
      :_FillValue = 1.0E37f; // float

    float refractoryDON(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_nitrogen meter-3";
      :long_name = "refractory dissolved organic nitrogen concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "refractoryDON, scalar, series";
      :_FillValue = 1.0E37f; // float

    float semilabileDOC(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_carbon meter-3";
      :long_name = "semilabile dissolved organic carbon concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "semilabileDOC, scalar, series";
      :_FillValue = 1.0E37f; // float

    float refractoryDOC(time3=6362, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_carbon meter-3";
      :long_name = "refractory dissolved organic carbon concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 s_rho lat_rho lon_rho ";
      :field = "refractoryDOC, scalar, series";
      :_FillValue = 1.0E37f; // float

    float PO4(time2=3369, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "millimole_PO4 meter-3";
      :long_name = "phosphate concentration";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time2_run time2 s_rho lat_rho lon_rho ";
      :field = "PO4, scalar, series";
      :_FillValue = 1.0E37f; // float

    float ISS1(time2=3369, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "milligram per litre";
      :long_name = "inorganic suspended solids class#1";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time2_run time2 s_rho lat_rho lon_rho ";
      :field = "ISS1, scalar, series";
      :_FillValue = 1.0E37f; // float

    float ISS2(time2=3369, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "milligram per litre";
      :long_name = "inorganic suspended solids class#2";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time2_run time2 s_rho lat_rho lon_rho ";
      :field = "ISS2, scalar, series";
      :_FillValue = 1.0E37f; // float

    float shflux(time3=6362, eta_rho=150, xi_rho=100);
      :units = "watt meter-2";
      :long_name = "surface net heat flux";
      :negative_value = "upward flux, cooling";
      :positive_value = "downward flux, heating";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "surface heat flux, scalar, series";
      :_FillValue = 1.0E37f; // float

    float ssflux(time3=6362, eta_rho=150, xi_rho=100);
      :units = "meter second-1";
      :long_name = "surface net salt flux, (E-P)*SALT";
      :negative_value = "upward flux, freshening (net precipitation)";
      :positive_value = "downward flux, salting (net evaporation)";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "surface net salt flux, scalar, series";
      :_FillValue = 1.0E37f; // float

    float latent(time3=6362, eta_rho=150, xi_rho=100);
      :units = "watt meter-2";
      :long_name = "net latent heat flux";
      :negative_value = "upward flux, cooling";
      :positive_value = "downward flux, heating";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "latent heat flux, scalar, series";
      :_FillValue = 1.0E37f; // float

    float sensible(time3=6362, eta_rho=150, xi_rho=100);
      :units = "watt meter-2";
      :long_name = "net sensible heat flux";
      :negative_value = "upward flux, cooling";
      :positive_value = "downward flux, heating";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "sensible heat flux, scalar, series";
      :_FillValue = 1.0E37f; // float

    float lwrad(time3=6362, eta_rho=150, xi_rho=100);
      :units = "watt meter-2";
      :long_name = "net longwave radiation flux";
      :negative_value = "upward flux, cooling";
      :positive_value = "downward flux, heating";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "longwave radiation, scalar, series";
      :_FillValue = 1.0E37f; // float

    float swrad(time3=6362, eta_rho=150, xi_rho=100);
      :units = "watt meter-2";
      :long_name = "solar shortwave radiation flux";
      :negative_value = "upward flux, cooling";
      :positive_value = "downward flux, heating";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time3_run time3 lat_rho lon_rho ";
      :field = "shortwave radiation, scalar, series";
      :_FillValue = 1.0E37f; // float

    float Aragonite(time=1981, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "n/a";
      :long_name = "aragonite saturation state";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time_run time s_rho lat_rho lon_rho ";
      :field = "Aragonite, scalar, series";
      :_FillValue = 1.0E37f; // float

    float pH(time=1981, s_rho=20, eta_rho=150, xi_rho=100);
      :units = "n/a";
      :long_name = "pH";
      :time = "ocean_time";
      :grid = "grid";
      :location = "face";
      :coordinates = "time_run time s_rho lat_rho lon_rho ";
      :field = "pH, scalar, series";
      :_FillValue = 1.0E37f; // float

    double s_rho(s_rho=20);
      :units = "";
      :long_name = "S-coordinate at RHO-points";
      :valid_min = -1.0; // double
      :valid_max = 0.0; // double
      :positive = "up";
      :standard_name = "ocean_s_coordinate_g2";
      :formula_terms = "s: s_rho C: Cs_r eta: zeta depth: h depth_c: hc";
      :field = "s_rho, scalar";
      :_CoordinateTransformType = "Vertical";
      :_CoordinateAxisType = "GeoZ";
      :_CoordinateZisPositive = "up";
      :_CoordinateAxes = "s_rho";

    double s_w(s_w=21);
      :units = "";
      :long_name = "S-coordinate at W-points";
      :valid_min = -1.0; // double
      :valid_max = 0.0; // double
      :positive = "up";
      :standard_name = "ocean_s_coordinate_g2";
      :formula_terms = "s: s_w C: Cs_w eta: zeta depth: h depth_c: hc";
      :field = "s_w, scalar";
      :_CoordinateTransformType = "Vertical";
      :_CoordinateAxisType = "GeoZ";
      :_CoordinateZisPositive = "up";
      :_CoordinateAxes = "s_w";

    double lon_rho(eta_rho=150, xi_rho=100);
      :units = "degrees_east";
      :long_name = "longitude of RHO-points";
      :standard_name = "longitude";
      :field = "lon_rho, scalar";
      :_CoordinateAxisType = "Lon";

    double lat_rho(eta_rho=150, xi_rho=100);
      :units = "degrees_north";
      :long_name = "latitude of RHO-points";
      :standard_name = "latitude";
      :field = "lat_rho, scalar";
      :_CoordinateAxisType = "Lat";

    double lon_u(eta_u=150, xi_u=99);
      :units = "degrees_east";
      :long_name = "longitude of U-points";
      :standard_name = "longitude";
      :field = "lon_u, scalar";
      :_CoordinateAxisType = "Lon";

    double lat_u(eta_u=150, xi_u=99);
      :units = "degrees_north";
      :long_name = "latitude of U-points";
      :standard_name = "latitude";
      :field = "lat_u, scalar";
      :_CoordinateAxisType = "Lat";

    double lon_v(eta_v=149, xi_v=100);
      :units = "degrees_east";
      :long_name = "longitude of V-points";
      :standard_name = "longitude";
      :field = "lon_v, scalar";
      :_CoordinateAxisType = "Lon";

    double lat_v(eta_v=149, xi_v=100);
      :units = "degrees_north";
      :long_name = "latitude of V-points";
      :standard_name = "latitude";
      :field = "lat_v, scalar";
      :_CoordinateAxisType = "Lat";

    double lon_psi(eta_psi=149, xi_psi=99);
      :units = "degrees_east";
      :long_name = "longitude of PSI-points";
      :standard_name = "longitude";
      :field = "lon_psi, scalar";
      :_CoordinateAxisType = "Lon";

    double lat_psi(eta_psi=149, xi_psi=99);
      :units = "degrees_north";
      :long_name = "latitude of PSI-points";
      :standard_name = "latitude";
      :field = "lat_psi, scalar";
      :_CoordinateAxisType = "Lat";

    double ocean_time(ocean_time=13);
      :units = "seconds since 2009-01-01 00:00:00";
      :long_name = "time since initialization";
      :calendar = "gregorian_proleptic";
      :field = "time, scalar, series";
      :_CoordinateAxisType = "Time";

    double time(time=1981);
      :long_name = "Forecast time for ForecastModelRunCollection";
      :standard_name = "time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "Time";

    double time_run(time=1981);
      :long_name = "run times for coordinate = time";
      :standard_name = "forecast_reference_time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "RunTime";

    double time_offset(time=1981);
      :long_name = "offset hour from start of run for coordinate = time";
      :standard_name = "forecast_period";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01T00:00:00Z";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "TimeOffset";

    double time1(time1=3002);
      :long_name = "Forecast time for ForecastModelRunCollection";
      :standard_name = "time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "Time";

    double time1_run(time1=3002);
      :long_name = "run times for coordinate = time1";
      :standard_name = "forecast_reference_time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "RunTime";

    double time1_offset(time1=3002);
      :long_name = "offset hour from start of run for coordinate = time1";
      :standard_name = "forecast_period";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01T00:00:00Z";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "TimeOffset";

    double time2(time2=3369);
      :long_name = "Forecast time for ForecastModelRunCollection";
      :standard_name = "time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "Time";

    double time2_run(time2=3369);
      :long_name = "run times for coordinate = time2";
      :standard_name = "forecast_reference_time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "RunTime";

    double time2_offset(time2=3369);
      :long_name = "offset hour from start of run for coordinate = time2";
      :standard_name = "forecast_period";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01T00:00:00Z";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "TimeOffset";

    double time3(time3=6362);
      :long_name = "Forecast time for ForecastModelRunCollection";
      :standard_name = "time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "Time";

    double time3_run(time3=6362);
      :long_name = "run times for coordinate = time3";
      :standard_name = "forecast_reference_time";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01 00:00:00.000 UTC";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "RunTime";

    double time3_offset(time3=6362);
      :long_name = "offset hour from start of run for coordinate = time3";
      :standard_name = "forecast_period";
      :calendar = "proleptic_gregorian";
      :units = "hours since 2017-03-01T00:00:00Z";
      :missing_value = NaN; // double
      :_CoordinateAxisType = "TimeOffset";

  // global attributes:
  :file = "chesroms_ECB_his.nc";
  :format = "netCDF-3 64bit offset file";
  :Conventions = "CF-1.4, SGRID-0.3";
  :type = "ROMS/TOMS history file";
  :title = "ChesROMS ECB Coupled Physical-Biogeochemical Model";
  :var_info = "varinfo_2019.dat";
  :rst_file = "chesroms_ECB_rst.nc";
  :his_file = "chesroms_ECB_his.nc";
  :avg_file = "chesroms_ECB_avg.nc";
  :sta_file = "chesroms_ECB_sta.nc";
  :grd_file = "/ches/data10/abever/WORK/FORECAST/ECB/GENERAL_FORCING/chesroms_grd_smoother_msl_mes.nc";
  :ini_file = "chesroms_ECB_ini.nc";
  :tide_file = "/ches/data10/abever/WORK/FORECAST/ECB/GENERAL_FORCING/ChesROMS_frc_tides_COMTEHYear3_2015_062016.nc";
  :frc_file_01 = "/ches/data10/abever/WORK/OCEAN_ACIDIFICATION/MODEL_INPUTS/FORECAST_SETUP/ChesROMS_atmo_nit_dep_forcing_2019_repeat_Pierre2017_AB.nc, /ches/data10/abever/WORK/OCEAN_ACIDIFICATION/MODEL_INPUTS/FORECAST_SETUP/ChesROMS_atmo_nit_dep_forcing_2020_repeat_Pierre2017_AB.nc, /ches/data10/abever/WORK/OCEAN_ACIDIFICATION/MODEL_INPUTS/FORECAST_SETUP/ChesROMS_atmo_nit_dep_forcing_2021_repeat_Pierre2017_AB.nc";
  :frc_file_02 = "chesroms_met_NF.nc";
  :bry_file = "chesroms_ECB_obc.nc";
  :script_file = "ocean_ECB_2019.in";
  :bpar_file = "bio_ECB_2019.in";
  :spos_file = "stations_ECB_011420.in";
  :NLM_LBC = "\nEDGE:           WEST   SOUTH  EAST   NORTH  \nzeta:           Clo    Clo    Cha    Clo    \nubar:           Clo    Clo    Fla    Clo    \nvbar:           Clo    Clo    Fla    Clo    \nu:              Clo    Clo    Rad    Clo    \nv:              Clo    Clo    Rad    Clo    \ntemp:           Clo    Clo    RadNud Clo    \nsalt:           Clo    Clo    RadNud Clo    \nNO3:            Clo    Clo    RadNud Clo    \nNH4:            Clo    Clo    RadNud Clo    \nchlorophyll:    Clo    Clo    Rad    Clo    \nphytoplankton:  Clo    Clo    Rad    Clo    \nzooplankton:    Clo    Clo    Rad    Clo    \nLdetritusN:     Clo    Clo    Rad    Clo    \nSdetritusN:     Clo    Clo    Rad    Clo    \nLdetritusC:     Clo    Clo    Rad    Clo    \nSdetritusC:     Clo    Clo    Rad    Clo    \nTIC:            Clo    Clo    RadNud Clo    \nalkalinity:     Clo    Clo    RadNud Clo    \noxygen:         Clo    Clo    RadNud Clo    \nsemilabileDON:  Clo    Clo    RadNud Clo    \nrefractoryDON:  Clo    Clo    RadNud Clo    \nsemilabileDOC:  Clo    Clo    RadNud Clo    \nrefractoryDOC:  Clo    Clo    RadNud Clo    \nPO4:            Clo    Clo    Rad    Clo    \nISS1:           Clo    Clo    Rad    Clo    \nISS2:           Clo    Clo    Rad    Clo    \ntke:            Clo    Clo    Rad    Clo";
  :svn_url = "https://www.myroms.org/svn/src/trunk";
  :svn_rev = "exported";
  :code_dir = "/ches/data10/abever/WORK/FORECAST/ECB/EXECUTABLE/VERSION_2019_JAN/CODE/SOURCE";
  :header_dir = "/ches/home00/abever/data10/WORK/FORECAST/ECB/EXECUTABLE/VERSION_2019_JAN/COMPILE";
  :header_file = "ecb_forecast.h";
  :os = "Linux";
  :cpu = "x86_64";
  :compiler_system = "pgi";
  :compiler_command = "/usr/local/amd64/seoul/pgi-14.3/mvapich2-1.9-ib/bin/mpif90";
  :compiler_flags = "-tp piledriver -m64 -Kieee -O3 -Mfree";
  :tiling = "002x012";
  :history = "Thu Jul 15 04:20:59 2021: ncks -A -v Aragonite,pH chesroms_ECB_his_tmp.nc chesroms_ECB_his.nc\nThu Jul  8 04:21:03 2021: ncks -A -v Aragonite,pH chesroms_ECB_his_tmp.nc chesroms_ECB_his.nc\nROMS/TOMS, Version 3.7, Thursday - July 8, 2021 -  3:51:35 AM ;\nFMRC Best Dataset";
  :ana_file = "ROMS/Functionals/ana_btflux.h, ROMS/Functionals/ana_rain.h, ROMS/Functionals/ana_stflux.h";
  :bio_file = "ROMS/Nonlinear/Biology/fennel.h";
  :CPP_options = "ChesROMS_ECB, ADD_FSOBC, ADD_M2OBC, ANA_BPFLUX, ANA_BSFLUX, ANA_BTFLUX, ANA_RAIN, ANA_SPFLUX, ANA_SSFLUX, ASSUMED_SHAPE, AVERAGES, BIO_FENNEL, BIO_SEDIMENT, !BOUNDARY_A BULK_FLUXES, CARBON, !COLLECT_ALL..., CURVGRID, DENITRIFICATION, DJ_GRADPS, DOUBLE_PRECISION, GLS_MIXING, LIMIT_STFLX_COOLING, KANTHA_CLAYSON, LONGWAVE_OUT, MASKING, MPI, NONLINEAR, NONLIN_EOS, N2S2_HORAVG, !OCMIP_OXYGEN_SC, OXYGEN, PERFECT_RESTART, POWER_LAW, PROFILE, K_GSCHEME, REDUCE_ALLGATHER, !RST_SINGLE, SALINITY, SOLAR_SOURCE, SOLVE3D, SSH_TIDES, STATIONS, TALK_NONCONSERV, TS_MPDATA, UV_ADV, UV_COR, UV_U3HADVECTION, UV_C4VADVECTION, UV_QDRAG, UV_TIDES, VAR_RHO_2D, WIND_MINUS_CURRENT";
  :NCO = "4.4.4";
  :_CoordSysBuilder = "ucar.nc2.dataset.conv.CF1Convention";
  :cdm_data_type = "GRID";
  :featureType = "GRID";
  :location = "Proto fmrc:ECB_ROMS";
  :protometadata_speech = "I\'d like to thank all the little people...";
}