Changeset 7771 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SOL

Timestamp:

2017-03-09T11:38:31+01:00 (7 years ago)

Author:

frrh

Message:

Apply optimisations to various areas of code replacing the use of
allocated pointers with straightforward direct ALLOCATE and DEALLOCATE
operations.

These optimisations largely have an impact in models featuring MEDUSA,
i.e. those with significant numbers of tracers, although they are
expected to have a small impact in all configurations.

Code developed and tested in NEMO branch branches/UKMO/dev_r5518_optim_GO6_alloc
Tested in stand-alone GO6-GSI8, GO6-GSI8-MEDUSA and UKESM coupled models.
NEMO ticket #1821 documents this change further.

File:

• branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SOL/solpcg.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SOL/solpcg.F90

@@ -93 +93 @@
       REAL(wp) ::   zgcad        ! temporary scalars
       REAL(wp), DIMENSION(2) ::   zsum
-      REAL(wp), POINTER, DIMENSION(:,:) ::   zgcr
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zgcr
       !!----------------------------------------------------------------------
       !
       IF( nn_timing == 1 )  CALL timing_start('sol_pcg')
       !
-      CALL wrk_alloc( jpi, jpj, zgcr )
+      ALLOCATE( zgcr(1:jpi,1:jpj) )
       !
       ! Initialization of the algorithm with standard PCG
@@ -210 +210 @@
       CALL lbc_lnk( gcx, c_solver_pt, 1. )      ! Output in gcx with lateral b.c. applied
       ! 
-      CALL wrk_dealloc( jpi, jpj, zgcr )
+      DEALLOCATE ( zgcr )
       !
       IF( nn_timing == 1 )  CALL timing_stop('sol_pcg')