source: branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcini_idtra.F90 @ 6164

Last change on this file since 6164 was 6164, checked in by jpalmier, 5 years ago

JPALM — 23-12-2015 — 1_ adapt CFC in MEDUSA branch - now Working properly 2_ add diagnostics to Ideal tracer

File size: 4.1 KB
Line 
1MODULE trcini_idtra
2   !!======================================================================
3   !!                         ***  MODULE trcini_idtra  ***
4   !! TOP :   initialisation of the IDEAL-TRACER tracers
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) from trcini.idtra.h90
7   !!----------------------------------------------------------------------
8#if defined key_idtra
9   !!----------------------------------------------------------------------
10   !!   'key_idtra'                                               IDEAL-TRACER tracers
11   !!----------------------------------------------------------------------
12   !! trc_ini_idtra      : IDEAL-TRACER model initialisation
13   !!----------------------------------------------------------------------
14   USE oce_trc         ! Ocean variables
15   USE par_trc         ! TOP parameters
16   USE trc             ! TOP variables
17   USE trcsms_idtra    ! IDEAL-TRACER sms trends
18   ! USE par_idtra       ! IDEAL-TRACER parameters
19   ! USE in_out_manager  ! I/O manager
20   ! USE lib_mpp
21   ! USE iom
22
23   IMPLICIT NONE
24   PRIVATE
25
26   PUBLIC   trc_ini_idtra   ! called by trcini.F90 module
27
28   INTEGER  ::   inum                   ! unit number
29
30   !!----------------------------------------------------------------------
31   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
32   !! $Id$
33   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
34   !!----------------------------------------------------------------------
35CONTAINS
36
37   SUBROUTINE trc_ini_idtra
38      !!----------------------------------------------------------------------
39      !!                     ***  trc_ini_idtra  ***
40      !!
41      !! ** Purpose :   initialization for idtra model
42      !!
43      !! ** Method  : - Read the namidtra namelist and check the parameter values
44      !!----------------------------------------------------------------------
45      INTEGER  ::    jn, jl 
46      !!----------------------------------------------------------------------
47
48      IF(lwp) WRITE(numout,*)
49      IF(lwp) WRITE(numout,*) ' trc_ini_idtra: initialisation of Ideal Tracers model'
50      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
51
52      IF( trc_sms_idtra_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'trc_ini_idtra:unable to allocate CFC arrays' )
53
54
55      ! Initialization of trn in case of  no restart
56      !----------------------------------------------
57      qtr_idtra(:,:,:) = 0._wp
58      inv_idtra(:,:,:) = 0._wp
59      IF( .NOT. ln_rsttr ) THEN
60         IF(lwp) THEN
61            WRITE(numout,*)
62            WRITE(numout,*) 'Initialization of id-tracers ; No restart : '
63            WRITE(numout,*) '                             ; Init field equal 1 at surface - zero elsewhere'
64            WRITE(numout,*) '                             ; qint idtra equal 0 '
65         ENDIF
66         qint_idtra(:,:,:) = 0._wp
67         DO jn = jp_idtra0, jp_idtra1
68             trn(:,:,:,jn) = 0.e0
69             trn(:,:,1,jn) = 1.0
70           IF(lwp) WRITE(numout,*) 'Idealise Tracer initialisation -- jn = ',jn
71         END DO
72      ENDIF
73
74
75      !   Ideal traceur do not need any atmospheric concentration.
76      ! We consider that sucface concentration is equal to 1,
77      ! that it is advectied within the water circulation,
78      ! and that it is regularly degraded as if it was a radiactive tracer (tricium for example)
79      ! But we can play with tha caracteristic time of
80      !--------------------------------------------------------------------
81
82
83
84      IF(lwp) WRITE(numout,*) 'Initialization of IDEAL-TRACER tracers done'
85      IF(lwp) WRITE(numout,*) ' '
86
87   END SUBROUTINE trc_ini_idtra
88
89#else
90   !!----------------------------------------------------------------------
91   !!   Dummy module                                         No IDEAL-TRACER tracers
92   !!----------------------------------------------------------------------
93CONTAINS
94   SUBROUTINE trc_ini_idtra             ! Empty routine
95
96
97   END SUBROUTINE trc_ini_idtra
98#endif
99
100   !!======================================================================
101END MODULE trcini_idtra
102                                               
103
104
105
Note: See TracBrowser for help on using the repository browser.