Changeset 3329


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

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

Location:
branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD
Files:
3 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 
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90

    r3326 r3329  
    217217      ENDIF 
    218218      ! 
    219       nkstp     = 0 
     219      nkstp     = nit000 - 1 
    220220      r1_2_rau0 = 0.5_wp / rau0 
    221221      ! 
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90

    r3327 r3329  
    142142      IF ( lk_vvl               )   CALL ctl_stop('trd_pen_init : PE trends not coded for variable volume') 
    143143      ! 
     144      nkstp     = nit000 - 1 
     145      ! 
    144146   END SUBROUTINE trd_pen_init 
    145147 
Note: See TracChangeset for help on using the changeset viewer.