Changeset 5038 for branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90
- Timestamp:
- 2015-01-20T15:26:13+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90
r4529 r5038 79 79 REAL(wp) :: zgrazp, zgrazm, zgrazsd 80 80 REAL(wp) :: zgrazmf, zgrazsf, zgrazpf 81 REAL(wp) :: zrfact2 82 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing 81 REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d 83 82 CHARACTER (len=25) :: charout 84 83 !!--------------------------------------------------------------------- … … 86 85 IF( nn_timing == 1 ) CALL timing_start('p4z_micro') 87 86 ! 88 IF( l n_diatrc .AND. lk_iomput ) CALL wrk_alloc( jpi, jpj, jpk, zgrazing )87 IF( lk_iomput ) CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 89 88 ! 90 89 DO jk = 1, jpkm1 … … 137 136 ! Various remineralization and excretion terms 138 137 ! -------------------------------------------- 139 zgrasrat = zgraztotf/ ( zgraztot + rtrn )140 zgrasratn = zgraztotn/ ( zgraztot + rtrn )138 zgrasrat = ( zgraztotf + rtrn ) / ( zgraztot + rtrn ) 139 zgrasratn = ( zgraztotn + rtrn ) / ( zgraztot + rtrn ) 141 140 zepshert = MIN( 1., zgrasratn, zgrasrat / ferat3) 142 141 zepsherv = zepshert * MIN( epsher, (1. - unass) * zgrasrat / ferat3, (1. - unass) * zgrasratn ) … … 191 190 END DO 192 191 ! 193 IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN 194 zrfact2 = 1.e3 * rfact2r 195 CALL iom_put( "GRAZ1" , zgrazing(:,:,:) * zrfact2 * tmask(:,:,:) ) ! Total grazing of phyto by zooplankton 192 IF( lk_iomput .AND. jnt == nrdttrc ) THEN 193 CALL wrk_alloc( jpi, jpj, jpk, zw3d ) 194 IF( iom_use( "GRAZ1" ) ) THEN 195 zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ! Total grazing of phyto by zooplankton 196 CALL iom_put( "GRAZ1", zw3d ) 197 ENDIF 198 CALL wrk_dealloc( jpi, jpj, jpk, zw3d ) 196 199 ENDIF 197 200 ! … … 202 205 ENDIF 203 206 ! 204 IF( l n_diatrc .AND. lk_iomput ) CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )207 IF( lk_iomput ) CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 205 208 ! 206 209 IF( nn_timing == 1 ) CALL timing_stop('p4z_micro') … … 235 238 READ ( numnatp_cfg, nampiszoo, IOSTAT = ios, ERR = 902 ) 236 239 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiszoo in configuration namelist', lwp ) 237 WRITE ( numonp, nampiszoo )240 IF(lwm) WRITE ( numonp, nampiszoo ) 238 241 239 242 IF(lwp) THEN ! control print
Note: See TracChangeset
for help on using the changeset viewer.