- Timestamp:
- 2017-05-09T12:14:45+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90
r7180 r8003 87 87 zcompaph = MAX( ( trb(ji,jj,jk,jpphy) - 1e-8 ), 0.e0 ) 88 88 zstep = xstep 89 # if defined key_degrad90 zstep = zstep * facvol(ji,jj,jk)91 # endif92 89 ! When highly limited by macronutrients, very small cells 93 90 ! dominate the community. As a consequence, aggregation … … 121 118 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca 122 119 tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca 123 #if defined key_kriest124 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortp125 tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + ztortp * xkr_dnano + zrespp * xkr_ddiat126 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zmortp * zfactfe127 #else128 120 tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zfracal * zmortp 129 121 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + ( 1. - zfracal ) * zmortp … … 132 124 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ( 1. - zfracal ) * zmortp * zfactfe 133 125 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zfracal * zmortp * zfactfe 134 #endif135 126 END DO 136 127 END DO … … 181 172 ! ------------------------------------------------------------ 182 173 zstep = xstep 183 # if defined key_degrad184 zstep = zstep * facvol(ji,jj,jk)185 # endif186 174 ! Phytoplankton respiration 187 175 ! ------------------------ … … 206 194 tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zmortp2 * zfactsi 207 195 tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zmortp2 * zfactsi 208 #if defined key_kriest209 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zmortp2210 tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + ztortp2 * xkr_ddiat + zrespp2 * xkr_daggr211 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zmortp2 * zfactfe212 #else213 196 tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zrespp2 + 0.5 * ztortp2 214 197 tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + 0.5 * ztortp2 … … 217 200 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + 0.5 * ztortp2 * zfactfe 218 201 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ( zrespp2 + 0.5 * ztortp2 ) * zfactfe 219 #endif220 202 END DO 221 203 END DO
Note: See TracChangeset
for help on using the changeset viewer.