Changed Smagorinsky length scale to harmonic mean of cell area, e.g., lenA2=(2*rA(i,j,bi,bj)/(dxF(I,J,bi,bj)+dyF(I,J,bi,bj)))**2
Changed Smagorinsky length scale to harmonic mean of cell area, e.g.,
lenA2=(2*rA(i,j,bi,bj)/(dxF(I,J,bi,bj)+dyF(I,J,bi,bj)))**2

Should be more robust for anisotropic grids.

altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h,altMITgcm/MITgcm66h