Changeset 3329 for branches/2012/dev_r3309_LOCEAN12_Ediag
- Timestamp:
- 2012-03-16T13:22:15+01:00 (12 years ago)
- 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 10 10 !! trd_init : initialization step 11 11 !!---------------------------------------------------------------------- 12 USE dom_oce ! ocean domain 12 13 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 17 19 USE in_out_manager ! I/O manager 18 20 USE lib_mpp ! MPP library … … 79 81 ! 80 82 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 91 85 !!gm : Potential BUG : 3D output only for vector invariant form! add a ctl_stop or code the flux form case 92 86 !!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 94 94 ! 95 95 END SUBROUTINE trd_init -
branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90
r3326 r3329 217 217 ENDIF 218 218 ! 219 nkstp = 0219 nkstp = nit000 - 1 220 220 r1_2_rau0 = 0.5_wp / rau0 221 221 ! -
branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90
r3327 r3329 142 142 IF ( lk_vvl ) CALL ctl_stop('trd_pen_init : PE trends not coded for variable volume') 143 143 ! 144 nkstp = nit000 - 1 145 ! 144 146 END SUBROUTINE trd_pen_init 145 147
Note: See TracChangeset
for help on using the changeset viewer.