Changeset 2334


Ignore:
Timestamp:
2010-10-29T08:15:27+02:00 (10 years ago)
Author:
gm
Message:

v3.3beta: Suppress of trabbc in diahsb

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

    r2287 r2334  
    22   !!====================================================================== 
    33   !!                       ***  MODULE  diahsb  *** 
    4    !! Ocean diagnostics: Heat salt and volume budgets 
     4   !! Ocean diagnostics: Heat, salt and volume budgets 
    55   !!====================================================================== 
    6    !! History :  NEMO 3.3  !  2010-09  (M. Leclair)  Original code  
    7    !!---------------------------------------------------------------------- 
    8    !! * Modules used 
     6   !! History :  3.3  ! 2010-09  (M. Leclair)  Original code  
     7   !!---------------------------------------------------------------------- 
     8 
     9   !!---------------------------------------------------------------------- 
    910   USE oce             ! ocean dynamics and tracers 
    1011   USE dom_oce         ! ocean space and time domain 
     
    2223   PRIVATE 
    2324 
    24    !! * Routine accessibility 
    25    PUBLIC dia_hsb        ! routine called by step.F90 
    26    PUBLIC dia_hsb_init   ! routine called by opa.F90 
    27  
    28    LOGICAL , PUBLIC ::   ln_diahsb  = .FALSE.   !: check the heat and salt budgets 
    29  
    30    !! * Module variables 
     25   PUBLIC   dia_hsb        ! routine called by step.F90 
     26   PUBLIC   dia_hsb_init   ! routine called by opa.F90 
     27 
     28   LOGICAL, PUBLIC ::   ln_diahsb  = .FALSE.   !: check the heat and salt budgets 
     29 
    3130   INTEGER                                 ::   numhsb                           ! 
    3231   REAL(dp)                                ::   surf_tot   , vol_tot             ! 
     
    4544   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    4645   !! $Id$ 
    47    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     46   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4847   !!---------------------------------------------------------------------- 
    4948 
     
    5453      !!                  ***  ROUTINE dia_hsb  *** 
    5554      !!      
    56       !! ** Purpose: Compute the ocean global heat content, salt content and volume 
    57       !!             non conservation 
     55      !! ** Purpose: Compute the ocean global heat content, salt content and volume conservation 
    5856      !!  
    5957      !! ** Method : - Compute the deviation of heat content, salt content and volume 
    60       !!        at the current time step from their values at nit000 
    61       !!      - Compute the contribution of forcing and remove it from these 
    62       !!                deviations 
     58      !!             at the current time step from their values at nit000 
     59      !!             - Compute the contribution of forcing and remove it from these deviations 
     60      !! 
    6361      !! ** Action : Write the results in the 'heat_salt_volume_budgets.txt' ASCII file 
    6462      !!--------------------------------------------------------------------------- 
     
    8280      z_frc_trd_s =           SUM( sbc_tsc(:,:,jp_sal) * surf(:,:) )     ! salt fluxes 
    8381      ! Add penetrative solar radiation 
    84       IF( ln_traqsr ) z_frc_trd_t = z_frc_trd_t + ro0cpr * SUM( qsr     (:,:) * surf(:,:) ) 
    85 #if defined key_trabbc 
     82      IF( ln_traqsr )   z_frc_trd_t = z_frc_trd_t + ro0cpr * SUM( qsr     (:,:) * surf(:,:) ) 
    8683      ! Add geothermal heat flux 
    87       IF( lk_trabbc ) z_frc_trd_t = z_frc_trd_t + ro0cpr * SUM( qgh_trd0(:,:) * surf(:,:) ) 
    88 #endif 
     84      IF( ln_trabbc )   z_frc_trd_t = z_frc_trd_t + ro0cpr * SUM( qgh_trd0(:,:) * surf(:,:) ) 
    8985      IF( lk_mpp ) THEN 
    9086         CALL mpp_sum( z_frc_trd_v ) 
     
    143139 
    1441409020  FORMAT(I5,11D15.7) 
    145  
     141      ! 
    146142   END SUBROUTINE dia_hsb 
     143 
    147144 
    148145   SUBROUTINE dia_hsb_init 
     
    161158      INTEGER            ::   jk       ! dummy loop indice 
    162159      INTEGER            ::   ierror   ! local integer 
    163  
     160      !! 
    164161      NAMELIST/namhsb/ ln_diahsb 
    165162      !!---------------------------------------------------------------------- 
    166  
     163      ! 
    167164      REWIND ( numnam )              ! Read Namelist namhsb  
    168165      READ   ( numnam, namhsb ) 
    169  
     166      ! 
    170167      IF(lwp) THEN                   ! Control print 
    171168         WRITE(numout,*) 
     
    175172         WRITE(numout,*) '      Switch for hsb diagnostic (T) or not (F)  ln_diahsb  = ', ln_diahsb 
    176173      ENDIF 
    177       !!--------------------------------------------------------------------------- 
    178174 
    179175      IF( .NOT. ln_diahsb )   RETURN 
     
    218214         vol_tot  = vol_tot + SUM( surf(:,:) * tmask(:,:,jk)     & 
    219215            &                      * fse3t_n(:,:,jk)         ) 
    220       ENDDO 
     216      END DO 
    221217      IF( lk_mpp ) THEN  
    222218         CALL mpp_sum( vol_tot ) 
Note: See TracChangeset for help on using the changeset viewer.