/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F

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

revision 1.16 by jahn, Fri Dec 27 17:29:00 2013 UTC revision 1.17 by jahn, Fri Sep 12 16:19:39 2014 UTC
# Line 1220  c ---- end steph's alternative Line 1220  c ---- end steph's alternative
1220                  dfetl=dfetl+fet_flx(i,j,k,bi,bj)                  dfetl=dfetl+fet_flx(i,j,k,bi,bj)
1221                  dsil=dsil+si_flx(i,j,k,bi,bj)                  dsil=dsil+si_flx(i,j,k,bi,bj)
1222  #endif  #endif
1223  c  
1224    #ifdef ALLOW_OBCS
1225              IF (useOBCS) THEN
1226                dpo4l  = dpo4l *maskInC(i,j,bi,bj)
1227                dno3l  = dno3l *maskInC(i,j,bi,bj)
1228                dfetl  = dfetl *maskInC(i,j,bi,bj)
1229                dsil   = dsil  *maskInC(i,j,bi,bj)
1230                ddopl  = ddopl *maskInC(i,j,bi,bj)
1231                ddonl  = ddonl *maskInC(i,j,bi,bj)
1232                ddofel = ddofel*maskInC(i,j,bi,bj)
1233                dpopl  = dpopl *maskInC(i,j,bi,bj)
1234                dponl  = dponl *maskInC(i,j,bi,bj)
1235                dpofel = dpofel*maskInC(i,j,bi,bj)
1236                dpsil  = dpsil *maskInC(i,j,bi,bj)
1237                dnh4l  = dnh4l *maskInC(i,j,bi,bj)
1238                dno2l  = dno2l *maskInC(i,j,bi,bj)
1239                DO nz = 1,nzmax
1240                 dzoop (nz) = dzoop (nz)*maskInC(i,j,bi,bj)
1241                 dzoon (nz) = dzoon (nz)*maskInC(i,j,bi,bj)
1242                 dzoofe(nz) = dzoofe(nz)*maskInC(i,j,bi,bj)
1243                 dzoosi(nz) = dzoosi(nz)*maskInC(i,j,bi,bj)
1244                ENDDO
1245                DO np = 1,npmax
1246                 dPhy(np) = dPhy(np)*maskInC(i,j,bi,bj)
1247    #ifdef GEIDER
1248    #ifdef DYNAMIC_CHL
1249                 dphychl(np) = dphychl(np)*maskInC(i,j,bi,bj)
1250    #endif
1251    #endif
1252                ENDDO
1253    #ifdef ALLOW_CDOM
1254                dcdoml = dcdoml*maskInC(i,j,bi,bj)
1255    #endif
1256    #ifdef ALLOW_CARBON
1257                ddicl = ddicl*maskInC(i,j,bi,bj)
1258                ddocl = ddocl*maskInC(i,j,bi,bj)
1259                dpocl = dpocl*maskInC(i,j,bi,bj)
1260                dpicl = dpicl*maskInC(i,j,bi,bj)
1261                dalkl = dalkl*maskInC(i,j,bi,bj)
1262                do2l  = do2l *maskInC(i,j,bi,bj)
1263                DO nz = 1,nzmax
1264                 dzoocl(nz) = dzoocl(nz)*maskInC(i,j,bi,bj)
1265                ENDDO
1266    #endif
1267              ENDIF
1268    #endif
1269    
1270  c now update main tracer arrays  c now update main tracer arrays
1271            dtplankton = PTRACERS_dTLev(k)/float(nsubtime)            dtplankton = PTRACERS_dTLev(k)/float(nsubtime)
1272            Ptr(i,j,k,bi,bj,iPO4 ) = Ptr(i,j,k,bi,bj,iPO4) +            Ptr(i,j,k,bi,bj,iPO4 ) = Ptr(i,j,k,bi,bj,iPO4) +
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22