- 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/sshwzv.F90
r7753 r7910 36 36 USE lbclnk ! ocean lateral boundary condition (or mpp link) 37 37 USE lib_mpp ! MPP library 38 USE wrk_nemo ! Memory Allocation39 38 USE timing ! Timing 40 39 USE wet_dry ! Wetting/Drying flux limting … … 74 73 INTEGER :: jk ! dummy loop indice 75 74 REAL(wp) :: z2dt, zcoef ! local scalars 76 REAL(wp), POINTER, DIMENSION(:,:) :: zhdiv ! 2D workspace75 REAL(wp), DIMENSION(jpi,jpj) :: zhdiv ! 2D workspace 77 76 !!---------------------------------------------------------------------- 78 77 ! 79 78 IF( nn_timing == 1 ) CALL timing_start('ssh_nxt') 80 79 ! 81 CALL wrk_alloc( jpi,jpj, zhdiv )82 80 ! 83 81 IF( kt == nit000 ) THEN … … 134 132 IF(ln_ctl) CALL prt_ctl( tab2d_1=ssha, clinfo1=' ssha - : ', mask1=tmask, ovlap=1 ) 135 133 ! 136 CALL wrk_dealloc( jpi, jpj, zhdiv )137 134 ! 138 135 IF( nn_timing == 1 ) CALL timing_stop('ssh_nxt') … … 161 158 REAL(wp) :: z1_2dt ! local scalars 162 159 REAL(wp), POINTER, DIMENSION(:,: ) :: z2d 163 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3d, zhdiv160 REAL(wp), DIMENSION(jpi,jpj,jpk) :: z3d, zhdiv 164 161 !!---------------------------------------------------------------------- 165 162 ! … … 180 177 ! 181 178 IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN ! z_tilde and layer cases 182 CALL wrk_alloc( jpi, jpj, jpk, zhdiv )183 179 ! 184 180 DO jk = 1, jpkm1 … … 200 196 END DO 201 197 ! IF( ln_vvl_layer ) wn(:,:,:) = 0.e0 202 CALL wrk_dealloc( jpi, jpj, jpk, zhdiv )203 198 ELSE ! z_star and linear free surface cases 204 199 DO jk = jpkm1, 1, -1 ! integrate from the bottom the hor. divergence
Note: See TracChangeset
for help on using the changeset viewer.