source: branches/NERC/dev_r5107_NOC_MEDUSA/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcnam_idtra.F90 @ 5707

Last change on this file since 5707 was 5707, checked in by acc, 5 years ago

JPALM —25-08-2015 — add MEDUSA in the branch. MEDUSA version already up-to-date with this trunk revision

File size: 4.7 KB
Line 
1MODULE trcnam_idtra
2   !!======================================================================
3   !!                         ***  MODULE trcnam_idtra  ***
4   !! TOP :   initialisation of some run parameters for IDEAL-TRACER chemical model
5   !!======================================================================
6   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) from trcnam.idtra.h90
7   !!----------------------------------------------------------------------
8#if defined key_idtra
9   !!----------------------------------------------------------------------
10   !!   'key_idtra'                                               IDEAL-TRACER tracers
11   !!----------------------------------------------------------------------
12   !! trc_nam_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 specific variable
18   USE in_out_manager  ! I/O manager
19
20   IMPLICIT NONE
21   PRIVATE
22
23   PUBLIC   trc_nam_idtra   ! called by trcnam.F90 module
24
25   !!----------------------------------------------------------------------
26   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
27   !! $Id: trcnam_idtra.F90 1581 2009-08-05 14:53:12Z smasson $
28   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
29   !!----------------------------------------------------------------------
30
31CONTAINS
32
33   SUBROUTINE trc_nam_idtra
34      !!-------------------------------------------------------------------
35      !!                  ***  ROUTINE trc_nam_idtra  ***
36      !!
37      !! ** Purpose :   Definition some run parameter for IDEAL-TRACER model
38      !!
39      !! ** Method  :   Read the namidtra namelist and check the parameter
40      !!       values called at the first timestep (nit000)
41      !!
42      !! ** input   :   Namelist namidtra
43      !!----------------------------------------------------------------------
44      REAL(wp) ::   tmp_decay          !! Years ; half time decay of our idealize tracer
45      REAL(wp) ::   TDECyr, TDEC   
46      CHARACTER(LEN=32)   ::   clname
47      !!
48      NAMELIST/namidtra/tmp_decay
49!! #if defined key_trc_diaadd
50!!      ! definition of additional diagnostic as a structure
51!!      INTEGER :: jl, jn
52!!
53!! #endif
54      !!
55
56!!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
57!! Jpalm -- 4-11-2014
58!! namelist for idealize tracer
59!! only thing in namelist is the chosen half time decay
60!! no atmospheric conditions, cause we do impose a surface concentration of 1,
61!! and no additionnal diagnostics,
62!! because the only thing we are interested in is the water mass concentration on this tracer.
63!!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
64
65
66      IF(lwp) WRITE(numout,*)
67      clname = 'namelist_idtra'
68      IF(lwp) WRITE(numout,*) ' trc_nam_idtra: read IDEAL-TRACER namelist'
69      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
70      !!
71      !! Open the namelist file :
72      !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~
73      CALL ctl_opn( numnatm, TRIM( clname ), 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
74
75      !! Read the namelists :
76      !!~~~~~~~~~~~~~~~~~~~~~~~
77      !! First namelist of our idealize tracer :
78      !! read the decay 1/2 time of our tracer, to define in the namelist.
79      !! tmp_decay = 1y ; 10y ; 100y or 1000y depending of which water mass you want to track
80      !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
81      READ(numnatm,namidtra)
82
83      IF(lwp) WRITE(numout,*) '   -  half time decay of our idealize tracer : ', tmp_decay
84
85      ! decroissance radioactive du traceur ideal
86      ! ---------------------------------------
87      ! TDECyr = 12.43/LOG(2.)             !! Tricium as example
88       TDECyr = tmp_decay/LOG(2.)          !! Idealise tracer -- with tmp_decay given in the idtracer namelist
89       TDEC = TDECyr*365.*24.*60.*60.      !! translate in second
90       FDEC = EXP( -rdt/TDEC )
91
92
93!! #if defined key_trc_diaadd  && ! defined key_iomput
94      !!
95      !!  -Here you can add tracers names to be read
96      !! in a namelist.
97      !!  -But this is not necessary with the iomput module
98      !! cause names are written in the Iodef file.
99      !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
100
101!! #endif
102
103   END SUBROUTINE trc_nam_idtra
104
105#else
106   !!----------------------------------------------------------------------
107   !!  Dummy module :                                                No IDEAL-TRACER
108   !!----------------------------------------------------------------------
109CONTAINS
110   SUBROUTINE trc_nam_idtra                      ! Empty routine
111   END  SUBROUTINE  trc_nam_idtra
112#endif
113
114   !!======================================================================
115END MODULE trcnam_idtra
116
117
118
119
120
121
Note: See TracBrowser for help on using the repository browser.