pkg/flt/flt_down.F

C $Header: /u/gcmpack/MITgcm/pkg/flt/flt_down.F,v 1.2 2008/12/03 01:32:35 jmc Exp $
C $Name:  $

#include "FLT_OPTIONS.h"

      SUBROUTINE FLT_DOWN (
     I                      myTime, myIter, myThid )

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

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