/[MITgcm]/MITgcm/pkg/dic/bio_export.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/bio_export.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.14 by dfer, Sun Aug 26 23:05:48 2007 UTC revision 1.15 by jmc, Tue Oct 9 00:01:42 2007 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
5  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
6    
# Line 5  CBOP Line 8  CBOP
8  C !ROUTINE: BIO_EXPORT  C !ROUTINE: BIO_EXPORT
9    
10  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
11        SUBROUTINE BIO_EXPORT( PTR_PO4 ,        SUBROUTINE BIO_EXPORT( PTR_PO4 ,
12  #ifdef ALLOW_FE  #ifdef ALLOW_FE
13       I           PTR_FE,       I           PTR_FE,
14  #endif    #endif
15       I           bioac,       I           bioac,
16       I           bi,bj,imin,imax,jmin,jmax,       I           bi,bj,imin,imax,jmin,jmax,
17       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
18    
19  c !DESCRIPTION:  c !DESCRIPTION:
20  C  Calculate biological activity and export                  C  Calculate biological activity and export
21    
22  C !USES: ===============================================================  C !USES: ===============================================================
23        IMPLICIT NONE        IMPLICIT NONE
# Line 103  C$TAF STORE lit = bio_export Line 106  C$TAF STORE lit = bio_export
106               thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)               thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
107               thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)               thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)
108               theps = 1.d-6               theps = 1.d-6
109               nutlimit= ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +               nutlimit= ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +
110       &                 ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2       &                 ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2
111  #else  #else
112               nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),               nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22