Bug in TOP_SRC/TRP/trdmld_trc.F90 Allocation problem.

[molines]

• From Aurelie Albert : bug in the routine trd_mld_bio_zint in module TOP_SRC/TRP/trdmld_trc.F90 : No deallocation at the end of the routine but a new allocation ==> memory leakage and eventual crash on memory fault...

.....
   SUBROUTINE trd_mld_bio_zint( ptrc_trdmld, ktrd )
      !!----------------------------------------------------------------------
      !!                  ***  ROUTINE trd_mld_bio_zint  ***
      !!
      !! ** Purpose :   Compute the vertical average of the 3D fields given as arguments
      !!                to the subroutine. This vertical average is performed from ocean
      !!                surface down to a chosen control surface.
      !!
      !! ** Method/usage :
      !!      The control surface can be either a mixed layer depth (time varying)
      !!      or a fixed surface (jk level or bowl).
      !!      Choose control surface with nctls in namelist NAMTRD :
      !!        nctls_trc = 0  : use mixed layer with density criterion
      !!        nctls_trc = 1  : read index from file 'ctlsurf_idx'
      !!        nctls_trc > 1  : use fixed level surface jk = nctls_trc
      !!      Note: in the remainder of the routine, the volume between the
      !!            surface and the control surface is called "mixed-layer"
      !!----------------------------------------------------------------------
      !!
      INTEGER                         , INTENT(in) ::   ktrd          ! bio trend index
      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) ::   ptrc_trdmld   ! passive trc trend
#if defined key_pisces_reduced
      !
      INTEGER ::   ji, jj, jk, isum
      REAL(wp), POINTER, DIMENSION(:,:) :: zvlmsk
      !!----------------------------------------------------------------------

      CALL wrk_alloc( jpi, jpj, zvlmsk )

..... core of the routine

      CALL wrk_alloc( jpi, jpj, zvlmsk )
#endif
      !
   END SUBROUTINE trd_mld_bio_zint

• found in 3.4, 3.5beta and 3.6...
• obvious fix ; --> CALL wrk_dealloc( jpi, jpj, zvlmsk )