Variables for NAURU99 experiments : date = year-month-day hour = decimal hour lat = latitude North ( - = South ) lon = longitude East ( - = West ) ws_s = true wind speed , ETL sonic ws_a = true wind speed, imet provane anemometer seaT_s = sea snake temperature, ETL 0.05 m depth seaT_t = tsg water temperature, 5 m depth airT_e = air temperature, ETL airT_i = imet air temperature qa_e = air specific humidity, ETL qa_i = imet air specific humidity swr = downward solar flux, ETL lwr = downward IR flux, ETL wdir_e = true wind direction, ETL sonic wdir_i = true wind direction, clockwise rel north, imet rain = rainrate, ETL STI optical rain gauge, uncorrected shf_c = sensible heat flux, covariance, ETL sonic anemometer shf_i = sensible heat flux, ID, ETL sonic anemometer shf_b = bulk sensible heat flux lhf_c = latent heat flux, covariance lhf_i = latent heat flux, ID lhf_b = bulk latent heat flux, includes Webb et al. correction hlwebb = correction to measured latent heat flux, Webb et al. head = ship heading, deg. clockwise rel north, SCS laser ring gyro reldir = relative wind direction, clockwise rel ship's bow, ETL sonic urel = relative wind speed, ETL sog = speed over ground, SCS gps cog = course over ground, SCS gps qse = sea surface specific humidity, from snake qs_tsg = imet bulk water specific humidity taucx = covariance streamwise stress, ETL sonic anemometer taucy = covariance cross-stream stress, ETL sonic anemometer tauib = ID streamwise stress, ETL sonic anemometer taub = bulk wind stress along mean wind sigoph = standard deviation of ophir fast hygrometer clear channel tiltx = flow tilt at ETL sonic anemometer, earth frame ushp = doppler log, SCS J = ship plume contamination index Jm = ship maneuver index hrain = rain heat flux, Gosnell et al. ct = ct^2 cq = cq^2 cu = cu^2 cw = cw^2 zu_etl = height of mean wind sensor, 17.7 m zt_etl = height of mean air temperature sensor, 15.5 m zq_etl = height of mean air humidity sensor, 15.5 m The first release of turbulent and bulk fluxes for Jasmine, Nauru99, Kwajex, and Moorings cruises. This document is the Readme for crus_hr.txt and crus_10.txt files where 'crus'= is described below. The _hr refers to hourly averages; the _10 to 10-minute averages. crus Name Jul. Day Dates Hours Nom. Lat/Lon naur nauru99 166-199 6/15-7/18 794 0.5 S 167 E The data files crus_hr.txt and contain measurements of turbulent and radiative fluxes plus bulk meteorological variables from the R/V Ronald H. Brown in various legs of the 'around the world cruise' of 1999. The files also contain bulk estimates of the turbulent fluxes computed using a recently updated version (2.6) of the COARE flux algorithm (see ftp://ftp.etl.noaa.gov/et7/users/cfairall/bulkalg/ for documentation). Most quantities have been subjected to one round of intercomparison/calibration scrutiny; however, there may be future modifications based on accounting for other sources of data and revised calibrations. Both direct (covariance) and inertial-dissipation (ID) turbulent flux calculations are included in this present data. The data in this file comes from three sources: The ETL motion-correct flux package [sonic anemometer acquired at 20.83 Hz, fast hygrometer acquired at 20 Hz, and 6-component motion measurements acquired at 10 Hz), the ships SCS system (acquired at 2 sec intervals), and the ETL mean measurement systems (sampled at 10 sec and averaged to 1 min). The sonic is 5 channels of data; the SCS file is 13 channels, and the ETL mean system is 16 channels. A series of programs are run that read these data files, decode them, compute covariance, variances, spectra, etc. at 10-min time resolution. A set of 10-min process files are written for each cruise. One particular file, crus_da.dat consists of 164 columns of 10-min data including the turbulent and mean variables used here. A final program reads the crus_da.dat files, applies various corrections, computes the ID and bulk fluxes and the data quality indices, and writes the crus_10.txt files. This program also contains routines that average the data to fixed one-hr time intervals, computes the ID and bulk fluxes from one-hr means, and writes the crus_hr.txt files. Note, that in the case of the turbulent variables (covariances, variances, structure function parameters) only 10-min values that pass the data indicator criteria are used in the one-hr average. Further experimental details are as follows: True wind speed is computed from the sonic anemometer using the ship's Laser ring gyro and the odec doppler log; thus, it is interpreted as the speed relative to the water. Some modest eyeball flow distortion corrections to the relative wind components have been used in an attempt to reduce the transitions when stopping for stations. Air temperature and humidity from the ETL (aspirated Vaisala HMP-235) were carefully compared with the CSIRO psychrometer values obtained during the JASMINE experiment and with handheld wet-dry bulk psychrometer readings after that. The ETL humidities were increased about 3% based on these intercomparisons; the air temperatures were left unchanged. Final values are good to about 1% RH and 0.2 C. At various times during some cruises the seasnake was removed from the water, in which case the ship's thermosalinograph, corrected for any warm layer effects, was substituted. Sea surface temperature measurements often showed contamination by the ship during stops and maneuvers at stations. The sensor was airborne a lot when the ship was underway in the stronger wind/rougher sea days and the accuracy under those conditions is uncertain. Longwave flux was obtained from 2 Eppley PIR units, logged and computed as per Fairall et al. Jtech, 1998. Shortwave flux was obtained from 2 Eppley PSP units. The original values were increased 5% based on intercomparisons with the BNL PRP. The rainrate was obtained from the ETL STI optical raingauge. Note, this is model org-815 - not the mini-org used on COARE. No correction for cosine response has been attempted. There was a +0.1 to 0.4 mm/hr bias that has been removed. Bulk estimates of air sea fluxes were computed using the COARE bulk algorithm version 2.6, which differs from the published version in only minor details. Turbulent fluxes were computed from motion-corrected time series of fast sensors. Covariance latent heat fluxes were obtained by cross-correlating the motion-corrected vertical velocity with fast humidity fluctuations from an OPHIR IR hygrometer. The turbulent fluctuations from the OPHIR were scaled by the ratio of the OPHIR mean humidity to the ETL Vaisala mean humidity. An additional scaling factor (4.5% increase) was used to account for the physical separation of the OPHIR and the sonic. Using a right-handed coordinate system with x boward, y to the port, and z up, the displacement vector from the OPHIR to the sonic is (1.0, 0.69, 1.57) m. We used Kristensen et al. (J. Atmos. Oceanic Tech., 14, 814-821, 1997) to estimate the correlation loss. On Jasmine and Nauru99 a closed-path Licor 6262 fast humidity/CO2 sensor was also used. The values from this unit were scaled to agree in with the corrected OPHIR values on a daily basis. This is necessary because some fraction (10-15%) of the correlation is lost in the Licor sampling tube. Because the IR hygrometers detect water vapor mass concentration (rho_v in kg/m^3), their water vapor - velocity correlations must be corrected as per Webb et al Hlatent = Le + hl_webb The values given for covariance and ID latent heat fluxes in the file are Le. Values for hl_webb are included in column 34. This should be applied to the covariance and ID values. It is already included in the bulk values given here. Sensible heat flux was computed from vertical velocity - sonic temperature covariance. The humidity contribution to sonic temperature was removed using the bulk latent heat flux. Simple data quality indicators have been use to edit the turbulence data. At each 10-min time step, values for ID turbulence variables were computed if the following criteria were met: jj=find((reldir<90 | reldir>270) & sig_h<8 &sig_u<2 &sp2<2 & org<5 & sqrt(ww)./ugw<.55+.002*U.^2 & sqrt(vv)./ugu<1.1); reldir is the relative wind direction sig_h the standard deviation of ship heading (deg) sig_u the standard deviation of the ship speed sp2 the standard deviation of cross-ship motion corrections org the rainrate ww the vertical velocity variance vv the cross stream velocity variance Otherwise, the ID variables were set to NaN (not a number). In processing the 10-min data to one-hr averages, only the jj rows were used in averaging the turbulence variables. If there were no valid values in the 1-hr interval, the turbulence variables were set to NaN. The criteria given above were subdivided to be approximately compatible with indices used in the past: J=ones(length(jdy),1); ii=find(sqrt(ww)./ugw<.55+.002*U.^2 & sqrt(vv)./ugu<1.1); J(ii)=0; jm=3*ones(length(jdy),1); ij=find(sig_h<8 & sig_u<2 & sp2<2); jm(ij)=sig_u(ij); *A value of J=0 implies no ship contamination. *A value of jm<3 implies no significant maneuver during the average. *Sigoph is an index of salt or rain contamination on the fast hygrometer (OPHIR) optics. Values for fluxes begin to be affected when sigoph exceeds 20 although a threshold of 50 gives acceptable data. *Turbulent fluxes are computed by converting the anemometer 3-component velocities to fixed earth coordinates, correcting the fast time series for ship motion, and re-setting the coordinate system normal to the 10-min mean flow through one rotation about the original vertical and one tilt. The variable tiltx gives the tilt used for the computation. Experience shows that tilts greater than about 10 deg give questionable fluxes. *Nauru99: Directly downwind of Nauru Island JD 185.12-188.00 and a few other periods.