pkg/dic/alk_surfforcing.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/alk_surfforcing.F,v 1.12 2008/04/07 20:31:16 dfer Exp $
C $Name:  $

#include "DIC_OPTIONS.h"

CBOP
C !ROUTINE: ALK_SURFFORCING

C !INTERFACE: ==========================================================
      SUBROUTINE ALK_SURFFORCING( PTR_ALK , GALK,
     I           bi,bj,imin,imax,jmin,jmax,
     I           myIter,myTime,myThid)

C !DESCRIPTION:
C Calculate the alkalinity change due to freshwater flux
C Update alkalinity tendency term

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

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
C  myIter               :: current timestep
C  myTime               :: current time
C  bi,bj                :: tile indices
C  PTR_ALK              :: alkalinity field
      INTEGER myIter, myThid
      _RL myTime
      _RL  PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
      INTEGER iMin,iMax,jMin,jMax, bi, bj

C !OUTPUT PARAMETERS: ==================================================
C  GALK                 :: tendency term of alkalinity
      _RL  GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy)

#ifdef ALLOW_PTRACERS
#ifdef DIC_BIOTIC

#ifdef ALLOW_OLD_VIRTUALFLUX
C !LOCAL VARIABLES: ====================================================
C  i,j                  :: loop indices
      INTEGER i,j
CEOP

        DO j=jmin,jmax
          DO i=imin,imax
            IF (maskC(i,j,1,bi,bj).NE.0. _d 0) THEN
c calculate virtual flux
c EminusPforV = dS/dt*(1/Sglob)
C NOTE: Be very careful with signs here!
C Positive EminusPforV => loss of water to atmos and increase
C in salinity. Thus, also increase in other surface tracers
C (i.e. positive virtual flux into surface layer)
              GALK(i,j)=gsm_ALK*surfaceForcingS(i,j,bi,bj)*
     &                  recip_drF(1)/gsm_S
c OR
c let virtual flux be zero
c             GALK(i,j)=0.d0
c
            ELSE
              GALK(i,j)=0. _d 0
            ENDIF
          ENDDO
         ENDDO
#endif /* ALLOW_OLD_VIRTUALFLUX */

#endif
#endif
        RETURN
        END
1	jmc	1.13	C $Header: /u/gcmpack/MITgcm/pkg/dic/alk_surfforcing.F,v 1.12 2008/04/07 20:31:16 dfer Exp $
2	jmc	1.9	C $Name: $
3
4	edhill	1.2	#include "DIC_OPTIONS.h"
5	stephd	1.1
6	stephd	1.3	CBOP
7			C !ROUTINE: ALK_SURFFORCING
8
9			C !INTERFACE: ==========================================================
10	jmc	1.9	SUBROUTINE ALK_SURFFORCING( PTR_ALK , GALK,
11	stephd	1.1	I bi,bj,imin,imax,jmin,jmax,
12			I myIter,myTime,myThid)
13
14	stephd	1.3	C !DESCRIPTION:
15			C Calculate the alkalinity change due to freshwater flux
16			C Update alkalinity tendency term
17
18			C !USES: ===============================================================
19	stephd	1.1	IMPLICIT NONE
20			#include "SIZE.h"
21			#include "DYNVARS.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "GRID.h"
25			#include "FFIELDS.h"
26	dfer	1.11	#include "DIC_VARS.h"
27	stephd	1.1
28	stephd	1.3	C !INPUT PARAMETERS: ===================================================
29			C myThid :: thread number
30			C myIter :: current timestep
31			C myTime :: current time
32			C bi,bj :: tile indices
33			C PTR_ALK :: alkalinity field
34	stephd	1.1	INTEGER myIter, myThid
35			_RL myTime
36			_RL PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
37	stephd	1.3	INTEGER iMin,iMax,jMin,jMax, bi, bj
38
39			C !OUTPUT PARAMETERS: ==================================================
40			C GALK :: tendency term of alkalinity
41	stephd	1.1	_RL GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42
43			#ifdef ALLOW_PTRACERS
44			#ifdef DIC_BIOTIC
45	stephd	1.3
46	jmc	1.13	#ifdef ALLOW_OLD_VIRTUALFLUX
47	stephd	1.3	C !LOCAL VARIABLES: ====================================================
48			C i,j :: loop indices
49			INTEGER i,j
50			CEOP
51	stephd	1.1
52	stephd	1.7	DO j=jmin,jmax
53			DO i=imin,imax
54	dfer	1.10	IF (maskC(i,j,1,bi,bj).NE.0. _d 0) THEN
55	stephd	1.1	c calculate virtual flux
56			c EminusPforV = dS/dt*(1/Sglob)
57			C NOTE: Be very careful with signs here!
58			C Positive EminusPforV => loss of water to atmos and increase
59			C in salinity. Thus, also increase in other surface tracers
60			C (i.e. positive virtual flux into surface layer)
61	stephd	1.5	GALK(i,j)=gsm_ALKsurfaceForcingS(i,j,bi,bj)
62			& recip_drF(1)/gsm_S
63	stephd	1.1	c OR
64			c let virtual flux be zero
65			c GALK(i,j)=0.d0
66			c
67			ELSE
68			GALK(i,j)=0. _d 0
69			ENDIF
70			ENDDO
71			ENDDO
72	dfer	1.8	#endif /* ALLOW_OLD_VIRTUALFLUX */
73	stephd	1.1
74			#endif
75			#endif
76			RETURN
77			END