Changeset 4468
- Timestamp:
- 2014-02-04T13:18:14+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r4462 r4468 773 773 USE wrk_nemo, ONLY: zmbk => wrk_2d_1 774 774 USE par_oce, ONLY: jpkf, jpkfm1 775 USE arpdebugging, ONLY: dump_array 775 776 !! 776 777 INTEGER :: ji, jj ! dummy loop indices … … 788 789 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~' 789 790 ! 790 !CALL dump_array(0, 'mbathy', mbathy, withHalos=.TRUE.)791 CALL dump_array(0, 'mbathy', mbathy, withHalos=.TRUE.) 791 792 792 793 mbkt(:,:) = MAX( mbathy(:,:) , 1 ) ! bottom k-index of T-level (=1 over land) … … 819 820 jpkorig = jpk 820 821 IF( domtrim_z ) THEN 821 !!$ DO jj = 1, jpj 822 !!$ DO ji = 1, jpi 823 !!$ mbkmax(ji,jj) = MIN(jpk, MAX(mbkt(ji,jj)+1, mbku(ji,jj), mbkv(ji,jj))) 824 !!$ ! write(*,*) narea, ': SMPDBG: ', ji, jj, mbkt(ji,jj), mbku(ji,jj), mbkv(ji,jj), mbkmax(ji,jj) 825 !!$ END DO 826 !!$ END DO 822 DO jj = 1, jpj 823 DO ji = 1, jpi 824 mbkmax(ji,jj) = MIN(jpk, MAX(mbkt(ji,jj)+1, mbku(ji,jj), mbkv(ji,jj))) 825 ! write(*,*) narea, ': SMPDBG: ', ji, jj, mbkt(ji,jj), mbku(ji,jj), mbkv(ji,jj), mbkmax(ji,jj) 826 !!$ IF(mbkmax(ji,jj) /= jpkf)THEN 827 !!$ WRITE (*,*) narea,': ARPDBG: mbkmax at ',ji,jj,' in {mpp_init3,zgr_bot_level} = ',mbkmax(ji,jj),jpkf 828 !!$ END IF 829 END DO 830 END DO 831 832 CALL dump_array(0, 'mbkmax', mbkmax) 833 827 834 ! mbkmax computed in mpp_init3 because needed before ANY halo 828 835 ! swaps can be performed … … 1301 1308 USE mapcomm_mod, ONLY: trimmed, cyclic_bc 1302 1309 USE mapcomm_mod, ONLY: nidx, eidx, sidx, widx 1303 1310 USE arpdebugging, ONLY: dump_array 1304 1311 !! DCSE_NEMO: wrk_nemo module variables renamed, need additional directives 1305 1312 !FTRANS gsigw3 :I :I :z … … 1477 1484 ENDIF 1478 1485 1479 !CALL dump_array(0, 'hbatt', hbatt, withHalos=.FALSE.)1486 CALL dump_array(0, 'hbatt', hbatt, withHalos=.FALSE.) 1480 1487 1481 1488 ! ! ==============================
Note: See TracChangeset
for help on using the changeset viewer.