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

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zint.F90

@@ -21 +21 @@
 
    PUBLIC   p4z_int  
+   PUBLIC   p4z_int_alloc
 
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &   !:
-      tgfunc,            &  !:  Temp. dependancy of various biological rates
-      tgfunc2               !:  Temp. dependancy of mesozooplankton rates 
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   tgfunc    !: Temp. dependancy of various biological rates
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   tgfunc2   !: Temp. dependancy of mesozooplankton rates
 
-   !! * Module variables
-   REAL(wp) :: &
-      xksilim = 16.5E-6        ! Half-saturation constant for the computation of the Si half-saturation constant
-
+   REAL(wp) ::   xksilim = 16.5e-6_wp   ! Half-saturation constant for the Si half-saturation constant computation
 
    !!----------------------------------------------------------------------
    !! 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
 
@@ -47 +43 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      !!
       INTEGER  ::   ji, jj
       REAL(wp) ::   zdum
@@ -54 +49 @@
       ! Computation of phyto and zoo metabolic rate
       ! -------------------------------------------
-
       tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) )
       tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) )
@@ -61 +55 @@
       ! constant for silica uptake
       ! ---------------------------------------------------
-
       DO ji = 1, jpi
          DO jj = 1, jpj
@@ -68 +61 @@
          END DO
       END DO
-
+      !
       IF( nday_year == nyear_len(1) ) THEN
          xksi    = xksimax
-         xksimax = 0.e0
+         xksimax = 0._wp
       ENDIF
       !
    END SUBROUTINE p4z_int
+
+
+   INTEGER FUNCTION p4z_int_alloc()
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_int_alloc  ***
+      !!----------------------------------------------------------------------
+      ALLOCATE( tgfunc(jpi,jpj,jpk), tgfunc2(jpi,jpj,jpk), STAT=p4z_int_alloc )
+      !
+      IF( p4z_int_alloc /= 0 )   CALL ctl_warn('p4z_int_alloc : failed to allocate arrays.')
+      !
+   END FUNCTION p4z_int_alloc
 
 #else