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

trunk/NEMOGCM/NEMO/OFF_SRC/domrea.F90

@@ -16 +16 @@
    USE dommsk          ! domain: masks
    USE lbclnk          ! lateral boundary condition - MPP exchanges
-   USE in_out_manager  ! I/O manager
+   USE lib_mpp 
+   USE in_out_manager
+   USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
 
    IMPLICIT NONE
@@ -53 +55 @@
       !!----------------------------------------------------------------------
       USE iom
+      USE wrk_nemo, ONLY: zmbk => wrk_2d_1, zprt => wrk_2d_2, zprw => wrk_2d_3
       !!
       INTEGER  ::   ji, jj, jk   ! dummy loop indices
       INTEGER  ::   ik, inum0 , inum1 , inum2 , inum3 , inum4   ! local integers
       REAL(wp) ::   zrefdep         ! local real
-      REAL(wp), DIMENSION(jpi,jpj) ::   zmbk, zprt, zprw   ! 2D workspace
       !!----------------------------------------------------------------------
 
@@ -63 +65 @@
       IF(lwp) WRITE(numout,*) 'dom_rea : read NetCDF mesh and mask information file(s)'
       IF(lwp) WRITE(numout,*) '~~~~~~~'
+
+      IF( wrk_in_use(2, 1,2,3)  ) THEN
+         CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.') ; RETURN
+      END IF
 
       zmbk(:,:) = 0._wp
@@ -141 +147 @@
          CALL iom_get( inum3, jpdom_data, 'e2u', e2u )
          CALL iom_get( inum3, jpdom_data, 'e2v', e2v )
+
+         e1e2t(:,:) = e1t(:,:) * e2t(:,:)                              ! surface at T-points
 
          CALL iom_get( inum3, jpdom_data, 'ff', ff )
@@ -314 +322 @@
       END SELECT
       !
+      IF( wrk_not_released(2, 1,2,3)  ) CALL ctl_stop('dom_rea:failed to release workspace arrays.')
+      !
    END SUBROUTINE dom_rea
 
@@ -327 +337 @@
       !! ** Action  : - update mbathy: level bathymetry (in level index)
       !!----------------------------------------------------------------------
+      USE wrk_nemo, ONLY: zmbk => wrk_2d_4
+      !
       INTEGER ::   ji, jj   ! dummy loop indices
-      REAL(wp), DIMENSION(jpi,jpj) ::   zmbk   ! 2D workspace 
-      !!----------------------------------------------------------------------
+      !!----------------------------------------------------------------------
+
       !
       IF(lwp) WRITE(numout,*)
       IF(lwp) WRITE(numout,*) '    zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels '
       IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~~~'
+      !
+      IF( wrk_in_use(2, 4) ) THEN
+         CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.')  ;  RETURN
+      END IF
       !
       mbkt(:,:) = MAX( mbathy(:,:) , 1 )    ! bottom k-index of T-level (=1 over land)
@@ -347 +363 @@
       zmbk(:,:) = REAL( mbkv(:,:), wp )   ;   CALL lbc_lnk(zmbk,'V',1.)   ;   mbkv  (:,:) = MAX( INT( zmbk(:,:) ), 1 )
       !
+      IF( wrk_not_released(2, 4) ) CALL ctl_stop('dom_rea:failed to release workspace arrays.')
+      !
    END SUBROUTINE zgr_bot_level