Changeset 12555 for branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
- Timestamp:
- 2020-03-13T17:33:02+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r6487 r12555 114 114 READ ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 ) 115 115 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp ) 116 IF(lwm ) WRITE ( numond, namzgr )116 IF(lwm .AND. nprint > 2) WRITE ( numond, namzgr ) 117 117 118 118 IF(lwp) THEN ! Control print … … 154 154 END IF 155 155 ! 156 IF( nprint == 1.AND. lwp ) THEN156 IF( nprint > 3 .AND. lwp ) THEN 157 157 WRITE(numout,*) ' MIN val mbathy ', MINVAL( mbathy(:,:) ), ' MAX ', MAXVAL( mbathy(:,:) ) 158 158 WRITE(numout,*) ' MIN val depth t ', MINVAL( gdept_0(:,:,:) ), & … … 170 170 & ' w ', MAXVAL( e3w_0(:,:,:) ) 171 171 ENDIF 172 ! 173 IF(lwp .AND. lflush) CALL flush(numout) 172 174 ! 173 175 IF( nn_timing == 1 ) CALL timing_stop('dom_zgr') … … 265 267 ENDIF 266 268 ENDIF 269 IF(lflush) CALL flush(numout) 267 270 ENDIF 268 271 … … 338 341 WRITE(numout, "(9x,' level gdept_1d gdepw_1d e3t_1d e3w_1d ')" ) 339 342 WRITE(numout, "(10x, i4, 4f9.2)" ) ( jk, gdept_1d(jk), gdepw_1d(jk), e3t_1d(jk), e3w_1d(jk), jk = 1, jpk ) 343 IF(lflush) CALL flush(numout) 340 344 ENDIF 341 345 DO jk = 1, jpk ! control positivity … … 598 602 ENDIF 599 603 ! 604 IF(lwp .AND. lflush) CALL flush(numout) 605 ! 600 606 IF( nn_timing == 1 ) CALL timing_stop('zgr_bat') 601 607 ! … … 651 657 ! 652 658 ENDIF 659 ! 660 IF(lwp .AND. lflush) CALL flush(numout) 653 661 ! 654 662 END SUBROUTINE zgr_bat_zoom … … 779 787 ENDIF 780 788 781 IF( lwp .AND. nprint == 1) THEN ! control print789 IF( lwp .AND. nprint > 3 ) THEN ! control print 782 790 WRITE(numout,*) 783 791 WRITE(numout,*) ' bathymetric field : number of non-zero T-levels ' … … 789 797 CALL wrk_dealloc( jpi, jpj, zbathy ) 790 798 ! 799 IF(lwp .AND. lflush) CALL flush(numout) 800 ! 791 801 IF( nn_timing == 1 ) CALL timing_stop('zgr_bat_ctl') 792 802 ! … … 818 828 IF(lwp) WRITE(numout,*) ' zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ' 819 829 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~' 830 IF(lwp .AND. lflush) CALL flush(numout) 820 831 ! 821 832 mbkt(:,:) = MAX( mbathy(:,:) , 1 ) ! bottom k-index of T-level (=1 over land) … … 862 873 IF(lwp) WRITE(numout,*) ' zgr_top_level : ocean top k-index of T-, U-, V- and W-levels ' 863 874 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~' 875 IF(lwp .AND. lflush) CALL flush(numout) 864 876 ! 865 877 mikt(:,:) = MAX( misfdep(:,:) , 1 ) ! top k-index of T-level (=1) … … 1246 1258 WRITE(numout,*) 1247 1259 WRITE(numout,*) 'domzgr gdep3w(mbathy)' ; CALL prihre(zprt(:,:,6),jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 1260 IF(lflush) CALL flush(numout) 1248 1261 ENDIF 1249 1262 ! … … 1273 1286 INTEGER :: id, jd, nprocd 1274 1287 INTEGER :: icompt, ibtest, ibtestim1, ibtestip1, ibtestjm1, ibtestjp1 ! (ISF) 1275 LOGICAL :: ll_print ! Allow control print for debugging1276 1288 REAL(wp) :: ze3tp , ze3wp ! Last ocean level thickness at T- and W-points 1277 1289 REAL(wp) :: zdepwp, zdepth ! Ajusted ocean depth to avoid too small e3t … … 1753 1765 ENDIF 1754 1766 1767 IF(lwp .AND. lflush) CALL flush(numout) 1768 1755 1769 CALL wrk_dealloc( jpi, jpj, zmask, zbathy, zrisfdep ) 1756 1770 CALL wrk_dealloc( jpi, jpj, zmisfdep, zmbathy ) … … 1829 1843 READ ( numnam_cfg, namzgr_sco, IOSTAT = ios, ERR = 902 ) 1830 1844 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr_sco in configuration namelist', lwp ) 1831 IF(lwm ) WRITE ( numond, namzgr_sco )1845 IF(lwm .AND. nprint > 2) WRITE ( numond, namzgr_sco ) 1832 1846 1833 1847 IF(lwp) THEN ! control print … … 1855 1869 WRITE(numout,*) ' Offset for Zb rn_zb_b = ',rn_zb_b 1856 1870 WRITE(numout,*) ' Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b' 1871 IF(lflush) CALL flush(numout) 1857 1872 ENDIF 1858 1873 … … 1984 1999 WRITE(numout,*) 1985 2000 CALL prihre( hbatt(1,1), jpi, jpj, 1, jpi, 1, 1, jpj, 1, 0._wp, numout ) 1986 IF( nprint == 1) THEN2001 IF( nprint > 3 ) THEN 1987 2002 WRITE(numout,*) ' bathy MAX ', MAXVAL( bathy(:,:) ), ' MIN ', MINVAL( bathy(:,:) ) 1988 2003 WRITE(numout,*) ' hbatt MAX ', MAXVAL( hbatt(:,:) ), ' MIN ', MINVAL( hbatt(:,:) ) … … 2045 2060 hiff(:,:) = MIN( hiff(:,:), hbatf(:,:) ) 2046 2061 2047 IF( nprint == 1.AND. lwp ) THEN2062 IF( nprint > 3 .AND. lwp ) THEN 2048 2063 WRITE(numout,*) ' MAX val hif t ', MAXVAL( hift (:,:) ), ' f ', MAXVAL( hiff (:,:) ), & 2049 2064 & ' u ', MAXVAL( hifu (:,:) ), ' v ', MAXVAL( hifv (:,:) ) … … 2139 2154 END DO 2140 2155 END DO 2141 IF( nprint == 1.AND. lwp ) WRITE(numout,*) ' MIN val mbathy h90 ', MINVAL( mbathy(:,:) ), &2156 IF( nprint > 3 .AND. lwp ) WRITE(numout,*) ' MIN val mbathy h90 ', MINVAL( mbathy(:,:) ), & 2142 2157 & ' MAX ', MAXVAL( mbathy(:,:) ) 2143 2158 2144 IF( nprint == 1.AND. lwp ) THEN ! min max values over the local domain2159 IF( nprint > 3 .AND. lwp ) THEN ! min max values over the local domain 2145 2160 WRITE(numout,*) ' MIN val mbathy ', MINVAL( mbathy(:,:) ), ' MAX ', MAXVAL( mbathy(:,:) ) 2146 2161 WRITE(numout,*) ' MIN val depth t ', MINVAL( gdept_0(:,:,:) ), & … … 2236 2251 END DO 2237 2252 ! 2253 IF(lwp .AND. lflush) CALL flush(numout) 2254 ! 2238 2255 CALL wrk_dealloc( jpi, jpj, zenv, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 ) 2239 2256 ! … … 2506 2523 z_gsigt(jk) = -fssig( REAL(jk,wp) ) 2507 2524 END DO 2508 IF( nprint == 1.AND. lwp ) WRITE(numout,*) 'z_gsigw 1 jpk ', z_gsigw(1), z_gsigw(jpk)2525 IF( nprint > 3 .AND. lwp ) WRITE(numout,*) 'z_gsigw 1 jpk ', z_gsigw(1), z_gsigw(jpk) 2509 2526 ! 2510 2527 ! Coefficients for vertical scale factors at w-, t- levels … … 2546 2563 END DO 2547 2564 END DO 2565 2566 IF(lwp .AND. lflush) CALL flush(numout) 2548 2567 2549 2568 CALL wrk_dealloc( jpk, z_gsigw, z_gsigt, z_gsi3w )
Note: See TracChangeset
for help on using the changeset viewer.