C **********************************************************************
C *                                                                    *
C *                                SOFTWARE LICENSING                  *
C *                                                                    *
C *                This program is free software; you can redistribute *
C *                it and/or modify it under the terms of the GNU      *
C *                General Public License as published by the Free     *
C *                Software Foundation, either Version 2 of the        *
C *                license, or (at your option) any later version.     *
C *                                                                    *
C *                This program is distributed in the hope that it     *
C *                will be useful, but without any warranty; without   *
C *                even the implied warranty of merchantability or     *
C *                fitness for a particular purpose. See the GNU       *
C *                General Public License for more details.            *
C *                                                                    *
C *                A copy of the GNU General Public License is         *
C *                available at http://www.gnu.org/copyleft/gpl.html   *
C *                or by writing to the Free Software Foundation, Inc.,*
C *                59 Temple Place - Suite 330, Boston, MA 02111, USA. *
C *                                                                    *
C **********************************************************************

      SUBROUTINE REORDR
C
C  CHANGES AGE OF LAYERS AT END OF EACH DAY AND RESETS TAUMAX
C
C  REVISION DATE :  JUNE 8, 1995
C
C****************************************************************
C
      INCLUDE 'comdeck'
C
C  24 HOUR REORDERING
C
      LMAXM=LAYMAX-1
      CRAT=FTIME2(LMAXM)/FTIME2(LAYMAX)
      DO 10 J=2,JM-1
        DO 10 I=2,IM-1
C
          IF (IBMSK(I,J).LT.0.OR.IBMSK(I,J).EQ.1.OR.
     +                              FSM(I,J).EQ.0.0) GOTO 10
C
C  ADD MATERIAL TO 7 DAY OLD LAYER
C
          IF (LAYER(LMAXM,I,J).EQ.1) THEN
            IF (KSED.GT.1) THEN
              IF (LAYER(LAYMAX,I,J).EQ.1) THEN        
                PSED1(LAYMAX,I,J)=(PSED1(LMAXM,I,J)*TSED(LMAXM,I,J)
     +          +PSED1(LAYMAX,I,J)*TSED(LAYMAX,I,J))/
     +          (TSED(LAYMAX,I,J)+TSED(LMAXM,I,J))
                PSED2(LAYMAX,I,J)=1.-PSED1(LAYMAX,I,J)
              ELSE
                PSED1(LAYMAX,I,J)=PSED1(LMAXM,I,J)
                PSED2(LAYMAX,I,J)=PSED2(LMAXM,I,J)
              ENDIF
            ELSE
              PSED1(LAYMAX,I,J)=1.
              PSED2(LAYMAX,I,J)=0.
            ENDIF
            LAYER(LAYMAX,I,J)=1
            TSED(LAYMAX,I,J)=TSED(LAYMAX,I,J)+TSED(LMAXM,I,J)
C
CPL ADD THE FOLLOWING LINES PRAVI 09/08/99
C
            EBMAX(LAYMAX,I,J)=CRAT*EBMAX(LMAXM,I,J)
            EBTOT(LAYMAX,I,J)=CRAT*EBTOT(LMAXM,I,J)
            EBCUR(LAYMAX,I,J)=CRAT*EBCUR(LMAXM,I,J)
            TAUMAX(LAYMAX,I,J)=TAUMAX(LMAXM,I,J)
            TAUCUR(LAYMAX,I,J)=TAUCUR(LMAXM,I,J)
CPL 09/08/99            EBTOT(LAYMAX,I,J)=EBTOT(LAYMAX,I,J)-TSED(LMAXM,I,J)
CPL 09/08/99            IF (EBTOT(LAYMAX,I,J).LT.0.0) EBTOT(LAYMAX,I,J)=0.0
C
            PSED1(LMAXM,I,J)=0.0
            PSED2(LMAXM,I,J)=0.0
            TSED(LMAXM,I,J)=0.0
            LAYER(LMAXM,I,J)=0
            TAUMAX(LMAXM,I,J)=0.0
            EBTOT(LMAXM,I,J)=0.0
            TAUCUR(LMAXM,I,J)=0.0
            EBMAX(LMAXM,I,J)=0.0
            EBCUR(LMAXM,I,J)=0.0
          ENDIF
C
C  ADD MATERIAL TO YOUNGER LAYERS 
C
         DO 16 LL=LAYMAX-1,2,-1
          IF (LAYER(LL-1,I,J).EQ.1) THEN
            IF (KSED.GT.1) THEN
              IF (LAYER(LL,I,J).EQ.1) THEN        
                PSED1(LL,I,J)=(PSED1(LL-1,I,J)
     +          *TSED(LL-1,I,J)+PSED1(LL,I,J)
     +          *TSED(LL,I,J))/
     +          (TSED(LL,I,J)+TSED(LL-1,I,J))
                PSED2(LL,I,J)=1.-PSED1(LL,I,J)
              ELSE
                PSED1(LL,I,J)=PSED1(LL-1,I,J)
                PSED2(LL,I,J)=PSED2(LL-1,I,J)
              ENDIF
            ELSE
              PSED1(LL,I,J)=1.
              PSED2(LL,I,J)=0.
            ENDIF
            LAYER(LL,I,J)=1
            TSED(LL,I,J)=TSED(LL,I,J)+TSED(LL-1,I,J)
C
            CRAT=FTIME2(LL-1)/FTIME2(LL)
C
            EBMAX(LL,I,J)=CRAT*EBMAX(LL-1,I,J)
            EBTOT(LL,I,J)=CRAT*EBTOT(LL-1,I,J)
            EBCUR(LL,I,J)=CRAT*EBCUR(LL-1,I,J)
            TAUMAX(LL,I,J)=TAUMAX(LL-1,I,J)
            TAUCUR(LL,I,J)=TAUCUR(LL-1,I,J)
C
            PSED1(LL-1,I,J)=0.0
            PSED2(LL-1,I,J)=0.0
            TSED(LL-1,I,J)=0.0
            LAYER(LL-1,I,J)=0
            TAUMAX(LL-1,I,J)=0.0
            EBTOT(LL-1,I,J)=0.0
            TAUCUR(LL-1,I,J)=0.0
            EBMAX(LL-1,I,J)=0.0
            EBCUR(LL-1,I,J)=0.0
          ENDIF
 16      CONTINUE
C
          PSED1(1,I,J)=0.0
          PSED2(1,I,J)=0.0
          TSED(1,I,J)=0.0
          LAYER(1,I,J)=0
          TAUMAX(1,I,J)=0.0
          EBTOT(1,I,J)=0.0
          TAUCUR(1,I,J)=0.0
          EBMAX(1,I,J)=0.0
          EBCUR(1,I,J)=0.0
 10   CONTINUE
C
      N24CNT=0
C
      RETURN
      END