/[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.27 by mlosch, Fri Jul 30 15:02:43 2010 UTC revision 1.28 by mlosch, Mon Aug 2 07:10:52 2010 UTC
# Line 75  CEOP Line 75  CEOP
75        CALL INSOL(myTime,sfac,bi,bj,myThid)        CALL INSOL(myTime,sfac,bi,bj,myThid)
76  #endif  #endif
77  #endif  #endif
       DO j=jmin,jmax  
        DO i=imin,imax  
78  cph: following init assumes nlev=nr  cph: following init assumes nlev=nr
79  cph: set in dic_init_fixed.F (ALLOW_AUTODIFF)  cph: set in dic_init_fixed.F (ALLOW_AUTODIFF)
80  C$TAF INIT bio_export = static, nr  C$TAF INIT bio_export = static, nsx*nsy*nr
81  C FORTRAN-77 dynamic memory uses adstore adresto  C FORTRAN-77 dynamic memory uses adstore adresto
82  CRG C$TAF INIT bio_export = memory  CRG C$TAF INIT bio_export = memory
83          DO j=jmin,jmax
84           DO i=imin,imax
85  #ifdef READ_PAR  #ifdef READ_PAR
86          lit(i,j)=PAR(i,j,bi,bj)          lit(i,j)=PAR(i,j,bi,bj)
87  #elif (defined USE_QSW)  #elif (defined USE_QSW)
# Line 99  c ice fraction needs to be taken into ac Line 99  c ice fraction needs to be taken into ac
99        ENDDO        ENDDO
100  c  c
101        DO k=1,nlev        DO k=1,nlev
102    C$TAF STORE lit = bio_export
103         DO j=jmin,jmax         DO j=jmin,jmax
104          DO i=imin,imax          DO i=imin,imax
 C$TAF STORE lit(i,j) = bio_export  
105           atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0)           atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0)
106           if (k.gt.1) atten=atten+(k0*drF(k-1)           if (k.gt.1) atten=atten+(k0*drF(k-1)
107       &                         *hFacC(i,j,k-1,bi,bj)*.5 _d 0)       &                         *hFacC(i,j,k-1,bi,bj)*.5 _d 0)
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22