Since 2006-06-06: ================= FNMOC High Resolution Ocean Analysis for GODAE These fileds are a FNMOC contribution to GODAE. Includes daily 3D temps, salts, geopotential and u,v current fields on a nearly global mercator grid. There are 2161 longitudes, 1051 latitudes and 34 levels in each file (2161x1051x34). See FILE FORMAT below The 34 grid levels are depths (m) and are as follows: z_lvl = 0., 2.5, 7.5, 12.5, 17.5, 25., 32.5, 40., 50., 62.5, 75., 100., 125., 150., 200., 300., 400., 500., 600., 700., 800., 900., 1000., 1100., 1200., 1300., 1400., 1500., 1750., 2000., 2500., 3000., 4000., 5000. NOTE: Prior to 2006-09-22: Only the 1st 31 levels were present Ascii format grid defintions may be found in the GrADS .ctl files. The grid point latitudes and longitudes are available in: grdlat_sfc_000000_000000_1o2161x1051_00000000_datafld grdlon_sfc_000000_000000_1o2161x1051_00000000_datafld These are the lat/lon at each grid point stored as described in FILE FORMAT below. Prior to 2006-06-06: ==================== There are 1441 longitudes, 701 latitudes and 32 levels in each file (1441x701x32). See FILE FORMAT below. The 32 grid levels are depths (m) and are as follows: z_lvl = 0., 2.5, 7.5, 12.5, 17.5, 25., 32.5, 40., 50., 62.5, 75., 100., 125., 150., 200., 300., 400., 500., 600., 700., 800., 900., 1000., 1100., 1200., 1300., 1400., 1500., 1750., 2000., 2500., 3000., Ascii format grid defintions may be found in the GrADS .ctl files. 32 bit direct access binary grid definitions are in: grdlat_sfc_000000_000000_1o1441x0701_00000000_datafld grdlon_sfc_000000_000000_1o1441x0701_00000000_datafld FILE FORMAT =========== The files are big-endian 32 bit IEEE floating point files. (In Fortan parlance, the files are direct access, unformatted files, 1 record per file.) Missing data are indicated by the value -999. The "*_sfc_*" files are 2-D dimensional data files: Each file contains the data for the M longitudes and N latitudes, at the surface; stored in that order -- M by N The "*_pre_*" files are 3-D dimensional data files: Each file contains the data for the M longitudes, N latitudes, and L levels; stored in that order -- M by N by L. Using Fortran as an example (compiled with the "big-endian" option; gfortran -fconvert=big-endian), to read the file: integer, parameter :: M = 2161 integer, parameter :: N = 1051 integer, parameter :: L = 34 integer, parameter :: pre_reclen = M * N * L * 4 !..record length in bytes integer, parameter :: sfc_reclen = M * N * 4 !..record length in bytes real :: data_3d(M, N, L) real :: data_2d(M, N) !..read 2-D data open (1, file="..._sfc_...", access='direct', form='unformatted', recl=sfc_reclen) read (1, rec=1) data_2d close (1) !..read 3-D data open (1, file="..._pre_...", access='direct', form='unformatted', recl=reclen) read (1, rec=1) data_3d close (1)