New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trctrp.F90 in branches/UKMO/r5936_CO6_CO5_shelfdiagnostic/NEMOGCM/NEMO/TOP_SRC/TRP – NEMO

source: branches/UKMO/r5936_CO6_CO5_shelfdiagnostic/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90 @ 7113

Last change on this file since 7113 was 7113, checked in by jcastill, 7 years ago

Remove again the svn keywords, as it did not work before

File size: 5.7 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 trabbl          ! bottom boundary layer               (trc_bbl routine)
18   USE trcbbl          ! bottom boundary layer               (trc_bbl routine)
19   USE trcdmp          ! internal damping                    (trc_dmp routine)
20   USE trcldf          ! lateral mixing                      (trc_ldf routine)
21   USE trcadv          ! advection                           (trc_adv routine)
22   USE trczdf          ! vertical diffusion                  (trc_zdf routine)
23   USE trcnxt          ! time-stepping                       (trc_nxt routine)
24   USE trcrad          ! positivity                          (trc_rad routine)
25   USE trcsbc          ! surface boundary condition          (trc_sbc routine)
26   USE zpshde          ! partial step: hor. derivative       (zps_hde routine)
27
28#if defined key_agrif
29   USE agrif_top_sponge ! tracers sponges
30   USE agrif_top_update ! tracers updates
31#endif
32
33   IMPLICIT NONE
34   PRIVATE
35
36   PUBLIC   trc_trp    ! called by trc_stp
37
38   !!----------------------------------------------------------------------
39   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
40   !! $Id$
41   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
42   !!----------------------------------------------------------------------
43
44CONTAINS
45
46   SUBROUTINE trc_trp( kt )
47      !!----------------------------------------------------------------------
48      !!                     ***  ROUTINE trc_trp  ***
49      !!                     
50      !! ** Purpose :   Management of passive tracers transport
51      !!
52      !! ** Method  : - Compute the passive tracers trends
53      !!              - Update the passive tracers
54      !!----------------------------------------------------------------------
55      INTEGER, INTENT( in ) ::  kt  ! ocean time-step index
56      !! ---------------------------------------------------------------------
57      !
58      IF( nn_timing == 1 )   CALL timing_start('trc_trp')
59      !
60      IF( .NOT. lk_c1d ) THEN
61         !
62                                CALL trc_sbc    ( kt )      ! surface boundary condition
63         IF( lk_trabbl )        CALL trc_bbl    ( kt )      ! advective (and/or diffusive) bottom boundary layer scheme
64         IF( ln_trcdmp )        CALL trc_dmp    ( kt )      ! internal damping trends
65         IF( ln_trcdmp_clo )    CALL trc_dmp_clo( kt )      ! internal damping trends on closed seas only
66                                CALL trc_adv    ( kt )      ! horizontal & vertical advection
67         !                                                         ! Partial top/bottom cell: GRADh( trb ) 
68         IF( ln_zps ) THEN
69           IF( ln_isfcav ) THEN ; CALL zps_hde_isf( kt, jptra, trb, pgtu=gtru, pgtv=gtrv, pgtui=gtrui, pgtvi=gtrvi )  ! both top & bottom
70           ELSE                 ; CALL zps_hde    ( kt, jptra, trb, gtru, gtrv )                                      !  only bottom
71           ENDIF
72         ENDIF
73         !                                                     
74                                CALL trc_ldf    ( kt )      ! lateral mixing
75#if defined key_agrif
76         IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_trc       ! tracers sponge
77#endif
78                                CALL trc_zdf    ( kt )      ! vertical mixing and after tracer fields
79                                CALL trc_nxt    ( kt )      ! tracer fields at next time step     
80         IF( ln_trcrad )        CALL trc_rad    ( kt )      ! Correct artificial negative concentrations
81
82#if defined key_agrif
83         IF( .NOT.Agrif_Root()) CALL Agrif_Update_Trc( kt ) ! Update tracer at AGRIF zoom boundaries : children only
84#endif
85         !
86      ELSE                                               ! 1D vertical configuration
87                                CALL trc_sbc( kt )            ! surface boundary condition
88         IF( ln_trcdmp )        CALL trc_dmp( kt )            ! internal damping trends
89                                CALL trc_zdf( kt )            ! vertical mixing and after tracer fields
90                                CALL trc_nxt( kt )            ! tracer fields at next time step     
91          IF( ln_trcrad )       CALL trc_rad( kt )            ! Correct artificial negative concentrations
92         !
93      END IF
94      !
95      IF( nn_timing == 1 )   CALL timing_stop('trc_trp')
96      !
97   END SUBROUTINE trc_trp
98
99#else
100   !!----------------------------------------------------------------------
101   !!   Dummy module :                                        No TOP models
102   !!----------------------------------------------------------------------
103CONTAINS
104   SUBROUTINE trc_trp( kt )              ! Empty routine
105      INTEGER, INTENT(in) ::   kt
106      WRITE(*,*) 'trc_trp: You should not have seen this print! error?', kt
107   END SUBROUTINE trc_trp
108#endif
109   
110   !!======================================================================
111END MODULE trctrp
Note: See TracBrowser for help on using the repository browser.