!!---------------------------------------------------------------------
!!
!!                        trcini.npzd.h
!!                        ------------- 
!!  purpose :
!!  ---------
!!     special initialisation for NPZD model
!!
!!   History :
!!   ---------
!!      original    : 99-09 (M. Levy) 
!!      additions   : 00-12 (E. Kestenare) add sediment computations
!!      03-2005 O.Aumont and A. El Moussaoui F90 
!! ==============================================================================  
 
!!---------------------------------------------------------------------
!!  TOP 1.0,  LOCEAN-IPSL (2005)
!! $Header$
!! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
!!---------------------------------------------------------------------
!! local variables
!! ---------------
      REAL ztest

!! 2. initialization of field for optical model
!! --------------------------------------------

          xze(:,:)=5.
          xpar(:,:,:)=0.

!! 3. initialization for passive tracer remineralisation-damping  array
!! -------------------------------------------------------------------------
 
      DO jn=1,jptra
          remdmp(:,jn)=tminr
      END DO

      IF(lwp) THEN
          WRITE(numout,*) ' '
          WRITE(numout,*) ' trcini: compute remineralisation-damping  '
          WRITE(numout,*) '         arrays for tracers'
      ENDIF

!! 4. sediments: no martin's remineralisation profile
!! -------------------------------------------------------------------------

      dminl = 0.
      dmin3 = 0. 

!! CALCUL DU MASK DE COTE
!!
        cmask=0.
        do ji=2,jpi-1
          do jj=2,jpj-1
            if (tmask(ji,jj,1) == 1) then
             ztest=tmask(ji+1,jj,1)*tmask(ji-1,jj,1)*tmask(ji,jj+1,1)
     .             *tmask(ji,jj-1,1)
             if (ztest == 0) cmask(ji,jj)=1.
             endif
          end do
        end do

        cmask(1,:)=cmask(jpi-1,:)
        cmask(jpi,:)=cmask(2,:)

!!     CALCUL DE LA SURFACE COTIERE
!!
         do ji=2,jpi-1
          do jj=2,jpj-1
          areacot=areacot+e1t(ji,jj)*e2t(ji,jj)*cmask(ji,jj)
          end do
         end do