- Timestamp:
- 2011-03-03T17:13:18+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OFF_SRC/domrea.F90
r2528 r2648 16 16 USE dommsk ! domain: masks 17 17 USE lbclnk ! lateral boundary condition - MPP exchanges 18 USE in_out_manager ! I/O manager 18 USE lib_mpp 19 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 19 20 20 21 IMPLICIT NONE … … 53 54 !!---------------------------------------------------------------------- 54 55 USE iom 56 USE wrk_nemo, ONLY: zmbk => wrk_2d_1, zprt => wrk_2d_2, zprw => wrk_2d_3 55 57 !! 56 58 INTEGER :: ji, jj, jk ! dummy loop indices 57 59 INTEGER :: ik, inum0 , inum1 , inum2 , inum3 , inum4 ! local integers 58 60 REAL(wp) :: zrefdep ! local real 59 REAL(wp), DIMENSION(jpi,jpj) :: zmbk, zprt, zprw ! 2D workspace60 61 !!---------------------------------------------------------------------- 61 62 … … 63 64 IF(lwp) WRITE(numout,*) 'dom_rea : read NetCDF mesh and mask information file(s)' 64 65 IF(lwp) WRITE(numout,*) '~~~~~~~' 66 67 IF( wrk_in_use(2, 1,2,3) ) THEN 68 CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.') ; RETURN 69 END IF 65 70 66 71 zmbk(:,:) = 0._wp … … 141 146 CALL iom_get( inum3, jpdom_data, 'e2u', e2u ) 142 147 CALL iom_get( inum3, jpdom_data, 'e2v', e2v ) 148 149 e1e2t(:,:) = e1t(:,:) * e2t(:,:) ! surface at T-points 143 150 144 151 CALL iom_get( inum3, jpdom_data, 'ff', ff ) … … 314 321 END SELECT 315 322 ! 323 IF( wrk_not_released(2, 1,2,3) ) CALL ctl_stop('dom_rea:failed to release workspace arrays.') 324 ! 316 325 END SUBROUTINE dom_rea 317 326 … … 327 336 !! ** Action : - update mbathy: level bathymetry (in level index) 328 337 !!---------------------------------------------------------------------- 338 USE wrk_nemo, ONLY: zmbk => wrk_2d_4 339 ! 329 340 INTEGER :: ji, jj ! dummy loop indices 330 REAL(wp), DIMENSION(jpi,jpj) :: zmbk ! 2D workspace331 !!---------------------------------------------------------------------- 341 !!---------------------------------------------------------------------- 342 332 343 ! 333 344 IF(lwp) WRITE(numout,*) 334 345 IF(lwp) WRITE(numout,*) ' zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ' 335 346 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~' 347 ! 348 IF( wrk_in_use(2, 4) ) THEN 349 CALL ctl_stop('dom_rea: ERROR: requested workspace arrays unavailable.') ; RETURN 350 END IF 336 351 ! 337 352 mbkt(:,:) = MAX( mbathy(:,:) , 1 ) ! bottom k-index of T-level (=1 over land) … … 347 362 zmbk(:,:) = REAL( mbkv(:,:), wp ) ; CALL lbc_lnk(zmbk,'V',1.) ; mbkv (:,:) = MAX( INT( zmbk(:,:) ), 1 ) 348 363 ! 364 IF( wrk_not_released(2, 4) ) CALL ctl_stop('dom_rea:failed to release workspace arrays.') 365 ! 349 366 END SUBROUTINE zgr_bot_level 350 367
Note: See TracChangeset
for help on using the changeset viewer.