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 8850 for branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90 – NEMO

Ignore:
Timestamp:
2017-11-30T09:30:44+01:00 (6 years ago)
Author:
gm
Message:

#1911 (ENHANCE-09): PART I.3 - phasing with trunk: tracer trends output (see #1877 trunk change from 86666 to 8698)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

    r7881 r8850  
    2828   USE oce_trc         ! ocean dynamics and tracers variables 
    2929   USE trc             ! ocean passive tracers variables 
    30    USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    31    USE prtctl_trc      ! Print control for debbuging 
    3230   USE trd_oce 
    3331   USE trdtra 
     
    3836   USE agrif_top_interp 
    3937# endif 
     38   ! 
     39   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     40   USE prtctl_trc      ! Print control for debbuging 
    4041 
    4142   IMPLICIT NONE 
    4243   PRIVATE 
    4344 
    44    PUBLIC   trc_nxt          ! routine called by step.F90 
     45   PUBLIC   trc_nxt   ! routine called by step.F90 
    4546 
    4647   REAL(wp) ::   rfact1, rfact2 
     
    8283      REAL(wp) ::   zfact            ! temporary scalar 
    8384      CHARACTER (len=22) :: charout 
    84       REAL(wp), POINTER, DIMENSION(:,:,:,:) ::  ztrdt  
     85      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) ::   ztrdt    ! 4D workspace 
    8586      !!---------------------------------------------------------------------- 
    8687      ! 
     
    102103 
    103104      IF( l_trdtrc )  THEN             ! trends: store now fields before the Asselin filter application 
    104          CALL wrk_alloc( jpi, jpj, jpk, jptra, ztrdt ) 
     105         ALLOCATE( ztrdt(jpi,jpj,jpk,jptra) ) 
    105106         ztrdt(:,:,:,:)  = trn(:,:,:,:) 
    106107      ENDIF 
     
    137138            END DO 
    138139         END DO 
    139          CALL wrk_dealloc( jpi, jpj, jpk, jptra, ztrdt )  
     140         DEALLOCATE( ztrdt )  
    140141      END IF 
    141142      ! 
Note: See TracChangeset for help on using the changeset viewer.