- Timestamp:
- 2019-12-16T14:36:58+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/TOP/PISCES/P4Z/p4zsed.F90
r12210 r12257 330 330 ENDIF 331 331 332 IF( lk_iomput ) THEN 333 IF( knt == nrdttrc ) THEN 334 zfact = 1.e+3 * rfact2r ! conversion from molC/l/kt to molN/m3/s 335 IF( iom_use("Nfix" ) ) CALL iom_put( "Nfix", nitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) ) ! nitrogen fixation 336 IF( iom_use("INTNFIX") ) THEN ! nitrogen fixation rate in ocean ( vertically integrated ) 337 zwork(:,:) = 0. 338 DO jk = 1, jpkm1 339 zwork(:,:) = zwork(:,:) + nitrpot(:,:,jk) * nitrfix * rno3 * zfact * e3t_n(:,:,jk) * tmask(:,:,jk) 340 ENDDO 341 CALL iom_put( "INTNFIX" , zwork ) 342 ENDIF 343 IF( iom_use("SedCal" ) ) CALL iom_put( "SedCal", zsedcal(:,:) * zfact ) 344 IF( iom_use("SedSi" ) ) CALL iom_put( "SedSi", zsedsi (:,:) * zfact ) 345 IF( iom_use("SedC" ) ) CALL iom_put( "SedC", zsedc (:,:) * zfact ) 346 IF( iom_use("Sdenit" ) ) CALL iom_put( "Sdenit", sdenit (:,:) * zfact * rno3 ) 347 ENDIF 332 IF( lk_iomput .AND. knt == nrdttrc ) THEN 333 zfact = 1.e+3 * rfact2r ! conversion from molC/l/kt to molN/m3/s 334 CALL iom_put( "Nfix", nitrpot(:,:,:) * nitrfix * rno3 * zfact * tmask(:,:,:) ) ! nitrogen fixation 335 CALL iom_put( "SedCal", zsedcal(:,:) * zfact ) 336 CALL iom_put( "SedSi" , zsedsi (:,:) * zfact ) 337 CALL iom_put( "SedC" , zsedc (:,:) * zfact ) 338 CALL iom_put( "Sdenit", sdenit (:,:) * zfact * rno3 ) 348 339 ENDIF 349 340 !
Note: See TracChangeset
for help on using the changeset viewer.