Changeset 3446 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
- Timestamp:
- 2012-08-10T13:13:55+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
r3443 r3446 81 81 REAL(wp) :: zgrazn, zgrazpoc, zgraznf, zgrazf 82 82 REAL(wp) :: zgrazfff, zgrazffe 83 REAL(wp) :: zrfact2 84 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing 83 85 CHARACTER (len=25) :: charout 84 REAL(wp) :: zrfact285 86 !!--------------------------------------------------------------------- 86 87 ! 87 88 IF( nn_timing == 1 ) CALL timing_start('p4z_meso') 88 89 ! 89 90 IF( ln_diatrc .AND. lk_iomput ) CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 91 ! 90 92 DO jk = 1, jpkm1 91 93 DO jj = 1, jpj … … 145 147 146 148 ! Total grazing ( grazing by microzoo is already computed in p4zmicro ) 147 grazing(ji,jj,jk) = grazing(ji,jj,jk) + zgraztot 149 IF( ln_diatrc .AND. lk_iomput ) zgrazing(ji,jj,jk) = zgraztot 150 148 151 ! Mesozooplankton efficiency 149 152 ! -------------------------- … … 213 216 END DO 214 217 ! 215 IF( ln_diatrc .AND. lk_iomput ) THEN218 IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN 216 219 zrfact2 = 1.e3 * rfact2r 217 grazing(:,:,:) = grazing(:,:,:) * zrfact2 * tmask(:,:,:) ! Total grazing of phyto by zoo 218 prodcal(:,:,:) = prodcal(:,:,:) * zrfact2 * tmask(:,:,:) ! Calcite production 219 IF( jnt == nrdttrc ) THEN 220 CALL iom_put( "GRAZ" , grazing ) ! Total grazing of phyto by zooplankton 221 CALL iom_put( "PCAL" , prodcal ) ! Calcite production 222 ENDIF 220 CALL iom_put( "GRAZ2", zgrazing(:,:,:) * zrfact2 * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 221 CALL iom_put( "PCAL" , prodcal(:,:,:) * zrfact2 * tmask(:,:,:) ) ! Calcite production 223 222 ENDIF 224 223 ! … … 228 227 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm) 229 228 ENDIF 229 ! 230 IF( ln_diatrc .AND. lk_iomput ) CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 230 231 ! 231 232 IF( nn_timing == 1 ) CALL timing_stop('p4z_meso')
Note: See TracChangeset
for help on using the changeset viewer.