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

trunk/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90

@@ -90 +90 @@
       !! ** Action :   Update pta arrays with the before rotated diffusion
       !!----------------------------------------------------------------------
-      USE oce         , zftu => ua   ! use ua as workspace
-      USE oce         , zftv => va   ! use va as workspace
-      !!
+      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
+      USE oce     , ONLY:   zftu => ua       , zftv  => va         ! (ua,va) used as workspace
+      USE wrk_nemo, ONLY:   zdkt => wrk_2d_1 , zdk1t => wrk_2d_2 , z2d  => wrk_2d_3   ! 2D workspace
+      USE wrk_nemo, ONLY:   zdit => wrk_3d_6 , zdjt  => wrk_3d_7 , ztfw => wrk_3d_8   ! 3D workspace
+      !
       INTEGER                              , INTENT(in   ) ::   kt         ! ocean time-step index
       CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator)
@@ -100 +102 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) ::   pta        ! tracer trend 
       REAL(wp)                             , INTENT(in   ) ::   pahtb0     ! background diffusion coef
-      !!
+      !
       INTEGER  ::  ji, jj, jk, jn   ! dummy loop indices
       REAL(wp) ::  zmsku, zabe1, zcof1, zcoef3   ! local scalars
       REAL(wp) ::  zmskv, zabe2, zcof2, zcoef4   !   -      -
       REAL(wp) ::  zcoef0, zbtr, ztra            !   -      -
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zdkt, zdk1t         ! 2D workspace
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdit, zdjt, ztfw    ! 3D workspace
 #if defined key_diaar5
-      REAL(wp), DIMENSION(jpi,jpj)     ::   z2d                 ! 2D workspace
       REAL(wp)                         ::   zztmp               ! local scalar
 #endif
       !!----------------------------------------------------------------------
+
+      IF( wrk_in_use(3, 6,7,8) .OR. wrk_in_use(2, 1,2,3) ) THEN
+          CALL ctl_stop('tra_ldf_iso : requested workspace array unavailable')   ;   RETURN
+      ENDIF
 
       IF( kt == nit000 )  THEN
@@ -288 +291 @@
       END DO
       !
+      IF( wrk_not_released(3, 6,7,8) .OR.   &
+          wrk_not_released(2, 1,2,3) )   CALL ctl_stop('tra_ldf_iso: failed to release workspace arrays')
+      !
    END SUBROUTINE tra_ldf_iso