To compute simple flux terms wT, wq from gridded data in NetCDF format we can use the CDO “mul” operator:

cdo  mul  W.nc  T.nc  WT.nc
cdo  mul  W.nc  Q.nc  WQ.nc

The resulting variable in the output file will have the name of the first variable “W”. So we can then rename the variables using “chname“:

cdo  chname,W,wT  WT.nc  wT.nc
rm WT.nc