1 |
adcroft |
1.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 |
|
|
|