pkg/flt/flt_down.F

C $Header: /u/gcmpack/MITgcm/pkg/flt/flt_down.F,v 1.3 2009/01/04 00:58:23 jmc Exp $
C $Name:  $

#include "FLT_OPTIONS.h"

      SUBROUTINE FLT_DOWN (
     I                      myTime, myIter, myThid )

C     ==================================================================
C     SUBROUTINE FLT_DOWN
C     ==================================================================
C     o This routine moves particles vertical from the surface to
C       the target depth
C     ==================================================================

C     !USES:
      IMPLICIT NONE

C     == global variables ==
#include "EEPARAMS.h"
#include "SIZE.h"
#include "FLT.h"

C     == routine arguments ==
      _RL myTime
      INTEGER myIter, myThid

C     == local variables ==
      INTEGER bi, bj
      INTEGER ip

C     == end of interface ==

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)

        DO ip=1,npart_tile(bi,bj)

         IF (
     &  (                             myTime.GE.tstart(ip,bi,bj))
     &   .AND.
     &  (  tend(ip,bi,bj).EQ.-1. .OR. myTime.LE.  tend(ip,bi,bj))
     &   .AND.
     &  (kpart(ip,bi,bj) .LE. flt_surf)
     &   .AND.
     &  (iup(ip,bi,bj)   .GT. 0.)
     &      ) THEN

c         IF (tstart(ip,bi,bj) .NE. -1.              .AND.
c     &      myTime    .GE. tstart(ip,bi,bj) .AND.
c     &      myTime    .LE.   tend(ip,bi,bj) .AND.
c     &      kpart(ip,bi,bj) .EQ. flt_surf .AND.
c     &      iup(ip,bi,bj)   .GT. 0.) THEN

C     move float to the target depth

           IF (mod(myTime,iup(ip,bi,bj)).EQ.itop(ip,bi,bj))
     &      kpart(ip,bi,bj) = kfloat(ip,bi,bj)

         ENDIF
        ENDDO

       ENDDO
      ENDDO

      RETURN
      END

1	jmc	1.4	C $Header: /u/gcmpack/MITgcm/pkg/flt/flt_down.F,v 1.3 2009/01/04 00:58:23 jmc Exp $
2	jmc	1.2	C $Name: $
3	adcroft	1.1
4	jmc	1.2	#include "FLT_OPTIONS.h"
5	adcroft	1.1
6	jmc	1.3	SUBROUTINE FLT_DOWN (
7			I myTime, myIter, myThid )
8	adcroft	1.1
9	jmc	1.3	C ==================================================================
10	jmc	1.4	C SUBROUTINE FLT_DOWN
11	jmc	1.3	C ==================================================================
12			C o This routine moves particles vertical from the surface to
13			C the target depth
14			C ==================================================================
15	jmc	1.4
16			C !USES:
17			IMPLICIT NONE
18	jmc	1.3
19			C == global variables ==
20	adcroft	1.1	#include "EEPARAMS.h"
21			#include "SIZE.h"
22			#include "FLT.h"
23
24	jmc	1.3	C == routine arguments ==
25			_RL myTime
26			INTEGER myIter, myThid
27	adcroft	1.1
28	jmc	1.3	C == local variables ==
29			INTEGER bi, bj
30			INTEGER ip
31	adcroft	1.1
32	jmc	1.3	C == end of interface ==
33	adcroft	1.1
34			DO bj=myByLo(myThid),myByHi(myThid)
35	jmc	1.4	DO bi=myBxLo(myThid),myBxHi(myThid)
36	adcroft	1.1
37	jmc	1.4	DO ip=1,npart_tile(bi,bj)
38	adcroft	1.1
39	jmc	1.3	IF (
40			& ( myTime.GE.tstart(ip,bi,bj))
41			& .AND.
42			& ( tend(ip,bi,bj).EQ.-1. .OR. myTime.LE. tend(ip,bi,bj))
43			& .AND.
44			& (kpart(ip,bi,bj) .LE. flt_surf)
45			& .AND.
46			& (iup(ip,bi,bj) .GT. 0.)
47			& ) THEN
48
49			c IF (tstart(ip,bi,bj) .NE. -1. .AND.
50			c & myTime .GE. tstart(ip,bi,bj) .AND.
51			c & myTime .LE. tend(ip,bi,bj) .AND.
52			c & kpart(ip,bi,bj) .EQ. flt_surf .AND.
53			c & iup(ip,bi,bj) .GT. 0.) THEN
54
55			C move float to the target depth
56
57	jmc	1.4	IF (mod(myTime,iup(ip,bi,bj)).EQ.itop(ip,bi,bj))
58	adcroft	1.1	& kpart(ip,bi,bj) = kfloat(ip,bi,bj)
59
60	jmc	1.3	ENDIF
61	jmc	1.4	ENDDO
62	adcroft	1.1
63	jmc	1.4	ENDDO
64	adcroft	1.1	ENDDO
65
66	jmc	1.3	RETURN
67			END
68	adcroft	1.1