August GDEM Climatology 6 IHFRMT : input hydrography file format (see below) 0 ICLIMA : option to compute mean fields (see below) 6 ICFRMT : input climatology file format (see below) 2 COORD : grid: [0] Cartesian, [1] spherical, [2] rotated spherical 159 IM : number of points in the zonal direction 183 JM : number of points in the meridional direction 0.0810 DX : zonal grid spacing (km or degrees) 0.0810 DY : meridional grid spacing (km or degrees) -126.183208 CLNG : domain centroid longitude (degrees) 34.7354214 CLAT : domain centroid latitude (degrees) 0.0 DELX : X-offset from transformation point to grid center. (km or deg) 0.0 DELY : Y-offset from transformation point to grid center. (km or deg) 31.5959681 THETAD : domain rotation (degrees, counterclockwise from EAST) 100 NNCE : number of influential points 556.0 RNCE : radius of influence (km) 1000.0 TNCE : influential time window (days) 0 RAVG : averaging radius for nearby observations (km) 0 TAVG : averaging time for nearby observations (days) 5000.0 ROUT : outlier radius (km), observations outside ROUT are rejected 1 ICORR : type of analytical correlation (see below) 70.0 XZERO : correlation, zonal zero crossing (km) 70.0 YZERO : correlation, meridional zero crossing (km) 27.5 XDCAY : zonal decorrelation (decay) scale (km) 27.5 YDCAY : meridional decorrelation (decay) scale (km) 1.0 TDCAY : temporal decorrelation (decay) scale (days) 0 UPHSE : correlation, zonal phase speed (km/day) 0 VPHSE : correlation, meridional phase speed (km/day) 0 CORANG : correlation, rotation angle (deg) counterclockwise from EAST 280.0 XZEROM : MEAN - correlation, zonal zero crossing (km) 280.0 YZEROM : MEAN - correlation, meridional zero crossing (km) 110.0 XDCAYM : MEAN - zonal decorrelation (decay) scale (km) 110.0 YDCAYM : MEAN - meridional decorrelation (decay) scale (km) 10.0 TDCAYM : MEAN - temporal decorrelation (decay) scale (days) 0 UPHSEM : MEAN - correlation, zonal phase speed (km/day) 0 VPHSEM : MEAN - correlation, meridional phase speed (km/day) 0 CORANGM: MEAN - correlation, rotation angle (deg) from EAST 0 ITRT : detrend option: [0] mean, [1] linear, [2] quadratic, [3] cubic 0.20 OBSERR : observation error (value to add to correlation matrix diagonal) 0.40 CLMERR : climatology error (value to add to correlation matrix diagonal) 2 NORD : order to of Shapiro filter to apply to each OA field 1 NFRQ : number of times to Shapiro filter the field (0 no filtering) 3 NFLD : number of fields to OA 2,3,4 IDFLD : field(s) type to OA, IDFLD(1),...,IDFLD(NFLD), (see ID below) 3 IDHOPT : dynamic height option: [0] no, [1,2,3,4] yes (see below) 2000.0 ZREF : level of no motion (in meters) for dynamic height computation 1 NDAY : number of days to OA. Line below, OADAY(1),...,OADAY(NDAY) 12867.0 31 NLEV : number of depth levels to OA. Line below, depth(s) in meters: 0.0 10.0 20.0 30.0 50.0 75.0 100.0 125.0 150.0 200.0 250.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0 1100.0 1200.0 1300.0 1400.0 1500.0 1750.0 2000.0 2500.0 3000.0 4000.0 5000.0 5500.0 1 IGRIDS : flag to read in GRIDS NetCDF file: [0] no, [1] yes 0 IWRT : flag to write out verbose output: [0] no, [1] yes oa_gdem_aug.nc oa_gdem_aug.echo /home/projects2/AOSN2/Data/Climo/GDEM/gdem_aug_big9km2.mods /home/projects2/AOSN2/Data/Climo/GDEM/gdem_aug_big9km2.mods /home/projects2/AOSN2/Grids/Big9km2/grids_AOSN2_big_nois.nc *** Above FILENAMES: 1st line: output (NetCDF) objectively analyzed field(s) 2nd line: output (ascii) objective analysis information 3nd line: input hydrographic data 4th line: input climatological data 5th line: input GRIDS NetCDF file *** IHFRMT, ICFRMT: input file types: [0] unformatted: depth, temperature, salinity, dynamic height [1] unformatted: depth, temperature and salinity (CTD,XCTD) [2] unformatted: depth and temperature (XBT,AXBT) [3] unformatted: depth and dynamic height [4] unformatted: depth and any oceanic field [5] unformatted: any field (a single value; NHPTS=NHVAR=1) [6] formatted: HOPS hydrography ASCII format NOTE: the read statement for options 0 to 5 is as follows: read (iunit) castid, hlat, hlng, jday, itime read (iunit) nhpts, ((fields(i,j),i=nhvar),j=1,nhpts) If IHFRMT=5, the following parameters are reset internally: NFLD=1, IDFLD(1)=0, NDAY=1, OADAY(1)=0, NLEV=1, Z(1)=-9. *** ICLIMA: mean field computation option [0] OA mean field from hydrographic observations only [1] OA mean field from climatologic observations only [2] OA mean field using both hydrographic and climatologic observations [3] Use global observation mean (a single value) NOTE: user must provide climatology file when using ICLIMA=1,2,3. For very sparse hydrographic data, it is recommended to use both hydrographic data and climatological data (e.g. Levitus, Bauer and Robinson, etc.) when computing mean fields. Climatology data is a more oceanographically correct field than any extrapolation of data to places of no observations. In dynamically active regions, the melding of hydrographic and climatological data to compute mean fields is like assuming non-stationary statistics for the mean. *** IDFLD: field types to objectively analyze: [0] generic oceanic field [1] pressure (depth) [2] temperature [3] salinity [4] dynamic height [5] density [6] sigma-t [7] sound speed *** ICORR: analytical correlation option: [1] Rotated or non-rotated, time displaced or not, Cartesian or spherical, isotropic (symmetric: circular) or anisotropic (asymmetric: elliptical) Gaussian correlation C(Xj,Yj,Tj,Xi,Yi,Ti) = ( 1 - a2 ) * EXP ( b2 + c2 ) where a2 = Rx^2 / xzero^2 + Ry^2 / yzero^2 b2 = - 0.5 * ( Rx^2 / xdcay^2 + Ry^2 / ydcay^2 ) c2 = - 0.5 * delt^2 / tdcay^2 Rx = xp * cos(phi) + yp * sin(phi) + delt * uphse Ry = yp * cos(phi) - xp * sin(phi) + delt * vphse phi = corang - thetad delt = Tj - Ti if Cartesian xp = Xj - Xi yp = Yj - Yi if spherical r = gcircle (Xj,Yj,Xi,Yi,dist,bearing) xp = dist*cos(bearing) yp = dist*sin(bearing) *** Correlation displacement speeds UPHSE and VPHSE are aligned with the correlation axes rotated by CORANG. *** XZEROM, YZEROM, XDCAYM, and YDCAYM correlation scales are only used to compute the mean fields. It is recommended to use values larger than the ones (XZERO,YZERO,XDCAY,YDCAY) used in the objective analysis interpolation. *** IDHOPT: Dynamic height computation option. This computation is only possible when pressure (or depth), temperature, and salinity are available at input (IHFRMT=0,1) and when requested with IDFLD=4. [0] DO NOT compute [1] integrate specific volume anomaly from ZREF to the surface [2] integrate specific volume anomaly from ZREF to the bottom [3] integrate in both [1] and [2] directions [4] integrate specific volume anomaly from the surface to the local bottom of the hydrographic profile. If dynamic height is provided at input (IHFRMT=3) and its interpolation is requested (IDFLD=4), set IDHOPT=0. *** If RAVG and TAVG are zero, then close observations are not averaged. *** NDAY can be used for statistical forecasting. If only one objective analysis is desired set NDAY=1 and if, in addition, it is desired to be centered at the mean observations time, set OADAY=0 (line below NDAY). *** The OADAY should be consisted with input data. For instance, Julian day or year day. *** The IWRT flag, when turned on, writes into the echo file a message every time that there are insufficient influential observations or every time that the interpolation is unsuccessful because a poor conditioned local covariance matrix. The mean field value is assigned at that particular interpolation point, if either of these two problems occur.