- Timestamp:
- 2017-04-13T16:21:08+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_no_wrk_alloc/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90
r7753 r7910 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 … … 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 ! … … 68 67 IF( nn_timing == 1 ) CALL timing_start('dyn_ldf_lap') 69 68 ! 70 CALL wrk_alloc( jpi, jpj, zcur, zdiv )71 69 ! 72 70 IF( kpass == 1 ) THEN ; zsign = 1._wp ! bilaplacian operator require a minus sign … … 107 105 END DO ! End of slab 108 106 ! ! =============== 109 CALL wrk_dealloc( jpi, jpj, zcur, zdiv )110 107 ! 111 108 IF( nn_timing == 1 ) CALL timing_stop('dyn_ldf_lap') … … 131 128 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pua, pva ! momentum trend 132 129 ! 133 REAL(wp), POINTER, DIMENSION(:,:,:) :: zulap, zvlap ! laplacian at u- and v-point130 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zulap, zvlap ! laplacian at u- and v-point 134 131 !!---------------------------------------------------------------------- 135 132 ! 136 133 IF( nn_timing == 1 ) CALL timing_start('dyn_ldf_blp') 137 134 ! 138 CALL wrk_alloc( jpi, jpj, jpk, zulap, zvlap )139 135 ! 140 136 IF( kt == nit000 ) THEN … … 154 150 CALL dyn_ldf_lap( kt, zulap, zvlap, pua, pva, 2 ) ! rotated laplacian applied to zlap (output in pta) 155 151 ! 156 CALL wrk_dealloc( jpi, jpj, jpk, zulap, zvlap )157 152 ! 158 153 IF( nn_timing == 1 ) CALL timing_stop('dyn_ldf_blp')
Note: See TracChangeset
for help on using the changeset viewer.