- Timestamp:
- 2017-02-01T11:22:41+01:00 (7 years ago)
- Location:
- branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90
r7617 r7627 64 64 ! 65 65 IF( nn_timing == 1 ) CALL timing_start('p4z_ligand') 66 ! 67 CALL wrk_alloc( jpi, jpj, jpk, zligrem, zligpr, zrligprod ) 66 68 ! 67 69 ! ------------------------------------------------------------------ -
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90
r7617 r7627 413 413 CALL iom_put( "REMINP" , zremipoc(:,:,:) * tmask(:,:,:) ) ! Remineralisation rate 414 414 CALL iom_put( "REMING" , zremigoc(:,:,:) * tmask(:,:,:) ) ! Remineralisation rate 415 CALL iom_put( "REMINF" , zfolimi(:,:,:) * tmask(:,:,:) * 1.e+ 3 * rfact2r) ! Remineralisation rate415 CALL iom_put( "REMINF" , zfolimi(:,:,:) * tmask(:,:,:) * 1.e+9 * zrfact2 ) ! Remineralisation rate 416 416 ENDIF 417 417 -
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90
r7617 r7627 79 79 CHARACTER (len=25) :: charout 80 80 REAL(wp), POINTER, DIMENSION(:,: ) :: ztempbac 81 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zolimi, zdepprod, zfacsi, zfacsib, znitr 81 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zdepprod, zdepeff, zfacsi, zfacsib 82 REAL(wp), POINTER, DIMENSION(:,:,:) :: znitr, zolimi, zfecbact 82 83 !!--------------------------------------------------------------------- 83 84 ! … … 86 87 ! Allocate temporary workspace 87 88 CALL wrk_alloc( jpi, jpj, ztempbac ) 88 CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, znitr, zfacsi, zfacsib ) 89 CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib ) 90 CALL wrk_alloc( jpi, jpj, jpk, zolimi, znitr, zfebact ) 89 91 90 92 ! Initialization of local variables … … 93 95 ! Initialisation of temprary arrys 94 96 zdepprod(:,:,:) = 1._wp 97 zdepeff (:,:,:) = 0.3_wp 95 98 ztempbac(:,:) = 0._wp 99 zfebact (:,:,:) = 0._wp 96 100 zfacsib(:,:,:) = xsilab / ( 1.0 - xsilab ) 97 101 zfacsi(:,:,:) = xsilab … … 113 117 zdepbac (ji,jj,jk) = zdepmin**0.683 * ztempbac(ji,jj) 114 118 zdepprod(ji,jj,jk) = zdepmin**0.273 119 zdepeff (ji,jj,jk) = 0.3 * zdepmin**0.3 115 120 ENDIF 116 121 END DO … … 208 213 zbactfer = feratb * rfact2 * prmax(ji,jj,jk) * xlimbacl(ji,jj,jk) & 209 214 & * trb(ji,jj,jk,jpfer) / ( xkferb + trb(ji,jj,jk,jpfer) ) & 210 & * zdepprod(ji,jj,jk) * zdep bac(ji,jj,jk)215 & * zdepprod(ji,jj,jk) * zdepeff(ji,jj,jk) * zdepbac(ji,jj,jk) 211 216 #if defined key_kriest 212 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.05 213 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.05 217 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.15 218 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.15 219 zfebact(ji,jj,jk) = zbactfer * 0.15 214 220 #else 215 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.16 216 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.12 217 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.04 221 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.39 222 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.3 223 tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.09 224 zfebact(ji,jj,jk) = zbactfer * 0.39 218 225 #endif 219 226 END DO … … 274 281 CALL iom_put( "NIT" , znitr(:,:,:) * rno3 * tmask(:,:,:) * zrfact2 ) ! 275 282 CALL iom_put( "BACT", zdepbac(:,:,:) * 1.E6 * tmask(:,:,:) ) ! Bacterial biomass 283 CALL iom_put( "FEBACT" , zfebact(:,:,:) * 1E9 * tmask(:,:,:) * zrfact2 ) ! Denitrification 276 284 277 285 ENDIF 278 286 ! 279 287 CALL wrk_dealloc( jpi, jpj, ztempbac ) 280 CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib, znitr ) 288 CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib ) 289 CALL wrk_dealloc( jpi, jpj, jpk, zolimi, znitr, zfebact ) 281 290 ! 282 291 IF( nn_timing == 1 ) CALL timing_stop('p4z_rem')
Note: See TracChangeset
for help on using the changeset viewer.