mean_grain_size.f

Go to the documentation of this file.

!                   **************************
                    SUBROUTINE mean_grain_size
!                   **************************
!
!
!***********************************************************************
! SISYPHE   V6P1                                   21/07/2011
!***********************************************************************
!
!brief    GEOMETRIC MEAN GRAIN SIZES OF ACTIVE-LAYER AND UNDER-LAYER.
!
!history  BUI MINH DUC
!+        **/11/2002
!+        V6P0
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE declarations_sisyphe
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER I,J
!
!-----------------------------------------------------------------------
!
!     UNLADM IS NEEDED FOR HUNZIKER
!
      DO j=1,npoin
        acladm%R(j) = 0.d0
        unladm%R(j) = 0.d0
        IF(nsicla.GT.1) THEN
          DO i=1,nsicla
            IF(avail(j,1,i).GT.0.d0) THEN
              acladm%R(j) = acladm%R(j) + fdm(i)*avail(j,1,i)
              unladm%R(j) = unladm%R(j) + fdm(i)*avail(j,2,i)
            ENDIF
          ENDDO
        ENDIF
        IF(acladm%R(j).LE.0.d0) acladm%R(j) = fdm(1)
        IF(unladm%R(j).LE.0.d0) unladm%R(j) = acladm%R(j)
      ENDDO

      CALL user_mean_grain_size
!
!-----------------------------------------------------------------------
!
      RETURN
      END