Changeset 5075 for branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
- Timestamp:
- 2015-02-11T11:50:34+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
r4624 r5075 83 83 REAL(wp) :: zgrazfffp, zgrazfffg, zgrazffep, zgrazffeg 84 84 CHARACTER (len=25) :: charout 85 REAL(wp) :: zrfact2 86 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing 85 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d 87 86 88 87 !!--------------------------------------------------------------------- … … 90 89 IF( nn_timing == 1 ) CALL timing_start('p4z_meso') 91 90 ! 92 IF( l n_diatrc .AND. lk_iomput ) THEN91 IF( lk_iomput ) THEN 93 92 CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 94 93 zgrazing(:,:,:) = 0._wp … … 163 162 zfrac = zproport * grazflux * zstep * wsbio4(ji,jj,jk) & 164 163 & * trn(ji,jj,jk,jpgoc) * trn(ji,jj,jk,jpmes) & 165 & * ( 0. 1 + 3.9* zratio2 / ( 1.**2 + zratio2 ) )164 & * ( 0.2 + 3.8 * zratio2 / ( 1.**2 + zratio2 ) ) 166 165 zfracfe = zfrac * trn(ji,jj,jk,jpbfe) / (trn(ji,jj,jk,jpgoc) + rtrn) 167 166 … … 186 185 187 186 ! Total grazing ( grazing by microzoo is already computed in p4zmicro ) 188 IF( l n_diatrc .AND. lk_iomput ) zgrazing(ji,jj,jk) = zgraztot187 IF( lk_iomput ) zgrazing(ji,jj,jk) = zgraztot 189 188 190 189 ! Mesozooplankton efficiency … … 249 248 END DO 250 249 ! 251 IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN 252 zrfact2 = 1.e3 * rfact2r 253 CALL iom_put( "GRAZ2", zgrazing(:,:,:) * zrfact2 * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 254 CALL iom_put( "PCAL" , prodcal(:,:,:) * zrfact2 * tmask(:,:,:) ) ! Calcite production 250 IF( lk_iomput .AND. jnt == nrdttrc ) THEN 251 CALL wrk_alloc( jpi, jpj, jpk, zw3d ) 252 IF( iom_use( "GRAZ2" ) ) THEN 253 zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Total grazing of phyto by zooplankton 254 CALL iom_put( "GRAZ2", zw3d ) 255 ENDIF 256 IF( iom_use( "PCAL" ) ) THEN 257 zw3d(:,:,:) = prodcal(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Calcite production 258 CALL iom_put( "PCAL", zw3d ) 259 ENDIF 260 CALL wrk_dealloc( jpi, jpj, jpk, zw3d ) 255 261 ENDIF 256 262 ! … … 261 267 ENDIF 262 268 ! 263 IF( l n_diatrc .AND. lk_iomput ) CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )269 IF( lk_iomput ) CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 264 270 ! 265 271 IF( nn_timing == 1 ) CALL timing_stop('p4z_meso')
Note: See TracChangeset
for help on using the changeset viewer.