UWM/COADS Monthly Climatology (1945-89) - North Atlantic, Grid #4
1 IOTYPE : type of output floating-point data: [0] single, [1] double
1 IHGRID : horizontal grid: [0] from GRID NetCDF, [1] specified below
80 IM : number of points in the zonal direction
55 JM : number of points in the meridional direction
1.0 DELLNG : longitude uniform grid spacing (degrees)
1.0 DELLAT : latitude uniform grid spacing (degrees)
-140.0 CLNGD : domain centroid longitude (degrees, west values are negative)
29.0 CLATD : domain centroid latitude (degrees, south values are negative)
0 ICLIMA : option to compute OA mean fields (see below)
40 NNCE : number of influential points
300 RNCE : radius of influence (km)
10 TNCE : influential time window (days)
0 RAVG : averaging radius for nearby observations (km)
0 TAVG : averaging time for nearby observations (days)
1000 ROUT : outlier radius (km), observations outside ROUT are rejected
1 ICORR : type of analytical correlation (see below)
300 XZERO : correlation, zonal zero crossing (km)
300 YZERO : correlation, meridional zero crossing (km)
300 XDCAY : zonal decorrelation (decay) scale (km)
300 YDCAY : meridional decorrelation (decay) scale (km)
10 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
300 XZEROM : MEAN - correlation, zonal zero crossing (km)
300 YZEROM : MEAN - correlation, meridional zero crossing (km)
300 XDCAYM : MEAN - zonal decorrelation (decay) scale (km)
300 YDCAYM : MEAN - meridional decorrelation (decay) scale (km)
0 ITRT : detrend option: [0] mean, [1] linear, [2] quadratic, [3] cubic
0.1 OBSERR : observation error (value to add to correlation matrix diagonal)
0.2 CLMERR : climatology error (value to add to correlation matrix diagonal)
0 IOAER : OA error type to output, [0] normalized 0 to 1, [1] dimensional
0 IWRT : flag to write out verbose output: [0] no, [1] yes
2 NORD : order to of Shapiro filter to apply to each OA field
0 NFRQ : number of times to Shapiro filter the field (0 no filtering)
1 NFLD : number of fields to OA
19 IDFLD : field(s) type to OA, IDFLD(1:NFLD), see ID below
0 IDHOPT : dynamic height option: [0] no, [1,2,3,4,5,6] yes (see below)
500 ZREF : level of no motion (in meters) for dynamic height computation
3 NDAY : number of days to OA. Line below, OADAY(1:NDAY)
5,15,25
1 NLEV : number of depth levels to OA. Line below, depth(s) in meters:
0
orig_hdat.hdat
orig_hdat.hdat
orig_hdat1.nc
orig_hdat.echo
!
!=====================================================================
! Copyright (c) 2000 Rutgers University ===
!=====================================================================
!
*** Above FILENAMES:
1st line: input hydrographic data in HDAT format.
2nd line: input climatological data in HDAT format.
3rd line: input GRID NetCDF file, if any.
4th line: output (NetCDF) objectively analyzed field(s)
5th line: output (ascii) objective analysis information
*** 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:
[ 1] pressure (depth)
[ 2] temperature
[ 3] salinity
[ 4] dynamic height
[ 5] potential temperature
[ 6] density
[ 7] sigma-t
[ 8] sound speed
[ 9] surface net heat flux
[10] surface freswater flux (E-P)
[11] precipitation rate
[12] evaporation rate
[13] solar shortwave radiation
[14] surface u-momentum stress
[15] surface v-momentum stress
[16] sea surface temperature (SST)
[17] surface net heat flux sensitivity to SST (dQ/dSST)
[18] bathymetry
[19] sea surface height
[20] incoming longwave radiation
[21] surface air temperature
[22] surface dew-point temperature
[23] surface u-wind
[24] surface v-wind
[25] relative humidity
[26] cloud fraction
[27] surface air pressure
[28] surface u-velocity
[29] surface v-velocity
*** ICORR: analytical correlation option:
[1] Rotated or non-rotated, time displaced or not,
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 local bottom of
the hydrographic profile to the surface
[5] integrate specific volume anomaly from the surface to the
local bottom of the hydrographic profile.
[6] same as [3] but if the local depth is shallower than
ZREF, only integrate from the local depth to the surface
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.