- Timestamp:
- 2019-12-14T14:57:23+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p4zmicro.F90
r12236 r12252 70 70 REAL(wp) :: zgrazp, zgrazm, zgrazsd 71 71 REAL(wp) :: zgrazmf, zgrazsf, zgrazpf 72 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo 73 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zw3d, zzligprod 72 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo, zzligprod 74 73 CHARACTER (len=25) :: charout 75 74 !!--------------------------------------------------------------------- 76 75 ! 77 76 IF( ln_timing ) CALL timing_start('p4z_micro') 78 !79 IF (ln_ligand) THEN80 ALLOCATE( zzligprod(jpi,jpj,jpk) )81 zzligprod(:,:,:) = 0._wp82 ENDIF83 77 ! 84 78 DO jk = 1, jpkm1 … … 187 181 END DO 188 182 ! 189 IF( lk_iomput ) THEN 190 IF( knt == nrdttrc ) THEN 191 ALLOCATE( zw3d(jpi,jpj,jpk) ) 192 IF( iom_use( "GRAZ1" ) ) THEN 193 zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Total grazing of phyto by zooplankton 194 CALL iom_put( "GRAZ1", zw3d ) 195 ENDIF 196 IF( iom_use( "FEZOO" ) ) THEN 197 zw3d(:,:,:) = zfezoo(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ! 198 CALL iom_put( "FEZOO", zw3d ) 199 ENDIF 200 IF( iom_use( "LPRODZ" ) .AND. ln_ligand ) THEN 201 zw3d(:,:,:) = zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) 202 CALL iom_put( "LPRODZ" , zw3d ) 203 ENDIF 204 DEALLOCATE( zw3d ) 183 IF( lk_iomput .AND. knt == nrdttrc ) THEN 184 zgrazing(:,:,jpk) = 0._wp 185 zfezoo (:,:,jpk) = 0._wp 186 CALL iom_put( "GRAZ1", zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 187 CALL iom_put( "FEZOO", zfezoo (:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 188 IF( ln_ligand ) THEN 189 zzligprod(:,:,jpk) = 0._wp 190 CALL iom_put( "LPRODZ", zzligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) ) 205 191 ENDIF 206 192 ENDIF 207 !208 IF (ln_ligand) DEALLOCATE( zzligprod )209 193 ! 210 194 IF(sn_cfctl%l_prttrc) THEN ! print mean trends (used for debugging)
Note: See TracChangeset
for help on using the changeset viewer.