- Timestamp:
- 2014-02-04T13:19:11+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
r4475 r4479 772 772 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 773 773 USE wrk_nemo, ONLY: zmbk => wrk_2d_1 774 USE par_oce, ONLY: jpkf, jpkfm1775 774 USE arpdebugging, ONLY: dump_array 776 775 !! … … 824 823 mbkmax(ji,jj) = MIN(jpk, MAX(mbkt(ji,jj)+1, mbku(ji,jj), mbkv(ji,jj))) 825 824 ! write(*,*) narea, ': SMPDBG: ', ji, jj, mbkt(ji,jj), mbku(ji,jj), mbkv(ji,jj), mbkmax(ji,jj) 826 !!$ IF(mbkmax(ji,jj) /= jpkf)THEN827 !!$ WRITE (*,*) narea,': ARPDBG: mbkmax at ',ji,jj,' in {mpp_init3,zgr_bot_level} = ',mbkmax(ji,jj),jpkf828 !!$ END IF829 825 END DO 830 826 END DO … … 834 830 ! mbkmax computed in mpp_init3 because needed before ANY halo 835 831 ! swaps can be performed 836 jpkf = MAXVAL( mbkmax(:,:) )837 WRITE(*,*) narea,': ARPDBG: shallowest pt and jpkf = ', &838 MINVAL(mbkmax(:,:)), jpkf839 832 ! Play rather fast and loose and just change the value of jpk 840 833 ! here... 841 jpk = jpkf 834 jpk = MAXVAL( mbkmax(:,:) ) 835 WRITE(*,*) narea,': ARPDBG: shallowest pt and new, sub-domain jpk = ', & 836 MINVAL(mbkmax(:,:)), jpk 842 837 jpkm1 = jpk - 1 843 838 ELSE 844 839 WRITE(*,*) narea,': ARPDBG: NOT trimming domain in z' 845 840 mbkmax(:,:) = jpk 846 jpkf = jpk847 841 END IF 848 jpkfm1 = jpkf - 1849 842 850 843 !
Note: See TracChangeset
for help on using the changeset viewer.