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

trunk/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90

@@ -27 +27 @@
       !!---------------------------------------------
 #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_3d_1, wrk_3d_2
+      USE wrk_nemo, ONLY: wrk_3d_3, wrk_3d_4
+      USE wrk_nemo, ONLY: wrk_3d_7, wrk_3d_6
+      USE wrk_nemo, ONLY: wrk_3d_8
+      !!
       INTEGER :: ji,jj,jk
       INTEGER :: spongearea
       REAL(wp) :: timecoeff
       REAL(wp) :: zta, zsa, zabe1, zabe2, zbtr
-      REAL(wp), DIMENSION(jpi,jpj) :: localviscsponge
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: tbdiff, sbdiff
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztu ,ztv, zsu ,zsv
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  ztab
+      REAL(wp), POINTER, DIMENSION(:,:) :: localviscsponge
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: tbdiff, sbdiff
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztu, zsu, ztv, zsv
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab
 
 #if defined SPONGE
+      localviscsponge => wrk_2d_1
+      tbdiff => wrk_3d_1 ;sbdiff => wrk_3d_2
+      ztu => wrk_3d_3 ; zsu => wrk_3d_4
+      ztv => wrk_3d_7 ; zsv => wrk_3d_6
+      ztab => wrk_3d_8
 
       timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot()
@@ -44 +55 @@
       Agrif_UseSpecialValue = .TRUE.
       ztab = 0.e0
-      CALL Agrif_Bc_Variable(ztab, ta,calledweight=timecoeff,procname=interptn)
+      CALL Agrif_Bc_Variable(ztab, ta_id,calledweight=timecoeff,procname=interptn)
       Agrif_UseSpecialValue = .FALSE.
 
@@ -52 +63 @@
       Agrif_SpecialValue=0.
       Agrif_UseSpecialValue = .TRUE.
-      CALL Agrif_Bc_Variable(ztab, sa,calledweight=timecoeff,procname=interpsn)
+      CALL Agrif_Bc_Variable(ztab, sa_id,calledweight=timecoeff,procname=interpsn)
       Agrif_UseSpecialValue = .FALSE.
 
       sbdiff(:,:,:) = sb(:,:,:) - ztab(:,:,:)
-
 
       spongearea = 2 + 2 * Agrif_irhox()
@@ -164 +174 @@
       !!---------------------------------------------
 #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_3d_1, wrk_3d_2
+      USE wrk_nemo, ONLY: wrk_3d_3, wrk_3d_4
+      USE wrk_nemo, ONLY: wrk_3d_5
+      !!
       INTEGER :: ji,jj,jk
       INTEGER :: spongearea
       REAL(wp) :: timecoeff
       REAL(wp) :: ze2u, ze1v, zua, zva, zbtr
-      REAL(wp), DIMENSION(jpi,jpj) :: localviscsponge
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztab, ubdiff, vbdiff,rotdiff,hdivdiff
+      REAL(wp), POINTER, DIMENSION(:,:) :: localviscsponge
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ubdiff, vbdiff
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: rotdiff, hdivdiff
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab
 
 #if defined SPONGE
+      localviscsponge => wrk_2d_1
+      ubdiff => wrk_3d_1 ; vbdiff => wrk_3d_2
+      rotdiff => wrk_3d_3 ; hdivdiff => wrk_3d_4
+      ztab => wrk_3d_5
 
       timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot()
@@ -179 +200 @@
       Agrif_UseSpecialValue = ln_spc_dyn
       ztab = 0.e0
-      CALL Agrif_Bc_Variable(ztab, ua,calledweight=timecoeff,procname=interpun)
+      CALL Agrif_Bc_Variable(ztab, ua_id,calledweight=timecoeff,procname=interpun)
       Agrif_UseSpecialValue = .FALSE.
 
@@ -187 +208 @@
       Agrif_SpecialValue=0.
       Agrif_UseSpecialValue = ln_spc_dyn
-      CALL Agrif_Bc_Variable(ztab, va,calledweight=timecoeff,procname=interpvn)
+      CALL Agrif_Bc_Variable(ztab, va_id,calledweight=timecoeff,procname=interpvn)
       Agrif_UseSpecialValue = .FALSE.
 
@@ -250 +271 @@
          spongedoneU = .TRUE.
     
-    spbtr3(:,:) = 1./( e1f(:,:) * e2f(:,:))
+     spbtr3(:,:) = 1./( e1f(:,:) * e2f(:,:))
       ENDIF