/[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.1 by stephd, Wed Jun 25 21:00:36 2003 UTC revision 1.4 by jmc, Tue Oct 28 22:19:53 2003 UTC
# Line 1  Line 1 
1  #include "CPP_OPTIONS.h"  #include "DIC_OPTIONS.h"
2  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
3    
4  CStartOfInterFace  CStartOfInterFace
# Line 46  c Line 46  c
46          call insol(myTime,sfac,bj)          call insol(myTime,sfac,bj)
47          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
48           DO i=1-OLx,sNx+OLx           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)             lit=sfac(j)
56             DO k=1,nlev             DO k=1,nlev
57               atten=(k0*delz(k)*.5)  C$TAF STORE lit = bio_export
58               if (k.gt.1) atten=(k0*delz(k-1)*.5)               atten=(k0*drF(k)*.5)
59                 if (k.gt.1) atten=(k0*drF(k-1)*.5)
60               lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj))               lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj))
61               if (lit.lt.0.d0.or.lit.gt.150) then               if (lit.lt.0.d0.or.lit.gt.150) then
62                   print*,'QQ lit', lit                   print*,'QQ lit', lit
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22