Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/C14b/trcini_c14b.F90

Timestamp:

2010-12-27T18:33:53+01:00 (14 years ago)

Author:

rblod

Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/C14b/trcini_c14b.F90 (modified) (3 diffs, 1 prop)

trunk/NEMOGCM/NEMO/TOP_SRC/C14b/trcini_c14b.F90

@@ -41 +41 @@
 
    !!---------------------------------------------------------------------
-   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
-   !! $Id: trcini_cfc.F90 1146 2008-06-25 11:42:56Z rblod $ 
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
+   !! $Id$ 
+   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -59 +59 @@
       !!----------------------------------------------------------------------
 
-      IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) ' trc_ini_cfc: initialisation of CFC chemical model'
+      !  Control consitency
+      CALL trc_ctl_c14b
+
+      IF(lwp) WRITE(numout,*) ''
+      IF(lwp) WRITE(numout,*) ' trc_ini_c14b: initialisation of Bomb C14 chemical model'
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
 
@@ -162 +165 @@
    END SUBROUTINE trc_ini_c14b
    
+   SUBROUTINE trc_ctl_c14b
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE trc_ctl_c14b  ***
+      !!
+      !! ** Purpose :   control the cpp options, namelist and files 
+      !!----------------------------------------------------------------------
+
+      IF(lwp) THEN
+          WRITE(numout,*) ' C14 bomb Model '
+          WRITE(numout,*) ' '
+      ENDIF
+
+      ! Check number of tracers
+      ! -----------------------   
+      IF( jp_c14b > 1) CALL ctl_stop( ' Change jp_c14b to be equal 1 in par_c14b.F90' )
+
+      ! Check tracer names
+      ! ------------------
+      IF ( ctrcnm(jpc14) /= 'C14B' ) THEN
+           ctrcnm(jpc14)  = 'C14B'
+           ctrcnl(jpc14)  = 'Bomb C14 concentration'
+      ENDIF
+
+      IF(lwp) THEN
+         CALL ctl_warn( ' we force tracer names' )
+         WRITE(numout,*) ' tracer nb: ',jpc14,' name = ',ctrcnm(jpc14), ctrcnl(jpc14)
+         WRITE(numout,*) ' '
+      ENDIF
+
+      ! Check tracer units
+      ! ------------------
+      IF( ctrcun(jpc14) /= 'ration' ) THEN
+          ctrcun(jpc14) = 'ration'
+          IF(lwp) THEN
+             CALL ctl_warn( ' we force tracer unit' )
+             WRITE(numout,*) ' tracer  ',ctrcnm(jpc14), 'UNIT= ',ctrcun(jpc14)
+             WRITE(numout,*) ' '
+          ENDIF
+       ENDIF
+      !
+   END SUBROUTINE trc_ctl_c14b
 #else
    !!----------------------------------------------------------------------