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

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

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


Revision 1.3 - (show annotations) (download)
Thu Oct 9 04:19:19 2003 UTC (20 years, 6 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 #include "DIC_OPTIONS.h"
2 #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 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 lit=sfac(j)
56 DO k=1,nlev
57 C$TAF STORE lit = bio_export
58 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