source: NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/MY_TRC/trcsms_my_trc.F90 @ 10420

Last change on this file since 10420 was 10420, checked in by smasson, 21 months ago

dev_r10164_HPC09_ESIWACE_PREP_MERGE: force STOP when fail to allocate array, see #2133

  • Property svn:keywords set to Id
File size: 3.4 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   USE trcbc, only : trc_bc
18
19   IMPLICIT NONE
20   PRIVATE
21
22   PUBLIC   trc_sms_my_trc       ! called by trcsms.F90 module
23   PUBLIC   trc_sms_my_trc_alloc ! called by trcini_my_trc.F90 module
24
25   ! Defined HERE the arrays specific to MY_TRC sms and ALLOCATE them in trc_sms_my_trc_alloc
26
27   !!----------------------------------------------------------------------
28   !! NEMO/TOP 4.0 , NEMO Consortium (2018)
29   !! $Id$
30   !! Software governed by the CeCILL license (see ./LICENSE)
31   !!----------------------------------------------------------------------
32CONTAINS
33
34   SUBROUTINE trc_sms_my_trc( kt )
35      !!----------------------------------------------------------------------
36      !!                     ***  trc_sms_my_trc  ***
37      !!
38      !! ** Purpose :   main routine of MY_TRC model
39      !!
40      !! ** Method  : -
41      !!----------------------------------------------------------------------
42      !
43      INTEGER, INTENT(in) ::   kt   ! ocean time-step index
44      INTEGER ::   jn   ! dummy loop index
45      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrmyt
46      !!----------------------------------------------------------------------
47      !
48      IF( ln_timing )   CALL timing_start('trc_sms_my_trc')
49      !
50      IF(lwp) WRITE(numout,*)
51      IF(lwp) WRITE(numout,*) ' trc_sms_my_trc:  MY_TRC model'
52      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
53
54      IF( l_trdtrc )  ALLOCATE( ztrmyt(jpi,jpj,jpk) )
55
56      CALL trc_bc ( kt )       ! tracers: surface and lateral Boundary Conditions
57
58      ! add here the call to BGC model
59
60      ! Save the trends in the mixed layer
61      IF( l_trdtrc ) THEN
62          DO jn = jp_myt0, jp_myt1
63            ztrmyt(:,:,:) = tra(:,:,:,jn)
64            CALL trd_trc( ztrmyt, jn, jptra_sms, kt )   ! save trends
65          END DO
66          DEALLOCATE( ztrmyt )
67      END IF
68      !
69      IF( ln_timing )   CALL timing_stop('trc_sms_my_trc')
70      !
71   END SUBROUTINE trc_sms_my_trc
72
73
74   INTEGER FUNCTION trc_sms_my_trc_alloc()
75      !!----------------------------------------------------------------------
76      !!              ***  ROUTINE trc_sms_my_trc_alloc  ***
77      !!----------------------------------------------------------------------
78      !
79      ! ALLOCATE here the arrays specific to MY_TRC
80      ! ALLOCATE( tab(...) , STAT=trc_sms_my_trc_alloc )
81      trc_sms_my_trc_alloc = 0      ! set to zero if no array to be allocated
82      !
83      IF( trc_sms_my_trc_alloc /= 0 ) CALL ctl_stop( 'STOP', 'trc_sms_my_trc_alloc : failed to allocate arrays' )
84      !
85   END FUNCTION trc_sms_my_trc_alloc
86
87   !!======================================================================
88END MODULE trcsms_my_trc
Note: See TracBrowser for help on using the repository browser.