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.
trcini_my_trc.F90 in branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC – NEMO

source: branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90 @ 6772

Last change on this file since 6772 was 6772, checked in by cbricaud, 8 years ago

clean in coarsening branch

  • Property svn:keywords set to Id
File size: 3.1 KB
Line 
1MODULE trcini_my_trc
2   !!======================================================================
3   !!                         ***  MODULE trcini_my_trc  ***
4   !! TOP :   initialisation of the MY_TRC tracers
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) Original code
7   !!----------------------------------------------------------------------
8#if defined key_my_trc
9   !!----------------------------------------------------------------------
10   !!   'key_my_trc'                                               CFC tracers
11   !!----------------------------------------------------------------------
12   !! trc_ini_my_trc   : MY_TRC model initialisation
13   !!----------------------------------------------------------------------
14   USE par_trc         ! TOP parameters
15   USE oce_trc
16   USE trc
17   USE trcsms_my_trc
18   USE dom_oce, ONLY : gdepw_1d,e3t_1d,nyear_len
19
20   IMPLICIT NONE
21   PRIVATE
22
23   PUBLIC   trc_ini_my_trc   ! called by trcini.F90 module
24
25   !!----------------------------------------------------------------------
26   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
27   !! $Id$
28   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
29   !!----------------------------------------------------------------------
30CONTAINS
31
32   SUBROUTINE trc_ini_my_trc
33      !!----------------------------------------------------------------------
34      !!                     ***  trc_ini_my_trc  *** 
35      !!
36      !! ** Purpose :   initialization for MY_TRC model
37      !!
38      !! ** Method  : - Read the namcfc namelist and check the parameter values
39      !!----------------------------------------------------------------------
40
41      !                       ! Allocate MY_TRC arrays
42      IF( trc_sms_my_trc_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'trc_ini_my_trc: unable to allocate MY_TRC arrays' )
43
44      IF(lwp) WRITE(numout,*)
45      IF(lwp) WRITE(numout,*) ' trc_ini_my_trc: passive tracer unit vector'
46      IF(lwp) WRITE(numout,*) ' To check conservation : '
47      IF(lwp) WRITE(numout,*) '   1 - No sea-ice model '
48      IF(lwp) WRITE(numout,*) '   2 - No runoff ' 
49      IF(lwp) WRITE(numout,*) '   3 - precipitation and evaporation equal to 1 : E=P=1 ' 
50      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
51
52      nlb_age = MINLOC( gdepw_1d, mask = gdepw_1d > age_depth, dim = 1  )
53      nl_age = nlb_age - 1
54      nla_age = nl_age - 1
55      frac_kill_age = (age_depth - gdepw_1d(nl_age))/e3t_1d(nl_age)
56      frac_add_age = 1._wp -  frac_kill_age
57      rryear = 1._wp / ( nyear_len(1) * rday )
58
59     
60      IF( .NOT. ln_rsttr ) trb(:,:,:,jp_myt0:jp_myt1) = 0._wp
61      IF( .NOT. ln_rsttr ) trn(:,:,:,jp_myt0:jp_myt1) = 0._wp
62      !
63   END SUBROUTINE trc_ini_my_trc
64
65#else
66   !!----------------------------------------------------------------------
67   !!   Dummy module                                        No MY_TRC model
68   !!----------------------------------------------------------------------
69CONTAINS
70   SUBROUTINE trc_ini_my_trc             ! Empty routine
71   END SUBROUTINE trc_ini_my_trc
72#endif
73
74   !!======================================================================
75END MODULE trcini_my_trc
Note: See TracBrowser for help on using the repository browser.