Changeset 2609 for branches

Timestamp:

2011-02-22T18:14:30+01:00 (13 years ago)

Author:

trackstand2

Message:

Replaced local workspace arrays with those in wrk_nemo module

Location:

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER

Files:

• trcbio.F90 (modified) (3 diffs)
• trcini_lobster.F90 (modified) (2 diffs)
• trcopt.F90 (modified) (3 diffs)
• trcsed.F90 (modified) (2 diffs)
• trcsms_lobster.F90 (modified) (2 diffs)

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcbio.F90

@@ -60 +60 @@
       !!              for passive tracers are saved for futher diagnostics.
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_use,  wrk_release
+      USE wrk_nemo, ONLY: wrk_3d_2, wrk_4d_1
+      !!
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
       !!
@@ -75 +78 @@
 #endif
 #if defined key_diatrc && defined key_iomput
-      REAL(wp), DIMENSION(jpi,jpj,17)     :: zw2d
-      REAL(wp), DIMENSION(jpi,jpj,jpk,3) :: zw3d
+      REAL(wp), POINTER,   DIMENSION(:,:,:) :: zw2d
+      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: zw3d
 #endif
       REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::   ztrbio
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
+
+#if defined key_diatrc && defined key_iomput
+      IF( (.NOT. wrk_use(3, 2)) .OR. (.NOT. wrk_use(4, 1)) )THEN
+         CALL ctl_stop('trc_bio : requested workspace arrays unavailable.')
+         RETURN
+      END IF
+      ! Set-up pointers into sub-arrays of workspaces
+      zw2d => wrk_3d_2(:,:,1:17)
+      zw3d => wrk_4d_1(:,:,:,1:3)
+#endif
 
       IF( kt == nit000 ) THEN
@@ -485 +498 @@
          CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
       ENDIF
-
+      !
+#if defined key_diatrc && defined key_iomput
+      IF( (.NOT. wrk_release(3, 2)) .OR. (.NOT. wrk_release(4, 1)) )THEN
+         CALL ctl_stop('trc_bio : failed to release workspace arrays.')
+      END IF
+#endif
+      !
    END SUBROUTINE trc_bio
 

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcini_lobster.F90

@@ -43 +43 @@
       !! ** purpose :   specific initialisation for LOBSTER bio-model
       !!----------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_use, wrk_release
+      USE wrk_nemo, ONLY: zrro => wrk_2d_1, zdm0 => wrk_3d_1
+      !!
       INTEGER  ::   ji, jj, jk, jn
       REAL(wp) ::   ztest, zfluo, zfluu
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zrro
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdm0
-      !!----------------------------------------------------------------------
+      !!----------------------------------------------------------------------
+
+      IF( (.NOT. wrk_use(2, 1)) .OR. (.NOT. wrk_use(3, 1)) )THEN
+         CALL ctl_stop('trc_ini_lobster : requested workspace arrays unavailable.')
+         RETURN
+      END IF
 
       !  Control consitency
@@ -253 +259 @@
       IF(lwp) WRITE(numout,*) ' '
 
+      IF( (.NOT. wrk_release(2, 1)) .OR. (.NOT. wrk_release(3, 1)) )THEN
+         CALL ctl_stop('trc_ini_lobster : failed to release workspace arrays.')
+      END IF
 
    END SUBROUTINE trc_ini_lobster

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcopt.F90

@@ -52 +52 @@
       !!                xze    ???
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_use, wrk_release
+      USE wrk_nemo, ONLY: zpar100 => wrk_2d_1, & ! irradiance at euphotic layer depth
+                          zpar0m  => wrk_2d_2    ! irradiance just below the surface
+      USE wrk_nemo, ONLY: zparr => wrk_3d_2, &   ! red and green compound of par
+                          zparg => wrk_3d_3
+      !!
       INTEGER, INTENT( in ) ::   kt   ! index of the time stepping
       !!
@@ -59 +65 @@
       REAL(wp) ::   zkr, zkg            ! total absorption coefficient in red and green
       REAL(wp) ::   zcoef               ! temporary scalar
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zpar100         ! irradiance at euphotic layer depth
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zpar0m          ! irradiance just below the surface
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zparr, zparg    ! red and green compound of par
 
       !!---------------------------------------------------------------------
+
+      IF( (.NOT. wrk_use(2, 1,2)) .OR. (.NOT. wrk_use(3, 2,3)) )THEN
+         CALL ctl_stop('trc_opt : requested workspace arrays unavailable.')
+         RETURN
+      END IF
 
       IF( kt == nit000 ) THEN
@@ -130 +138 @@
       ENDIF
       !
+      IF( (.NOT. wrk_release(2, 1,2)) .OR. (.NOT. wrk_release(3, 2,3)) )THEN
+         CALL ctl_stop('trc_opt : failed to release workspace arrays.')
+      END IF
+      !
    END SUBROUTINE trc_opt
 

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsed.F90

@@ -56 +56 @@
       !!              trend of passive tracers is saved for futher diagnostics.
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_use, wrk_release
+      USE wrk_nemo, ONLY: zwork => wrk_3d_2
+      USE wrk_nemo, ONLY: zw2d  => wrk_2d_1 ! only used (if defined 
+                                            ! key_diatrc && defined key_iomput)
+      !!
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
       !!
       INTEGER  ::   ji, jj, jk, jl
       REAL(wp) ::   ztra
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwork
-#if defined key_diatrc && defined key_iomput
-      REAL(wp), DIMENSION(jpi,jpj) ::  zw2d
-#endif
       REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrbio
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
+
+      IF( (.NOT. wrk_use(3,2)) .OR. (.NOT. wrk_use(2,1)) )THEN
+         CALL ctl_stop('trc_sed : requested workspace arrays unavailable.')
+         RETURN
+      END IF
 
       IF( kt == nit000 ) THEN
@@ -144 +150 @@
       ENDIF
 
+      IF( (.NOT. wrk_release(3,2)) .OR. (.NOT. wrk_release(2,1)) )THEN
+         CALL ctl_stop('trc_sed : failed to release workspace arrays.')
+      END IF
+
    END SUBROUTINE trc_sed
 

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsms_lobster.F90

@@ -45 +45 @@
       !!
       !! ** Method  : - ???
-      !! -------------------------------------------------------------------------------------
+      !! --------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_use, wrk_release
+      USE wrk_nemo, ONLY: ztrlob => wrk_3d_1   ! used for lobster sms trends
+      !!
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
       INTEGER :: jn
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   ztrlob   ! used for lobster sms trends
-      !!
+      !! --------------------------------------------------------------------
+
+      IF(.NOT. wrk_use(3,1))THEN
+         CALL ctl_stop('trc_sms_lobster : requested workspace array unavailable.')
+         RETURN
+      END IF
 
       CALL trc_opt( kt )      ! optical model
@@ -64 +71 @@
 
       IF( lk_trdmld_trc )  CALL trd_mld_bio( kt )   ! trends: Mixed-layer
+
+      IF(.NOT. wrk_release(3,1))THEN
+         CALL ctl_stop('trc_sms_lobster : failed to release workspace array.')
+      END IF
 
    END SUBROUTINE trc_sms_lobster