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 3329 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90 – NEMO

Ignore:
Timestamp:
2012-03-16T13:22:15+01:00 (12 years ago)
Author:
gm
Message:

Ediag branche: #927 add the missing trdpen initialisation in trdini

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90

    r3325 r3329  
    1010   !!   trd_init      : initialization step 
    1111   !!---------------------------------------------------------------------- 
     12   USE dom_oce        ! ocean domain 
    1213   USE trd_oce        ! trends: ocean variables 
    13    USE trdken         ! trends: 3D kinetic energy 
    14    USE trdglo         ! trends: global domain integral constraints properties 
    15    USE trdmld         ! ocean active mixed layer tracers trends  
    16    USE trdvor         ! ocean vorticity trends  
     14   USE trdken         ! trends: 3D kinetic   energy 
     15   USE trdpen         ! trends: 3D potential energy 
     16   USE trdglo         ! trends: global domain averaged tracers and dynamics 
     17   USE trdmld         ! trends: mixed layer averaged trends (tracer only) 
     18   USE trdvor         ! trends: vertical averaged vorticity  
    1719   USE in_out_manager ! I/O manager 
    1820   USE lib_mpp        ! MPP library 
     
    7981      ! 
    8082       
    81       IF( ln_PE_trd .OR. ln_dyn_mld )   & 
    82          CALL ctl_stop( 'PE, or ML on momentum are not yet coded we stop' ) 
    83  
    84       ! 
    85       IF( ln_glo_trd )   CALL trd_glo_init      ! integral constraints trends 
    86       IF( ln_tra_mld )   CALL trd_mld_init      ! mixed-layer trends (active  tracers)   
    87       IF( ln_vor_trd )   CALL trd_vor_init      ! vorticity trends 
    88       IF( ln_KE_trd  )   CALL trd_ken_init      ! 3D Kinetic energy diagnostics 
    89  
    90       ! 
     83      IF( lk_vvl )   CALL ctl_stop( 'trend diagnostics with variable volume not validated' ) 
     84       
    9185!!gm  : Potential BUG : 3D output only for vector invariant form!  add a ctl_stop or code the flux form case 
    9286!!gm  : bug/pb for vertical advection of tracer in vvl case: add T.dt[eta] in the output...  
    93       ! 
     87 
     88      !                             ! diagnostic initialization   
     89      IF( ln_glo_trd )   CALL trd_glo_init      ! global domain averaged trends 
     90      IF( ln_tra_mld )   CALL trd_mld_init      ! mixed-layer          trends   
     91      IF( ln_vor_trd )   CALL trd_vor_init      ! barotropic vorticity trends 
     92      IF( ln_KE_trd  )   CALL trd_ken_init      ! 3D Kinetic    energy trends 
     93      IF( ln_PE_trd  )   CALL trd_pen_init      ! 3D Potential  energy trends 
    9494      ! 
    9595   END SUBROUTINE trd_init 
Note: See TracChangeset for help on using the changeset viewer.