Changeset 8568 for branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90
- Timestamp:
- 2017-09-27T16:29:24+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90
r7753 r8568 19 19 USE in_out_manager ! I/O manager 20 20 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 21 USE wrk_nemo ! Memory Allocation22 21 USE timing ! Timing 23 22 … … 31 30 # include "vectopt_loop_substitute.h90" 32 31 !!---------------------------------------------------------------------- 33 !! NEMO/OPA 3.7 , NEMO Consortium (2014)32 !! NEMO/OPA 4.0 , NEMO Consortium (2017) 34 33 !! $Id$ 35 34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 57 56 REAL(wp) :: zsign ! local scalars 58 57 REAL(wp) :: zua, zva ! local scalars 59 REAL(wp), POINTER, DIMENSION(:,:) ::zcur, zdiv58 REAL(wp), DIMENSION(jpi,jpj) :: zcur, zdiv 60 59 !!---------------------------------------------------------------------- 61 60 ! … … 66 65 ENDIF 67 66 ! 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') 71 68 ! 72 69 IF( kpass == 1 ) THEN ; zsign = 1._wp ! bilaplacian operator require a minus sign … … 107 104 END DO ! End of slab 108 105 ! ! =============== 109 CALL wrk_dealloc( jpi, jpj, zcur, zdiv )110 106 ! 111 IF( nn_timing == 1 )CALL timing_stop('dyn_ldf_lap')107 IF( ln_timing ) CALL timing_stop('dyn_ldf_lap') 112 108 ! 113 109 END SUBROUTINE dyn_ldf_lap … … 131 127 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pua, pva ! momentum trend 132 128 ! 133 REAL(wp), POINTER, DIMENSION(:,:,:) :: zulap, zvlap ! laplacian at u- and v-point129 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zulap, zvlap ! laplacian at u- and v-point 134 130 !!---------------------------------------------------------------------- 135 131 ! 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') 139 133 ! 140 134 IF( kt == nit000 ) THEN … … 154 148 CALL dyn_ldf_lap( kt, zulap, zvlap, pua, pva, 2 ) ! rotated laplacian applied to zlap (output in pta) 155 149 ! 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') 159 151 ! 160 152 END SUBROUTINE dyn_ldf_blp
Note: See TracChangeset
for help on using the changeset viewer.