Hi kwerner,

Thanks for your reply. I looked at the code in this file, but I think what I want is the standard deviation at each eta level for the horizontal wind speed, and w velocity. It appears the above climate diagnostics uses 2d arrays only (ij) since it is at a specific level?

My goal is to output a value for the standard deviation of the horizontal wind speed, and standard deviation of the w velocity for each grid cell at all levels. (So instead of getting ('u', 'v', 'w',...) WRF output at each model level, I would get ('u', 'v', 'w', 'standard deviation of horizontal wind speed', 'standard deviation of w velocity',...).

My understanding is that the standard deviation would be based on all the individual time step values within each WRF history interval. Would/should I edit a different Fortran file to do this? Thanks.