source: NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/TOP/TRP/trctrp.F90 @ 12110

Last change on this file since 12110 was 12110, checked in by cetlod, 8 months ago

merge dev_r11219_TOP-01_cethe_PISCES_LBC onto dev_r12072_TOP-01_ENHANCE-11_CEthe

  • Property svn:keywords set to Id
File size: 5.8 KB
Line 
1MODULE trctrp
2   !!======================================================================
3   !!                       ***  MODULE trctrp  ***
4   !! Ocean Physics    : manage the passive tracer transport
5   !!======================================================================
6   !! History :   1.0  !  2004-03 (C. Ethe) Original code
7   !!             3.3  !  2010-07 (C. Ethe) Merge TRA-TRC
8   !!----------------------------------------------------------------------
9#if defined key_top
10   !!----------------------------------------------------------------------
11   !!   'key_top'                                                TOP models
12   !!----------------------------------------------------------------------
13   !!   trc_trp        : passive tracer transport
14   !!----------------------------------------------------------------------
15   USE oce_trc         ! ocean dynamics and active tracers variables
16   USE trc             ! ocean passive tracers variables
17   USE trcbbl          ! bottom boundary layer               (trc_bbl routine)
18   USE trcdmp          ! internal damping                    (trc_dmp routine)
19   USE trcldf          ! lateral mixing                      (trc_ldf routine)
20   USE trcadv          ! advection                           (trc_adv routine)
21   USE trczdf          ! vertical diffusion                  (trc_zdf routine)
22   USE trcnxt          ! time-stepping                       (trc_nxt routine)
23   USE trcrad          ! positivity                          (trc_rad routine)
24   USE trcsbc          ! surface boundary condition          (trc_sbc routine)
25   USE trcbc           ! Tracers boundary condtions          ( trc_bc routine)
26   USE zpshde          ! partial step: hor. derivative       (zps_hde routine)
27   USE bdy_oce   , ONLY: ln_bdy
28   USE trcbdy          ! BDY open boundaries
29
30#if defined key_agrif
31   USE agrif_top_sponge ! tracers sponges
32#endif
33
34   IMPLICIT NONE
35   PRIVATE
36
37   PUBLIC   trc_trp    ! called by trc_stp
38
39   !!----------------------------------------------------------------------
40   !! NEMO/TOP 4.0 , NEMO Consortium (2018)
41   !! $Id$
42   !! Software governed by the CeCILL license (see ./LICENSE)
43   !!----------------------------------------------------------------------
44
45CONTAINS
46
47   SUBROUTINE trc_trp( kt )
48      !!----------------------------------------------------------------------
49      !!                     ***  ROUTINE trc_trp  ***
50      !!                     
51      !! ** Purpose :   Management of passive tracers transport
52      !!
53      !! ** Method  : - Compute the passive tracers trends
54      !!              - Update the passive tracers
55      !!----------------------------------------------------------------------
56      INTEGER, INTENT( in ) ::  kt  ! ocean time-step index
57      !! ---------------------------------------------------------------------
58      !
59      IF( ln_timing )   CALL timing_start('trc_trp')
60      !
61      IF( .NOT. lk_c1d ) THEN
62         !
63                                CALL trc_sbc    ( kt )      ! surface boundary condition
64         IF( ln_trcbc .AND. lltrcbc .AND. kt /= nit000 )  &
65                                CALL trc_bc     ( kt )      ! tracers: surface and lateral Boundary Conditions
66         IF( ln_trabbl )        CALL trc_bbl    ( kt )      ! advective (and/or diffusive) bottom boundary layer scheme
67         IF( ln_trcdmp )        CALL trc_dmp    ( kt )      ! internal damping trends
68         IF( ln_bdy )           CALL trc_bdy_dmp( kt )      ! BDY damping trends
69                                CALL trc_adv    ( kt )      ! horizontal & vertical advection
70         !                                                         ! Partial top/bottom cell: GRADh( trb ) 
71         IF( ln_zps ) THEN
72           IF( ln_isfcav ) THEN ; CALL zps_hde_isf( kt, jptra, trb, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi )  ! both top & bottom
73           ELSE                 ; CALL zps_hde    ( kt, jptra, trb, gtru, gtrv )                                      !  only bottom
74           ENDIF
75         ENDIF
76         !                                                     
77                                CALL trc_ldf    ( kt )      ! lateral mixing
78#if defined key_agrif
79         IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_trc       ! tracers sponge
80#endif
81                                CALL trc_zdf    ( kt )      ! vertical mixing and after tracer fields
82                                CALL trc_nxt    ( kt )      ! tracer fields at next time step     
83         IF( ln_trcrad )        CALL trc_rad    ( kt )      ! Correct artificial negative concentrations
84         IF( ln_trcdmp_clo )    CALL trc_dmp_clo( kt )      ! internal damping trends on closed seas only
85
86         !
87      ELSE                                               ! 1D vertical configuration
88                                CALL trc_sbc( kt )            ! surface boundary condition
89         IF( ln_trcdmp )        CALL trc_dmp( kt )            ! internal damping trends
90                                CALL trc_zdf( kt )            ! vertical mixing and after tracer fields
91                                CALL trc_nxt( kt )            ! tracer fields at next time step     
92          IF( ln_trcrad )       CALL trc_rad( kt )            ! Correct artificial negative concentrations
93         !
94      END IF
95      !
96      IF( ln_timing )   CALL timing_stop('trc_trp')
97      !
98   END SUBROUTINE trc_trp
99
100#else
101   !!----------------------------------------------------------------------
102   !!   Dummy module :                                        No TOP models
103   !!----------------------------------------------------------------------
104CONTAINS
105   SUBROUTINE trc_trp( kt )              ! Empty routine
106      INTEGER, INTENT(in) ::   kt
107      WRITE(*,*) 'trc_trp: You should not have seen this print! error?', kt
108   END SUBROUTINE trc_trp
109#endif
110   
111   !!======================================================================
112END MODULE trctrp
Note: See TracBrowser for help on using the repository browser.