- 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/LIM_SRC_3/limadv_umx.F90
r7753 r7910 21 21 USE lbclnk ! lateral boundary conditions -- MPP exchanges 22 22 USE lib_mpp ! MPP library 23 USE wrk_nemo ! work arrays24 23 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 25 24 USE timing ! Timing … … 65 64 REAL(wp) :: zfp_ui, zfp_vj ! - - 66 65 REAL(wp) :: zfm_ui, zfm_vj ! - - 67 REAL(wp), POINTER, DIMENSION(:,:) :: zt_ups, zfu_ups, zfv_ups, ztrd, zfu_ho, zfv_ho, zt_u, zt_v66 REAL(wp), DIMENSION(jpi,jpj) :: zt_ups, zfu_ups, zfv_ups, ztrd, zfu_ho, zfv_ho, zt_u, zt_v 68 67 !!---------------------------------------------------------------------- 69 68 ! 70 69 IF( nn_timing == 1 ) CALL timing_start('lim_adv_umx') 71 70 ! 72 CALL wrk_alloc( jpi,jpj, zt_ups, zfu_ups, zfv_ups, ztrd, zfu_ho, zfv_ho, zt_u, zt_v )73 71 ! 74 72 ! … … 146 144 ! 147 145 ! 148 CALL wrk_dealloc( jpi,jpj, zt_ups, zfu_ups, zfv_ups, ztrd, zfu_ho, zfv_ho, zt_u, zt_v )149 146 ! 150 147 IF( nn_timing == 1 ) CALL timing_stop('lim_adv_umx') … … 174 171 INTEGER :: ji, jj ! dummy loop indices 175 172 REAL(wp) :: zc_box ! - - 176 REAL(wp), POINTER, DIMENSION(:,:) :: zzt173 REAL(wp), DIMENSION(jpi,jpj) :: zzt 177 174 !!---------------------------------------------------------------------- 178 175 ! 179 176 IF( nn_timing == 1 ) CALL timing_start('macho') 180 177 ! 181 CALL wrk_alloc( jpi,jpj, zzt )182 178 ! 183 179 IF( MOD( (kt - 1) / nn_fsbc , 2 ) == 0 ) THEN !== odd ice time step: adv_x then adv_y ==! … … 219 215 ENDIF 220 216 ! 221 CALL wrk_dealloc( jpi,jpj, zzt )222 217 ! 223 218 IF( nn_timing == 1 ) CALL timing_stop('macho') … … 245 240 INTEGER :: ji, jj ! dummy loop indices 246 241 REAL(wp) :: zcu, zdx2, zdx4 ! - - 247 REAL(wp), POINTER, DIMENSION(:,:) :: ztu1, ztu2, ztu3, ztu4242 REAL(wp), DIMENSION(jpi,jpj) :: ztu1, ztu2, ztu3, ztu4 248 243 !!---------------------------------------------------------------------- 249 244 ! 250 245 IF( nn_timing == 1 ) CALL timing_start('ultimate_x') 251 246 ! 252 CALL wrk_alloc( jpi,jpj, ztu1, ztu2, ztu3, ztu4 )253 247 ! 254 248 ! !-- Laplacian in i-direction --! … … 346 340 END SELECT 347 341 ! 348 CALL wrk_dealloc( jpi,jpj, ztu1, ztu2, ztu3, ztu4 )349 342 ! 350 343 IF( nn_timing == 1 ) CALL timing_stop('ultimate_x') … … 372 365 INTEGER :: ji, jj ! dummy loop indices 373 366 REAL(wp) :: zcv, zdy2, zdy4 ! - - 374 REAL(wp), POINTER, DIMENSION(:,:) :: ztv1, ztv2, ztv3, ztv4367 REAL(wp), DIMENSION(jpi,jpj) :: ztv1, ztv2, ztv3, ztv4 375 368 !!---------------------------------------------------------------------- 376 369 ! 377 370 IF( nn_timing == 1 ) CALL timing_start('ultimate_y') 378 371 ! 379 CALL wrk_alloc( jpi,jpj, ztv1, ztv2, ztv3, ztv4 )380 372 ! 381 373 ! !-- Laplacian in j-direction --! … … 474 466 END SELECT 475 467 ! 476 CALL wrk_dealloc( jpi,jpj, ztv1, ztv2, ztv3, ztv4 )477 468 ! 478 469 IF( nn_timing == 1 ) CALL timing_stop('ultimate_y') … … 502 493 REAL(wp) :: zpos, zneg, zbt, za, zb, zc, zbig, zsml, z1_dt ! local scalars 503 494 REAL(wp) :: zau, zbu, zcu, zav, zbv, zcv, zup, zdo ! - - 504 REAL(wp), POINTER, DIMENSION(:,:) :: zbetup, zbetdo, zbup, zbdo, zmsk, zdiv495 REAL(wp), DIMENSION(jpi,jpj) :: zbetup, zbetdo, zbup, zbdo, zmsk, zdiv 505 496 !!---------------------------------------------------------------------- 506 497 ! 507 498 IF( nn_timing == 1 ) CALL timing_start('nonosc_2d') 508 499 ! 509 CALL wrk_alloc( jpi,jpj, zbetup, zbetdo, zbup, zbdo, zmsk, zdiv )510 500 ! 511 501 zbig = 1.e+40_wp … … 578 568 CALL lbc_lnk_multi( paa, 'U', -1., pbb, 'V', -1. ) ! lateral boundary condition (changed sign) 579 569 ! 580 CALL wrk_dealloc( jpi,jpj, zbetup, zbetdo, zbup, zbdo, zmsk, zdiv )581 570 ! 582 571 IF( nn_timing == 1 ) CALL timing_stop('nonosc_2d')
Note: See TracChangeset
for help on using the changeset viewer.