Ignore:
Timestamp:
2014-12-22T12:56:38+01:00 (6 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/p4zmicro.F90

    r4800 r4996  
    7979      REAL(wp) :: zgrazp, zgrazm, zgrazsd 
    8080      REAL(wp) :: zgrazmf, zgrazsf, zgrazpf 
    81       REAL(wp) :: zrfact2 
    82       REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing 
     81      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d 
    8382      CHARACTER (len=25) :: charout 
    8483      !!--------------------------------------------------------------------- 
     
    8685      IF( nn_timing == 1 )  CALL timing_start('p4z_micro') 
    8786      ! 
    88       IF( ln_diatrc .AND. lk_iomput )  CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
     87      IF( lk_iomput )  CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
    8988      ! 
    9089      DO jk = 1, jpkm1 
     
    191190      END DO 
    192191      ! 
    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 ) 
    196199      ENDIF 
    197200      ! 
     
    202205      ENDIF 
    203206      ! 
    204       IF( ln_diatrc .AND. lk_iomput )  CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
     207      IF( lk_iomput )  CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
    205208      ! 
    206209      IF( nn_timing == 1 )  CALL timing_stop('p4z_micro') 
Note: See TracChangeset for help on using the changeset viewer.