Does reanalysis use a constant drag coefficient to compute the surface wind stress?
[12/11/96] Drag coefficient is a function of bulk Richardson number and is not constant. • How does NCEP/NCAR reanalaysis compute the potential vorticity on isentropic surfaces?[12/06/96] First the static stability N**2=g/T*(dT/dz+g/cp) is computed on model levels. Then the winds, temperature and static stability are interpolated to isentropic surfaces linearly in log(theta). (Outside the model domain, the fields are held constant for now and will later be compressed out of the final product.) Then the absolute vorticity zeta is spectrally computed on the isentropic surfaces, where the shortest wavelength in the spectral domain is about 4 grid lengths. (The vorticity is computed in the T36 spectral domain for the 2.5×2.5 degree grid.) The density rho=(T/theta)**(cp/R)*p0/(R*T) is also computed at this time directly on the isentropic surfaces. Finally, the NCEP potential vorticity is computed as zeta*N**2/(g*rho). Thus the units of NCEP PV are m**2/s/kg, which is different from the usua