Changeset 2677 for branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

Timestamp:

2011-03-09T15:39:40+01:00 (13 years ago)

Author:

rblod

Message:

Commit in NST_SRC for agrif and dynamic memory

File:

• branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90 (modified) (4 diffs)

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

@@ -9 +9 @@
    USE agrif_oce
    USE trc
+   USE lib_mpp
 
    IMPLICIT NONE
@@ -28 +29 @@
       !!---------------------------------------------
 #include "domzgr_substitute.h90"
-
+      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY: wrk_2d_1
+      USE wrk_nemo, ONLY: wrk_4d_1, wrk_4d_2, wrk_4d_3, wrk_4d_4
+      !! 
       INTEGER :: ji,jj,jk,jl
       INTEGER :: spongearea
       REAL(wp) :: timecoeff
       REAL(wp) :: ztra, zabe1, zabe2, zbtr
-      REAL(wp), DIMENSION(jpi,jpj) :: localviscsponge
-      REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) :: trbdiff
-      REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) :: ztru ,ztrv
-      REAL(wp), DIMENSION(jpi,jpj,jpk,jptra) ::  ztab
+      REAL(wp), POINTER, DIMENSION(:,:) :: localviscsponge
+      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: trbdiff, ztru, ztrv, ztab
 
 #if defined SPONGE_TOP
+      IF( wrk_in_use(4, 1,2,3,4) .OR. wrk_in_use(2, 1) ) THEN
+         CALL ctl_stop('Agrif_Sponge_trc : requested workspace arrays unavailable')
+         RETURN
+      ENDIF
+      localviscsponge => wrk_2d_1
+      trbdiff(:,:,:,:) => wrk_4d_1(:,:,:,1:jptra)
+      ztru   (:,:,:,:) => wrk_4d_2(:,:,:,1:jptra)
+      ztrv   (:,:,:,:) => wrk_4d_3(:,:,:,1:jptra)
+      ztab   (:,:,:,:) => wrk_4d_4(:,:,:,1:jptra)
 
       timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot()
@@ -45 +56 @@
       Agrif_UseSpecialValue = .TRUE.
       ztab = 0.e0
-      CALL Agrif_Bc_Variable(ztab, tra,calledweight=timecoeff,procname=interptrn)
+      CALL Agrif_Bc_Variable(ztab, tra_id,calledweight=timecoeff,procname=interptrn)
       Agrif_UseSpecialValue = .FALSE.
 
@@ -143 +154 @@
       ENDDO
       ENDDO
+ 
+      IF( wrk_not_released(4, 1,2,3,4) .OR. wrk_not_released(2, 1) ) THEN
+         CALL ctl_stop('Agrif_Sponge_trc : failed to release workspace arrays.')
+         RETURN
+      ENDIF
 
 #endif