source: branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsms_medusa.F90 @ 8521

Last change on this file since 8521 was 8521, checked in by marc, 3 years ago

gmed:#346 Merging branches/NERC/dev_r5518_GO6_MEDUSA_conserv

File size: 4.1 KB
Line 
1MODULE trcsms_medusa
2   !!======================================================================
3   !!                         ***  MODULE trcsms_medusa  ***
4   !! TOP :   Main module of the MEDUSA tracers
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) Original code
7   !!              -   !  2008-08  (K. Popova) adaptation for MEDUSA
8   !!              -   !  2008-11  (A. Yool) continuing adaptation for MEDUSA
9   !!              -   !  2010-03  (A. Yool) updated for branch inclusion
10   !!              -   !  2017-08  (A. Yool) amend for slow detritus bug
11   !!----------------------------------------------------------------------
12#if defined key_medusa
13   !!----------------------------------------------------------------------
14   !!   'key_medusa'                                       bio tracers
15   !!----------------------------------------------------------------------
16   !! trc_sms_medusa   : MEDUSA_TRC model main routine
17   !!----------------------------------------------------------------------
18   USE par_trc         ! TOP parameters
19   USE oce_trc
20   USE trc
21   USE trcbio_medusa
22   USE trcopt_medusa
23   USE trcsed_medusa
24   USE trcavg_medusa
25
26
27   IMPLICIT NONE
28   PRIVATE
29
30   PUBLIC   trc_sms_medusa   ! called by trcsms.F90 module
31
32   !!----------------------------------------------------------------------
33   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
34   !! $Id$
35   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
36   !!----------------------------------------------------------------------
37
38CONTAINS
39
40   SUBROUTINE trc_sms_medusa( kt )
41      !!----------------------------------------------------------------------
42      !!                     ***  trc_sms_medusa  *** 
43      !!
44      !! ** Purpose :   main routine of MEDUSA_TRC model
45      !!
46      !! ** Method  : -
47      !!----------------------------------------------------------------------
48      INTEGER, INTENT(in) :: kt   ! ocean time-step index
49
50# if defined key_debug_medusa
51         IF(lwp) WRITE(numout,*) ' MEDUSA inside trc_sms_medusa'
52         CALL flush(numout)
53# endif
54
55      IF( kt == nittrc000 ) THEN
56       IF(lwp) WRITE(numout,*)
57       IF(lwp) WRITE(numout,*) ' trc_sms_medusa:  MEDUSA model'
58       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
59      ENDIF
60
61      CALL trc_avg_medusa( kt ) ! rolling average module
62# if defined key_debug_medusa
63      IF(lwp) WRITE(numout,*) ' MEDUSA done trc_avg_medusa'
64      CALL flush(numout)
65# endif
66     
67      CALL trc_opt_medusa( kt ) ! optical model
68# if defined key_debug_medusa
69      IF(lwp) WRITE(numout,*) ' MEDUSA done trc_opt_medusa'
70      CALL flush(numout)
71# endif
72
73      !! AXY & JPALM (28/02/17): call dust before trc_bio_medusa (because of coupling)
74      CALL trc_sed_medusa_dust( kt ) ! dust submodel
75# if defined key_debug_medusa
76      IF(lwp) WRITE(numout,*) ' MEDUSA done trc_sed_medusa_dust'
77      CALL flush(numout)
78# endif
79
80# if defined key_kill_medusa
81      !! MEDUSA skipped
82      IF(lwp) WRITE(numout,*) ' MEDUSA killed at kt =', kt
83      CALL flush(numout)
84# else
85      CALL trc_bio_medusa( kt ) ! biological model
86#  if defined key_debug_medusa
87      IF(lwp) WRITE(numout,*) ' MEDUSA done trc_bio_medusa'
88      CALL flush(numout)
89#  endif
90     
91!! AXY (08/08/2017): remove call to buggy subroutine (now handled by detritus.F90)
92!!       CALL trc_sed_medusa( kt ) ! sedimentation model
93!! #  if defined key_debug_medusa
94!!       IF(lwp) WRITE(numout,*) ' MEDUSA done trc_sed_medusa'
95!!       CALL flush(numout)
96!! #  endif
97# endif
98
99   END SUBROUTINE trc_sms_medusa
100   
101#else
102   !!----------------------------------------------------------------------
103   !!   Dummy module                                        No MEDUSA model
104   !!----------------------------------------------------------------------
105CONTAINS
106   SUBROUTINE trc_sms_medusa( kt )             ! Empty routine
107      INTEGER, INTENT( in ) ::   kt
108      WRITE(*,*) 'trc_sms_medusa: You should not have seen this print! error?', kt
109   END SUBROUTINE trc_sms_medusa
110#endif
111
112   !!======================================================================
113END MODULE trcsms_medusa
114
Note: See TracBrowser for help on using the repository browser.