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 8568 for branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90 – NEMO

Ignore:
Timestamp:
2017-09-27T16:29:24+02:00 (7 years ago)
Author:
gm
Message:

#1911 (ENHANCE-09): PART I.2 - _NONE option + remove zts + see associated wiki page

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90

    r7753 r8568  
    1919   USE in_out_manager ! I/O manager 
    2020   USE lbclnk         ! ocean lateral boundary conditions (or mpp link) 
    21    USE wrk_nemo       ! Memory Allocation 
    2221   USE timing         ! Timing 
    2322 
     
    3130#  include "vectopt_loop_substitute.h90" 
    3231   !!---------------------------------------------------------------------- 
    33    !! NEMO/OPA 3.7 , NEMO Consortium (2014) 
     32   !! NEMO/OPA 4.0 , NEMO Consortium (2017) 
    3433   !! $Id$  
    3534   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     
    5756      REAL(wp) ::   zsign        ! local scalars 
    5857      REAL(wp) ::   zua, zva     ! local scalars 
    59       REAL(wp), POINTER, DIMENSION(:,:) ::  zcur, zdiv 
     58      REAL(wp), DIMENSION(jpi,jpj) ::   zcur, zdiv 
    6059      !!---------------------------------------------------------------------- 
    6160      ! 
     
    6665      ENDIF 
    6766      ! 
    68       IF( nn_timing == 1 )   CALL timing_start('dyn_ldf_lap') 
    69       ! 
    70       CALL wrk_alloc( jpi, jpj, zcur, zdiv )  
     67      IF( ln_timing )   CALL timing_start('dyn_ldf_lap') 
    7168      ! 
    7269      IF( kpass == 1 ) THEN   ;   zsign =  1._wp      ! bilaplacian operator require a minus sign 
     
    107104      END DO                                           !   End of slab 
    108105      !                                                ! =============== 
    109       CALL wrk_dealloc( jpi, jpj, zcur, zdiv )  
    110106      ! 
    111       IF( nn_timing == 1 )  CALL timing_stop('dyn_ldf_lap') 
     107      IF( ln_timing )   CALL timing_stop('dyn_ldf_lap') 
    112108      ! 
    113109   END SUBROUTINE dyn_ldf_lap 
     
    131127      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pua, pva   ! momentum trend 
    132128      ! 
    133       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zulap, zvlap   ! laplacian at u- and v-point 
     129      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zulap, zvlap   ! laplacian at u- and v-point 
    134130      !!---------------------------------------------------------------------- 
    135131      ! 
    136       IF( nn_timing == 1 )  CALL timing_start('dyn_ldf_blp') 
    137       ! 
    138       CALL wrk_alloc( jpi, jpj, jpk, zulap, zvlap )  
     132      IF( ln_timing )   CALL timing_start('dyn_ldf_blp') 
    139133      ! 
    140134      IF( kt == nit000 )  THEN 
     
    154148      CALL dyn_ldf_lap( kt, zulap, zvlap, pua, pva, 2 )   ! rotated laplacian applied to zlap (output in pta) 
    155149      ! 
    156       CALL wrk_dealloc( jpi, jpj, jpk, zulap, zvlap )  
    157       ! 
    158       IF( nn_timing == 1 )  CALL timing_stop('dyn_ldf_blp') 
     150      IF( ln_timing )   CALL timing_stop('dyn_ldf_blp') 
    159151      ! 
    160152   END SUBROUTINE dyn_ldf_blp 
Note: See TracChangeset for help on using the changeset viewer.