Changeset 11467 for NEMO/branches/UKMO/dev_r10037_GPU/src/OCE/DYN/dynnxt.F90

Timestamp:

2019-08-22T11:49:08+02:00 (5 years ago)

Author:

andmirek

Message:

Ticket #2197 allocate arrays at the beggining of the run

File:

• NEMO/branches/UKMO/dev_r10037_GPU/src/OCE/DYN/dynnxt.F90 (modified) (5 diffs)

NEMO/branches/UKMO/dev_r10037_GPU/src/OCE/DYN/dynnxt.F90

@@ -92 +92 @@
       !!               un,vn   now horizontal velocity of next time-step
       !!----------------------------------------------------------------------
+      USE scoce, ONLY : zue => scr2D1, zve => scr2D2, &
+                       ze3u_f => scr1, ze3v_f => scr2, zua => scr3, zva => scr4
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index
       !
@@ -98 +100 @@
       REAL(wp) ::   zue3a, zue3n, zue3b, zuf, zcoef    ! local scalars
       REAL(wp) ::   zve3a, zve3n, zve3b, zvf, z1_2dt   !   -      -
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zue, zve
-      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   ze3u_f, ze3v_f, zua, zva 
       !!----------------------------------------------------------------------
       !
       IF( ln_timing    )   CALL timing_start('dyn_nxt')
-      IF( ln_dynspg_ts )   ALLOCATE( zue(jpi,jpj)     , zve(jpi,jpj)     )
-      IF( l_trddyn     )   ALLOCATE( zua(jpi,jpj,jpk) , zva(jpi,jpj,jpk) )
       !
       IF( kt == nit000 ) THEN
@@ -275 +273 @@
             ELSE                          ! Asselin filter applied on thickness weighted velocity
                !
-               ALLOCATE( ze3u_f(jpi,jpj,jpk) , ze3v_f(jpi,jpj,jpk) )
                ! Before filtered scale factor at (u/v)-points stored in ze3u_f, ze3v_f
                CALL dom_vvl_interpol( e3t_b(:,:,:), ze3u_f, 'U' )
@@ -302 +299 @@
                e3v_b(:,:,1:jpkm1) = ze3v_f(:,:,1:jpkm1)
                !
-               DEALLOCATE( ze3u_f , ze3v_f )
             ENDIF
             !
@@ -368 +364 @@
          &                       tab3d_2=vn, clinfo2=' Vn: '       , mask2=vmask )
       ! 
-      IF( ln_dynspg_ts )   DEALLOCATE( zue, zve )
-      IF( l_trddyn     )   DEALLOCATE( zua, zva )
       IF( ln_timing    )   CALL timing_stop('dyn_nxt')
       !