- Timestamp:
- 2017-02-01T11:22:41+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zrem.F90
r7617 r7627 80 80 CHARACTER (len=25) :: charout 81 81 REAL(wp), POINTER, DIMENSION(:,: ) :: ztempbac 82 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zwork1, zdepprod, zfacsi, zfacsib, znitr 82 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zdepprod, zdepeff, zfacsi, zfacsib 83 REAL(wp), POINTER, DIMENSION(:,:,:) :: znitr, zfebact, zwork1 83 84 !!--------------------------------------------------------------------- 84 85 ! … … 87 88 ! Allocate temporary workspace 88 89 CALL wrk_alloc( jpi, jpj, ztempbac ) 89 CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zwork1, zfacsi, zfacsib, znitr ) 90 CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib) 91 CALL wrk_alloc( jpi, jpj, jpk, znitr, zfebact, zwork1 ) 90 92 91 93 ! Initialisation of temprary arrys 92 94 zdepprod(:,:,:) = 1._wp 95 zdepeff (:,:,:) = 0.3_wp 93 96 ztempbac(:,:) = 0._wp 97 zfebact (:,:,:) = 0._wp 94 98 zfacsib(:,:,:) = xsilab / ( 1.0 - xsilab ) 95 99 zfacsi(:,:,:) = xsilab … … 111 115 zdepbac (ji,jj,jk) = zdepmin**0.683 * ztempbac(ji,jj) 112 116 zdepprod(ji,jj,jk) = zdepmin**0.273 117 zdepeff (ji,jj,jk) = 0.3 * zdepmin**0.3 113 118 ENDIF 114 119 END DO … … 220 225 zbactfer = feratb * rfact2 * prmaxp(ji,jj,jk) * xlimbacl(ji,jj,jk) & 221 226 & * biron(ji,jj,jk) / ( xkferb + biron(ji,jj,jk) ) & 222 & * zdepprod(ji,jj,jk) * zdep bac(ji,jj,jk)227 & * zdepprod(ji,jj,jk) * zdepeff(ji,jj,jk) * zdepbac(ji,jj,jk) 223 228 #if defined key_kriest 224 229 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.05 225 230 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.05 231 zfebact(ji,jj,jk) = zbactfer * 0.15 226 232 #else 227 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.12 228 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.09 229 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.03 233 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.39 234 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.3 235 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.09 236 zfebact(ji,jj,jk) = zbactfer * 0.39 230 237 #endif 231 238 END DO … … 284 291 CALL iom_put( "DENIT" , denitrc(:,:,:) * rdenit * rno3 * tmask(:,:,:) * zrfact2 ) ! Denitrification 285 292 CALL iom_put( "NIT" , znitr(:,:,:) * rno3 * tmask(:,:,:) * zrfact2 ) ! 286 CALL iom_put( "BACT" , zdepbac(:,:,:) * 1.E6 * tmask(:,:,:) ) ! Bacterial biomass 293 CALL iom_put( "BACT", zdepbac(:,:,:) * 1.E6 * tmask(:,:,:) ) ! Bacterial biomass 294 CALL iom_put( "FEBACT" , zfebact(:,:,:) * 1E9 * tmask(:,:,:) * zrfact2 ) ! Denitrification 287 295 ENDIF 288 296 ! 289 297 CALL wrk_dealloc( jpi, jpj, ztempbac ) 290 CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zwork1, zfacsi, zfacsib, znitr ) 298 CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib ) 299 CALL wrk_dealloc( jpi, jpj, jpk, zwork1, zfebact, znitr ) 291 300 ! 292 301 IF( nn_timing == 1 ) CALL timing_stop('p5z_rem')
Note: See TracChangeset
for help on using the changeset viewer.