Hi,

I run the composet "CIAF" in the CESM1.2. This composet has the interannual forcing COREv2 datasets.

I found the wind stress variables "TAUX/TAUY" are contained in the POP2 output history file. So I want to know how the wind stress be calculated from the wind speed. At first, I calculate the wind stress with empirical formula：

zonal windstress = rho*Cd*sqrt(u*u+v*v)*u

meridional windstress = rho*Cd*sqrt(u*u+v*v)*v

where rho is the air density, Cd is the drag coefficient, u and v are zonal and meridional wind speed, respectively. Here, rho and Cd are constants.

However, the result is not the same as that in output file. Can anyone knows the code calculating the wind stress?

Thanks.

Best regards.