1 |
C $Header: $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "PTRACERS_OPTIONS.h" |
5 |
|
6 |
CBOP |
7 |
C !ROUTINE: PTRACERS_DEBUG |
8 |
|
9 |
C !INTERFACE: ========================================================== |
10 |
SUBROUTINE PTRACERS_DEBUG( |
11 |
I myThid ) |
12 |
|
13 |
C !DESCRIPTION: |
14 |
C Calculates tendancy for passive tracers and integrates forward |
15 |
C in time. |
16 |
|
17 |
C !USES: =============================================================== |
18 |
IMPLICIT NONE |
19 |
#include "SIZE.h" |
20 |
#include "EEPARAMS.h" |
21 |
#include "PTRACERS.h" |
22 |
|
23 |
C !INPUT PARAMETERS: =================================================== |
24 |
C myThid :: thread number |
25 |
INTEGER myThid |
26 |
|
27 |
C !OUTPUT PARAMETERS: ================================================== |
28 |
C none |
29 |
|
30 |
#ifdef ALLOW_PTRACERS |
31 |
|
32 |
C !LOCAL VARIABLES: ==================================================== |
33 |
C i,j,k,bi,bj,iTracer :: loop indices |
34 |
INTEGER i,j,iTracer |
35 |
CHARACTER*(9) strng |
36 |
CEOP |
37 |
|
38 |
C Loop over tracers |
39 |
DO iTracer=1,PTRACERS_num |
40 |
|
41 |
WRITE(strng(1:9),'(a7,i2.2)') 'PTRACER',iTracer |
42 |
CALL DEBUG_STATS_RL(Nr, |
43 |
& pTracer(1-Olx,1-Oly,1,1,1,iTracer),strng,myThid) |
44 |
WRITE(strng(1:9),'(a7,i2.2)') 'GTRACER',iTracer |
45 |
CALL DEBUG_STATS_RL(Nr, |
46 |
& gPtr(1-Olx,1-Oly,1,1,1,iTracer),strng,myThid) |
47 |
WRITE(strng(1:9),'(a7,i2.2)') 'GNm1TRC',iTracer |
48 |
CALL DEBUG_STATS_RL(Nr, |
49 |
& gPtrNm1(1-Olx,1-Oly,1,1,1,iTracer),strng,myThid) |
50 |
|
51 |
C end of tracer loop |
52 |
ENDDO |
53 |
|
54 |
#endif /* ALLOW_PTRACERS */ |
55 |
|
56 |
RETURN |
57 |
END |