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 10251 for branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90 – NEMO

Ignore:
Timestamp:
2018-10-29T15:20:26+01:00 (5 years ago)
Author:
kingr
Message:

Rolled back to r10247 - i.e., undid merge of pkg br and 3.6_stable br

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90

    r10249 r10251  
    1919 
    2020   PUBLIC trc_wri_my_trc  
    21 #if defined key_tracer_budget 
    22    REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:), SAVE :: trb_temp ! slwa 
    23 #endif 
    24  
    2521 
    2622#  include "top_substitute.h90" 
    2723CONTAINS 
    2824 
    29 #if defined key_tracer_budget 
    30    SUBROUTINE trc_wri_my_trc (kt, fl) ! slwa 
    31 #else 
    3225   SUBROUTINE trc_wri_my_trc 
    33 #endif 
    3426      !!--------------------------------------------------------------------- 
    3527      !!                     ***  ROUTINE trc_wri_trc  *** 
     
    3729      !! ** Purpose :   output passive tracers fields  
    3830      !!--------------------------------------------------------------------- 
    39 #if defined key_tracer_budget 
    40       INTEGER, INTENT( in ), OPTIONAL     :: fl  
    41       INTEGER, INTENT( in )               :: kt 
    42       REAL(wp), DIMENSION(jpi,jpj,jpk)    :: trpool !tracer pool temporary output 
    43 #else 
    44       INTEGER, INTENT( in )               :: kt 
    45 #endif 
    4631      CHARACTER (len=20)   :: cltra 
    47       INTEGER              :: jn,jk ! JC TODO jk defined here but may not be used 
     32      INTEGER              :: jn 
    4833      !!--------------------------------------------------------------------- 
    4934  
    5035      ! write the tracer concentrations in the file 
    5136      ! --------------------------------------- 
    52  
    53  
    54 #if defined key_tracer_budget 
    55       IF( PRESENT(fl)) THEN 
    56 ! depth integrated 
    57 ! for strict budgetting write this out at end of timestep as an average between 'now' and 'after' at kt 
    58          DO jn = jp_myt0, jp_myt1  
    59           IF(ln_trdtrc (jn))THEN 
    60             trpool(:,:,:) = 0.5 * ( trn(:,:,:,jn) * fse3t_a(:,:,:) +  & 
    61                                         trb_temp(:,:,:,jn) * fse3t(:,:,:) ) 
    62             cltra = TRIM( ctrcnm(jn) )//"e3t"     ! depth integrated output 
    63             IF( kt == nittrc000 ) write(6,*)'output pool ',cltra 
    64             DO jk = 1, jpk 
    65                trpool(:,:,jk) = trpool(:,:,jk) 
    66             END DO 
    67             CALL iom_put( cltra, trpool) 
    68  
    69           END IF 
    70          END DO 
    71  
    72       ELSE 
    73  
    74          IF( kt == nittrc000 ) THEN 
    75            ALLOCATE(trb_temp(jpi,jpj,jpk,jp_my_trc))  ! slwa 
    76          ENDIF 
    77          trb_temp(:,:,:,:)=trn(:,:,:,:) ! slwa save for tracer budget (unfiltered trn) 
    78  
    79  
    80       END IF 
    81 #else 
    8237      DO jn = jp_myt0, jp_myt1 
    8338         cltra = TRIM( ctrcnm(jn) )                  ! short title for tracer 
    8439         CALL iom_put( cltra, trn(:,:,:,jn) ) 
    8540      END DO 
    86 #endif 
    8741      ! 
    8842   END SUBROUTINE trc_wri_my_trc 
     
    9448   PUBLIC trc_wri_my_trc 
    9549CONTAINS 
    96 #if defined key_tracer_budget 
    97    SUBROUTINE trc_wri_my_trc (kt, fl) ! slwa 
    98       INTEGER, INTENT( in ), OPTIONAL     :: fl  
    99       INTEGER, INTENT( in )               :: kt 
    100 #else 
    101    ! JC TODO Subroutine arguments (kt) inconsistent with earlier definition 
    102    SUBROUTINE trc_wri_my_trc (kt) 
    103       INTEGER, INTENT( in )               :: kt 
    104 #endif 
     50   SUBROUTINE trc_wri_my_trc                     ! Empty routine   
    10551   END SUBROUTINE trc_wri_my_trc 
    10652#endif 
Note: See TracChangeset for help on using the changeset viewer.