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 32 for trunk/NEMO/OPA_SRC/DIA/diagap.F90 – NEMO

Ignore:
Timestamp:
2004-02-17T10:20:15+01:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE001 : First major NEMO update

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DIA/diagap.F90

    r3 r32  
    2121 
    2222   IMPLICIT NONE 
     23   PRIVATE 
     24 
     25   !! * Routine accessibility 
     26   PUBLIC dia_gap     ! called in step.F90 module 
    2327 
    2428   !! * Shared module variables 
    25    LOGICAL, PUBLIC, PARAMETER ::   & 
    26       lk_diagap = .TRUE.         ! model-data diagnostics flag 
     29   LOGICAL, PUBLIC, PARAMETER ::   lk_diagap = .TRUE.   !: model-data diagnostics flag 
    2730 
    2831   !! * Module variables 
    2932   INTEGER ::                 & 
    30 !!!   numgap,                 &  ! logical unit for differences diagnostic 
     33!???  numgap,                 &  ! logical unit for differences diagnostic 
    3134      ngap  ,                 &  ! time step frequency 
    3235      nprg                       ! switch for control print 
     
    125128            END DO 
    126129         END DO 
    127 #if defined key_mpp 
    128          CALL mpp_sum( volk, jpk ) 
    129 #endif 
     130         IF( lk_mpp )   CALL mpp_sum( volk, jpk )   ! sum over the global domain 
     131 
    130132         volkr(:) = 0.e0 
    131133         DO jk = 1, jpk 
     
    233235            smodg(jpk) = smodg(jpk) + smodg(jk) * volk(jk) / vol 
    234236         END DO   
    235  
    236 #if defined key_mpp 
    237           CALL mpp_sum( tdtag, jpk ) 
    238           CALL mpp_sum( sdtag, jpk ) 
    239           CALL mpp_sum( tmodg, jpk ) 
    240           CALL mpp_sum( smodg, jpk ) 
    241 #endif 
     237          IF( lk_mpp)   CALL mpp_sum( tdtag, jpk )   ! sum over the global domain 
     238          IF( lk_mpp)   CALL mpp_sum( sdtag, jpk ) 
     239          IF( lk_mpp)   CALL mpp_sum( tmodg, jpk ) 
     240          IF( lk_mpp)   CALL mpp_sum( smodg, jpk ) 
    242241 
    243242          ! 3.  Averaged output in file numgap 
     
    291290#else 
    292291   !!---------------------------------------------------------------------- 
    293    !!   Default option :                                       Empty module 
    294    !!---------------------------------------------------------------------- 
    295    LOGICAL, PUBLIC, PARAMETER ::   lk_diagap = .FALSE.   ! 'key_diagap' flag 
     292   !!   Default option :                                       Dummy module 
     293   !!---------------------------------------------------------------------- 
     294   LOGICAL, PUBLIC, PARAMETER ::   lk_diagap = .FALSE.    !: diagap flag 
    296295CONTAINS 
    297    SUBROUTINE dia_gap( kt )           ! Empty routine 
    298       WRITE(*,*) kt                       ! no warning in compilation phase 
     296   SUBROUTINE dia_gap( kt )           ! Dummy routine 
     297      WRITE(*,*) 'dia_gap: You should not have seen this print! error?', kt 
    299298   END SUBROUTINE dia_gap 
    300299#endif 
Note: See TracChangeset for help on using the changeset viewer.