Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/trc_oce.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/OPA_SRC/trc_oce.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

@@ -13 +13 @@
    USE in_out_manager  ! I/O manager
    USE dom_oce         ! ocean space and time domain
+   USE lib_mpp         ! MPP library
 
    IMPLICIT NONE
@@ -20 +21 @@
    PUBLIC   trc_oce_rgb_read   ! routine called by traqsr.F90
    PUBLIC   trc_oce_ext_lev    ! function called by traqsr.F90 at least
- 
-   REAL(wp), PUBLIC                          ::   r_si2   !: largest depth of extinction (blue & 0.01 mg.m-3)  (RGB)
-   REAL(wp), PUBLIC , DIMENSION(jpi,jpj,jpk) ::   etot3   !: light absortion coefficient
+   PUBLIC   trc_oce_alloc      ! function called by nemogcm.F90
+
+   REAL(wp), PUBLIC                                      ::   r_si2   !: largest depth of extinction (blue & 0.01 mg.m-3)  (RGB)
+   REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) ::   etot3   !: light absortion coefficient
 
 #if defined key_top && defined key_pisces
@@ -53 +55 @@
    !! NEMO/OPA 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
+
+   INTEGER FUNCTION trc_oce_alloc()
+      !!----------------------------------------------------------------------
+      !!                  ***  trc_oce_alloc  ***
+      !!----------------------------------------------------------------------
+      ALLOCATE( etot3(jpi,jpj,jpk)   , STAT= trc_oce_alloc )
+      !
+      IF( trc_oce_alloc /= 0 )   CALL ctl_warn('trc_oce_alloc: failed to allocate etot3 array')
+   END FUNCTION trc_oce_alloc
+
 
    SUBROUTINE trc_oce_rgb( prgb )
@@ -239 +250 @@
       !
       ! It is not necessary to compute anything bellow the following depth
-      zhext = prldex * ( LOG(10.e0) * zprec + LOG(pqsr_frc) )
-     
+      zhext = prldex * ( LOG(10._wp) * zprec + LOG(pqsr_frc) )
+      !
       ! Level of light extinction
       pjl = jpkm1
@@ -250 +261 @@
    END FUNCTION trc_oce_ext_lev
 
-
    !!======================================================================
 END MODULE trc_oce