bug in the outputs from dia_ar5_hst

in diaar5.F90 (dia_ar5_hst), it seems that the surface advective and diffusive fluxes are counted twice. Here:

      z2d(:,:) = puflx(:,:,1)
      DO_3D( 0, 0, 0, 0, 1, jpkm1 )
         z2d(ji,jj) = z2d(ji,jj) + puflx(ji,jj,jk)
      END_3D

Should be:

      z2d(:,:) = 0._wp
      DO_3D( 0, 0, 0, 0, 1, jpkm1 )
         z2d(ji,jj) = z2d(ji,jj) + puflx(ji,jj,jk)
      END_3D

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information

Admin message