- Timestamp:
- 2011-11-17T17:38:56+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OFF_SRC/domrea.F90
r2977 r3149 19 19 USE lib_mpp 20 20 USE in_out_manager 21 USE wrk_nemo , ONLY: wrk_in_use, wrk_not_released21 USE wrk_nemo_2 22 22 23 23 IMPLICIT NONE … … 56 56 !!---------------------------------------------------------------------- 57 57 USE iom 58 USE wrk_nemo, ONLY: zmbk => wrk_2d_1, zprt => wrk_2d_2, zprw => wrk_2d_359 58 !! 60 59 INTEGER :: ji, jj, jk ! dummy loop indices 61 60 INTEGER :: ik, inum0 , inum1 , inum2 , inum3 , inum4 ! local integers 62 61 REAL(wp) :: zrefdep ! local real 62 REAL(wp), POINTER, DIMENSION(:,:) :: zmbk, zprt, zprw 63 63 !!---------------------------------------------------------------------- 64 64 … … 67 67 IF(lwp) WRITE(numout,*) '~~~~~~~' 68 68 69 IF( wrk_in_use(2, 1,2,3) ) THEN 70 CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.') ; RETURN 71 END IF 69 CALL wrk_alloc( jpi, jpj, zmbk, zprt, zprw ) 72 70 73 71 zmbk(:,:) = 0._wp … … 331 329 END SELECT 332 330 ! 333 IF( wrk_not_released(2, 1,2,3) ) CALL ctl_stop('dom_rea:failed to release workspace arrays.')331 CALL wrk_dealloc( jpi, jpj, zmbk, zprt, zprw ) 334 332 ! 335 333 END SUBROUTINE dom_rea … … 346 344 !! ** Action : - update mbathy: level bathymetry (in level index) 347 345 !!---------------------------------------------------------------------- 348 USE wrk_nemo, ONLY: zmbk => wrk_2d_4349 346 ! 350 347 INTEGER :: ji, jj ! dummy loop indices 348 REAL(wp), POINTER, DIMENSION(:,:) :: zmbk 351 349 !!---------------------------------------------------------------------- 352 350 … … 356 354 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~' 357 355 ! 358 IF( wrk_in_use(2, 4) ) THEN 359 CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.') ; RETURN 360 END IF 356 CALL wrk_alloc( jpi, jpj, zmbk ) 361 357 ! 362 358 mbkt(:,:) = MAX( mbathy(:,:) , 1 ) ! bottom k-index of T-level (=1 over land) … … 372 368 zmbk(:,:) = REAL( mbkv(:,:), wp ) ; CALL lbc_lnk(zmbk,'V',1.) ; mbkv (:,:) = MAX( INT( zmbk(:,:) ), 1 ) 373 369 ! 374 IF( wrk_not_released(2, 4) ) CALL ctl_stop('dom_rea:failed to release workspace arrays.')370 CALL wrk_dealloc( jpi, jpj, zmbk ) 375 371 ! 376 372 END SUBROUTINE zgr_bot_level
Note: See TracChangeset
for help on using the changeset viewer.