source: NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/TOP/MY_TRC/trcsms_my_trc.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: 3.3 KB
Line 
1MODULE trcsms_my_trc
2   !!======================================================================
3   !!                         ***  MODULE trcsms_my_trc  ***
4   !! TOP :   Main module of the MY_TRC tracers
5   !!======================================================================
6   !! History :      !  2007  (C. Ethe, G. Madec)  Original code
7   !!                !  2016  (C. Ethe, T. Lovato) Revised architecture
8   !!----------------------------------------------------------------------
9   !! trc_sms_my_trc       : MY_TRC model main routine
10   !! trc_sms_my_trc_alloc : allocate arrays specific to MY_TRC sms
11   !!----------------------------------------------------------------------
12   USE par_trc         ! TOP parameters
13   USE oce_trc         ! Ocean variables
14   USE trc             ! TOP variables
15   USE trd_oce
16   USE trdtrc
17
18   IMPLICIT NONE
19   PRIVATE
20
21   PUBLIC   trc_sms_my_trc       ! called by trcsms.F90 module
22   PUBLIC   trc_sms_my_trc_alloc ! called by trcini_my_trc.F90 module
23
24   ! Defined HERE the arrays specific to MY_TRC sms and ALLOCATE them in trc_sms_my_trc_alloc
25
26   !!----------------------------------------------------------------------
27   !! NEMO/TOP 4.0 , NEMO Consortium (2018)
28   !! $Id$
29   !! Software governed by the CeCILL license (see ./LICENSE)
30   !!----------------------------------------------------------------------
31CONTAINS
32
33   SUBROUTINE trc_sms_my_trc( kt )
34      !!----------------------------------------------------------------------
35      !!                     ***  trc_sms_my_trc  ***
36      !!
37      !! ** Purpose :   main routine of MY_TRC model
38      !!
39      !! ** Method  : -
40      !!----------------------------------------------------------------------
41      !
42      INTEGER, INTENT(in) ::   kt   ! ocean time-step index
43      INTEGER ::   jn   ! dummy loop index
44      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrmyt
45      !!----------------------------------------------------------------------
46      !
47      IF( ln_timing )   CALL timing_start('trc_sms_my_trc')
48      !
49      IF(lwp) WRITE(numout,*)
50      IF(lwp) WRITE(numout,*) ' trc_sms_my_trc:  MY_TRC model'
51      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
52
53      IF( l_trdtrc )  ALLOCATE( ztrmyt(jpi,jpj,jpk) )
54
55      ! add here the call to BGC model
56
57      ! Save the trends in the mixed layer
58      IF( l_trdtrc ) THEN
59          DO jn = jp_myt0, jp_myt1
60            ztrmyt(:,:,:) = tra(:,:,:,jn)
61            CALL trd_trc( ztrmyt, jn, jptra_sms, kt )   ! save trends
62          END DO
63          DEALLOCATE( ztrmyt )
64      END IF
65      !
66      IF( ln_timing )   CALL timing_stop('trc_sms_my_trc')
67      !
68   END SUBROUTINE trc_sms_my_trc
69
70
71   INTEGER FUNCTION trc_sms_my_trc_alloc()
72      !!----------------------------------------------------------------------
73      !!              ***  ROUTINE trc_sms_my_trc_alloc  ***
74      !!----------------------------------------------------------------------
75      !
76      ! ALLOCATE here the arrays specific to MY_TRC
77      ! ALLOCATE( tab(...) , STAT=trc_sms_my_trc_alloc )
78      trc_sms_my_trc_alloc = 0      ! set to zero if no array to be allocated
79      !
80      IF( trc_sms_my_trc_alloc /= 0 ) CALL ctl_stop( 'STOP', 'trc_sms_my_trc_alloc : failed to allocate arrays' )
81      !
82   END FUNCTION trc_sms_my_trc_alloc
83
84   !!======================================================================
85END MODULE trcsms_my_trc
Note: See TracBrowser for help on using the repository browser.