New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7545 for branches/UKMO/MEDUSA_optim_MG_MS_RH/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90 – NEMO

Ignore:
Timestamp:
2017-01-11T12:27:34+01:00 (7 years ago)
Author:
frrh
Message:

Apply Maff's changes with the exception of
NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
which will clash with changes made to the source code since
the optimisation work started.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/MEDUSA_optim_MG_MS_RH/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r5733 r7545  
    7878      ! 
    7979      INTEGER ::   jk   ! dummy loop index 
    80       REAL(wp), POINTER, DIMENSION(:,:,:) :: zun, zvn, zwn 
     80      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zun, zvn, zwn 
     81 
    8182      !!---------------------------------------------------------------------- 
    8283      ! 
    8384      IF( nn_timing == 1 )  CALL timing_start('tra_adv') 
    8485      ! 
    85       CALL wrk_alloc( jpi, jpj, jpk, zun, zvn, zwn ) 
     86 
     87      ALLOCATE( zun(1:jpi, 1:jpj, 1:jpk) ) 
     88      ALLOCATE( zvn(1:jpi, 1:jpj, 1:jpk) ) 
     89      ALLOCATE( zwn(1:jpi, 1:jpj, 1:jpk) ) 
     90 
    8691      !                                          ! set time step 
    8792      IF( neuler == 0 .AND. kt == nit000 ) THEN     ! at nit000 
     
    108113      zvn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom 
    109114      zwn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom 
     115 
    110116      ! 
    111117      IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif )   & 
     
    155161         &                       tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    156162      ! 
     163      DEALLOCATE ( zun, zvn, zwn ) 
    157164      IF( nn_timing == 1 )  CALL timing_stop( 'tra_adv' ) 
    158165      ! 
    159       CALL wrk_dealloc( jpi, jpj, jpk, zun, zvn, zwn ) 
    160166      !                                           
    161167   END SUBROUTINE tra_adv 
Note: See TracChangeset for help on using the changeset viewer.