pkg/ptracers/ptracers_write_state.F

C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_write_state.F,v 1.1 2002/03/04 19:01:29 adcroft Exp $
C $Name:  $

#include "PTRACERS_OPTIONS.h"

CBOP
C !ROUTINE: PTRACERS_WRITE_STATE

C !INTERFACE: ==========================================================
      SUBROUTINE PTRACERS_WRITE_STATE( myIter,myTime,myThid )

C !DESCRIPTION:
C     Calculates tendancy for passive tracers and integrates forward
C     in time.

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PTRACERS.h"

C !INPUT PARAMETERS: ===================================================
C  myIter               :: time-step number
C  myTime               :: model time
C  myThid               :: thread number
      INTEGER myIter
      _RL myTime
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  none

#ifdef ALLOW_PTRACERS

C !LOCAL VARIABLES: ====================================================
C  iTracer              :: loop indices
C  suff                 :: character buffer for creating suffix in filename
C  iRec                 :: record number
      INTEGER iTracer
      CHARACTER*(MAX_LEN_MBUF) suff
      INTEGER iRec
CEOP

C Loop over tracers
      DO iTracer=1,PTRACERS_numInUse

#ifdef MULTIPLE_RECORD_STATE_FILES
C Write each snap-shot as a record in a time-history file
        IF (PTRACERS_useRecords) THEN
C        Write each passive tracer as a record
         WRITE(suff,'(A)') 'PTRACERS'
         iRec = int ( (myCurrentTime-startTime) / dumpFreq +1.5 )
         iRec = (iRec-1)*PTRACERS_num + iTracer + 1
        ELSE
C        Write each passive tracer in its own file
         WRITE(pref,'(A7,I2.2)') 'PTRACER',iTracer
         iRec = int ( (myCurrentTime-startTime) / dumpFreq +1.5 )
        ENDIF
#else
C Write each snap-shot as a new file with a 10-digit suffix to
C indicate time-step number
        IF (PTRACERS_useRecords) THEN
C        Write each passive tracer as a record
         WRITE(suff,'(A,I10.10)') 'PTRACERS',myIter
         iRec=iTracer
        ELSE
C        Write each passive tracer in its own file
         WRITE(suff,'(A7,I2.2,A1,I10.10)') 'PTRACER',iTracer,'.',myIter
         iRec=1
        ENDIF
#endif /* MULTIPLE_RECORD_STATE_FILES */

        CALL WRITE_REC_XYZ_RL(suff,
     &      pTracer(1-Olx,1-Oly,1,1,1,iTracer),iRec,myIter,myThid)

C End of tracer loop
      ENDDO

#endif /* ALLOW_PTRACERS */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_write_state.F,v 1.1 2002/03/04 19:01:29 adcroft Exp $
2	C $Name: $
3
4	#include "PTRACERS_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: PTRACERS_WRITE_STATE
8
9	C !INTERFACE: ==========================================================
10	SUBROUTINE PTRACERS_WRITE_STATE( myIter,myTime,myThid )
11
12	C !DESCRIPTION:
13	C Calculates tendancy for passive tracers and integrates forward
14	C in time.
15
16	C !USES: ===============================================================
17	IMPLICIT NONE
18	#include "SIZE.h"
19	#include "EEPARAMS.h"
20	#include "PTRACERS.h"
21
22	C !INPUT PARAMETERS: ===================================================
23	C myIter :: time-step number
24	C myTime :: model time
25	C myThid :: thread number
26	INTEGER myIter
27	_RL myTime
28	INTEGER myThid
29
30	C !OUTPUT PARAMETERS: ==================================================
31	C none
32
33	#ifdef ALLOW_PTRACERS
34
35	C !LOCAL VARIABLES: ====================================================
36	C iTracer :: loop indices
37	C suff :: character buffer for creating suffix in filename
38	C iRec :: record number
39	INTEGER iTracer
40	CHARACTER*(MAX_LEN_MBUF) suff
41	INTEGER iRec
42	CEOP
43
44	C Loop over tracers
45	DO iTracer=1,PTRACERS_numInUse
46
47	#ifdef MULTIPLE_RECORD_STATE_FILES
48	C Write each snap-shot as a record in a time-history file
49	IF (PTRACERS_useRecords) THEN
50	C Write each passive tracer as a record
51	WRITE(suff,'(A)') 'PTRACERS'
52	iRec = int ( (myCurrentTime-startTime) / dumpFreq +1.5 )
53	iRec = (iRec-1)*PTRACERS_num + iTracer + 1
54	ELSE
55	C Write each passive tracer in its own file
56	WRITE(pref,'(A7,I2.2)') 'PTRACER',iTracer
57	iRec = int ( (myCurrentTime-startTime) / dumpFreq +1.5 )
58	ENDIF
59	#else
60	C Write each snap-shot as a new file with a 10-digit suffix to
61	C indicate time-step number
62	IF (PTRACERS_useRecords) THEN
63	C Write each passive tracer as a record
64	WRITE(suff,'(A,I10.10)') 'PTRACERS',myIter
65	iRec=iTracer
66	ELSE
67	C Write each passive tracer in its own file
68	WRITE(suff,'(A7,I2.2,A1,I10.10)') 'PTRACER',iTracer,'.',myIter
69	iRec=1
70	ENDIF
71	#endif /* MULTIPLE_RECORD_STATE_FILES */
72
73	CALL WRITE_REC_XYZ_RL(suff,
74	& pTracer(1-Olx,1-Oly,1,1,1,iTracer),iRec,myIter,myThid)
75
76	C End of tracer loop
77	ENDDO
78
79	#endif /* ALLOW_PTRACERS */
80
81	RETURN
82	END