1 |
C $Header: $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "FLT_CPPOPTIONS.h" |
5 |
|
6 |
subroutine flt_down ( |
7 |
I myCurrentIter, |
8 |
I myCurrentTime, |
9 |
I myThid |
10 |
& ) |
11 |
|
12 |
c ================================================================== |
13 |
c SUBROUTINE flt_down |
14 |
c ================================================================== |
15 |
c |
16 |
c o This routine moves particles vertical from the surface to |
17 |
c the target depth |
18 |
c |
19 |
c ================================================================== |
20 |
c SUBROUTINE flt_down |
21 |
c ================================================================== |
22 |
|
23 |
c == global variables == |
24 |
|
25 |
#include "EEPARAMS.h" |
26 |
#include "SIZE.h" |
27 |
#include "FLT.h" |
28 |
|
29 |
c == routine arguments == |
30 |
|
31 |
INTEGER myCurrentIter, myThid |
32 |
_RL myCurrentTime |
33 |
INTEGER bi, bj |
34 |
|
35 |
c == local variables == |
36 |
|
37 |
integer ip |
38 |
|
39 |
c == end of interface == |
40 |
|
41 |
|
42 |
DO bj=myByLo(myThid),myByHi(myThid) |
43 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
44 |
|
45 |
do ip=1,npart_tile(bi,bj) |
46 |
|
47 |
if( |
48 |
& ( myCurrentTime.ge.tstart(ip,bi,bj)) |
49 |
& .and. |
50 |
& ( tend(ip,bi,bj).eq.-1. .or. myCurrentTime.le. tend(ip,bi,bj)) |
51 |
& .and. |
52 |
& (kpart(ip,bi,bj) .le. flt_surf) |
53 |
& .and. |
54 |
& (iup(ip,bi,bj) .gt. 0.) |
55 |
& ) then |
56 |
|
57 |
c if(tstart(ip,bi,bj) .ne. -1. .and. |
58 |
c & myCurrentTime .ge. tstart(ip,bi,bj) .and. |
59 |
c & myCurrentTime .le. tend(ip,bi,bj) .and. |
60 |
c & kpart(ip,bi,bj) .eq. flt_surf .and. |
61 |
c & iup(ip,bi,bj) .gt. 0.) then |
62 |
|
63 |
c move float to the target depth |
64 |
c |
65 |
if(mod(myCurrentTime,iup(ip,bi,bj)).eq.itop(ip,bi,bj)) |
66 |
& kpart(ip,bi,bj) = kfloat(ip,bi,bj) |
67 |
|
68 |
endif |
69 |
enddo |
70 |
|
71 |
ENDDO |
72 |
ENDDO |
73 |
|
74 |
return |
75 |
end |
76 |
|