Changeset 2648 for branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OFF_SRC/dommsk.F90

Timestamp:

2011-03-03T17:13:18+01:00 (13 years ago)

Author:

cetlod

Message:

Changed OFF_SRC component to use dynamic memory

File:

• branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OFF_SRC/dommsk.F90 (modified) (5 diffs)

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OFF_SRC/dommsk.F90

@@ -12 +12 @@
    USE oce             ! ocean dynamics and tracers
    USE dom_oce         ! ocean space and time domain
-   USE in_out_manager  ! I/O manager
+   USE lib_mpp
 
    IMPLICIT NONE
@@ -19 +19 @@
    PUBLIC   dom_msk    ! routine called by inidom.F90
 
-#if defined key_degrad
-   !! ------------------------------------------------
-   !! Degradation method
-   !! --------------------------------------------------
-   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jpk) ::   facvol  !! volume for degraded regions
-#endif
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   facvol  !! volume for degraded regions
 
    !! * Substitutions
@@ -34 +29 @@
    !!----------------------------------------------------------------------
 CONTAINS
+
    
    SUBROUTINE dom_msk
@@ -49 +45 @@
       !!               tpol     : ???
       !!----------------------------------------------------------------------
+      USE wrk_nemo, ONLY: iwrk_in_use, iwrk_not_released
+      USE wrk_nemo, ONLY: imsk => iwrk_2d_1
       INTEGER  ::   ji, jk                   ! dummy loop indices
       INTEGER  ::   iif, iil, ijf, ijl       ! local integers
-      INTEGER, DIMENSION(jpi,jpj) ::  imsk   ! 2D workspace
       !!---------------------------------------------------------------------
       !
+      IF( iwrk_in_use(2, 1) ) THEN
+         CALL ctl_stop('dom_msk: ERROR: requested workspace arrays unavailable.')  ;  RETURN
+      END IF
+      !
+      CALL dom_msk_alloc
+
       ! Interior domain mask (used for global sum)
       ! --------------------
@@ -95 +98 @@
       ENDIF
       !
+      IF( iwrk_not_released(2, 1) ) CALL ctl_stop('dom_msk: failed to release workspace arrays.')
+      !
    END SUBROUTINE dom_msk
+
+   SUBROUTINE dom_msk_alloc()
+      !!---------------------------------------------------------------------
+      !!                 ***  ROUTINE dom_msk_alloc  ***
+      !!---------------------------------------------------------------------
+#if defined key_degrad
+      INTEGER :: ierr 
+
+      ALLOCATE( facvol(jpi,jpj,jpk), STAT=ierr )
+      IF( ierr /= 0 )  &
+        &           CALL ctl_stop('STOP', 'dom_msk : unable to allocate facvol array')
+#endif
+
+   END SUBROUTINE dom_msk_alloc
 
    !!======================================================================