Changeset 3168 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

Timestamp:

2011-11-22T10:48:38+01:00 (13 years ago)

Author:

cbricaud

Message:

add timing, change dynamical allocation and correct small bugs

File:

• branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90 (modified) (7 diffs)

branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

@@ -19 +19 @@
    USE lib_mpp        ! distribued memory computing library
    USE iom            ! I/O manager library
+   USE timing         ! preformance summary
+   USE wrk_nemo_2     ! working arrays
 
    IMPLICIT NONE
@@ -65 +67 @@
       !! ** Purpose :   compute and output some AR5 diagnostics
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   zarea_ssh => wrk_2d_1 , zbotpres => wrk_2d_2   ! 2D workspace
-      USE wrk_nemo, ONLY:   zrhd      => wrk_3d_1 , zrhop    => wrk_3d_2   ! 3D      -
-      USE wrk_nemo, ONLY:   ztsn      => wrk_4d_1                          ! 4D      -
       !
       INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
@@ -74 +72 @@
       INTEGER  ::   ji, jj, jk                      ! dummy loop arguments
       REAL(wp) ::   zvolssh, zvol, zssh_steric, zztmp, zarho, ztemp, zsal, zmass
+      !
+      REAL(wp), POINTER, DIMENSION(:,:)     :: zarea_ssh , zbotpres       ! 2D workspace 
+      REAL(wp), POINTER, DIMENSION(:,:,:)   :: zrhd , zrhop               ! 3D workspace
+      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztsn                       ! 4D workspace
       !!--------------------------------------------------------------------
-
-      IF( wrk_in_use(2, 1,2) .OR.   &
-          wrk_in_use(3, 1,2) .OR.   &
-          wrk_in_use(4, 1)   ) THEN
-         CALL ctl_stop('dia_ar5: requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
+      IF( nn_timing == 1 )   CALL timing_start('dia_ar5')
+ 
+      CALL wrk_alloc( jpi , jpj              , zarea_ssh , zbotpres )
+      CALL wrk_alloc( jpi , jpj , jpk        , zrhd      , zrhop    )
+      CALL wrk_alloc( jpi , jpj , jpk , jpts , ztsn                 )
 
       CALL iom_put( 'cellthc', fse3t(:,:,:) )
@@ -160 +161 @@
       CALL iom_put( 'saltot' , zsal  )
       !
-      IF( wrk_not_released(2, 1,2) .OR.   &
-          wrk_not_released(3, 1,2) .OR.   &
-          wrk_not_released(4, 1)   )   CALL ctl_stop('dia_ar5: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi , jpj              , zarea_ssh , zbotpres )
+      CALL wrk_dealloc( jpi , jpj , jpk        , zrhd      , zrhop    )
+      CALL wrk_dealloc( jpi , jpj , jpk , jpts , ztsn                 )
+      !
+      IF( nn_timing == 1 )   CALL timing_stop('dia_ar5')
       !
    END SUBROUTINE dia_ar5
@@ -173 +176 @@
       !! ** Purpose :   initialization for AR5 diagnostic computation
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   wrk_4d_1      ! 4D workspace
-      !
       INTEGER  ::   inum
       INTEGER  ::   ik
@@ -183 +183 @@
       !!----------------------------------------------------------------------
       !
-      IF(wrk_in_use(4, 1) ) THEN
-         CALL ctl_stop('dia_ar5_init: requested workspace array unavailable.')   ;   RETURN
-      ENDIF
-      zsaldta => wrk_4d_1(:,:,:,1:2)
-
+      IF( nn_timing == 1 )   CALL timing_start('dia_ar5_init')
+      !
+      CALL wrk_alloc( jpi , jpj , jpk, jpts, zsaldta )
       !                                      ! allocate dia_ar5 arrays
       IF( dia_ar5_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dia_ar5_init : unable to allocate arrays' )
@@ -221 +219 @@
       ENDIF
       !
-      IF( wrk_not_released(4, 1) )   CALL ctl_stop('dia_ar5_init: failed to release workspace array')
+      CALL wrk_dealloc( jpi , jpj , jpk, jpts, zsaldta )
+      !
+      IF( nn_timing == 1 )   CALL timing_stop('dia_ar5_init')
       !
    END SUBROUTINE dia_ar5_init