pkg/dic/alk_surfforcing.F

#include "DIC_OPTIONS.h"
#include "GCHEM_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_ABIOTIC.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

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

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