matrix_example/code/ptracers_forcing_surf.F

C $Header: /u/gcmpack/MITgcm/verification/matrix_example/code/ptracers_forcing_surf.F,v 1.1 2005/04/18 20:14:55 spk Exp $
C $Name:  $

#include "PTRACERS_OPTIONS.h"

CBOP
C !ROUTINE: PTRACERS_FORCING_SURF

C !INTERFACE: ==========================================================
      SUBROUTINE PTRACERS_FORCING_SURF(
     I                            bi, bj, iMin, iMax, jMin, jMax,
     I                            myTime,myIter,myThid )

C !DESCRIPTION:
C     Precomputes surface forcing term for pkg/ptracers.
C     Precomputation is needed because of non-local KPP transport term,
C     routine KPP_TRANSPORT_PTR.

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PTRACERS_SIZE.h"
#include "PTRACERS.h"
#include "PARAMS.h"
#include "FFIELDS.h"
c #include "DYNVARS.h"
#include "GRID.h"

C !INPUT PARAMETERS: ===================================================
C  bi,bj                :: tile indices
C  myTime               :: model time
C  myIter               :: time-step number
C  myThid               :: thread number
      INTEGER bi, bj, iMin, iMax, jMin, jMax
      _RL myTime
      INTEGER myIter
      INTEGER myThid

#ifdef ALLOW_PTRACERS

C !LOCAL VARIABLES: ====================================================
C  i,j                  :: loop indices
      INTEGER i, j, iTracer
CEOP

      IF (.NOT.useMATRIX) THEN
C Example of how to add forcing at the surface
       DO iTracer=1,PTRACERS_numInUse
          DO j = jMin, jMax
             DO i = iMin, iMax
                surfaceForcingPtr(i,j,bi,bj,iTracer) =
     &               surfaceForcingS(i,j,bi,bj)
             ENDDO
          ENDDO
       ENDDO
      ENDIF

#endif /* ALLOW_PTRACERS */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/verification/matrix_example/code/ptracers_forcing_surf.F,v 1.1 2005/04/18 20:14:55 spk Exp $
2	C $Name: $
3
4	#include "PTRACERS_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: PTRACERS_FORCING_SURF
8
9	C !INTERFACE: ==========================================================
10	SUBROUTINE PTRACERS_FORCING_SURF(
11	I bi, bj, iMin, iMax, jMin, jMax,
12	I myTime,myIter,myThid )
13
14	C !DESCRIPTION:
15	C Precomputes surface forcing term for pkg/ptracers.
16	C Precomputation is needed because of non-local KPP transport term,
17	C routine KPP_TRANSPORT_PTR.
18
19	C !USES: ===============================================================
20	IMPLICIT NONE
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PTRACERS_SIZE.h"
24	#include "PTRACERS.h"
25	#include "PARAMS.h"
26	#include "FFIELDS.h"
27	c #include "DYNVARS.h"
28	#include "GRID.h"
29
30	C !INPUT PARAMETERS: ===================================================
31	C bi,bj :: tile indices
32	C myTime :: model time
33	C myIter :: time-step number
34	C myThid :: thread number
35	INTEGER bi, bj, iMin, iMax, jMin, jMax
36	_RL myTime
37	INTEGER myIter
38	INTEGER myThid
39
40	#ifdef ALLOW_PTRACERS
41
42	C !LOCAL VARIABLES: ====================================================
43	C i,j :: loop indices
44	INTEGER i, j, iTracer
45	CEOP
46
47	IF (.NOT.useMATRIX) THEN
48	C Example of how to add forcing at the surface
49	DO iTracer=1,PTRACERS_numInUse
50	DO j = jMin, jMax
51	DO i = iMin, iMax
52	surfaceForcingPtr(i,j,bi,bj,iTracer) =
53	& surfaceForcingS(i,j,bi,bj)
54	ENDDO
55	ENDDO
56	ENDDO
57	ENDIF
58
59	#endif /* ALLOW_PTRACERS */
60
61	RETURN
62	END