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

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

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

Revision 1.3 - (hide annotations) (download)
Thu Oct 9 04:19:19 2003 UTC (20 years, 8 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51o_pre, checkpoint51n_pre, checkpoint51l_post, checkpoint51o_post, checkpoint51i_post, checkpoint51l_pre, checkpoint51n_post, checkpoint51m_post
Branch point for: tg2-branch, checkpoint51n_branch
Changes since 1.2: +1 -1 lines 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 edhill 1.3 #include "DIC_OPTIONS.h"
2 stephd 1.1 #include "GCHEM_OPTIONS.h"
3    
4     CStartOfInterFace
5     SUBROUTINE BIO_EXPORT( PTR_PO4 ,
6     #ifdef ALLOW_FE
7     I PTR_FE,
8     #endif
9     I bioac,
10     I bi,bj,imin,imax,jmin,jmax,
11     I myIter,myTime,myThid)
12    
13     C /==========================================================\
14     C | SUBROUTINE BIO_EXPORT |
15     C | o Calculate biological activity and export |
16     C |==========================================================|
17     IMPLICIT NONE
18    
19     C == GLobal variables ==
20     #include "SIZE.h"
21     #include "DYNVARS.h"
22     #include "EEPARAMS.h"
23     #include "PARAMS.h"
24     #include "GRID.h"
25     #include "DIC_ABIOTIC.h"
26     #include "DIC_BIOTIC.h"
27    
28     C == Routine arguments ==
29     INTEGER myIter
30     _RL myTime
31     INTEGER myThid
32     _RL PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
33     #ifdef ALLOW_FE
34     _RL PTR_FE(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
35     #endif
36     _RL bioac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR)
37     _RL sfac(1-OLy:sNy+OLy)
38     _RL lit, atten
39     INTEGER imin, imax, jmin, jmax, bi, bj
40    
41     #ifdef ALLOW_PTRACERS
42     #ifdef DIC_BIOTIC
43     C == Local variables ==
44     INTEGER I,J,k
45     c
46     call insol(myTime,sfac,bj)
47     DO j=1-OLy,sNy+OLy
48     DO i=1-OLx,sNx+OLx
49 stephd 1.2 C Fortran-90
50     CRG C$TAF INIT bio_export = static, nlev
51     C FORTRAN-77 dynamic memory uses adstore adresto
52     CRG C$TAF INIT bio_export = memory
53     C FORTRAN-77 with know max of nlev
54     C$TAF INIT bio_export = static, 10
55 stephd 1.1 lit=sfac(j)
56     DO k=1,nlev
57 stephd 1.2 C$TAF STORE lit = bio_export
58 stephd 1.1 atten=(k0*delz(k)*.5)
59     if (k.gt.1) atten=(k0*delz(k-1)*.5)
60     lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj))
61     if (lit.lt.0.d0.or.lit.gt.150) then
62     print*,'QQ lit', lit
63     endif
64     bioac(i,j,k)=alpha(i,j,bi,bj)*
65     & PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)*
66     #ifdef ALLOW_FE
67     & PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)*
68     #endif
69     & lit/(lit+lit0)*maskC(i,j,k,bi,bj)
70     ENDDO
71     ENDDO
72     ENDDO
73     c
74     #endif
75     #endif
76     RETURN
77     END
  ViewVC Help
Powered by ViewVC 1.1.22