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

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

@@ -29 +29 @@
    PUBLIC   p4z_rem         ! called in p4zbio.F90
    PUBLIC   p4z_rem_init    ! called in trcsms_pisces.F90
-
-   !! * Shared module variables
+   PUBLIC   p4z_rem_alloc
+
    REAL(wp), PUBLIC ::   &
      xremik  = 0.3_wp      ,  & !:
@@ -39 +39 @@
      oxymin  = 1.e-6_wp         !:
 
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::    & !:
-     &                   denitr                     !: denitrification array
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   denitr   !: denitrification array
 
 
@@ -48 +47 @@
    !! 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
 
@@ -61 +59 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY:   ztempbac => wrk_2d_1
+      USE wrk_nemo, ONLY:   zdepbac  => wrk_3d_2 , zfesatur => wrk_3d_2 , zolimi => wrk_3d_4
+      !
       INTEGER, INTENT(in) ::   kt ! ocean time step
+      !
       INTEGER  ::   ji, jj, jk
       REAL(wp) ::   zremip, zremik , zlam1b
@@ -72 +75 @@
 #endif
       REAL(wp) ::   zlamfac, zonitr, zstep
-      REAL(wp), DIMENSION(jpi,jpj)     ::   ztempbac
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdepbac, zfesatur, zolimi
       CHARACTER (len=25) :: charout
-
       !!---------------------------------------------------------------------
 
+      IF(  wrk_in_use(2, 1)  .OR.  wrk_in_use(3, 2,3,4)  ) THEN
+         CALL ctl_stop('p4z_rem: requested workspace arrays unavailable')   ;   RETURN
+      ENDIF
 
        ! Initialisation of temprary arrys
-       zdepbac (:,:,:) = 0.0
-       zfesatur(:,:,:) = 0.0
-       zolimi  (:,:,:) = 0.0
-       ztempbac(:,:)   = 0.0
+       zdepbac (:,:,:) = 0._wp
+       zfesatur(:,:,:) = 0._wp
+       zolimi  (:,:,:) = 0._wp
+       ztempbac(:,:)   = 0._wp
 
       !  Computation of the mean phytoplankton concentration as
       !  a crude estimate of the bacterial biomass
       !   --------------------------------------------------
-
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -362 +364 @@
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zscave * zdenom2
 #endif
-
-            END DO
-         END DO
-      END DO
-      !
-
-       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
+            END DO
+         END DO
+      END DO
+      !
+
+      IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
          WRITE(charout, FMT="('rem5')")
          CALL prt_ctl_trc_info(charout)
          CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
-       ENDIF
-
-       !     Update the arrays TRA which contain the biological sources and sinks
-       !     --------------------------------------------------------------------
+      ENDIF
+
+      !     Update the arrays TRA which contain the biological sources and sinks
+      !     --------------------------------------------------------------------
 
       DO jk = 1, jpkm1
@@ -385 +386 @@
          tra(:,:,jk,jpdic) = tra(:,:,jk,jpdic) + zolimi(:,:,jk) + denitr(:,:,jk)
          tra(:,:,jk,jptal) = tra(:,:,jk,jptal) + denitr(:,:,jk) * rno3 * rdenit
-     END DO
-
-       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
+      END DO
+
+      IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
          WRITE(charout, FMT="('rem6')")
          CALL prt_ctl_trc_info(charout)
          CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
-       ENDIF
-
+      ENDIF
+      !
+      IF(  wrk_not_released(2, 1)     .OR.   &
+           wrk_not_released(3, 2,3,4)  )   CALL ctl_stop('p4z_rem: failed to release workspace arrays')
+      !
    END SUBROUTINE p4z_rem
 
+
    SUBROUTINE p4z_rem_init
-
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE p4z_rem_init  ***
@@ -408 +412 @@
       !!
       !!----------------------------------------------------------------------
-
       NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xlam1, oxymin
+      !!----------------------------------------------------------------------
 
       REWIND( numnat )                     ! read numnat
@@ -425 +429 @@
          WRITE(numout,*) '    halk saturation constant for anoxia       oxymin    =', oxymin
       ENDIF
-
-      nitrfac(:,:,:) = 0.0
-      denitr (:,:,:) = 0.0  
-
+      !
+      nitrfac(:,:,:) = 0._wp
+      denitr (:,:,:) = 0._wp
+      !
    END SUBROUTINE p4z_rem_init
+
+
+   INTEGER FUNCTION p4z_rem_alloc()
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_rem_alloc  ***
+      !!----------------------------------------------------------------------
+      ALLOCATE( denitr(jpi,jpj,jpk), STAT=p4z_rem_alloc )
+      !
+      IF( p4z_rem_alloc /= 0 )   CALL ctl_warn('p4z_rem_alloc: failed to allocate arrays')
+      !
+   END FUNCTION p4z_rem_alloc
 
 #else