pkg/dic/bio_export.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/bio_export.F,v 1.18 2007/11/30 17:47:41 stephd Exp $
C $Name:  $

#include "DIC_OPTIONS.h"
#include "GCHEM_OPTIONS.h"

CBOP
C !ROUTINE: BIO_EXPORT

C !INTERFACE: ==========================================================
      SUBROUTINE BIO_EXPORT( PTR_PO4 ,
#ifdef ALLOW_FE
     I           PTR_FE,
#endif
     I           bioac,
     I           bi,bj,imin,imax,jmin,jmax,
     I           myIter,myTime,myThid)

c !DESCRIPTION:
C  Calculate biological activity and export

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DIC_ABIOTIC.h"
#include "DIC_BIOTIC.h"
#ifdef USE_QSW
#include "FFIELDS.h"
#endif

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
C  myIter               :: current timestep
C  myTime               :: current time
C  PTR_PO4              :: phosphate tracer field
C  PTR_FE               :: iron tracer field
      INTEGER myIter
      _RL myTime
      INTEGER myThid
      _RL  PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
#ifdef ALLOW_FE
      _RL  PTR_FE(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
#endif
      INTEGER imin, imax, jmin, jmax, bi, bj


C !OUTPUT PARAMETERS: ==================================================
C  bioac               :: biological productivity (will be split
C                         between export and dissolved pool)
      _RL bioac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR)

#ifdef ALLOW_PTRACERS
#ifdef DIC_BIOTIC

C !LOCAL VARIABLES: ====================================================
C  i,j,k                  :: loop indices
       INTEGER I,J,k
      _RL sfac(1-OLy:sNy+OLy)
      _RL lit, atten
      _RL nutlimit
#ifdef AD_SAFE
      _RL thx, thy, theps
#endif
CEOP


#ifndef READ_PAR
#ifndef USE_QSW
        call insol(myTime,sfac,bj)
#endif
#endif
        DO j=jmin,jmax
         DO i=imin,imax
C Fortran-90
CRG C$TAF INIT bio_export = static, nlev
C FORTRAN-77 dynamic memory uses adstore adresto
CRG C$TAF INIT bio_export = memory
C FORTRAN-77 with know max of nlev
C$TAF INIT bio_export = static, 10
#ifdef READ_PAR
           lit=PAR(i,j,bi,bj)
#else
#ifdef USE_QSW
           lit=-parfrac*Qsw(i,j,bi,bj)*maskC(i,j,1,bi,bj)
#else
           lit=sfac(j)
#endif
#endif

           if (QSW_underice) then
c if using Qsw and seaice/thsice or coupled, then ice 
c fraction is already taken into account
              lit=lit
           else
              lit=lit*(1. _d 0 - FIce(i,j,bi,bj))
           endif
c
           DO k=1,nlev
C$TAF STORE lit = bio_export
             atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0)
             if (k.gt.1) atten=atten+(k0*drF(k-1)
     &                         *hFacC(i,j,k-1,bi,bj)*.5 _d 0)
             lit=lit*exp(-atten)
             IF (lit.LT.0. _d 0.OR.lit.GT.350. _d 0) THEN
                 print*,'QQ lit', lit
             ENDIF

#ifdef ALLOW_FE
#ifdef AD_SAFE
             thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
             thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)
             theps = 1. _d -6
             nutlimit= ( 1. _d 0 - tanh((thx-thy)/theps) ) * thx/2. _d 0
     &                +( 1. _d 0 + tanh((thx-thy)/theps) ) * thy/2. _d 0
#else
             nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),
     &                        PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE) )
#endif
#else
             nutlimit=    PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
#endif

             bioac(i,j,k)=alpha(i,j,bi,bj)*
     &                    lit/(lit+lit0)*maskC(i,j,k,bi,bj)*
     &                    nutlimit
            ENDDO
          ENDDO
        ENDDO
c
#endif
#endif
       RETURN
       END
1	stephd	1.19	C $Header: /u/gcmpack/MITgcm/pkg/dic/bio_export.F,v 1.18 2007/11/30 17:47:41 stephd Exp $
2	jmc	1.15	C $Name: $
3
4	edhill	1.3	#include "DIC_OPTIONS.h"
5	stephd	1.1	#include "GCHEM_OPTIONS.h"
6
7	stephd	1.5	CBOP
8			C !ROUTINE: BIO_EXPORT
9
10			C !INTERFACE: ==========================================================
11	jmc	1.15	SUBROUTINE BIO_EXPORT( PTR_PO4 ,
12	stephd	1.1	#ifdef ALLOW_FE
13	jmc	1.15	I PTR_FE,
14			#endif
15	stephd	1.1	I bioac,
16			I bi,bj,imin,imax,jmin,jmax,
17			I myIter,myTime,myThid)
18
19	stephd	1.5	c !DESCRIPTION:
20	jmc	1.15	C Calculate biological activity and export
21	stephd	1.5
22			C !USES: ===============================================================
23	stephd	1.1	IMPLICIT NONE
24			#include "SIZE.h"
25			#include "DYNVARS.h"
26			#include "EEPARAMS.h"
27			#include "PARAMS.h"
28			#include "GRID.h"
29			#include "DIC_ABIOTIC.h"
30			#include "DIC_BIOTIC.h"
31	stephd	1.12	#ifdef USE_QSW
32	stephd	1.11	#include "FFIELDS.h"
33			#endif
34	stephd	1.1
35	stephd	1.5	C !INPUT PARAMETERS: ===================================================
36			C myThid :: thread number
37			C myIter :: current timestep
38			C myTime :: current time
39			C PTR_PO4 :: phosphate tracer field
40			C PTR_FE :: iron tracer field
41	stephd	1.1	INTEGER myIter
42			_RL myTime
43			INTEGER myThid
44			_RL PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
45			#ifdef ALLOW_FE
46			_RL PTR_FE(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
47			#endif
48	stephd	1.5	INTEGER imin, imax, jmin, jmax, bi, bj
49
50
51			C !OUTPUT PARAMETERS: ==================================================
52			C bioac :: biological productivity (will be split
53			C between export and dissolved pool)
54	stephd	1.1	_RL bioac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR)
55
56			#ifdef ALLOW_PTRACERS
57			#ifdef DIC_BIOTIC
58	stephd	1.5
59			C !LOCAL VARIABLES: ====================================================
60			C i,j,k :: loop indices
61	stephd	1.1	INTEGER I,J,k
62	stephd	1.5	_RL sfac(1-OLy:sNy+OLy)
63			_RL lit, atten
64	stephd	1.10	_RL nutlimit
65			#ifdef AD_SAFE
66			_RL thx, thy, theps
67			#endif
68	stephd	1.5	CEOP
69	stephd	1.10
70
71	stephd	1.8	#ifndef READ_PAR
72	stephd	1.12	#ifndef USE_QSW
73	stephd	1.1	call insol(myTime,sfac,bj)
74	stephd	1.8	#endif
75	stephd	1.11	#endif
76	stephd	1.9	DO j=jmin,jmax
77			DO i=imin,imax
78	stephd	1.2	C Fortran-90
79			CRG C$TAF INIT bio_export = static, nlev
80			C FORTRAN-77 dynamic memory uses adstore adresto
81			CRG C$TAF INIT bio_export = memory
82			C FORTRAN-77 with know max of nlev
83			C$TAF INIT bio_export = static, 10
84	stephd	1.8	#ifdef READ_PAR
85			lit=PAR(i,j,bi,bj)
86			#else
87	stephd	1.12	#ifdef USE_QSW
88	stephd	1.13	lit=-parfracQsw(i,j,bi,bj)maskC(i,j,1,bi,bj)
89	stephd	1.11	#else
90	stephd	1.1	lit=sfac(j)
91	stephd	1.8	#endif
92	stephd	1.11	#endif
93	stephd	1.18
94	stephd	1.19	if (QSW_underice) then
95			c if using Qsw and seaice/thsice or coupled, then ice
96			c fraction is already taken into account
97			lit=lit
98			else
99			lit=lit*(1. _d 0 - FIce(i,j,bi,bj))
100			endif
101			c
102	stephd	1.1	DO k=1,nlev
103	stephd	1.2	C$TAF STORE lit = bio_export
104	dfer	1.17	atten=(k0drF(k)hFacC(i,j,k,bi,bj)*.5 _d 0)
105	stephd	1.7	if (k.gt.1) atten=atten+(k0*drF(k-1)
106	dfer	1.17	& hFacC(i,j,k-1,bi,bj).5 _d 0)
107	dfer	1.14	lit=lit*exp(-atten)
108	dfer	1.17	IF (lit.LT.0. _d 0.OR.lit.GT.350. _d 0) THEN
109	stephd	1.1	print*,'QQ lit', lit
110	dfer	1.17	ENDIF
111	stephd	1.10
112	stephd	1.1	#ifdef ALLOW_FE
113	stephd	1.10	#ifdef AD_SAFE
114			thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
115			thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)
116	dfer	1.17	theps = 1. _d -6
117			nutlimit= ( 1. _d 0 - tanh((thx-thy)/theps) ) * thx/2. _d 0
118			& +( 1. _d 0 + tanh((thx-thy)/theps) ) * thy/2. _d 0
119	stephd	1.5	#else
120	stephd	1.10	nutlimit= min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),
121			& PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE) )
122	stephd	1.1	#endif
123	stephd	1.10	#else
124			nutlimit= PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
125			#endif
126
127			bioac(i,j,k)=alpha(i,j,bi,bj)*
128			& lit/(lit+lit0)maskC(i,j,k,bi,bj)
129			& nutlimit
130	stephd	1.1	ENDDO
131			ENDDO
132			ENDDO
133			c
134			#endif
135			#endif
136			RETURN
137			END