/[MITgcm]/MITgcm/pkg/bling/bling_remineralization.F
ViewVC logotype

Diff of /MITgcm/pkg/bling/bling_remineralization.F

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

revision 1.7 by mmazloff, Thu Oct 27 17:55:31 2016 UTC revision 1.8 by mmazloff, Fri Jan 27 23:11:31 2017 UTC
# Line 163  C          we are not yet at the bottom Line 163  C          we are not yet at the bottom
163  C  Sinking speed is evaluated at the bottom of the cell  C  Sinking speed is evaluated at the bottom of the cell
164            depth_l=-rF(k+1)            depth_l=-rF(k+1)
165            IF (depth_l .LE. wsink0z)  THEN            IF (depth_l .LE. wsink0z)  THEN
166             wsink = wsink0             wsink = wsink0_2d(i,j,bi,bj)
167            ELSE            ELSE
168             wsink = wsinkacc * (depth_l - wsink0z) + wsink0             wsink = wsinkacc * (depth_l - wsink0z) + wsink0_2d(i,j,bi,bj)
169            ENDIF            ENDIF
170    
171  C  Nutrient remineralization lengthscale  C  Nutrient remineralization lengthscale
172  C  Not an e-folding scale: this term increases with remineralization.  C  Not an e-folding scale: this term increases with remineralization.
173            zremin = gamma_POM * ( PTR_O2(i,j,k)**2 /            zremin = gamma_POM2d(i,j,bi,bj) * ( PTR_O2(i,j,k)**2 /
174       &               (k_O2**2 + PTR_O2(i,j,k)**2) * (1-remin_min)       &               (k_O2**2 + PTR_O2(i,j,k)**2) * (1-remin_min)
175       &               + remin_min )/(wsink + epsln)       &               + remin_min )/(wsink + epsln)
176    
# Line 351  c     &           *CtoP/NUTfac*12.01/wsi Line 351  c     &           *CtoP/NUTfac*12.01/wsi
351       &             * MasstoN)**(0.58)*FreeFe       &             * MasstoN)**(0.58)*FreeFe
352  #else  #else
353              Fe_ads_org(i,j,k) =              Fe_ads_org(i,j,k) =
354       &           kFE_org*(PONflux_l/(epsln + wsink0)       &           kFE_org*(PONflux_l/(epsln + wsink0_2d(i,j,bi,bj))
355       &             * MasstoN)**(0.58)*FreeFe       &             * MasstoN)**(0.58)*FreeFe
356  #endif  #endif
357             ENDIF             ENDIF
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22