Changeset 2104 for branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD
- Timestamp:
- 2010-09-17T14:35:46+02:00 (14 years ago)
- Location:
- branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD/trdicp.F90
r2082 r2104 4 4 !! Ocean diagnostics: ocean tracers and dynamic trends 5 5 !!===================================================================== 6 !! History : ! 91-12 (G. Madec) 7 !! ! 92-06 (M. Imbard) add time step frequency 8 !! ! 96-01 (G. Madec) terrain following coordinates 9 !! 8.5 ! 02-06 (G. Madec) F90: Free form and module 10 !! 9.0 ! 04-08 (C. Talandier) New trends organization 6 !! History : 1.0 ! 2004-08 (C. Talandier) New trends organization 11 7 !!---------------------------------------------------------------------- 12 8 #if defined key_trdtra || defined key_trddyn || defined key_esopa … … 14 10 !! 'key_trdtra' or active tracers trends diagnostics 15 11 !! 'key_trddyn' momentum trends diagnostics 16 !!----------------------------------------------------------------------17 12 !!---------------------------------------------------------------------- 18 13 !! trd_icp : compute the basin averaged properties for tra/dyn … … 48 43 # include "vectopt_loop_substitute.h90" 49 44 !!---------------------------------------------------------------------- 50 !! OPA 9.0 , LOCEAN-IPSL (2005)51 !! $Id$ 45 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 46 !! $Id$ 52 47 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 53 48 !!---------------------------------------------------------------------- -
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD/trdmod.F90
r2026 r2104 4 4 !! Ocean diagnostics: ocean tracers and dynamic trends 5 5 !!===================================================================== 6 !! History : 9.0 ! 04-08 (C. Talandier) Original code 7 !! ! 05-04 (C. Deltel) Add Asselin trend in the ML budget 6 !! History : 1.0 ! 2004-08 (C. Talandier) Original code 7 !! - ! 2005-04 (C. Deltel) Add Asselin trend in the ML budget 8 !! 3.3 ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase 8 9 !!---------------------------------------------------------------------- 9 10 #if defined key_trdtra || defined key_trddyn || defined key_trdmld || defined key_trdvor || defined key_esopa … … 36 37 # include "vectopt_loop_substitute.h90" 37 38 !!---------------------------------------------------------------------- 38 !! OPA 9.0 , LOCEAN-IPSL (2005)39 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 39 40 !! $Id$ 40 41 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 216 217 217 218 ENDIF 218 219 ! 219 220 END SUBROUTINE trd_mod 220 221 221 # 222 #else 222 223 !!---------------------------------------------------------------------- 223 224 !! Default case : Empty module … … 230 231 CONTAINS 231 232 SUBROUTINE trd_mod(ptrd3dx, ptrd3dy, ktrd , ctype, kt ) ! Empty routine 232 REAL, DIMENSION(:,:,:), INTENT( in ) :: & 233 ptrd3dx, & ! Temperature or U trend 234 ptrd3dy ! Salinity or V trend 235 INTEGER, INTENT( in ) :: ktrd ! momentum or tracer trend index 236 INTEGER, INTENT( in ) :: kt ! Time step 237 CHARACTER(len=3), INTENT( in ) :: ctype ! momentum or tracers trends type 238 WRITE(*,*) 'trd_3d: You should not have seen this print! error ?', ptrd3dx(1,1,1) 239 WRITE(*,*) ' " ": You should not have seen this print! error ?', ptrd3dy(1,1,1) 240 WRITE(*,*) ' " ": You should not have seen this print! error ?', ktrd 241 WRITE(*,*) ' " ": You should not have seen this print! error ?', ctype 242 WRITE(*,*) ' " ": You should not have seen this print! error ?', kt 233 REAL :: ptrd3dx(:,:,:), ptrd3dy(:,:,:) 234 INTEGER :: ktrd, kt 235 CHARACTER(len=3) :: ctype 236 WRITE(*,*) 'trd_3d: You should not have seen this print! error ?', ptrd3dx(1,1,1), ptrd3dy(1,1,1) 237 WRITE(*,*) ' " ": You should not have seen this print! error ?', ktrd, ctype, kt 243 238 END SUBROUTINE trd_mod 244 # 239 #endif 245 240 246 241 SUBROUTINE trd_mod_init … … 251 246 !!---------------------------------------------------------------------- 252 247 USE in_out_manager ! I/O manager 253 248 !! 254 249 NAMELIST/namtrd/ nn_trd, nn_ctls, cn_trdrst_in, cn_trdrst_out, ln_trdmld_restart, rn_ucf, ln_trdmld_instant 255 250 !!---------------------------------------------------------------------- -
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD/trdmod_oce.F90
r2026 r2104 4 4 !! Ocean trends : set tracer and momentum trend variables 5 5 !!====================================================================== 6 !! History : 9.0 !04-08 (C. Talandier) Original code6 !! History : 1.0 ! 2004-08 (C. Talandier) Original code 7 7 !!---------------------------------------------------------------------- 8 8 USE trdicp_oce ! ocean momentum/tracers bassin properties trends variables … … 38 38 LOGICAL , PUBLIC :: l_trdtrc = .FALSE. !: tracers trend flag 39 39 # endif 40 ! !!!Active tracers trends indexes41 INTEGER, PUBLIC, PARAMETER :: jptra_trd_xad = 1 !: x- horizontal advection42 INTEGER, PUBLIC, PARAMETER :: jptra_trd_yad = 2 !: y- horizontal advection43 INTEGER, PUBLIC, PARAMETER :: jptra_trd_zad = 3 !: z- vertical advection44 INTEGER, PUBLIC, PARAMETER :: jptra_trd_ldf = 4 !: lateral diffusion45 INTEGER, PUBLIC, PARAMETER :: jptra_trd_zdf = 5 !: vertical diffusion (Kz)46 INTEGER, PUBLIC, PARAMETER :: jptra_trd_bbc = 6 !: Bottom Boundary Condition (geoth. flux)47 INTEGER, PUBLIC, PARAMETER :: jptra_trd_bbl = 7 !: Bottom Boundary Layer (diffusive/convective)48 INTEGER, PUBLIC, PARAMETER :: jptra_trd_npc = 8 !: static instability mixing49 INTEGER, PUBLIC, PARAMETER :: jptra_trd_dmp = 9 !: damping50 INTEGER, PUBLIC, PARAMETER :: jptra_trd_qsr = 10 !: penetrative solar radiation51 INTEGER, PUBLIC, PARAMETER :: jptra_trd_nsr = 11 !: non solar radiation52 INTEGER, PUBLIC, PARAMETER :: jptra_trd_atf = 12 !: Asselin correction40 ! !!!* Active tracers trends indexes 41 INTEGER, PUBLIC, PARAMETER :: jptra_trd_xad = 1 !: x- horizontal advection 42 INTEGER, PUBLIC, PARAMETER :: jptra_trd_yad = 2 !: y- horizontal advection 43 INTEGER, PUBLIC, PARAMETER :: jptra_trd_zad = 3 !: z- vertical advection 44 INTEGER, PUBLIC, PARAMETER :: jptra_trd_ldf = 4 !: lateral diffusion 45 INTEGER, PUBLIC, PARAMETER :: jptra_trd_zdf = 5 !: vertical diffusion (Kz) 46 INTEGER, PUBLIC, PARAMETER :: jptra_trd_bbc = 6 !: Bottom Boundary Condition (geoth. flux) 47 INTEGER, PUBLIC, PARAMETER :: jptra_trd_bbl = 7 !: Bottom Boundary Layer (diffusive/convective) 48 INTEGER, PUBLIC, PARAMETER :: jptra_trd_npc = 8 !: static instability mixing 49 INTEGER, PUBLIC, PARAMETER :: jptra_trd_dmp = 9 !: damping 50 INTEGER, PUBLIC, PARAMETER :: jptra_trd_qsr = 10 !: penetrative solar radiation 51 INTEGER, PUBLIC, PARAMETER :: jptra_trd_nsr = 11 !: non solar radiation 52 INTEGER, PUBLIC, PARAMETER :: jptra_trd_atf = 12 !: Asselin correction 53 53 #if defined key_top 54 ! !* Passive tracers trends indexes55 INTEGER, PUBLIC, PARAMETER :: jptra_trd_sms = 13 !: sources m. sinks56 INTEGER, PUBLIC, PARAMETER :: jptra_trd_radn = 14 !: corr. trn<0 in trcrad57 INTEGER, PUBLIC, PARAMETER :: jptra_trd_radb = 15 !: corr. trb<0 in trcrad (like atf)54 ! !!!* Passive tracers trends indexes 55 INTEGER, PUBLIC, PARAMETER :: jptra_trd_sms = 13 !: sources m. sinks 56 INTEGER, PUBLIC, PARAMETER :: jptra_trd_radn = 14 !: corr. trn<0 in trcrad 57 INTEGER, PUBLIC, PARAMETER :: jptra_trd_radb = 15 !: corr. trb<0 in trcrad (like atf) 58 58 #endif 59 59 60 ! !!!Momentum trends indexes61 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_hpg = 1 !: hydrostatic pressure gradient62 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_keg = 2 !: kinetic energy gradient63 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_rvo = 3 !: relative vorticity64 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_pvo = 4 !: planetary vorticity65 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_ldf = 5 !: lateral diffusion66 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_had = 6 !: horizontal advection67 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_zad = 7 !: vertical advection68 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_zdf = 8 !: vertical diffusion69 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_spg = 9 !: surface pressure gradient70 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_dat = 10 !: damping term71 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_swf = 11 !: surface wind forcing72 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_bfr = 12 !: bottom friction60 ! !!!* Momentum trends indexes 61 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_hpg = 1 !: hydrostatic pressure gradient 62 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_keg = 2 !: kinetic energy gradient 63 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_rvo = 3 !: relative vorticity 64 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_pvo = 4 !: planetary vorticity 65 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_ldf = 5 !: lateral diffusion 66 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_had = 6 !: horizontal advection 67 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_zad = 7 !: vertical advection 68 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_zdf = 8 !: vertical diffusion 69 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_spg = 9 !: surface pressure gradient 70 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_dat = 10 !: damping term 71 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_swf = 11 !: surface wind forcing 72 INTEGER, PUBLIC, PARAMETER :: jpdyn_trd_bfr = 12 !: bottom friction 73 73 74 74 !!---------------------------------------------------------------------- 75 !! NEMO/OPA 3. 2 , LOCEAN-IPSL (2009)76 !! $Id$ 75 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 76 !! $Id$ 77 77 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 78 78 !!====================================================================== -
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRD/trdmod_trc.F90
r2026 r2104 5 5 !!====================================================================== 6 6 !!---------------------------------------------------------------------- 7 !! Dummy module NO TOP use7 !! Dummy module NO TOP use 8 8 !!---------------------------------------------------------------------- 9 9 CONTAINS 10 10 11 11 SUBROUTINE trd_mod_trc( ptrtrd, kjn, ktrd, kt ) 12 INTEGER , INTENT( in ) :: kt ! time step 13 INTEGER , INTENT( in ) :: kjn ! tracer index 14 INTEGER , INTENT( in ) :: ktrd ! tracer trend index 15 REAL, DIMENSION(:,:,:), INTENT( in ) :: ptrtrd ! Temperature or U trend 12 INTEGER :: kt, kjn, ktrd 13 REAL :: ptrtrd(:,:,:) 16 14 WRITE(*,*) 'trd_mod_trc_trp : You should not have seen this print! error?', ptrtrd(1,1,1) 17 WRITE(*,*) ' " " : You should not have seen this print! error?', kjn 18 WRITE(*,*) ' " " : You should not have seen this print! error?', ktrd 19 WRITE(*,*) ' " " : You should not have seen this print! error?', kt 15 WRITE(*,*) ' " " : You should not have seen this print! error?', kjn, ktrd, kt 20 16 END SUBROUTINE trd_mod_trc 21 17 18 !!====================================================================== 22 19 END MODULE trdmod_trc
Note: See TracChangeset
for help on using the changeset viewer.