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

Timestamp:

2011-03-30T17:58:35+02:00 (13 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/C14b/trcini_c14b.F90 (modified) (10 diffs)

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

@@ -4 +4 @@
    !! TOP :   initialisation of the C14 bomb tracer
    !!======================================================================
-   !! History : Original ! 2005-10  (Z. Lachkar) 
-   !!               2.0  ! 2007-12  (C. Ethe ) 
+   !! History :  1.0  ! 2005-10  (Z. Lachkar) Original code
+   !!            2.0  ! 2007-12  (C. Ethe) 
    !!----------------------------------------------------------------------
 #if defined key_c14b
@@ -17 +17 @@
    USE trc             ! TOP variables
    USE trcsms_c14b     ! C14 sms trends
-   USE in_out_manager  ! I/O manager
 
    IMPLICIT NONE
@@ -24 +23 @@
    PUBLIC   trc_ini_c14b   ! called by trcini.F90 module
 
-   INTEGER  ::   &     ! With respect to data file !!
-     jpybeg = 1765 , & !: starting year for C14
-     jpyend = 2002     !: ending year for C14
-
-   INTEGER  ::   &   
-      nrec   ,  & ! number of year in CO2 Concentrations file
-      nmaxrec 
-
-   INTEGER  ::   inum1, inum2               ! unit number
-
-   REAL(wp) ::     &
-     ys40 = -40. ,    &             ! 40 degrees south
-     ys20 = -20. ,    &             ! 20 degrees south
-     yn20 =  20. ,    &             ! 20 degrees north
-     yn40 =  40.                    ! 40 degrees north
-
-   !!---------------------------------------------------------------------
+   !                             ! With respect to data file !!
+   INTEGER  ::   jpybeg = 1765   ! starting year for C14
+   INTEGER  ::   jpyend = 2002   ! ending year for C14
+   INTEGER  ::   nrec            ! number of year in CO2 Concentrations file
+   INTEGER  ::   nmaxrec 
+   INTEGER  ::   inum1, inum2    ! unit number
+
+   REAL(wp) ::   ys40 = -40.     ! 40 degrees south
+   REAL(wp) ::   ys20 = -20.     ! 20 degrees south
+   REAL(wp) ::   yn20 =  20.     ! 20 degrees north
+   REAL(wp) ::   yn40 =  40.     ! 40 degrees north
+
+   !!----------------------------------------------------------------------
    !! NEMO/TOP 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
-   !!----------------------------------------------------------------------
-
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
+   !!----------------------------------------------------------------------
 CONTAINS
 
@@ -59 +53 @@
       !!----------------------------------------------------------------------
 
-      !  Control consitency
-      CALL trc_ctl_c14b
+      !                     ! Allocate C14b arrays
+      IF( trc_sms_c14b_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'trc_ini_c14b : unable to allocate C14b arrays' )
+
+      CALL trc_ctl_c14b     !  Control consitency
 
       IF(lwp) WRITE(numout,*) ''
@@ -69 +65 @@
       ! Initialization of boundaries conditions
       ! --------------------------------------- 
-      qtr_c14(:,:) = 0.e0
+      qtr_c14(:,:) = 0._wp
       
       ! Initialization of qint in case of  no restart 
@@ -78 +74 @@
             WRITE(numout,*) 'Initialization de qint ; No restart : qint equal zero '
          ENDIF
-         trn     (:,:,:,jpc14) = 0.e0
-         qint_c14(:,:        ) = 0.e0
+         trn     (:,:,:,jpc14) = 0._wp
+         qint_c14(:,:        ) = 0._wp
       ENDIF
 
@@ -156 +152 @@
                  fareaz(ji,jj,3) = 0.
             ENDIF
-          END DO
-        END DO
-
+         END DO
+      END DO
       !
       IF(lwp) WRITE(numout,*) 'Initialization of C14 bomb tracer done'
       IF(lwp) WRITE(numout,*) ' '
-
+      !
    END SUBROUTINE trc_ini_c14b
-   
+
+
    SUBROUTINE trc_ctl_c14b
       !!----------------------------------------------------------------------
@@ -179 +175 @@
       ! Check number of tracers
       ! -----------------------   
-      IF( jp_c14b > 1) CALL ctl_stop( ' Change jp_c14b to be equal 1 in par_c14b.F90' )
+      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'
+      IF( ctrcnm(jpc14) /= 'C14B' ) THEN
+          ctrcnm(jpc14)  = 'C14B'
+          ctrcnl(jpc14)  = 'Bomb C14 concentration'
       ENDIF
 
@@ -197 +193 @@
       ! ------------------
       IF( ctrcun(jpc14) /= 'ration' ) THEN
-          ctrcun(jpc14) = 'ration'
+          ctrcun(jpc14)  = 'ration'
           IF(lwp) THEN
              CALL ctl_warn( ' we force tracer unit' )
@@ -206 +202 @@
       !
    END SUBROUTINE trc_ctl_c14b
+   
 #else
    !!----------------------------------------------------------------------