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