c1440_llc2160 spans 20200119_2100z to 20210326_1800z MItgcm timesteps 0000000000 to 0000829200 MITgcm deltaT = 45 s Available hourly output is: Eta.*.data sea surface height (m) KPPhbl.*.data mixing layer depth (m) PhiBot.*.data bottom pressure (m^2/s^2) PhiBot = g * ( etaN + \int (rho/rhoConst - 1) dr ) convert PhiBot to m of water: divide by g = 9.81 m/s^2 convert PhiBot to pressure in Pa: multiply by rhoConst = 1027.5 kg/m^3 PhiBot = g * ( etaN + \int (rho/rhoConst - 1) dr ) that is, PhiBot is the anomaly relative to Depth * rhoConst * g so that absolute bottom pressure in Pa is: Depth * rhoConst * g + PHIBOT * rhoConst SIarea.*.data fractional ice-covered area for 5 categories [0 to 1] SIheff.*.data effective ice thickness for 5 categories (m) SIhsnow.*.data effective snow thickness for 5 categories (m) SItice ice surface temperature for 5 categories (deg K) SIuice.*.data zonal (relative to grid) ice velocity, >0 from West to East (m/s) SIvice.*.data merid. (relative to grid) ice velocity, >0 from South to North (m/s) Salt.*.data salinity (g/kg) Theta.*.data potential temperature (deg C) U.*.data zonal (relative to grid) velocity, >0 from West to East (m/s) specified on Southwest C-grid U point V.*.data merid. (relative to grid) velocity, >0 from South to North (m/s) specified on Southwest C-grid V point W.*.data vertical velocity (m/s) oceFWflx.*.data net upward freshwater flux, >0 increases salinity (kg/m^2/s) oceQnet.*.data net upward surface heat flux (including shortwave), >0 decreases theta (W/m^2) oceQsw.*.data net upward shortwave radiation, >0 decreases theta (W/m^2) oceSflux.*.data net upward salt flux, >0 decreases salinity (g/m^2/s) oceTAUX.*.data zonal (relative to grid) surface wind stress, >0 increases uVel (N/m^2) specified on Southwest C-grid U point oceTAUY.*.data meridional (relative to grid) surf. wind stress, >0 increases vVel (N/m^2) specified on Southwest C-grid V point Please note that U, V, oceTAUX, oceTAUY, SIuice, and SIvice are aligned relative to model grid, not geographical coordinates, and that they are specified at the SouthWest C-grid velocity points. All other scalar fields are specified at the tracer point, i.e., the center of each grid box. Sea Surface Height (SSH) components: SSH = Eta + sIceLoad / rhoConst sIceLoad = SIheff * SEAICE_rhoIce + SIhsnow * SEAICE_rhoSnow SEAICE_rhoIce = 910 SEAICE_rhoSnow = 330 rhoConst = 1027.5 gravity = 9.81 https://github.com/MITgcm/MITgcm/files/7521065/steric.pdf Steric = SSH - PhiBot / gravity + (AtmosPres + TidalPotential) / rhoConst / gravity - AtmosPres - TidalPotential ================== atmospheric fields ================== PS/ U10M/ V10M/ from /nobackupp17/dmenemen/DYAMOND/c1440_llc2160/holding/inst_15mn_2d_asm_Mx/ C1440 grid interpolated to LLC2160 grid first file (ie, *.0000000000.data) from "21:15" instead of "21:00" because "21:00" is unavailable. Unlike oceTAUX/oceTAUY which is exactly aligned w/ llc-grid, U10M/V10M is Eastward/Northward on llc-grid, similar to " UE_VEL_C|UMR MR|m/s |Eastward Velocity (m/s) (cell center) VN_VEL_C|VMR MR|m/s |Northward Velocity (m/s) (cell center) " from 9.1.4.4. MITgcm kernel available diagnostics list at https://mitgcm.readthedocs.io/en/latest/outp_pkgs/outp_pkgs.html ========== SSH fields ========== SSH/ SSH_notides/ SSH_noIB/ SSH_steric/ SSH = Eta + sIceLoad / rhoConst where sIceLoad = SIheff * SEAICE_rhoIce + SIhsnow * SEAICE_rhoSnow SEAICE_rhoIce = 910 SEAICE_rhoSnow = 330 rhoConst = 1027.5 gravity = 9.81 SSH_notides = SSH - tides (by "t_tides" package) SSH_noIB = SSH_notides + (PS - PS_globalmean)/rhoConst/gravity SSH_steric = SSH - PhiBot/gravity + PS/rhoConst/gravity where PhiBot = PhiBot - PhiBot_timemean (over last 12 months) PS = PS - surf_pRef (=101325) refer to https://github.com/MITgcm/MITgcm/files/7521065/steric.pdf =========== Additional integration details are here: https://github.com/MITgcm-contrib/llc_hires/tree/master/c1440_llc2160 and in subdirectory global_SSH, obtained from: https://github.com/MITgcm-contrib/llc_hires/tree/master/c1440_llc2160/global_SSH