New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4996 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90 – NEMO

Ignore:
Timestamp:
2014-12-22T12:56:38+01:00 (9 years ago)
Author:
cetlod
Message:

trunk : New diagnostics in PISCES model, see ticket #1424

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

    r4624 r4996  
    7878      REAL(wp) ::   zofer2 
    7979#endif 
    80       REAL(wp) ::   zonitr, zstep, zrfact2 
     80      REAL(wp) ::   zonitr, zstep, zfact 
    8181      CHARACTER (len=25) :: charout 
    82       REAL(wp), POINTER, DIMENSION(:,:  ) :: ztempbac  
    83       REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zolimi, zdepprod 
     82      REAL(wp), POINTER, DIMENSION(:,:  ) :: ztempbac 
     83      REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zolimi, zdepprod, zw3d 
    8484      !!--------------------------------------------------------------------- 
    8585      ! 
     
    315315      END DO 
    316316 
    317       IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN 
    318           zrfact2 = 1.e3 * rfact2r 
    319           CALL iom_put( "REMIN" , zolimi(:,:,:) * tmask(:,:,:) * zrfact2 )  ! Remineralisation rate 
    320           CALL iom_put( "DENIT" , denitr(:,:,:) * rdenit * rno3 * tmask(:,:,:) * zrfact2  )  ! Denitrification 
    321       ENDIF 
     317      IF( jnt == nrdttrc ) THEN 
     318          CALL wrk_alloc( jpi, jpj, jpk, zw3d ) 
     319          zfact = 1.e+3 * rfact2r  !  conversion from mol/l/kt to  mol/m3/s 
     320          ! 
     321          IF( iom_use( "REMIN" ) )  THEN 
     322              zw3d(:,:,:) = zolimi(:,:,:) * tmask(:,:,:) * zfact !  Remineralisation rate 
     323              CALL iom_put( "REMIN"  , zw3d ) 
     324          ENDIF 
     325          IF( iom_use( "DENIT" ) )  THEN 
     326              zw3d(:,:,:) = denitr(:,:,:) * rdenit * rno3 * tmask(:,:,:) * zfact ! Denitrification 
     327              CALL iom_put( "DENIT"  , zw3d ) 
     328          ENDIF 
     329          ! 
     330          CALL wrk_dealloc( jpi, jpj, jpk, zw3d ) 
     331       ENDIF 
    322332 
    323333      IF(ln_ctl)   THEN  ! print mean trends (used for debugging) 
Note: See TracChangeset for help on using the changeset viewer.