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 7571 for branches/UKMO/dev_r5518_MEDUSA_optim_MG_MS_RH/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90 – NEMO

Ignore:
Timestamp:
2017-01-17T11:33:49+01:00 (7 years ago)
Author:
frrh
Message:

Commit first phase of changes - these work OK in a UKESM vn0.5
suite employing the GO6 package branch at:

branches/UKMO/dev_r5518_GO6_package@7206

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_MEDUSA_optim_MG_MS_RH/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r7570 r7571  
    8282      REAL(wp) ::   zv, z0v, zzwy, z0w        !   -      - 
    8383      REAL(wp) ::   ztra, zbtr, zdt, zalpha   !   -      - 
    84       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zslpx, zslpy   ! 3D workspace 
    85       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zwx  , zwy     ! -      -  
     84      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zslpx, zslpy   ! 3D workspace 
     85      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zwx  , zwy     ! -      -  
     86 
    8687      !!---------------------------------------------------------------------- 
    8788      ! 
    8889      IF( nn_timing == 1 )  CALL timing_start('tra_adv_muscl') 
    8990      ! 
    90       CALL wrk_alloc( jpi, jpj, jpk, zslpx, zslpy, zwx, zwy ) 
     91      ALLOCATE( zslpx(1:jpi, 1:jpj, 1:jpk) ) 
     92      ALLOCATE( zslpy(1:jpi, 1:jpj, 1:jpk) ) 
     93      ALLOCATE( zwx  (1:jpi, 1:jpj, 1:jpk) ) 
     94      ALLOCATE( zwy  (1:jpi, 1:jpj, 1:jpk) ) 
    9195      ! 
    9296      IF( kt == kit000 )  THEN 
     
    291295      END DO 
    292296      ! 
    293       CALL wrk_dealloc( jpi, jpj, jpk, zslpx, zslpy, zwx, zwy ) 
     297      DEALLOCATE( zslpx ) 
     298      DEALLOCATE( zslpy ) 
     299      DEALLOCATE( zwx   ) 
     300      DEALLOCATE( zwy   ) 
    294301      ! 
    295302      IF( nn_timing == 1 )  CALL timing_stop('tra_adv_muscl') 
Note: See TracChangeset for help on using the changeset viewer.