Changeset 3325 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90

Timestamp:

2012-03-12T15:44:43+01:00 (12 years ago)

Author:

gm

Message:

Ediag branche: #927 add Kinetic Energy trend diagnostics (trdken.F90)

File:

• branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90 (modified) (2 diffs)

branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90

@@ -11 +11 @@
    !!----------------------------------------------------------------------
    USE trd_oce        ! trends: ocean variables
-!   USE ldftra_oce     ! ocean active tracers lateral physics
-   USE trdglo         ! ocean bassin integral constraints properties
+   USE trdken         ! trends: 3D kinetic energy
+   USE trdglo         ! trends: global domain integral constraints properties
    USE trdmld         ! ocean active mixed layer tracers trends 
    USE trdvor         ! ocean vorticity trends 
@@ -79 +79 @@
       !
       
-      IF( ln_KE_trd .OR. ln_PE_trd .OR. ln_dyn_mld )   &
-         CALL ctl_stop( 'KE, PE, aur ML on momentum are not yet coded we stop' )
+      IF( ln_PE_trd .OR. ln_dyn_mld )   &
+         CALL ctl_stop( 'PE, or ML on momentum are not yet coded we stop' )
 
       !
-      IF( ln_glo_trd )    CALL trd_glo_init       ! integral constraints trends
-      IF( ln_tra_mld )    CALL trd_mld_init       ! mixed-layer trends (active  tracers)  
-      IF( ln_vor_trd )    CALL trd_vor_init       ! vorticity trends        
+      IF( ln_glo_trd )   CALL trd_glo_init      ! integral constraints trends
+      IF( ln_tra_mld )   CALL trd_mld_init      ! mixed-layer trends (active  tracers)  
+      IF( ln_vor_trd )   CALL trd_vor_init      ! vorticity trends
+      IF( ln_KE_trd  )   CALL trd_ken_init      ! 3D Kinetic energy diagnostics
+
       !
 !!gm  : Potential BUG : 3D output only for vector invariant form!  add a ctl_stop or code the flux form case