Changeset 11101
- Timestamp:
- 2019-06-11T16:10:28+02:00 (4 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package/NEMOGCM
- Files:
-
- 159 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/CONFIG/SHARED/namelist_ref
r10149 r11101 1080 1080 &namctl ! Control prints & Benchmark 1081 1081 !----------------------------------------------------------------------- 1082 ln_ctl = .false. ! trends control print (expensive!) 1083 nn_print = 0 ! level of print (0 no extra print) 1084 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1085 nn_ictle = 0 ! end i indice of control sum multi processor runs 1086 nn_jctls = 0 ! start j indice of control over a subdomain) 1087 nn_jctle = 0 ! end j indice of control 1088 nn_isplt = 1 ! number of processors in i-direction 1089 nn_jsplt = 1 ! number of processors in j-direction 1090 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 1091 ! (no physical validity of the results) 1092 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1082 ln_ctl = .false. ! trends control print (expensive!) 1083 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1084 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1085 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1086 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1087 sn_cfctl%l_layout = .FALSE. ! 1088 sn_cfctl%l_mppout = .FALSE. ! 1089 sn_cfctl%l_mpptop = .FALSE. ! 1090 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1091 sn_cfctl%procmax = 0 ! Maximum area number for reporting [default:0] 1092 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 1093 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info/solver.stat file 1094 nn_print = 0 ! level of print (0 no extra print) 1095 ln_flush = .FALSE. ! flush numout after each print. NOT recommended for production run 1096 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1097 nn_ictle = 0 ! end i indice of control sum multi processor runs 1098 nn_jctls = 0 ! start j indice of control over a subdomain) 1099 nn_jctle = 0 ! end j indice of control 1100 nn_isplt = 1 ! number of processors in i-direction 1101 nn_jsplt = 1 ! number of processors in j-direction 1102 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 1103 ! (no physical validity of the results) 1104 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1093 1105 / 1094 1106 !----------------------------------------------------------------------- -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90
r10477 r11101 767 767 & TRIM( c_asmbal ) // ' at timestep = ', kt 768 768 ENDIF 769 770 IF(lwp .AND. lflush) CALL flush(numout) 769 771 770 772 END SUBROUTINE asm_bgc_bal_wri … … 1140 1142 ENDIF 1141 1143 ! 1144 IF(lwp .AND. lflush) CALL flush(numout) 1142 1145 END SUBROUTINE phyto2d_asm_inc 1143 1146 … … 1327 1330 ENDIF 1328 1331 ! 1332 IF(lwp .AND. lflush) CALL flush(numout) 1329 1333 END SUBROUTINE phyto3d_asm_inc 1330 1334 … … 1598 1602 ENDIF 1599 1603 ! 1604 IF(lwp .AND. lflush) CALL flush(numout) 1600 1605 END SUBROUTINE pco2_asm_inc 1601 1606 … … 1821 1826 #endif 1822 1827 ! 1828 IF(lwp .AND. lflush) CALL flush(numout) 1829 1823 1830 END SUBROUTINE ph_asm_inc 1824 1831 … … 2186 2193 ENDIF 2187 2194 ! 2195 IF(lwp .AND. lflush) CALL flush(numout) 2188 2196 END SUBROUTINE bgc3d_asm_inc 2189 2197 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90
r10149 r11101 183 183 ENDIF 184 184 ! 185 IF(lwp .AND. lflush) CALL flush(numout) 185 186 END SUBROUTINE asm_bkg_wri 186 187 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r10149 r11101 164 164 READ ( numnam_cfg, nam_asminc, IOSTAT = ios, ERR = 902 ) 165 165 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_asminc in configuration namelist', lwp ) 166 IF(lwm ) WRITE ( numond, nam_asminc )166 IF(lwm .AND. nprint > 2) WRITE ( numond, nam_asminc ) 167 167 168 168 ! Control print … … 207 207 WRITE(numout,*) ' Choice of MLD for BGC assimilation mld_choice_bgc = ', mld_choice_bgc 208 208 WRITE(numout,*) ' Maximum absolute chlorophyll increment (<=0 = off) rn_maxchlinc = ', rn_maxchlinc 209 IF(lflush) CALL flush(numout) 209 210 ENDIF 210 211 … … 242 243 WRITE(numout,*) ' iitiaustr_date = ', iitiaustr_date 243 244 WRITE(numout,*) ' iitiaufin_date = ', iitiaufin_date 245 IF(lflush) CALL flush(numout) 244 246 ENDIF 245 247 … … 369 371 WRITE(numout,*) ' Time-integrated weight = ', ztotwgt 370 372 WRITE(numout,*) ' ===================================' 373 IF(lflush) CALL flush(numout) 371 374 ENDIF 372 375 … … 590 593 ENDIF 591 594 ! 595 IF(lflush) CALL flush(numout) 596 592 597 END SUBROUTINE asm_inc_init 593 598 … … 735 740 WRITE(numout,*) 'tra_asm_inc : Tracer IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 736 741 WRITE(numout,*) '~~~~~~~~~~~~' 742 IF(lflush) CALL flush(numout) 737 743 ENDIF 738 744 … … 864 870 & kt,' with IAU weight = ', wgtiau(it) 865 871 WRITE(numout,*) '~~~~~~~~~~~~' 872 IF(lflush) CALL flush(numout) 866 873 ENDIF 867 874 … … 940 947 & kt,' with IAU weight = ', wgtiau(it) 941 948 WRITE(numout,*) '~~~~~~~~~~~~' 949 IF(lflush) CALL flush(numout) 942 950 ENDIF 943 951 … … 1030 1038 & kt,' with IAU weight = ', wgtiau(it) 1031 1039 WRITE(numout,*) '~~~~~~~~~~~~' 1040 IF(lflush) CALL flush(numout) 1032 1041 ENDIF 1033 1042 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90
r8447 r11101 255 255 READ ( numnam_cfg, nambias, IOSTAT = ios, ERR = 902 ) 256 256 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambias in configuration namelist', lwp ) 257 IF(lwm ) WRITE ( numond, nambias )257 IF(lwm .AND. nprint > 2) WRITE ( numond, nambias ) 258 258 259 259 … … 320 320 321 321 WRITE(numout,*) ' time step is = ',rdt,'you choose to write pcbias at nn_bias_itwrt = ',nn_bias_itwrt,'and end of iau is rday/rdt=',rday/rdt 322 IF(lflush) call flush(numout) 322 323 ENDIF 323 324 IF( .NOT. ln_bias ) RETURN … … 561 562 ENDWHERE 562 563 ENDIF 563 564 565 IF(lwp .AND. lflush) call flush(numout) 564 566 565 567 END SUBROUTINE bias_init … … 895 897 ENDIF 896 898 ENDIF 897 899 IF(lwp .AND. lflush) call flush(numout) 898 900 END SUBROUTINE tra_bias 899 901 … … 971 973 & rhd_pc, gru_pc , grv_pc ) 972 974 ENDIF 975 976 IF(lwp .AND. lflush) call flush(numout) 973 977 974 978 END SUBROUTINE dyn_bias … … 1005 1009 ENDIF 1006 1010 ! 1011 IF(lwp .AND. lflush) call flush(numout) 1012 1007 1013 END SUBROUTINE bias_opn 1008 1014 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90
r6487 r11101 532 532 READ ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 ) 533 533 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp ) 534 IF(lwm ) WRITE ( numond, nambdy_dta )534 IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy_dta ) 535 535 536 536 cn_dir_array(ib_bdy) = cn_dir … … 893 893 894 894 ENDDO ! ib_bdy 895 896 IF(lwp .AND. lflush) CALL flush(numout) 895 897 896 898 IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init') -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90
r6487 r11101 110 110 IF( nn_timing == 1 ) CALL timing_start('bdy_init') 111 111 112 IF(lwp) WRITE(numout,*) 113 IF(lwp) WRITE(numout,*) 'bdy_init : initialization of open boundaries' 114 IF(lwp) WRITE(numout,*) '~~~~~~~~' 112 IF(lwp) THEN 113 WRITE(numout,*) 114 WRITE(numout,*) 'bdy_init : initialization of open boundaries' 115 WRITE(numout,*) '~~~~~~~~' 116 IF(lflush) CALL flush(numout) 117 ENDIF 115 118 ! 116 119 … … 131 134 READ ( numnam_cfg, nambdy, IOSTAT = ios, ERR = 902 ) 132 135 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in configuration namelist', lwp ) 133 IF(lwm ) WRITE ( numond, nambdy )136 IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy ) 134 137 135 138 ! ----------------------------------------- … … 420 423 READ ( numnam_cfg, nambdy_index, IOSTAT = ios, ERR = 904 ) 421 424 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_index in configuration namelist', lwp ) 422 IF(lwm ) WRITE ( numond, nambdy_index )425 IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy_index ) 423 426 424 427 SELECT CASE ( TRIM(ctypebdy) ) … … 1348 1351 ENDIF 1349 1352 1353 IF(lwp .AND. lflush) CALL flush(numout) 1354 1350 1355 CALL wrk_dealloc(jpi,jpj,zfmask) 1351 1356 … … 1722 1727 IF ((itest==0).AND.(lwp)) WRITE(numout,*) 'NO open boundary corner found' 1723 1728 ! 1729 IF(lwp. AND. lflush) CALL flush(numout) 1730 ! 1724 1731 ! Other tests TBD: 1725 1732 ! segments completly on land … … 1757 1764 IF(lwp) WRITE(numout,*) ' ========== have different open bdy schemes' 1758 1765 IF(lwp) WRITE(numout,*) 1766 IF(lflush) CALL flush(numout) 1759 1767 nstop = nstop + 1 1760 1768 ENDIF -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90
r9583 r11101 99 99 IF(lwp) WRITE(numout,*) 'bdytide_init : initialization of tidal harmonic forcing at open boundaries' 100 100 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 101 IF(lflush) CALL flush(numout) 101 102 ENDIF 102 103 … … 118 119 READ ( numnam_cfg, nambdy_tide, IOSTAT = ios, ERR = 902 ) 119 120 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_tide in configuration namelist', lwp ) 120 IF(lwm ) WRITE ( numond, nambdy_tide )121 IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy_tide ) 121 122 ! ! Parameter control and print 122 IF(lwp) WRITE(numout,*) ' '123 IF(lwp) WRITE(numout,*) ' Namelist nambdy_tide : tidal harmonic forcing at open boundaries'124 IF(lwp) WRITE(numout,*) ' read tidal data in 2d files: ', ln_bdytide_2ddta125 IF(lwp) WRITE(numout,*) ' assume complex conjugate : ', ln_bdytide_conj126 IF(lwp) WRITE(numout,*) ' Number of tidal components to read: ', nb_harmo127 IF(lwp) THEN128 123 IF(lwp) THEN 124 WRITE(numout,*) ' ' 125 WRITE(numout,*) ' Namelist nambdy_tide : tidal harmonic forcing at open boundaries' 126 WRITE(numout,*) ' read tidal data in 2d files: ', ln_bdytide_2ddta 127 WRITE(numout,*) ' assume complex conjugate : ', ln_bdytide_conj 128 WRITE(numout,*) ' Number of tidal components to read: ', nb_harmo 129 WRITE(numout,*) ' Tidal components: ' 129 130 DO itide = 1, nb_harmo 130 131 WRITE(numout,*) ' ', Wave(ntide(itide))%cname_tide 131 132 END DO 132 ENDIF 133 IF(lwp) WRITE(numout,*) ' ' 133 WRITE(numout,*) ' ' 134 IF(lflush) CALL flush(numout) 135 ENDIF 134 136 135 137 ! Allocate space for tidal harmonics data - get size from OBC data arrays … … 333 335 WRITE(numout,*) 'bdytide_update : (re)Initialization of the tidal bdy forcing at kt=',kt 334 336 WRITE(numout,*) '~~~~~~~~~~~~~~ ' 337 IF(lflush) CALL flush(numout) 335 338 ENDIF 336 339 ! … … 448 451 WRITE(numout,*) 'bdy_tide_dta : Refresh nodal factors for tidal open bdy data at kt=',kt 449 452 WRITE(numout,*) '~~~~~~~~~~~~~~ ' 453 IF(lflush) CALL flush(numout) 450 454 ENDIF 451 455 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/BDY/bdyvol.F90
r9583 r11101 87 87 IF(lwp) WRITE(numout,*)'bdy_vol : Correction of velocities along unstructured OBC' 88 88 IF(lwp) WRITE(numout,*)'~~~~~~~' 89 IF(lflush) CALL flush(numout) 89 90 END IF 90 91 … … 164 165 IF(lwp) WRITE(numout,*)' correction velocity zubtpecor =', zubtpecor , '(m/s)' 165 166 IF(lwp) WRITE(numout,*)' cumulated transport ztranst =', ztranst , '(m3/s)' 167 IF(lflush) CALL flush(numout) 166 168 END IF 167 169 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/C1D/c1d.F90
r6486 r11101 56 56 READ ( numnam_cfg, namc1d, IOSTAT = ios, ERR = 902 ) 57 57 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc1d in configuration namelist', lwp ) 58 IF(lwm ) WRITE ( numond, namc1d )58 IF(lwm .AND. nprint > 2) WRITE ( numond, namc1d ) 59 59 60 60 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/C1D/dtauvd.F90
r6486 r11101 73 73 READ ( numnam_cfg, namc1d_uvd, IOSTAT = ios, ERR = 902 ) 74 74 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc1d_uvd in configuration namelist', lwp ) 75 IF(lwm ) WRITE ( numond, namc1d_uvd )75 IF(lwm .AND. nprint > 2) WRITE ( numond, namc1d_uvd ) 76 76 77 77 ! ! force the initialization when dyndmp is used -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/C1D/dyndmp.F90
r6486 r11101 87 87 READ ( numnam_cfg, namc1d_dyndmp, IOSTAT = ios, ERR = 902 ) 88 88 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc1d_dyndmp in configuration namelist', lwp ) 89 IF(lwm ) WRITE ( numond, namc1d_dyndmp )89 IF(lwm .AND. nprint > 2) WRITE ( numond, namc1d_dyndmp ) 90 90 91 91 IF(lwp) THEN ! control print … … 100 100 WRITE(numout,*) ' Damping file name cn_resto = ', cn_resto 101 101 WRITE(numout,*) 102 IF(lflush) CALL flush(numout) 102 103 ENDIF 103 104 … … 128 129 CALL iom_close( imask ) 129 130 ENDIF 131 ! 132 IF(lwp .AND. lflush) CALL flush(numout) 130 133 ! 131 134 END SUBROUTINE dyn_dmp_init -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90
r6486 r11101 93 93 READ ( numnam_cfg, namcrs, IOSTAT = ios, ERR = 902 ) 94 94 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcrs in configuration namelist', lwp ) 95 IF(lwm ) WRITE ( numond, namcrs )95 IF(lwm .AND. nprint > 2) WRITE ( numond, namcrs ) 96 96 97 97 IF(lwp) THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r9583 r11101 151 151 READ ( numnam_cfg, namdct, IOSTAT = ios, ERR = 902 ) 152 152 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdct in configuration namelist', lwp ) 153 IF(lwm ) WRITE ( numond, namdct )153 IF(lwm .AND. nprint >2) WRITE ( numond, namdct ) 154 154 155 155 IF( lwp ) THEN … … 185 185 transports_3d(:,:,:,:)=0.0 186 186 transports_2d(:,:,:) =0.0 187 187 ! 188 IF(lwp .AND. lflush) CALL flush(numout) 189 ! 188 190 IF( nn_timing == 1 ) CALL timing_stop('dia_dct_init') 189 191 ! … … 301 303 CALL wrk_dealloc( nb_sec_max,nb_type_class,nb_class_max , zsum ) 302 304 ENDIF 303 305 ! 306 IF(lwp .AND. lflush) CALL flush(numout) 307 ! 304 308 IF( nn_timing == 1 ) CALL timing_stop('dia_dct') 305 309 ! … … 503 507 504 508 nb_sec = jsec-1 !number of section read in the file 505 509 ! 510 IF(lwp .AND. lflush) CALL flush(numout) 511 506 512 CALL wrk_dealloc( nb_point_max, directemp ) 507 513 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90
r6487 r11101 95 95 READ ( numnam_cfg, nam_diaharm, IOSTAT = ios, ERR = 902 ) 96 96 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diaharm in configuration namelist', lwp ) 97 IF(lwm ) WRITE ( numond, nam_diaharm )97 IF(lwm .AND. nprint > 2) WRITE ( numond, nam_diaharm ) 98 98 ! 99 99 IF(lwp) THEN … … 160 160 ALLOCATE( ana_temp(jpi,jpj,2*nb_ana,3) ) 161 161 ana_temp(:,:,:,:) = 0._wp 162 163 IF(lwp .AND. lflush) CALL flush(numout) 162 164 163 165 END SUBROUTINE dia_harm_init … … 231 233 CALL wrk_alloc( jpi , jpj , jpmax_harmo , 2 , ana_amp ) 232 234 233 IF(lwp) WRITE(numout,*) 234 IF(lwp) WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 235 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 235 IF(lwp) THEN 236 WRITE(numout,*) 237 WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 238 WRITE(numout,*) '~~~~~~~~~~~~' 239 IF(lflush) CALL flush(numout) 240 ENDIF 236 241 237 242 ztime_ini = nit000_han*rdt ! Initial time in seconds at the beginning of analysis … … 400 405 #endif 401 406 IF(lwp) WRITE(numout,*) ' ' 407 IF(lwp .AND. lflush) CALL flush(numout) 402 408 403 409 ! A) Elevation -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
r9321 r11101 361 361 READ ( numnam_cfg, namhsb, IOSTAT = ios, ERR = 902 ) 362 362 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namhsb in configuration namelist', lwp ) 363 IF(lwm ) WRITE ( numond, namhsb )363 IF(lwm .AND. nprint > 2) WRITE ( numond, namhsb ) 364 364 365 365 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90
r9583 r11101 131 131 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 132 132 IF(lwp) WRITE(numout,*) 133 IF(lwp .AND. lflush) CALL flush(numout) 133 134 ENDIF 134 135 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90
r9830 r11101 469 469 READ ( numnam_cfg, namptr, IOSTAT = ios, ERR = 902 ) 470 470 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namptr in configuration namelist', lwp ) 471 IF(lwm ) WRITE ( numond, namptr )471 IF(lwm .AND. nprint > 2) WRITE ( numond, namptr ) 472 472 473 473 IF(lwp) THEN ! Control print … … 478 478 WRITE(numout,*) ' Poleward heat & salt transport (T) or not (F) ln_diaptr = ', ln_diaptr 479 479 WRITE(numout,*) ' Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = ', ln_subbas 480 IF(lflush) CALL flush(numout) 480 481 ENDIF 481 482 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r10020 r11101 462 462 INTEGER, INTENT( in ) :: kt ! ocean time-step index 463 463 !! 464 LOGICAL :: ll_print = .FALSE. ! =T print and flush numout465 464 CHARACTER (len=40) :: clhstnam, clop, clmx ! local names 466 465 INTEGER :: inum = 11 ! temporary logical unit … … 489 488 ! ----------------- 490 489 491 ! local variable for debugging492 ll_print = .FALSE.493 ll_print = ll_print .AND. lwp494 495 490 ! Define frequency of output and means 496 491 zdt = rdt … … 532 527 IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma, & 533 528 ' limit storage in depth = ', ipk 534 529 535 530 ! WRITE root name in date.file for use by postpro 536 531 IF(lwp) THEN … … 806 801 CALL histend( nid_W, snc4chunks=snc4set ) 807 802 808 IF(lwp) WRITE(numout,*) 809 IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization' 810 IF(ll_print) CALL FLUSH(numout ) 803 IF(lwp) THEN 804 WRITE(numout,*) 805 WRITE(numout,*) 'End of NetCDF Initialization' 806 IF(lflush) CALL flush(numout) 807 ENDIF 811 808 812 809 ENDIF … … 987 984 ! 988 985 IF( nn_timing == 1 ) CALL timing_stop('dia_wri') 986 ! 987 IF(lwp .AND. lflush) CALL flush(numout) 989 988 ! 990 989 END SUBROUTINE dia_wri -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90
r6491 r11101 78 78 IF(lwp) WRITE(numout,*)'dom_clo : closed seas ' 79 79 IF(lwp) WRITE(numout,*)'~~~~~~~' 80 IF(lwp .AND. lflush) CALL flush(numout) 80 81 81 82 ! initial values … … 403 404 CALL lbc_lnk( emp , 'T', 1._wp ) 404 405 ! 406 IF(lwp .AND. lflush) CALL flush(numout) 407 ! 405 408 IF( nn_timing == 1 ) CALL timing_stop('sbc_clo') 406 409 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90
r9321 r11101 130 130 IF(lwp) WRITE(numout,'(a,i6,a,i2,a,i2,a,i8,a,i8)')' =======>> 1/2 time step before the start of the run DATE Y/M/D = ', & 131 131 & nyear, '/', nmonth, '/', nday, ' nsec_day:', nsec_day, ' nsec_week:', nsec_week 132 IF(lwp .AND. lflush) CALL flush(numout) 132 133 133 134 ! Up to now, calendar parameters are related to the end of previous run (nit000-1) … … 253 254 CALL ymds2ju( nyear, 01, 01, 0.0, fjulstartyear ) 254 255 ! 255 IF(lwp) WRITE(numout,'(a,i8,a,i4.4,a,i2.2,a,i2.2,a,i3.3)') '======>> time-step =', kt, & 256 IF(lwp .AND. nprint > 0) THEN 257 WRITE(numout,'(a,i8,a,i4.4,a,i2.2,a,i2.2,a,i3.3)') '======>> time-step =', kt, & 256 258 & ' New day, DATE Y/M/D = ', nyear, '/', nmonth, '/', nday, ' nday_year = ', nday_year 257 IF(lwp)WRITE(numout,'(a,i8,a,i7,a,i5)') ' nsec_year = ', nsec_year, &259 WRITE(numout,'(a,i8,a,i7,a,i5)') ' nsec_year = ', nsec_year, & 258 260 & ' nsec_month = ', nsec_month, ' nsec_day = ', nsec_day, ' nsec_week = ', nsec_week 261 IF(lflush) CALL flush(numout) 262 ENDIF 259 263 ENDIF 260 264 … … 321 325 CASE ( 2 ) ; WRITE(numout,*) ' nrstdt = 2 : calendar parameters read in restart' 322 326 END SELECT 323 WRITE(numout,*) 327 WRITE(numout,*) ' ' 328 IF(lflush) CALL flush(numout) 324 329 ENDIF 325 330 ! Control of date … … 353 358 WRITE(numout,*) ' number of elapsed days since the begining of run : ', adatrj 354 359 WRITE(numout,*) 360 IF(lflush) CALL flush(numout) 355 361 ENDIF 356 362 ! 357 363 ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN 358 364 ! 359 IF( kt == nitrst ) THEN 360 IF(lwp) WRITE(numout,*) 361 IF(lwp) WRITE(numout,*) 'rst_write : write oce restart file kt =', kt 362 IF(lwp) WRITE(numout,*) '~~~~~~~' 365 IF( kt == nitrst .AND. lwp .AND. nprint > 0) THEN 366 WRITE(numout,*) 367 WRITE(numout,*) 'rst_write : write oce restart file kt =', kt 368 WRITE(numout,*) '~~~~~~~' 369 IF(lflush) CALL flush(numout) 363 370 ENDIF 364 371 ! calendar control -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90
r6491 r11101 79 79 WRITE(numout,*) 'dom_init : domain initialization' 80 80 WRITE(numout,*) '~~~~~~~~' 81 IF(lflush) CALL flush(numout) 81 82 ENDIF 82 83 ! … … 160 161 READ ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 ) 161 162 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp ) 162 IF(lwm ) WRITE ( numond, namrun )163 IF(lwm .AND. nprint > 2) WRITE ( numond, namrun ) 163 164 ! 164 165 IF(lwp) THEN ! control print … … 193 194 WRITE(numout,*) ' overwrite an existing file ln_clobber = ', ln_clobber 194 195 WRITE(numout,*) ' NetCDF chunksize (bytes) nn_chunksz = ', nn_chunksz 196 IF(lflush) CALL flush(numout) 195 197 ENDIF 196 198 … … 251 253 READ ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 ) 252 254 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp ) 253 IF(lwm ) WRITE ( numond, namdom )255 IF(lwm .AND. nprint > 2) WRITE ( numond, namdom ) 254 256 255 257 IF(lwp) THEN … … 293 295 WRITE(numout,*) ' ppkth2 = ', ppkth2 294 296 WRITE(numout,*) ' ppacr2 = ', ppacr2 297 IF(lflush) CALL flush(numout) 295 298 ENDIF 296 299 … … 313 316 READ ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 314 317 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 315 IF(lwm ) WRITE( numond, namcla )318 IF(lwm .AND. nprint > 2) WRITE( numond, namcla ) 316 319 317 320 IF(lwp) THEN … … 319 322 WRITE(numout,*) ' Namelist namcla' 320 323 WRITE(numout,*) ' cross land advection nn_cla = ', nn_cla 324 IF(lflush) CALL flush(numout) 321 325 ENDIF 322 326 IF ( nn_cla .EQ. 1 ) THEN … … 337 341 READ ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 ) 338 342 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in configuration namelist', lwp ) 339 IF(lwm ) WRITE( numond, namnc4 )343 IF(lwm .AND. nprint > 2) WRITE( numond, namnc4 ) 340 344 341 345 IF(lwp) THEN ! control print … … 346 350 WRITE(numout,*) ' number of chunks in k-dimension nn_nchunks_k = ', nn_nchunks_k 347 351 WRITE(numout,*) ' apply netcdf4/hdf5 chunking & compression ln_nc4zip = ', ln_nc4zip 352 IF(lflush) CALL flush(numout) 348 353 ENDIF 349 354 … … 406 411 WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, iima2, ijma2 407 412 WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2 413 IF(lflush) CALL flush(numout) 408 414 ENDIF 409 415 ! … … 467 473 WRITE(numout,*) 'dom_stiff : maximum grid stiffness ratio: ', zrxmax 468 474 WRITE(numout,*) '~~~~~~~~~' 475 IF(lflush) CALL flush(numout) 469 476 ENDIF 470 477 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90
r6486 r11101 56 56 IF( jperio == 5 ) WRITE(numout,*) ' jperio= 5, north fold with F-point pivot' 57 57 IF( jperio == 6 ) WRITE(numout,*) ' jperio= 6, cyclic east-west and north fold with F-point pivot' 58 IF(lflush) CALL flush(numout) 58 59 ENDIF 59 60 ! … … 105 106 mj1(jj) = MAX( 0, MIN( jj - jpjzoom + 1 - njmpp + 1, jpj ) ) 106 107 END DO 107 IF(lwp ) THEN ! control print108 IF(lwp .AND. nprint > 0) THEN ! control print 108 109 WRITE(numout,*) 109 110 WRITE(numout,*) 'dom_glo : domain: data / local ' … … 136 137 WRITE(numout,*) ' ending index' 137 138 WRITE(numout,25) (mj1(jj),jj = 1,jpjdta) 139 IF(lflush) CALL flush(numout) 138 140 ENDIF 139 141 25 FORMAT( 100(10x,19i4,/) ) … … 195 197 ENDIF 196 198 ! 199 IF(lwp .AND. lflush) CALL flush(numout) 200 ! 197 201 END SUBROUTINE dom_glo 198 202 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90
r6487 r11101 421 421 gphi0 = zphi1 + zsin_alpha * ze1deg * FLOAT( jpjglo-2 ) 422 422 423 IF( nprint ==1.AND. lwp ) THEN423 IF( nprint > 3 .AND. lwp ) THEN 424 424 WRITE(numout,*) ' ze1', ze1, 'cosalpha', zcos_alpha, 'sinalpha', zsin_alpha 425 425 WRITE(numout,*) ' ze1deg', ze1deg, 'glam0', glam0, 'gphi0', gphi0 … … 507 507 508 508 509 IF( nprint == 1.AND. lwp ) THEN509 IF( nprint > 3 .AND. lwp ) THEN 510 510 WRITE(numout,*) ' e1u e2u ' 511 511 CALL prihre( e1u,jpi,jpj,jpi-5,jpi,1,jpj-5,jpj,1,0.,numout ) … … 600 600 ENDIF 601 601 ! 602 IF(lwp .AND. lflush) CALL flush(numout) 603 ! 602 604 IF( nn_timing == 1 ) CALL timing_stop('dom_hgr') 603 605 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
r9830 r11101 157 157 READ ( numnam_cfg, namlbc, IOSTAT = ios, ERR = 902 ) 158 158 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlbc in configuration namelist', lwp ) 159 IF(lwm ) WRITE ( numond, namlbc )159 IF(lwm .AND. nprint > 2) WRITE ( numond, namlbc ) 160 160 161 161 IF(lwp) THEN ! control print … … 514 514 ! CAUTION : The fmask may be further modified in dyn_vor_init ( dynvor.F90 ) 515 515 516 IF( nprint == 1.AND. lwp ) THEN ! Control print516 IF( nprint > 3 .AND. lwp ) THEN ! Control print 517 517 imsk(:,:) = INT( tmask_i(:,:) ) 518 518 WRITE(numout,*) ' tmask_i : ' … … 561 561 CALL wrk_dealloc( jpi, jpj, imsk ) 562 562 CALL wrk_dealloc( jpi, jpj, zwf ) 563 ! 564 IF(lwp .AND. lflush) CALL flush(numout) 563 565 ! 564 566 IF( nn_timing == 1 ) CALL timing_stop('dom_msk') … … 640 642 nicoa(inw,1,jk) = ji 641 643 njcoa(inw,1,jk) = jj 642 IF( nprint == 1) WRITE(numout,*) ' west : ', jk, inw, ji, jj644 IF( nprint > 3 ) WRITE(numout,*) ' west : ', jk, inw, ji, jj 643 645 ENDIF 644 646 zaa = tmask(ji,jj,jk) + tmask(ji,jj+1,jk) … … 647 649 nicoa(ine,2,jk) = ji 648 650 njcoa(ine,2,jk) = jj 649 IF( nprint == 1) WRITE(numout,*) ' east : ', jk, ine, ji, jj651 IF( nprint > 3 ) WRITE(numout,*) ' east : ', jk, ine, ji, jj 650 652 ENDIF 651 653 END DO … … 667 669 nicoa(ins,3,jk) = ji 668 670 njcoa(ins,3,jk) = jj 669 IF( nprint == 1) WRITE(numout,*) ' south : ', jk, ins, ji, jj671 IF( nprint > 3 ) WRITE(numout,*) ' south : ', jk, ins, ji, jj 670 672 ENDIF 671 673 zaa = tmask(ji+1,jj,jk) + tmask(ji,jj,jk) … … 674 676 nicoa(inn,4,jk) = ji 675 677 njcoa(inn,4,jk) = jj 676 IF( nprint == 1) WRITE(numout,*) ' north : ', jk, inn, ji, jj678 IF( nprint > 3 ) WRITE(numout,*) ' north : ', jk, inn, ji, jj 677 679 ENDIF 678 680 END DO -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90
r6486 r11101 105 105 106 106 END SELECT 107 108 IF(lwp .AND. lflush) CALL flush(numout) 107 109 108 110 END SUBROUTINE dom_stp -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90
r9321 r11101 133 133 IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 134 134 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 135 IF(lwp .AND. lflush) CALL flush(numout) 135 136 136 137 ! choose vertical coordinate (z_star, z_tilde or layer) … … 292 293 IF(lwp) WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors' 293 294 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 295 IF(lflush) CALL flush(numout) 294 296 ENDIF 295 297 … … 483 485 IF( ln_vvl_dbg .AND. .NOT. ll_do_bclinic ) THEN ! - ML - test: control prints for debuging 484 486 ! 485 IF( lwp ) WRITE(numout, *) 'kt =', kt487 IF( lwp .AND. nprint > 2) WRITE(numout, *) 'kt =', kt 486 488 IF ( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 487 489 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( zht(:,:) ) ) … … 525 527 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain 526 528 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 529 IF(lflush) CALL flush(numout) 527 530 END IF 528 531 … … 593 596 IF(lwp) WRITE(numout,*) 'dom_vvl_sf_swp : - time filter and swap of scale factors' 594 597 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ - interpolate scale factors and compute depths for next time step' 598 IF(lflush) CALL flush(numout) 595 599 ENDIF 596 600 … … 866 870 neuler = 0 867 871 ENDIF 872 ! 873 IF(lwp .AND. lflush) CALL flush(numout) 868 874 ! ! ----------- ! 869 875 IF( ln_vvl_zstar ) THEN ! z_star case ! … … 907 913 ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file 908 914 ! ! =================== 909 IF(lwp ) WRITE(numout,*) '---- dom_vvl_rst ----'915 IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- dom_vvl_rst ----' 910 916 ! ! --------- ! 911 917 ! ! all cases ! … … 955 961 READ ( numnam_cfg, nam_vvl, IOSTAT = ios, ERR = 902 ) 956 962 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_vvl in configuration namelist', lwp ) 957 IF(lwm ) WRITE ( numond, nam_vvl )963 IF(lwm .AND. nprint > 2) WRITE ( numond, nam_vvl ) 958 964 959 965 IF(lwp) THEN ! Namelist print … … 1008 1014 ! IF( ln_vvl_kepe ) WRITE(numout,*) ' kinetic to potential energy transfer : option used' 1009 1015 ! IF( .NOT. ln_vvl_kepe ) WRITE(numout,*) ' kinetic to potential energy transfer : option not used' 1016 IF(lflush) CALL flush(numout) 1010 1017 ENDIF 1011 1018 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90
r9830 r11101 90 90 IF(lwp) WRITE(numout,*) 'dom_wri : create NetCDF mesh and mask information file(s)' 91 91 IF(lwp) WRITE(numout,*) '~~~~~~~' 92 IF(lwp .AND. lflush) CALL flush(numout) 92 93 93 94 clnam0 = 'mesh_mask' ! filename (mesh and mask informations) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r6487 r11101 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 ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90
r6486 r11101 77 77 READ ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 ) 78 78 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtsd in configuration namelist', lwp ) 79 IF(lwm ) WRITE ( numond, namtsd )79 IF(lwm .AND. nprint > 2) WRITE ( numond, namtsd ) 80 80 81 81 IF( PRESENT( ld_tradmp ) ) ln_tsd_tradmp = .TRUE. ! forces the initialization when tradmp is used … … 93 93 WRITE(numout,*) ' T & S data not used' 94 94 ENDIF 95 IF(lflush) CALL flush(numout) 95 96 ENDIF 96 97 ! … … 216 217 WRITE(numout,*) 217 218 WRITE(numout,*) 'dta_tsd: interpolates T & S data onto the s- or mixed s-z-coordinate mesh' 219 IF(lflush) CALL flush(numout) 218 220 ENDIF 219 221 ! … … 275 277 ENDIF 276 278 ! 277 IF( lwp .AND. kt == nit000 ) THEN279 IF( lwp .AND. kt == nit000 .AND. nprint > 1) THEN 278 280 WRITE(numout,*) ' temperature Levitus ' 279 281 WRITE(numout,*) … … 294 296 CALL prihre( ptsd(:,:,jpkm1,jp_sal), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout ) 295 297 WRITE(numout,*) 298 IF(lflush) CALL flush(numout) 296 299 ENDIF 297 300 ! … … 299 302 ! (data used only for initialisation) 300 303 IF(lwp) WRITE(numout,*) 'dta_tsd: deallocte T & S arrays as they are only use to initialize the run' 304 IF(lwp .AND. lflush) CALL flush(numout) 301 305 DEALLOCATE( sf_tsd(jp_tem)%fnow ) ! T arrays in the structure 302 306 IF( sf_tsd(jp_tem)%ln_tint ) DEALLOCATE( sf_tsd(jp_tem)%fdta ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90
r6486 r11101 79 79 IF(lwp) WRITE(numout,*) 'istate_ini : Initialization of the dynamics and tracers' 80 80 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 81 IF(lwp .AND. lflush) CALL flush(numout) 81 82 82 83 CALL dta_tsd_init ! Initialisation of T & S input data … … 209 210 IF(lwp) WRITE(numout,*) 'istate_t_s : Philander s initial temperature profile' 210 211 IF(lwp) WRITE(numout,*) '~~~~~~~~~~ and constant salinity (',zsal,' psu)' 212 IF(lwp .AND. lflush) CALL flush(numout) 211 213 ! 212 214 DO jk = 1, jpk … … 359 361 ! 360 362 END SELECT 363 ! 364 IF(lwp .AND. lflush) CALL flush(numout) 361 365 ! 362 366 END SUBROUTINE istate_eel … … 443 447 ENDIF 444 448 ! 449 IF(lwp .AND. lflush) CALL flush(numout) 450 ! 445 451 END SUBROUTINE istate_gyre 446 452 … … 566 572 CALL wrk_dealloc( jpi, jpj, jpk, zprn) 567 573 ! 574 IF(lwp .AND. lflush) CALL flush(numout) 575 ! 568 576 END SUBROUTINE istate_uvg 569 577 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90
r6488 r11101 207 207 WRITE(numout,*) 208 208 WRITE(numout,*) ' smallest real computer value rsmall = ', rsmall 209 IF(lflush) CALL flush(numout) 209 210 ENDIF 210 211 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90
r6487 r11101 104 104 IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and relative vorticity' 105 105 IF(lwp) WRITE(numout,*) '~~~~~~~ NOT optimal for auto-tasking case' 106 IF(lwp .AND. lflush) CALL flush(numout) 106 107 ENDIF 107 108 … … 285 286 IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and' 286 287 IF(lwp) WRITE(numout,*) '~~~~~~~ relative vorticity' 288 IF(lwp .AND. lflush) CALL flush(numout) 287 289 ENDIF 288 290 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90
r6486 r11101 108 108 READ ( numnam_cfg, namdyn_adv, IOSTAT = ios, ERR = 902 ) 109 109 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_adv in configuration namelist', lwp ) 110 IF(lwm ) WRITE ( numond, namdyn_adv )110 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_adv ) 111 111 112 112 IF(lwp) THEN ! Namelist print … … 152 152 IF( nadv == 3 ) WRITE(numout,*) ' flux form : UBS scheme is used' 153 153 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection formulation' 154 IF(lflush) CALL flush(numout) 154 155 ENDIF 155 156 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r6486 r11101 66 66 WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection' 67 67 WRITE(numout,*) '~~~~~~~~~~~~' 68 IF(lflush) CALL flush(numout) 68 69 ENDIF 69 70 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
r6486 r11101 90 90 IF(lwp) WRITE(numout,*) 'dyn_adv_ubs : UBS flux form momentum advection' 91 91 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 92 IF(lwp .AND. lflush) CALL flush(numout) 92 93 ENDIF 93 94 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r8447 r11101 136 136 IF(lwp) THEN 137 137 WRITE(numout,*) " ! restore original density" 138 IF(lflush) CALL flush(numout) 138 139 ENDIF 139 140 rhd(:,:,:) = z_rhd_st(:,:,:) ! restore original density … … 176 177 READ ( numnam_cfg, namdyn_hpg, IOSTAT = ios, ERR = 902 ) 177 178 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_hpg in configuration namelist', lwp ) 178 IF(lwm ) WRITE ( numond, namdyn_hpg )179 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_hpg ) 179 180 ! 180 181 IF(lwp) THEN ! Control print … … 190 191 WRITE(numout,*) ' s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = ', ln_hpg_prj 191 192 WRITE(numout,*) ' time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = ', ln_dynhpg_imp 193 IF(lflush) CALL flush(numout) 192 194 ENDIF 193 195 ! … … 260 262 IF(lwp) WRITE(numout,*) 'dyn:hpg_zco : hydrostatic pressure gradient trend' 261 263 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ z-coordinate case ' 264 IF(lwp .AND. lflush) CALL flush(numout) 262 265 ENDIF 263 266 … … 325 328 IF(lwp) WRITE(numout,*) 'dyn:hpg_zps : hydrostatic pressure gradient trend' 326 329 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ z-coordinate with partial steps - vector optimization' 330 IF(lwp .AND. lflush) CALL flush(numout) 327 331 ENDIF 328 332 … … 423 427 IF(lwp) WRITE(numout,*) 'dyn:hpg_sco : hydrostatic pressure gradient trend' 424 428 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, OPA original scheme used' 429 IF(lwp .AND. lflush) CALL flush(numout) 425 430 ENDIF 426 431 … … 514 519 IF(lwp) WRITE(numout,*) 'dyn:hpg_isf : hydrostatic pressure gradient trend for ice shelf' 515 520 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, OPA original scheme used' 521 IF(lwp .AND. lflush) CALL flush(numout) 516 522 ENDIF 517 523 … … 768 774 IF(lwp) WRITE(numout,*) 'dyn:hpg_djc : hydrostatic pressure gradient trend' 769 775 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, density Jacobian with cubic polynomial scheme' 776 IF(lwp .AND. lflush) CALL flush(numout) 770 777 ENDIF 771 778 … … 1001 1008 IF(lwp) WRITE(numout,*) 'dyn:hpg_prj : hydrostatic pressure gradient trend' 1002 1009 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, cubic spline pressure Jacobian' 1010 IF(lwp .AND. lflush) CALL flush(numout) 1003 1011 ENDIF 1004 1012 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90
r6486 r11101 88 88 IF(lwp) WRITE(numout,*) 'dyn_keg : kinetic energy gradient trend, scheme number=', kscheme 89 89 IF(lwp) WRITE(numout,*) '~~~~~~~' 90 IF(lflush) CALL flush(numout) 90 91 ENDIF 91 92 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90
r6486 r11101 100 100 IF(lwp) WRITE(numout,*) 'dyn_ldf : no lateral diffusion on momentum setup' 101 101 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 102 IF(lflush) CALL flush(numout) 102 103 ENDIF 103 104 END SELECT … … 139 140 WRITE(numout,*) ' horizontal (geopotential) ln_dynldf_hor = ', ln_dynldf_hor 140 141 WRITE(numout,*) ' iso-neutral ln_dynldf_iso = ', ln_dynldf_iso 142 IF(lflush) CALL flush(numout) 141 143 ENDIF 142 144 … … 229 231 IF( nldf == 4 ) WRITE(numout,*) ' laplacian and bilaplacian operators' 230 232 IF( nldf == 5 ) WRITE(numout,*) ' rotated laplacian and bilaplacian operators' 233 IF(lflush) CALL flush(numout) 231 234 ENDIF 232 235 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r6486 r11101 90 90 WRITE(numout,*) 'dyn_ldf_bilap : iso-level bilaplacian operator' 91 91 WRITE(numout,*) '~~~~~~~~~~~~~' 92 IF(lwp .AND. lflush) CALL flush(numout) 92 93 ENDIF 93 94 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90
r8280 r11101 98 98 IF(lwp) WRITE(numout,*) 'dyn_ldf_bilapg : horizontal biharmonic operator in s-coordinate' 99 99 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 100 IF(lwp .AND. lflush) CALL flush(numout) 100 101 ! ! allocate dyn_ldf_bilapg arrays 101 102 IF( dyn_ldf_bilapg_alloc() /= 0 ) CALL ctl_stop('STOP', 'dyn_ldf_bilapg: failed to allocate arrays') … … 123 124 IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)), & 124 125 & ' wi', sqrt(MAXVAL(wslpi)) , ' wj', sqrt(MAXVAL(wslpj)) 126 IF(lwp .AND. lflush) CALL flush(numout) 125 127 endif 126 128 !!end -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_iso.F90
r6486 r11101 127 127 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate horizontal diffusive operator' 128 128 ! ! allocate dyn_ldf_bilap arrays 129 IF(lflush) CALL flush(numout) 129 130 IF( dyn_ldf_iso_alloc() /= 0 ) CALL ctl_stop('STOP', 'dyn_ldf_iso: failed to allocate arrays') 130 131 ENDIF … … 151 152 IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)), & 152 153 & ' wi', sqrt(MAXVAL(wslpi)) , ' wj', sqrt(MAXVAL(wslpj)) 154 IF(lwp .AND. lflush) CALL flush(numout) 153 155 endif 154 156 !!end -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90
r6486 r11101 70 70 IF(lwp) WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator' 71 71 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 72 IF(lflush) CALL flush(numout) 72 73 ENDIF 73 74 ! ! =============== -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90
r6486 r11101 126 126 READ ( numnam_cfg, namdyn_nept, IOSTAT = ios, ERR = 902 ) 127 127 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_nept in configuration namelist', lwp ) 128 IF(lwm ) WRITE ( numond, namdyn_nept )128 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_nept ) 129 129 130 130 IF(lwp) THEN ! Control print … … 145 145 WRITE(numout,*) 146 146 ENDIF 147 IF(lflush) CALL flush(numout) 147 148 ENDIF 148 149 ! … … 304 305 CALL lbc_lnk( zunep, 'U', -1.0_wp ) 305 306 CALL lbc_lnk( zvnep, 'V', -1.0_wp ) 306 WRITE(numout,*) ' zunep: min, max = ', unemin,unemax 307 WRITE(numout,*) ' zvnep: min, max = ', vnemin,vnemax 308 WRITE(numout,*) 307 IF(lwp .AND. nprint > 0) THEN 308 WRITE(numout,*) ' zunep: min, max = ', unemin,unemax 309 WRITE(numout,*) ' zvnep: min, max = ', vnemin,vnemax 310 WRITE(numout,*) 311 ENDIF 309 312 310 313 !! Compute, once and for all, the horizontal divergence (zhdivnep) … … 328 331 END DO 329 332 END DO 330 WRITE(numout,*) ' zhdivnep: min, max = ', zhdivmin,zhdivmax 331 WRITE(numout,*) ' zmrotnep: min, max = ', zmrotmin,zmrotmax 332 WRITE(numout,*) 333 334 IF(lwp .AND. nprint > 0) THEN 335 WRITE(numout,*) ' zhdivnep: min, max = ', zhdivmin,zhdivmax 336 WRITE(numout,*) ' zmrotnep: min, max = ', zmrotmin,zmrotmax 337 WRITE(numout,*) 338 IF(lflush) CALL flush(numout) 339 ENDIF 333 340 334 341 !! Deallocate temporary workspace arrays, which are all local to … … 384 391 IF(lwp) WRITE(numout,*) 'horizontal velocity divergence and' 385 392 IF(lwp) WRITE(numout,*) 'relative vorticity of Neptune flow' 393 IF(lwp .AND. lflush) CALL flush(numout) 386 394 #if defined key_noslip_accurate 387 395 !!---------------------------------------------------------------------- -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90
r6487 r11101 119 119 IF(lwp) WRITE(numout,*) 'dyn_nxt : time stepping' 120 120 IF(lwp) WRITE(numout,*) '~~~~~~~' 121 IF(lflush) CALL flush(numout) 121 122 ENDIF 122 123 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90
r6486 r11101 231 231 WRITE(numout,*) ' Free surface with time splitting lk_dynspg_ts = ', lk_dynspg_ts 232 232 WRITE(numout,*) ' Filtered free surface cst volume lk_dynspg_flt = ', lk_dynspg_flt 233 IF(lflush) CALL flush(numout) 233 234 ENDIF 234 235 … … 266 267 IF( nspg == 1 ) WRITE(numout,*) ' free surface with time splitting scheme' 267 268 IF( nspg == 2 ) WRITE(numout,*) ' filtered free surface' 269 IF(lflush) CALL flush(numout) 268 270 ENDIF 269 271 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_exp.F90
r9583 r11101 75 75 ! 76 76 IF( lk_vvl .AND. lwp ) WRITE(numout,*) ' lk_vvl=T : spg is included in dynhpg' 77 IF(lflush) CALL flush(numout) 77 78 ENDIF 78 79 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90
r9321 r11101 120 120 IF(lwp) WRITE(numout,*) 'dyn_spg_flt : surface pressure gradient trend' 121 121 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ (free surface constant volume case)' 122 IF(lwp .AND. lflush) CALL flush(numout) 122 123 123 124 ! set to zero free surface specific arrays -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r9321 r11101 195 195 IF(lwp) WRITE(numout,*) '~~~~~~~~~~ free surface with time splitting' 196 196 IF(lwp) WRITE(numout,*) 197 IF(lflush) CALL flush(numout) 197 198 ! 198 199 IF (neuler==0) ll_init=.TRUE. … … 1088 1089 READ ( numnam_cfg, namsplit, IOSTAT = ios, ERR = 902 ) 1089 1090 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsplit in configuration namelist', lwp ) 1090 IF(lwm ) WRITE ( numond, namsplit )1091 IF(lwm .AND. nprint > 2) WRITE ( numond, namsplit ) 1091 1092 ! 1092 1093 ! ! Max courant number for ext. grav. waves … … 1161 1162 IF(lwp) WRITE(numout,*) ' Barotropic time step [s] is :', rdtbt 1162 1163 IF(lwp) WRITE(numout,*) ' Maximum Courant number is :', zcmax 1164 IF(lwp .AND. lflush) CALL flush(numout) 1163 1165 ! 1164 1166 IF ((.NOT.ln_bt_av).AND.(.NOT.ln_bt_fw)) THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynvor.F90
r8280 r11101 234 234 IF(lwp) WRITE(numout,*) 'dyn:vor_ene : vorticity term: energy conserving scheme' 235 235 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 236 IF(lflush) CALL flush(numout) 236 237 ENDIF 237 238 … … 346 347 IF(lwp) WRITE(numout,*) 'dyn:vor_mix : vorticity term: mixed energy/enstrophy conserving scheme' 347 348 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 349 IF(lflush) CALL flush(numout) 348 350 ENDIF 349 351 … … 463 465 IF(lwp) WRITE(numout,*) 'dyn:vor_ens : vorticity term: enstrophy conserving scheme' 464 466 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 467 IF(lflush) CALL flush(numout) 465 468 ENDIF 466 469 … … 587 590 IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 588 591 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 592 IF(lflush) CALL flush(numout) 589 593 #if ! defined key_vvl 590 594 IF( .NOT.ALLOCATED(ze3f) ) THEN … … 731 735 READ ( numnam_cfg, namdyn_vor, IOSTAT = ios, ERR = 902 ) 732 736 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_vor in configuration namelist', lwp ) 733 IF(lwm ) WRITE ( numond, namdyn_vor )737 IF(lwm .AND. nprint >2) WRITE ( numond, namdyn_vor ) 734 738 735 739 IF(lwp) THEN ! Namelist print … … 743 747 WRITE(numout,*) ' enstrophy and energy conserving scheme ln_dynvor_een = ', ln_dynvor_een 744 748 WRITE(numout,*) ' enstrophy and energy conserving scheme (old) ln_dynvor_een_old= ', ln_dynvor_een_old 749 IF(lflush) CALL flush(numout) 745 750 ENDIF 746 751 … … 798 803 IF( nvor == 3 ) WRITE(numout,*) ' vorticity scheme : energy and enstrophy conserving scheme' 799 804 IF( nvor == -1 ) WRITE(numout,*) ' esopa test: use all lateral physics options' 805 IF(lflush) CALL flush(numout) 800 806 ENDIF 801 807 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90
r6486 r11101 74 74 IF(lwp)WRITE(numout,*) 75 75 IF(lwp)WRITE(numout,*) 'dyn_zad : arakawa advection scheme' 76 IF(lwp .AND. lflush) CALL flush(numout) 76 77 ENDIF 77 78 … … 188 189 IF(lwp)WRITE(numout,*) 189 190 IF(lwp)WRITE(numout,*) 'dyn_zad_zts : arakawa advection scheme with sub-timesteps' 191 IF(lwp .AND. lflush) CALL flush(numout) 190 192 ENDIF 191 193 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90
r6486 r11101 136 136 IF( nzdf == 0 ) WRITE(numout,*) ' Explicit time-splitting scheme' 137 137 IF( nzdf == 1 ) WRITE(numout,*) ' Implicit (euler backward) scheme' 138 IF(lflush) CALL flush(numout) 138 139 ENDIF 139 140 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_exp.F90
r6486 r11101 73 73 WRITE(numout,*) 'dyn_zdf_exp : vertical momentum diffusion - explicit operator' 74 74 WRITE(numout,*) '~~~~~~~~~~~ ' 75 IF(lflush) CALL flush(numout) 75 76 ENDIF 76 77 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
r6795 r11101 81 81 IF(lwp) WRITE(numout,*) 'dyn_zdf_imp : vertical momentum diffusion implicit operator' 82 82 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 83 IF(lwp .AND. lflush) CALL flush(numout) 83 84 ! 84 85 IF( lk_vvl ) THEN ; r_vvl = 1._wp ! Variable volume indicator -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r8400 r11101 87 87 IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height' 88 88 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 89 IF(lwp .AND. lflush) CALL flush(numout) 89 90 ! 90 91 ENDIF … … 184 185 IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity ' 185 186 IF(lwp) WRITE(numout,*) '~~~~~ ' 187 IF(lwp .AND. lflush) CALL flush(numout) 186 188 ! 187 189 wn(:,:,jpk) = 0._wp ! bottom boundary condition: w=0 (set once for all) … … 265 267 IF(lwp) WRITE(numout,*) 'ssh_swp : Asselin time filter and swap of sea surface height' 266 268 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 269 IF(lwp .AND. lflush) CALL flush(numout) 267 270 ENDIF 268 271 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/FLO/floats.F90
r9583 r11101 96 96 READ ( numnam_cfg, namflo, IOSTAT = ios, ERR = 902 ) 97 97 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namflo in configuration namelist', lwp ) 98 IF(lwm ) WRITE ( numond, namflo )98 IF(lwm .AND. nprint > 2) WRITE ( numond, namflo ) 99 99 ! 100 100 IF(lwp) THEN ! control print -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90
r6486 r11101 65 65 IF(lwp) WRITE(numout,*) '~~~~~~~~' 66 66 IF(lwp) WRITE(numout,*) ' jpnfl = ',jpnfl 67 IF(lwp .AND. lflush) CALL flush(numout) 67 68 68 69 !-------------------------! … … 72 73 73 74 IF(lwp) WRITE(numout,*) ' float restart file read' 75 IF(lwp .AND. lflush) CALL flush(numout) 74 76 75 77 ! open the restart file … … 93 95 94 96 IF(lwp) WRITE(numout,*) ' add new floats' 97 IF(lwp .AND. lflush) CALL flush(numout) 95 98 96 99 IF( ln_ariane )THEN !Add new floats with ariane convention … … 148 151 DO jfl = kfl_start,kfl_end 149 152 READ(inum,*) flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash 150 if(lwp)write(numout,*)'read:',jfl,flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash ; call flush(numout) 153 IF (lwp) THEN 154 write(numout,*)'read:',jfl,flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash 155 IF(lflush) CALL flush(numout) 156 ENDIF 151 157 END DO 152 158 CLOSE(inum) … … 275 281 IF(lwp) WRITE(numout,*)'geographical position',flxx(jfl),flyy(jfl),flzz(jfl) 276 282 IF(lwp) WRITE(numout,*)'index position',tpifl(jfl),tpjfl(jfl),tpkfl(jfl) 283 IF(lwp .AND. lflush) CALL flush(numout) 277 284 ENDIF 278 285 END DO -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90
r8046 r11101 56 56 INTEGER :: imx ! temporary integer for max berg class 57 57 LOGICAL, SAVE :: ll_first_call = .TRUE. 58 LOGICAL :: ll_write 58 59 !!---------------------------------------------------------------------- 59 60 ! … … 68 69 ! nn_coupled_iceshelf_fluxes uninitialised unless lk_oasis=true 69 70 IF( nn_coupled_iceshelf_fluxes .gt. 0 ) THEN 70 71 ll_write = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) .AND. lwp 71 72 ! Adjust total calving rates so that sum of iceberg calving and iceshelf melting in the northern 72 73 ! and southern hemispheres equals rate of increase of mass of greenland and antarctic ice sheets … … 80 81 81 82 ! check 82 IF(l wp) WRITE(numout, *) 'Greenland iceberg calving climatology (kg/s) : ',zgreenland_calving_sum83 IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving climatology (kg/s) : ',zgreenland_calving_sum 83 84 zgreenland_calving_sum = SUM( berg_grid%calving(:,:) * greenland_icesheet_mask(:,:) ) 84 85 IF( lk_mpp ) CALL mpp_sum( zgreenland_calving_sum ) 85 IF(l wp) WRITE(numout, *) 'Greenland iceberg calving adjusted value (kg/s) : ',zgreenland_calving_sum86 IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving adjusted value (kg/s) : ',zgreenland_calving_sum 86 87 87 88 zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) ) … … 92 93 93 94 ! check 94 IF(l wp) WRITE(numout, *) 'Antarctica iceberg calving climatology (kg/s) : ',zantarctica_calving_sum95 IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving climatology (kg/s) : ',zantarctica_calving_sum 95 96 zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) ) 96 97 IF( lk_mpp ) CALL mpp_sum( zantarctica_calving_sum ) 97 IF(l wp) WRITE(numout, *) 'Antarctica iceberg calving adjusted value (kg/s) : ',zantarctica_calving_sum98 IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving adjusted value (kg/s) : ',zantarctica_calving_sum 98 99 99 100 ENDIF … … 218 219 CALL lbc_lnk( berg_grid%stored_heat, 'T', 1._wp ) 219 220 ! 220 IF( nn_verbose_level > 0 .AND. icntmax > 1 ) WRITE(numicb,*) 'icb_clv: icnt=', icnt,' on', narea221 IF( nn_verbose_level > 0 .AND. icntmax > 1 .AND. numicb.NE.-1) WRITE(numicb,*) 'icb_clv: icnt=', icnt,' on', narea 221 222 ! 222 223 END SUBROUTINE icb_clv -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbdia.F90
r6491 r11101 311 311 'returned',calving_ret_net) 312 312 ENDIF 313 WRITE( numicb, '("calved by class = ",i6,20(",",i6))') (nbergs_calved_by_class(ik),ik=1,nclasses) 314 IF ( nspeeding_tickets > 0 ) WRITE( numicb, '("speeding tickets issued = ",i6)') nspeeding_tickets 313 IF(numicb.NE.-1) THEN 314 WRITE( numicb, '("calved by class = ",i6,20(",",i6))') (nbergs_calved_by_class(ik),ik=1,nclasses) 315 IF ( nspeeding_tickets > 0 ) WRITE( numicb, '("speeding tickets issued = ",i6)') nspeeding_tickets 316 ENDIF 315 317 316 318 nbergs_start = nbergs_end … … 418 420 stored_start = SUM( berg_grid%stored_ice(:,:,:) ) 419 421 IF( lk_mpp ) CALL mpp_sum( stored_start ) 420 WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored mass=',stored_start,' kg'422 IF(numicb.NE.-1) WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored mass=',stored_start,' kg' 421 423 ! 422 424 stored_heat_start = SUM( berg_grid%stored_heat(:,:) ) 423 425 IF( lk_mpp ) CALL mpp_sum( stored_heat_start ) 424 WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored heat=',stored_heat_start,' J'426 IF(numicb.NE.-1) WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored heat=',stored_heat_start,' J' 425 427 ENDIF 426 428 ! … … 495 497 ! 496 498 IF ( PRESENT(kbergs) ) THEN 497 WRITE(numicb,100) cd_budgetstr // ' state:', & 499 IF(numicb.NE.-1) & 500 WRITE(numicb,100) cd_budgetstr // ' state:', & 498 501 cd_startstr // ' start', pstartval, cd_budgetunits, & 499 502 cd_endstr // ' end', pendval, cd_budgetunits, & … … 501 504 '# of bergs', kbergs 502 505 ELSE 503 WRITE(numicb,100) cd_budgetstr // ' state:', & 506 IF(numicb.NE.-1) & 507 WRITE(numicb,100) cd_budgetstr // ' state:', & 504 508 cd_startstr // ' start', pstartval, cd_budgetunits, & 505 509 cd_endstr // ' end', pendval, cd_budgetunits, & … … 517 521 !!---------------------------------------------------------------------- 518 522 ! 519 WRITE(numicb,200) cd_budgetstr // ' check:', & 523 IF(numicb.NE.-1) & 524 WRITE(numicb,200) cd_budgetstr // ' check:', & 520 525 cd_startstr, pstartval, cd_budgetunits, & 521 526 cd_endstr, pendval, cd_budgetunits, & … … 538 543 & MAX( 1.e-30, MAX( abs( pendval - pstartval ) , ABS( pinval - poutval ) ) ) 539 544 540 WRITE(numicb,200) cd_budgetstr // ' budget:', & 545 IF(numicb.NE.-1) & 546 & WRITE(numicb,200) cd_budgetstr // ' budget:', & 541 547 & cd_instr // ' in', pinval, cd_budgetunits, & 542 548 & cd_outstr // ' out', poutval, cd_budgetunits, & … … 554 560 INTEGER, INTENT(in) :: pstartval, pendval 555 561 ! 556 WRITE(numicb,100) cd_budgetstr // ' state:', & 562 IF(numicb.NE.-1) & 563 & WRITE(numicb,100) cd_budgetstr // ' state:', & 557 564 & cd_startstr // ' start', pstartval, & 558 565 & cd_endstr // ' end', pendval, & … … 571 578 !!---------------------------------------------------------------------- 572 579 ! 573 WRITE(numicb,200) cd_budgetstr // ' budget:', & 580 IF(numicb.NE.-1) & 581 & WRITE(numicb,200) cd_budgetstr // ' budget:', & 574 582 cd_instr // ' in', pinval, & 575 583 cd_outstr // ' out', poutval, & -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbdyn.F90
r6486 r11101 371 371 ENDIF 372 372 ! ! check the speed and acceleration limits 373 IF( ABS( zuveln ) > pp_vel_lim .OR. ABS( zvveln ) > pp_vel_lim ) & 374 WRITE(numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive velocity' 375 IF( ABS( pax ) > pp_accel_lim .OR. ABS( pay ) > pp_accel_lim ) & 376 WRITE(numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive acceleration' 373 IF( ABS( zuveln ) > pp_vel_lim .OR. ABS( zvveln ) > pp_vel_lim ) THEN 374 IF(numicb.NE.-1) THEN 375 WRITE(numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive velocity' 376 CALL flush(numicb) 377 ELSE 378 WRITE(*,'("pe=",i3,x,a)') narea,'Dump triggered by excessive velocity' 379 CALL flush(6) 380 ENDIF 381 ENDIF 382 IF( ABS( pax ) > pp_accel_lim .OR. ABS( pay ) > pp_accel_lim ) THEN 383 IF(numicb.NE.-1) THEN 384 WRITE(numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive acceleration' 385 CALL flush(numicb) 386 ELSE 387 WRITE(*,'("pe=",i3,x,a)') narea,'Dump triggered by excessive acceleration' 388 CALL flush(6) 389 ENDIF 390 ENDIF 377 391 ! 378 392 END SUBROUTINE icb_accel -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90
r6823 r11101 78 78 ! ! note that we choose to do this on all processors since we cannot 79 79 ! ! predict where icebergs will be ahead of time 80 CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 80 numicb=-1 81 IF(nn_verbose_level>0 .AND. nprint>0) THEN 82 CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 83 ELSE 84 IF(lwp) CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 85 ENDIF 81 86 82 87 ! set parameters (mostly from namelist) … … 187 192 ENDIF 188 193 ! 189 IF( nn_verbose_level > 0 ) THEN194 IF( nn_verbose_level > 0 .AND. numicb.NE.-1) THEN 190 195 WRITE(numicb,*) 'processor ', narea 191 196 WRITE(numicb,*) 'jpi, jpj ', jpi, jpj … … 241 246 CALL iom_close( inum ) ! close file 242 247 ! 243 WRITE(numicb,*) 244 WRITE(numicb,*) ' calving read in a file' 248 IF(numicb.NE.-1) THEN 249 WRITE(numicb,*) 250 WRITE(numicb,*) ' calving read in a file' 251 ENDIF 245 252 ALLOCATE( sf_icb(1), STAT=istat1 ) ! Create sf_icb structure (calving) 246 253 ALLOCATE( sf_icb(1)%fnow(jpi,jpj,1), STAT=istat2 ) … … 336 343 ibergs = icb_utl_count() 337 344 IF( lk_mpp ) CALL mpp_sum(ibergs) 338 WRITE(numicb,'(a,i6,a)') 'diamonds, icb_ini_gen: ',ibergs,' were generated'345 IF(numicb.NE.-1) WRITE(numicb,'(a,i6,a)') 'diamonds, icb_ini_gen: ',ibergs,' were generated' 339 346 ! 340 347 END SUBROUTINE icb_ini_gen … … 368 375 READ ( numnam_cfg, namberg, IOSTAT = ios, ERR = 902 ) 369 376 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namberg in configuration namelist', lwp ) 370 IF(lwm ) WRITE ( numond, namberg )377 IF(lwm .AND. nprint > 2) WRITE ( numond, namberg ) 371 378 #else 372 379 IF(lwp) THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icblbc.F90
r6486 r11101 197 197 ENDIF 198 198 199 IF( nn_verbose_level >= 2 ) THEN199 IF( nn_verbose_level >= 2 .AND. numicb .NE. -1) THEN 200 200 WRITE(numicb,*) 'processor west : ', ipe_W 201 201 WRITE(numicb,*) 'processor east : ', ipe_E … … 240 240 this => this%next 241 241 ibergs_to_send_e = ibergs_to_send_e + 1 242 IF( nn_verbose_level >= 4 ) THEN242 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 243 243 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to east' 244 244 CALL flush( numicb ) … … 253 253 this => this%next 254 254 ibergs_to_send_w = ibergs_to_send_w + 1 255 IF( nn_verbose_level >= 4 ) THEN255 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 256 256 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to west' 257 257 CALL flush( numicb ) … … 267 267 END DO 268 268 ENDIF 269 IF( nn_verbose_level >= 3 ) THEN269 IF( nn_verbose_level >= 3 .AND. numicb .NE. -1) THEN 270 270 WRITE(numicb,*) 'bergstep ',nktberg,' send ew: ', ibergs_to_send_e, ibergs_to_send_w 271 271 CALL flush(numicb) … … 301 301 ibergs_rcvd_from_w = INT( zwebergs(2) ) 302 302 END SELECT 303 IF( nn_verbose_level >= 3 ) THEN303 IF( nn_verbose_level >= 3 .AND. numicb .NE. -1) THEN 304 304 WRITE(numicb,*) 'bergstep ',nktberg,' recv ew: ', ibergs_rcvd_from_w, ibergs_rcvd_from_e 305 305 CALL flush(numicb) … … 315 315 IF( ibergs_to_send_e > 0 .AND. l_isend ) CALL mpi_wait( iml_req1, iml_stat, iml_err ) 316 316 DO i = 1, ibergs_rcvd_from_e 317 IF( nn_verbose_level >= 4 ) THEN317 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 318 318 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_e%data(16,i)),' from east' 319 319 CALL flush( numicb ) … … 335 335 IF( ibergs_to_send_e > 0 .AND. l_isend ) CALL mpi_wait( iml_req3, iml_stat, iml_err ) 336 336 DO i = 1, ibergs_rcvd_from_e 337 IF( nn_verbose_level >= 4 ) THEN337 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 338 338 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_e%data(16,i)),' from east' 339 339 CALL flush( numicb ) … … 342 342 END DO 343 343 DO i = 1, ibergs_rcvd_from_w 344 IF( nn_verbose_level >= 4 ) THEN344 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 345 345 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_w%data(16,i)),' from west' 346 346 CALL flush( numicb ) … … 356 356 IF( ibergs_to_send_w > 0 .AND. l_isend ) CALL mpi_wait( iml_req4, iml_stat, iml_err ) 357 357 DO i = 1, ibergs_rcvd_from_w 358 IF( nn_verbose_level >= 4 ) THEN358 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 359 359 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_w%data(16,i)),' from west' 360 360 CALL flush( numicb ) … … 378 378 this => this%next 379 379 ibergs_to_send_n = ibergs_to_send_n + 1 380 IF( nn_verbose_level >= 4 ) THEN380 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 381 381 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to north' 382 382 CALL flush( numicb ) … … 388 388 this => this%next 389 389 ibergs_to_send_s = ibergs_to_send_s + 1 390 IF( nn_verbose_level >= 4 ) THEN390 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 391 391 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to south' 392 392 CALL flush( numicb ) … … 399 399 END DO 400 400 ENDIF 401 if( nn_verbose_level >= 3 ) then401 if( nn_verbose_level >= 3 .AND. numicb .NE. -1) then 402 402 write(numicb,*) 'bergstep ',nktberg,' send ns: ', ibergs_to_send_n, ibergs_to_send_s 403 403 call flush(numicb) … … 432 432 ibergs_rcvd_from_s = INT( zsnbergs(2) ) 433 433 END SELECT 434 if( nn_verbose_level >= 3 ) then434 if( nn_verbose_level >= 3 .AND. numicb .NE. -1) then 435 435 write(numicb,*) 'bergstep ',nktberg,' recv ns: ', ibergs_rcvd_from_s, ibergs_rcvd_from_n 436 436 call flush(numicb) … … 446 446 IF( ibergs_to_send_n > 0 .AND. l_isend ) CALL mpi_wait( iml_req1, iml_stat, iml_err ) 447 447 DO i = 1, ibergs_rcvd_from_n 448 IF( nn_verbose_level >= 4 ) THEN448 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 449 449 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_n%data(16,i)),' from north' 450 450 CALL flush( numicb ) … … 466 466 IF( ibergs_to_send_n > 0 .AND. l_isend ) CALL mpi_wait( iml_req3, iml_stat, iml_err ) 467 467 DO i = 1, ibergs_rcvd_from_n 468 IF( nn_verbose_level >= 4 ) THEN468 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 469 469 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_n%data(16,i)),' from north' 470 470 CALL flush( numicb ) … … 473 473 END DO 474 474 DO i = 1, ibergs_rcvd_from_s 475 IF( nn_verbose_level >= 4 ) THEN475 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 476 476 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_s%data(16,i)),' from south' 477 477 CALL flush( numicb ) … … 487 487 IF( ibergs_to_send_s > 0 .AND. l_isend ) CALL mpi_wait( iml_req4, iml_stat, iml_err ) 488 488 DO i = 1, ibergs_rcvd_from_s 489 IF( nn_verbose_level >= 4 ) THEN489 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 490 490 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_s%data(16,i)),' from south' 491 491 CALL flush( numicb ) … … 495 495 END SELECT 496 496 497 IF( nn_verbose_level > 0 ) THEN497 IF( nn_verbose_level > 0 .AND. numicb .NE. -1) THEN 498 498 ! compare the number of icebergs on this processor from the start to the end 499 499 ibergs_end = icb_utl_count() … … 535 535 IF( npolj /= 0 .AND. jpni > 1 ) CALL icb_lbc_mpp_nfld( ) 536 536 537 IF( nn_verbose_level > 0 537 IF( nn_verbose_level > 0) THEN 538 538 i = 0 539 539 this => first_berg … … 547 547 ijne .GT. mjg(nicbej)) THEN 548 548 i = i + 1 549 WRITE(numicb,*) 'berg lost in halo: ', this%number(:),iine,ijne 550 WRITE(numicb,*) ' ', nimpp, njmpp 551 WRITE(numicb,*) ' ', nicbdi, nicbei, nicbdj, nicbej 552 CALL flush( numicb ) 549 IF(numicb .NE. -1) THEN 550 WRITE(numicb,*) 'berg lost in halo: ', this%number(:),iine,ijne 551 WRITE(numicb,*) ' ', nimpp, njmpp 552 WRITE(numicb,*) ' ', nicbdi, nicbei, nicbdj, nicbej 553 CALL flush( numicb ) 554 ENDIF 553 555 ENDIF 554 556 this => this%next … … 659 661 IF( nicbfldproc(jjn) .eq. INT(znbergs(1)) ) EXIT 660 662 END DO 661 IF( jjn .GT. jpni ) write( numicb,*) 'ICB ERROR'663 IF( jjn .GT. jpni ) write(*,*) 'ICB ERROR' 662 664 nicbfldexpect(jjn) = INT( znbergs(2) ) 663 665 !IF ( nicbfldexpect(jjn) .GT. 0) write(numicb,*) 'ICB expecting ',nicbfldexpect(jjn),' from ', nicbfldproc(jjn) … … 714 716 tmpberg => this 715 717 ibergs_to_send = ibergs_to_send + 1 716 IF( nn_verbose_level >= 4 ) THEN718 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 717 719 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for north fold' 718 720 CALL flush( numicb ) … … 727 729 END DO 728 730 ENDIF 729 if( nn_verbose_level >= 3 ) then731 if( nn_verbose_level >= 3 .AND. numicb.NE.-1) then 730 732 write(numicb,*) 'bergstep ',nktberg,' send nfld: ', ibergs_to_send 731 733 call flush(numicb) … … 758 760 ! 759 761 DO jk = 1, ibergs_to_rcv 760 IF( nn_verbose_level >= 4 ) THEN762 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 761 763 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_f%data(16,jk)),' from north fold' 762 764 CALL flush( numicb ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbrst.F90
r9321 r11101 134 134 ! Sanity check 135 135 jn = icb_utl_count() 136 IF (n n_verbose_level >= 0) &136 IF (nprint > 0 .AND. lwp) & 137 137 WRITE(numout,'(2(a,i5))') 'icebergs, read_restart_bergs: # bergs =',jn,' on PE',narea-1 138 138 IF( lk_mpp ) THEN … … 141 141 CALL mpp_sum(jn) 142 142 ENDIF 143 IF(lwp) WRITE(numout,'(a,i5,a,i5,a)') 'icebergs, read_restart_bergs: there were',ibergs_in_file, & 143 IF(lwp .AND. nprint > 1) & 144 & WRITE(numout,'(a,i5,a,i5,a)') 'icebergs, read_restart_bergs: there were',ibergs_in_file, & 144 145 & ' bergs in the restart file and', jn,' bergs have been read' 145 146 ! … … 148 149 ! 149 150 IF(nn_timing == 2) CALL timing_stop('iom_rstget') 150 IF( lwp .and. n n_verbose_level>= 0) WRITE(numout,'(a)') 'icebergs, read_restart_bergs: completed'151 IF( lwp .and. nprint >= 0) WRITE(numout,'(a)') 'icebergs, read_restart_bergs: completed' 151 152 ! 152 153 END SUBROUTINE icb_rst_read … … 190 191 WRITE(cl_filename,'(A,"_icebergs_",A,"_restart.nc")') TRIM(cexper), TRIM(ADJUSTL(clkt)) 191 192 ENDIF 192 IF (nn_verbose_level >= 0) WRITE(numout,'(2a)') 'icebergs, write_restart: creating ',TRIM(cl_path)//TRIM(cl_filename) 193 IF (nprint > 0 .AND. lwp) & 194 WRITE(numout,'(2a)') 'icebergs, write_restart: creating ',TRIM(cl_path)//TRIM(cl_filename) 193 195 194 196 nret = NF90_CREATE(TRIM(cl_path)//TRIM(cl_filename), NF90_CLOBBER, ncid) … … 324 326 ENDIF 325 327 ENDDO 326 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_ice written'328 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_ice written' 327 329 328 330 nret = NF90_PUT_VAR( ncid, nkountid, num_bergs(:) ) … … 331 333 nret = NF90_PUT_VAR( ncid, nsheatid, berg_grid%stored_heat(:,:) ) 332 334 IF (nret .ne. NF90_NOERR) CALL ctl_stop('icebergs, write_restart: nf_put_var stored_heat failed') 333 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_heat written'335 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_heat written' 334 336 335 337 nret = NF90_PUT_VAR( ncid, ncalvid , src_calving(:,:) ) … … 337 339 nret = NF90_PUT_VAR( ncid, ncalvhid, src_calving_hflx(:,:) ) 338 340 IF (nret .ne. NF90_NOERR) CALL ctl_stop('icebergs, write_restart: nf_put_var calving_hflx failed') 339 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: calving written'341 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: calving written' 340 342 341 343 IF ( ASSOCIATED(first_berg) ) THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbstp.F90
r6486 r11101 88 88 89 89 ! write out time 90 IF( ll_verbose ) WRITE(numicb,9100) nktberg, ndastp, nsec_day90 IF( ll_verbose .AND. numicb.NE.-1) WRITE(numicb,9100) nktberg, ndastp, nsec_day 91 91 9100 FORMAT('kt= ',i8, ' day= ',i8,' secs=',i8) 92 92 … … 166 166 IF( nn_sample_rate .GT. 0 ) CALL icb_trj_end() 167 167 168 IF(lwp) WRITE(numout,'(a,i6)') 'icebergs: icb_end complete', narea 169 CALL flush( numicb ) 170 CLOSE( numicb ) 168 IF(lwp .AND. nprint>3) WRITE(numout,'(a,i6)') 'icebergs: icb_end complete', narea 169 IF(numicb.NE.-1) THEN 170 CALL flush( numicb ) 171 CLOSE( numicb ) 172 ENDIF 171 173 ! 172 174 END SUBROUTINE icb_end -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbtrj.F90
r7733 r11101 83 83 ELSE ; WRITE(cl_filename,'("trajectory_icebergs_",A,"-",A ,".nc")') TRIM(ADJUSTL(cldate_ini)), TRIM(ADJUSTL(cldate_end)) 84 84 ENDIF 85 IF ( lwp .AND. n n_verbose_level >= 0) WRITE(numout,'(2a)') 'icebergs, icb_trj_init: creating ',TRIM(cl_filename)85 IF ( lwp .AND. nprint > 2) WRITE(numout,'(2a)') 'icebergs, icb_trj_init: creating ',TRIM(cl_filename) 86 86 87 87 iret = NF90_CREATE(TRIM(cl_filename), NF90_CLOBBER, ntrajid) … … 246 246 this=>this%next 247 247 END DO 248 IF( lwp . and. nn_verbose_level> 0 ) WRITE(numout,*) 'trajectory write to frame ', jn248 IF( lwp .AND. nprint > 0 ) WRITE(numout,*) 'trajectory write to frame ', jn 249 249 num_traj = jn 250 250 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/ICB/icbutl.F90
r6486 r11101 661 661 ! 662 662 pt => berg%current_point 663 IF(numicb.NE.-1) THEN 663 664 WRITE(numicb, 9200) kt, berg%number(1), & 664 665 pt%xi, pt%yj, pt%lon, pt%lat, pt%uvel, pt%vvel, & 665 666 pt%uo, pt%vo, pt%ua, pt%va, pt%ui, pt%vi 666 667 CALL flush( numicb ) 668 ENDIF 667 669 9200 FORMAT(5x,i5,2x,i10,6(2x,2f10.4)) 668 670 ! … … 685 687 ! 686 688 this => first_berg 687 IF( ASSOCIATED(this) ) THEN 688 WRITE(numicb,'(a," pe=(",i3,")")' ) cd_label, narea 689 WRITE(numicb,'(a8,4x,a6,12x,a5,15x,a7,19x,a3,17x,a5,17x,a5,17x,a5)' ) & 689 IF(numicb.NE.-1) THEN 690 IF( ASSOCIATED(this) ) THEN 691 WRITE(numicb,'(a," pe=(",i3,")")' ) cd_label, narea 692 WRITE(numicb,'(a8,4x,a6,12x,a5,15x,a7,19x,a3,17x,a5,17x,a5,17x,a5)' ) & 690 693 & 'timestep', 'number', 'xi,yj','lon,lat','u,v','uo,vo','ua,va','ui,vi' 694 ENDIF 691 695 ENDIF 692 696 DO WHILE( ASSOCIATED(this) ) … … 697 701 inbergs = ibergs 698 702 IF( lk_mpp ) CALL mpp_sum(inbergs) 699 IF( ibergs > 0 ) WRITE(numicb,'(a," there are",i5," bergs out of",i6," on PE ",i4)') &703 IF( ibergs > 0 .AND. numicb.NE.-1) WRITE(numicb,'(a," there are",i5," bergs out of",i6," on PE ",i4)') & 700 704 & cd_label, ibergs, inbergs, narea 701 705 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90
r6491 r11101 99 99 !!---------------------------------------------------------------------- 100 100 LOGICAL :: ln_ctl !: run control for debugging 101 TYPE :: sn_ctl !: optional use structure for finer control over output selection 102 LOGICAL :: l_config = .FALSE. !: activate/deactivate finer control 103 ! Note if l_config is True then ln_ctl is ignored. 104 ! Otherwise setting ln_ctl True is equivalent to setting 105 ! all the following logicals in this structure True 106 LOGICAL :: l_runstat = .FALSE. !: Produce/do not produce run.stat file (T/F) 107 LOGICAL :: l_trcstat = .FALSE. !: Produce/do not produce tracer.stat file (T/F) 108 LOGICAL :: l_oceout = .FALSE. !: Produce all ocean.outputs (T) or just one (F) 109 LOGICAL :: l_layout = .FALSE. !: Produce all layout.dat files (T) or just one (F) 110 LOGICAL :: l_mppout = .FALSE. !: Produce/do not produce mpp.output_XXXX files (T/F) 111 LOGICAL :: l_mpptop = .FALSE. !: Produce/do not produce mpp.top.output_XXXX files (T/F) 112 ! Optional subsetting of processor report files 113 ! Default settings of 0/1000000/1 should ensure all areas report. 114 ! Set to a more restrictive range to select specific areas 115 INTEGER :: procmin = 0 !: Minimum narea to output 116 INTEGER :: procmax = 1000000 !: Maximum narea to output 117 INTEGER :: procincr = 1 !: narea increment to output 118 INTEGER :: ptimincr = 1 !: timestep increment to output (time.step and run.stat) 119 END TYPE 120 TYPE(sn_ctl) :: sn_cfctl !: run control structure for selective output 101 121 INTEGER :: nn_timing !: run control for timing 102 122 INTEGER :: nn_print !: level of print (0 no print) … … 109 129 INTEGER :: nn_bench !: benchmark parameter (0/1) 110 130 INTEGER :: nn_bit_cmp = 0 !: bit reproducibility (0/1) 111 131 LOGICAL :: ln_flush !: namelist flush numout after write statement 112 132 ! 133 LOGICAL :: lflush !: flush numout after write statement 113 134 INTEGER :: nprint, nictls, nictle, njctls, njctle, isplt, jsplt, nbench !: OLD namelist names 114 135 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90
r9830 r11101 484 484 iom_file(jf)%nfid = 0 ! free the id 485 485 IF( PRESENT(kiomid) ) kiomid = 0 ! return 0 as id to specify that the file was closed 486 IF(lwp ) WRITE(numout,*) TRIM(clinfo)//' close file: '//TRIM(iom_file(jf)%name)//' ok'486 IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' close file: '//TRIM(iom_file(jf)%name)//' ok' 487 487 ELSEIF( PRESENT(kiomid) ) THEN 488 488 WRITE(ctmp1,*) '--->', kiomid … … 924 924 925 925 IF( istop == nstop ) THEN ! no additional errors until this point... 926 IF(lwp) WRITE(numout,"(10x,' read ',a,' (rec: ',i6,') in ',a,' ok')") TRIM(cdvar), itime, TRIM(iom_file(kiomid)%name) 926 IF(lwp .AND. nprint > 1) & 927 WRITE(numout,"(10x,' read ',a,' (rec: ',i6,') in ',a,' ok')") TRIM(cdvar), itime, TRIM(iom_file(kiomid)%name) 927 928 928 929 !--- overlap areas and extra hallows (mpp) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/iom_nf90.F90
r6491 r11101 84 84 ! ! ============= 85 85 IF( ldwrt ) THEN ! ... in write mode 86 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in WRITE mode' 86 IF(lwp .AND. nprint > 0) THEN 87 WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in WRITE mode' 88 IF(lflush) CALL flush(numout) 89 ENDIF 87 90 IF( snc4set%luse ) THEN 88 91 CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_WRITE , if90id ), clinfo) … … 92 95 CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy ), clinfo) 93 96 ELSE ! ... in read mode 94 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in READ mode' 97 IF(lwp .AND. nprint > 0) THEN 98 WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in READ mode' 99 IF(lflush) CALL flush(numout) 100 ENDIF 95 101 CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_NOWRITE, if90id, chunksize = ichunk ), clinfo) 96 102 ENDIF … … 103 109 cdname = TRIM(cltmp) 104 110 ENDIF 105 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' create new file: '//TRIM(cdname)//' in WRITE mode' 111 IF(lwp .AND. nprint > 0) THEN 112 WRITE(numout,*) TRIM(clinfo)//' create new file: '//TRIM(cdname)//' in WRITE mode' 113 IF(lflush) CALL flush(numout) 114 ENDIF 106 115 107 116 IF( llclobber ) THEN ; imode = IOR( NF90_64BIT_OFFSET, NF90_CLOBBER ) … … 109 118 ENDIF 110 119 IF( snc4set%luse ) THEN 111 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' creating file: '//TRIM(cdname)//' in hdf5 (netcdf4) mode' 120 IF(lwp .AND. nprint > 0) THEN 121 WRITE(numout,*) TRIM(clinfo)//' creating file: '//TRIM(cdname)//' in hdf5 (netcdf4) mode' 122 IF(lflush) CALL flush(numout) 123 ENDIF 112 124 CALL GET_NF90_SYMBOL("NF90_HDF5", ihdf5) 113 125 IF( llclobber ) THEN ; imode = IOR(ihdf5, NF90_CLOBBER) … … 158 170 & len = iom_file(kiomid)%lenuld ), clinfo ) 159 171 ENDIF 160 IF(lwp) WRITE(numout,*) ' ---> '//TRIM(cdname)//' OK' 172 IF(lwp .AND. nprint > 0) THEN 173 WRITE(numout,*) ' ---> '//TRIM(cdname)//' OK' 174 IF(lflush) CALL flush(numout) 175 ENDIF 161 176 ELSE 162 177 kiomid = 0 ! return error flag … … 431 446 ! trick: defined to 0 to say that dimension variables are defined but not yet written 432 447 iom_file(kiomid)%dimsz(1, 1) = 0 433 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' define dimension variables done' 448 IF(lwp .AND. nprint > 1) THEN 449 WRITE(numout,*) TRIM(clinfo)//' define dimension variables done' 450 IF(lflush) CALL flush(numout) 451 ENDIF 434 452 ENDIF 435 453 ! define the data if it is not already done … … 498 516 CALL iom_nf90_check(SET_NF90_DEF_VAR_CHUNKING(if90id, idvar, ichunkalg, ichunksz), clinfo) 499 517 CALL iom_nf90_check(SET_NF90_DEF_VAR_DEFLATE(if90id, idvar, ishuffle, ideflate, ideflate_level), clinfo) 500 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' chunked ok. Chunks sizes: ', ichunksz 501 ENDIF 502 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' defined ok' 518 IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' chunked ok. Chunks sizes: ', ichunksz 519 ENDIF 520 IF(lwp .AND. nprint > 1) THEN 521 WRITE(numout,*) TRIM(clinfo)//' defined ok' 522 IF(lflush) CALL flush(numout) 523 ENDIF 503 524 ELSE 504 525 idvar = kvid … … 543 564 CALL iom_nf90_check(NF90_INQUIRE_DIMENSION( if90id, 3, len = iom_file(kiomid)%dimsz(1,3) ), clinfo) 544 565 iom_file(kiomid)%dimsz(1 , 4) = 1 ! unlimited dimension 545 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' write dimension variables done' 566 IF(lwp .AND. nprint > 1) THEN 567 WRITE(numout,*) TRIM(clinfo)//' write dimension variables done' 568 IF(lflush) CALL flush(numout) 569 ENDIF 546 570 ENDIF 547 571 ENDIF … … 561 585 IF( iom_file(kiomid)%luld(idvar) ) iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar) & 562 586 & = iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar) + 1 563 IF(lwp) WRITE(numout,*) TRIM(clinfo)//' written ok' 587 IF(lwp .AND. nprint > 0) THEN 588 WRITE(numout,*) TRIM(clinfo)//' written ok' 589 IF(lflush) CALL flush(numout) 590 ENDIF 564 591 ENDIF 565 592 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/prtctl.F90
r6486 r11101 369 369 WRITE(j_id,*) 370 370 WRITE(j_id,*) 371 IF(lflush) CALL flush(j_id) 371 372 372 373 9000 FORMAT(a41,i4.4,a14) … … 515 516 IF(lwp) WRITE(numout,*)' sum ilcitl(1,j) = ', zjdom, ' jpjglo = ', jpjglo 516 517 IF(lwp) WRITE(numout,*) 518 IF(lwp .AND. lflush) CALL flush(numout) 517 519 518 520 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r9321 r11101 110 110 ' open ocean restart NetCDF file: ',TRIM(clpath)//clname 111 111 END SELECT 112 IF ( snc4set%luse ) WRITE(numout,*) ' opened for NetCDF4 chunking and compression' 113 IF( kt == nitrst - 1 ) THEN ; WRITE(numout,*) ' kt = nitrst - 1 = ', kt 114 ELSE ; WRITE(numout,*) ' kt = ' , kt 112 IF(nprint > 1) THEN 113 IF ( snc4set%luse ) WRITE(numout,*) ' opened for NetCDF4 chunking and compression' 114 IF( kt == nitrst - 1) THEN 115 WRITE(numout,*) ' kt = nitrst - 1 = ', kt 116 ELSE 117 WRITE(numout,*) ' kt = ' , kt 118 ENDIF 115 119 ENDIF 120 IF(lflush) CALL flush(numout) 116 121 ENDIF 117 122 ! … … 209 214 IF ( snc4set%luse ) WRITE(numout,*) 'rst_read : configured with NetCDF4 support' 210 215 WRITE(numout,*) '~~~~~~~~' 216 IF(lflush) CALL flush(numout) 211 217 ENDIF 212 218 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LBC/cla.F90
r6486 r11101 83 83 IF(lwp) WRITE(numout,*) 'div_cla : cross land advection on hdiv ' 84 84 IF(lwp) WRITE(numout,*) '~~~~~~~~' 85 IF(lflush) CALL flush(numout) 85 86 ! 86 87 IF( nbab == 1 ) CALL cla_bab_el_mandeb('ini') ! Bab el Mandeb ( Red Sea - Indian ocean ) … … 120 121 IF(lwp) WRITE(numout,*) 'tra_cla : cross land advection on tracers ' 121 122 IF(lwp) WRITE(numout,*) '~~~~~~~~' 123 IF(lflush) CALL flush(numout) 122 124 ENDIF 123 125 ! … … 147 149 IF(lwp) WRITE(numout,*) 'cla_dynspg : cross land advection on (ua,va) ' 148 150 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 151 IF(lflush) CALL flush(numout) 149 152 ENDIF 150 153 ! … … 181 184 IF(lwp) WRITE(numout,*) 'cla_init : cross land advection initialisation ' 182 185 IF(lwp) WRITE(numout,*) '~~~~~~~~~' 186 IF(lwp .AND. lflush) CALL flush(numout) 183 187 ! 184 188 ! ! Allocate arrays for this module -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90
r8356 r11101 319 319 mynode = mpprank 320 320 321 IF( mynode == 0 ) THEN321 IF( mynode == 0 .AND. nprint > 2) THEN 322 322 CALL ctl_opn( kumond, TRIM(ldname), 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. , 1 ) 323 323 WRITE(kumond, nammpp) … … 3989 3989 IF( iost == 0 ) THEN 3990 3990 IF(ldwp) THEN 3991 WRITE(kout,*) ' file : ', clfile,' open ok' 3992 WRITE(kout,*) ' unit = ', knum 3993 WRITE(kout,*) ' status = ', cdstat 3994 WRITE(kout,*) ' form = ', cdform 3995 WRITE(kout,*) ' access = ', cdacce 3991 IF(nprint > 0) WRITE(kout,*) ' file : ', clfile,' open ok' 3992 IF(nprint > 2) THEN 3993 WRITE(kout,*) ' unit = ', knum 3994 WRITE(kout,*) ' status = ', cdstat 3995 WRITE(kout,*) ' form = ', cdform 3996 WRITE(kout,*) ' access = ', cdacce 3997 ENDIF 3996 3998 WRITE(kout,*) 3997 3999 ENDIF -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90
r6498 r11101 73 73 WRITE(numout,*) ' nimpp = ', nimpp 74 74 WRITE(numout,*) ' njmpp = ', njmpp 75 IF(lflush) CALL flush(numout) 75 76 ENDIF 76 77 … … 137 138 IF(lwp) WRITE(numout,*) 'mpp_init : Message Passing MPI' 138 139 IF(lwp) WRITE(numout,*) '~~~~~~~~' 140 IF(lflush) CALL flush(numout) 139 141 140 142 … … 287 289 IF(lwp) WRITE(numout,*)' sum ilcit(1,j)=', zjdom, ' jpjglo=', jpjglo 288 290 IF(lwp) WRITE(numout,*) 289 290 IF(lwp) THEN 291 IF(lwp .AND. lflush) CALL flush(numout) 292 293 IF(lwp .AND. nprint>1) THEN 291 294 ifreq = 4 292 295 il1 = 1 … … 305 308 il1 = il1+ifreq 306 309 END DO 310 IF(lflush) CALL flush(numout) 311 ENDIF 307 312 9200 FORMAT(' ***',20('*************',a3)) 308 313 9203 FORMAT(' * ',20(' * ',a3)) … … 310 315 9202 FORMAT(' ',i3,' * ',20(i3,' x',i3,' * ')) 311 316 9204 FORMAT(' * ',20(' ',i3,' * ')) 312 ENDIF313 317 314 318 ! 5. From global to local … … 364 368 365 369 ! Save processor layout in layout.dat file 366 IF (lwp ) THEN370 IF (lwp .AND. nprint > 1) THEN 367 371 CALL ctl_opn( inum, 'layout.dat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE., narea ) 368 372 WRITE(inum,'(a)') ' jpnij jpi jpj jpk jpiglo jpjglo' … … 430 434 ! Periodicity : no corner if nbondi = 2 and nperio != 1 431 435 432 IF(lwp ) THEN436 IF(lwp .and. nprint > 0) THEN 433 437 WRITE(numout,*) ' nproc = ', nproc 434 438 WRITE(numout,*) ' nowe = ', nowe , ' noea = ', noea … … 447 451 WRITE(numout,*) ' jprecj = ', jprecj , ' npnw = ', npnw 448 452 WRITE(numout,*) 453 IF(lflush) CALL flush(numout) 449 454 ENDIF 450 455 … … 456 461 CALL mpp_ini_north 457 462 IF(lwp) WRITE(numout,*) ' mpp_init : North fold boundary prepared for jpni >1' 463 IF(lwp .AND. lflush) CALL flush(numout) 458 464 ENDIF 459 465 … … 505 511 idid(2) = 2 506 512 507 IF(lwp ) THEN513 IF(lwp .AND. nprint > 0) THEN 508 514 WRITE(numout,*) 509 515 WRITE(numout,*) 'mpp_init_ioipsl : iloc = ', iloc (1), iloc (2) … … 511 517 WRITE(numout,*) ' ihals = ', ihals(1), ihals(2) 512 518 WRITE(numout,*) ' ihale = ', ihale(1), ihale(2) 519 IF(lflush) CALL flush(numout) 513 520 ENDIF 514 521 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90
r6498 r11101 87 87 READ ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 ) 88 88 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp ) 89 IF(lwm ) WRITE ( numond, namzgr )89 IF(lwm .AND. nprint > 2) WRITE ( numond, namzgr ) 90 90 91 91 IF(lwp)WRITE(numout,*) … … 93 93 IF(lwp)WRITE(numout,*) '~~~~~~~~' 94 94 IF(lwp)WRITE(numout,*) ' ' 95 IF(lwp .AND. lflush) CALL flush(numout) 95 96 96 97 IF( jpni*jpnj < jpnij ) CALL ctl_stop( ' jpnij > jpni x jpnj impossible' ) … … 184 185 IF(lwp) WRITE(numout,*) ' sum ilcj(1,j)=',zjdom,' jpjglo=',jpjglo 185 186 IF(lwp) WRITE(numout,*) 187 IF(lwp .AND. lflush) CALL flush(numout) 186 188 187 189 … … 364 366 il1 = il1+ifreq 365 367 END DO 368 IF(lflush) CALL flush(numout) 366 369 9400 FORMAT(' ***',20('*************',a3)) 367 370 9403 FORMAT(' * ',20(' * ',a3)) … … 572 575 WRITE(numout,*) ' jprecj = ', jprecj , ' npnw = ', npnw 573 576 WRITE(numout,*) 577 IF(lflush) CALL flush(numout) 574 578 ENDIF 575 579 … … 581 585 CALL mpp_ini_north 582 586 IF(lwp) WRITE(numout,*) ' mpp_init2 : North fold boundary prepared for jpni >1' 587 IF(lwp .AND. lflush) CALL flush(numout) 583 588 ENDIF 584 589 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90
r6486 r11101 81 81 READ ( numnam_cfg, namdyn_ldf, IOSTAT = ios, ERR = 902 ) 82 82 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_ldf in configuration namelist', lwp ) 83 IF(lwm ) WRITE ( numond, namdyn_ldf )83 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_ldf ) 84 84 85 85 IF(lwp) THEN ! Parameter print … … 159 159 nkahm_smag = 1 160 160 #endif 161 161 ! 162 IF(lwp .AND. lflush) CALL flush(numout) 162 163 ! 163 164 END SUBROUTINE ldf_dyn_init … … 210 211 ENDIF 211 212 ! 213 IF(lwp .AND. lflush) CALL flush(numout) 214 ! 212 215 END SUBROUTINE ldf_zpf_1d 213 216 … … 252 255 ENDIF 253 256 ! 257 IF(lwp .AND. lflush) CALL flush(numout) 258 ! 254 259 END SUBROUTINE ldf_zpf_1d_3d 255 260 … … 293 298 ENDIF 294 299 ! 300 IF(lwp .AND. lflush) CALL flush(numout) 301 ! 295 302 END SUBROUTINE ldf_zpf_3d 296 303 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90
r6486 r11101 41 41 IF(lwp) WRITE(numout,*) 'ldf_dyn_c2d : 2d lateral eddy viscosity coefficient' 42 42 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 43 IF(lwp .AND. lflush) CALL flush(numout) 43 44 44 45 ! harmonic operator (ahm1, ahm2) : ( T- and F- points) (used for laplacian operators … … 85 86 CALL prihre(ahm2,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 86 87 ENDIF 88 IF(lwp .AND. lflush) CALL flush(numout) 87 89 ENDIF 88 90 … … 119 121 CALL prihre(ahm4,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 120 122 ENDIF 123 IF(lwp .AND. lflush) CALL flush(numout) 121 124 ENDIF 122 125 ! … … 160 163 IF(lwp) WRITE(numout,*) '~~~~~~ --' 161 164 IF(lwp) WRITE(numout,*) ' orca ocean configuration' 165 IF(lwp .AND. lflush) CALL flush(numout) 162 166 163 167 IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN … … 316 320 ENDIF 317 321 ! 322 IF(lwp .AND. lflush) CALL flush(numout) 323 ! 318 324 CALL wrk_dealloc( jpi , jpj , icof ) 319 325 ! … … 358 364 IF(lwp) WRITE(numout,*) '~~~~~~ --' 359 365 IF(lwp) WRITE(numout,*) ' orca_r1 configuration' 366 IF(lwp .AND. lflush) CALL flush(numout) 360 367 361 368 IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN … … 519 526 ENDIF 520 527 ! 528 IF(lwp .AND. lflush) CALL flush(numout) 529 ! 521 530 CALL wrk_dealloc( jpi , jpj , icof ) 522 531 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90
r9975 r11101 179 179 ENDIF 180 180 ENDIF 181 ! 182 IF(lwp .AND. lflush) CALL flush(numout) 181 183 ! 182 184 CALL wrk_dealloc( jpk, zcoef ) … … 432 434 ENDIF 433 435 ! 436 IF(lwp .AND. lflush) CALL flush(numout) 437 ! 434 438 CALL wrk_dealloc( jpi , jpj , icof ) 435 439 CALL wrk_dealloc( jpk , zcoef ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_smag.F90
r9583 r11101 106 106 107 107 108 IF(lwp) WRITE(numout,*) 109 IF(lwp) WRITE(numout,*) 'ldf_dyn_smag : 3D lateral eddy viscosity coefficient' 110 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 111 108 IF(lwp) THEN 109 WRITE(numout,*) 110 WRITE(numout,*) 'ldf_dyn_smag : 3D lateral eddy viscosity coefficient' 111 WRITE(numout,*) '~~~~~~~~~~~' 112 IF(lflush) CALL flush(numout) 113 ENDIF 114 112 115 ENDIF 113 116 … … 172 175 173 176 ENDDO ! jpk 174 ahm1(:,:,jpk) = ahm1(:,:,jpkm1) 175 ahm2(:,:,jpk) = ahm2(:,:,jpkm1) 176 IF(lwp.and.kt==nit000) WRITE(numout,'(36x," ahm ", 7x)') 177 178 ahm1(:,:,jpk) = ahm1(:,:,jpkm1) 179 ahm2(:,:,jpk) = ahm2(:,:,jpkm1) 180 181 IF(lwp.and.kt==nit000) THEN 182 WRITE(numout,'(36x," ahm ", 7x)') 183 IF(lflush) CALL flush(numout) 184 ENDIF 185 186 IF(lwp.and.kt==nit000) THEN 177 187 DO jk = 1, jpk 178 179 IF(lwp.and.kt==nit000) WRITE(numout,'(30x,E10.2,8x,i3)') ahm1(jpi/2,jpj/2,jk), jk 188 WRITE(numout,'(30x,E10.2,8x,i3)') ahm1(jpi/2,jpj/2,jk), jk 180 189 END DO 190 IF(lflush) CALL flush(numout) 191 ENDIF 192 181 193 CALL lbc_lnk( ahm1, 'T', 1. ) ! Lateral boundary conditions on ( ahtt ) 182 194 CALL lbc_lnk( ahm2, 'F', 1. ) ! Lateral boundary conditions on ( ahtt ) … … 245 257 ahm4(:,:,jpk) = ahm4(:,:,jpkm1) 246 258 247 DO jk = 1, jpk 248 IF( kt == nit000 ) THEN 249 250 IF(lwp) WRITE(numout,'(30x,E10.2,8x,i3)') ahm3(jpi/2,jpj/2,jk), jk 259 IF( kt == nit000 .AND. lwp) THEN 260 DO jk = 1, jpk 261 WRITE(numout,'(30x,E10.2,8x,i3)') ahm3(jpi/2,jpj/2,jk), jk 262 IF(lflush) CALL flush(numout) 263 END DO 251 264 ENDIF 252 END DO 265 253 266 CALL lbc_lnk( ahm3, 'U', 1. ) ! Lateral boundary conditions 254 267 CALL lbc_lnk( ahm4, 'V', 1. ) … … 269 282 IF( lk_mpp ) CALL mpp_max( zdeltat ) ! max over the global domain 270 283 ! 271 IF( MOD( kt, nwrite ) == 1 .AND. lwp ) WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap: abs(ahm) max: ', zdeltat 284 IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 285 WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap: abs(ahm) max: ', zdeltat 286 IF(lflush) CALL flush(numout) 287 ENDIF 272 288 ENDIF 273 289 If(ln_dynldf_bilap)THEN … … 282 298 IF( lk_mpp ) CALL mpp_max( zdeltat ) ! max over the global domain 283 299 ! 284 IF( MOD( kt, nwrite ) == 1 .AND. lwp ) WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 300 IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 301 WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 302 IF(lflush) CALL flush(numout) 303 ENDIF 285 304 ! 286 305 ENDIF -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90
r6498 r11101 414 414 IF(lwp) THEN 415 415 WRITE(numout,*) ' Horizontal mixing in s-coordinate: slope = slope of s-surfaces' 416 IF(lflush) CALL flush(numout) 416 417 ENDIF 417 418 … … 451 452 IF(lwp) WRITE(numout,*) ' max slop: u',SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)), & 452 453 & ' wi', sqrt(MAXVAL(wslpi)), ' wj', sqrt(MAXVAL(wslpj)) 454 IF(lwp .AND. lflush) CALL flush(numout) 453 455 endif 454 456 … … 842 844 WRITE(numout,*) 'ldf_slp_init : direction of lateral mixing' 843 845 WRITE(numout,*) '~~~~~~~~~~~~' 846 IF(lflush) CALL flush(numout) 844 847 ENDIF 845 848 … … 865 868 IF(ln_sco .AND. (ln_traldf_hor .OR. ln_dynldf_hor )) THEN 866 869 IF(lwp) WRITE(numout,*) ' Horizontal mixing in s-coordinate: slope = slope of s-surfaces' 870 IF(lwp .AND. lflush) CALL flush(numout) 867 871 868 872 ! geopotential diffusion in s-coordinates on tracers and/or momentum -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90
r6486 r11101 85 85 READ ( numnam_cfg, namtra_ldf, IOSTAT = ios, ERR = 902 ) 86 86 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_ldf in configuration namelist', lwp ) 87 IF(lwm ) WRITE ( numond, namtra_ldf )87 IF(lwm .AND. nprint > 2) WRITE ( numond, namtra_ldf ) 88 88 89 89 IF(lwp) THEN ! control print … … 175 175 IF(lwp.and.ln_traldf_bilap) CALL ctl_stop(' SMAGORINSKY + BILAPLACIAN - UNSTABLE OR NON_CONSERVATIVE' ) 176 176 #endif 177 177 ! 178 IF(lwp .AND. lflush) CALL flush(numout) 178 179 ! 179 180 END SUBROUTINE ldf_tra_init -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c2d.h90
r6486 r11101 36 36 IF(lwp) WRITE(numout,*) ' ldf_tra_c2d : 2D eddy diffusivity and eddy' 37 37 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~ -- induced velocity coefficients' 38 IF(lwp .AND. lflush) CALL flush(numout) 38 39 ELSE 39 40 IF(lwp) WRITE(numout,*) 40 41 IF(lwp) WRITE(numout,*) ' ldf_tra2d : 2D eddy diffusivity coefficient' 41 42 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~ --' 43 IF(lwp .AND. lflush) CALL flush(numout) 42 44 ENDIF 43 45 … … 75 77 IF(lwp) WRITE(numout,*) ' Constant values used for eddy diffusivity coefficients' 76 78 IF(lwp) WRITE(numout,*) ' Variation lat/lon only for eddy induced velocity coefficients' 79 IF(lwp .AND. lflush) CALL flush(numout) 77 80 ENDIF 78 81 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c3d.h90
r6486 r11101 46 46 IF(lwp) WRITE(numout,*) ' Coefficients set to constant' 47 47 ENDIF 48 48 49 49 50 50 ! biharmonic operator (T-point) … … 121 121 122 122 # endif 123 124 IF(lwp .AND. lflush) CALL flush(numout) 125 123 126 END SUBROUTINE ldf_tra_c3d -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_smag.F90
r9583 r11101 99 99 IF(lwp) WRITE(numout,*) 100 100 IF(lwp) WRITE(numout,*) 101 IF(lwp .AND. lflush) CALL flush(numout) 101 102 ENDIF 102 103 … … 108 109 ahtt(:,:,:) = rn_aht_0 109 110 IF( ln_traldf_bilap ) THEN 110 IF( lwp .AND. kt == nit000) WRITE(numout,* )'ldf_tra_smag :no bilaplacian Smagorinsky diffusivity' 111 IF( lwp .AND. kt == nit000) WRITE(numout,* )'ldf_tra_smag :bilaplacian diffusivity set to constant' 111 IF( lwp .AND. kt == nit000) THEN 112 WRITE(numout,* )'ldf_tra_smag :no bilaplacian Smagorinsky diffusivity' 113 WRITE(numout,* )'ldf_tra_smag :bilaplacian diffusivity set to constant' 114 IF(lflush) CALL flush(numout) 115 ENDIF 112 116 ENDIF 113 117 … … 190 194 CALL prihre( ahtw(:,:,1), jpi, jpj, 1, jpi, 1, & 191 195 & 1, jpj, 1, 1.e-1, numout ) 196 IF(lflush) CALL flush(numout) 192 197 ENDIF 193 198 ENDIF -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r6486 r11101 244 244 READ ( numnam_cfg, namobs, IOSTAT = ios, ERR = 902 ) 245 245 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namobs in configuration namelist', lwp ) 246 IF(lwm ) WRITE ( numond, namobs )246 IF(lwm .AND. nprint > 2) WRITE ( numond, namobs ) 247 247 248 248 ! Count number of files for each type -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90
r6498 r11101 260 260 READ ( numnam_cfg, namsbc_alb, IOSTAT = ios, ERR = 902 ) 261 261 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_alb in configuration namelist', lwp ) 262 IF(lwm ) WRITE ( numond, namsbc_alb )262 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_alb ) 263 263 ! 264 264 IF(lwp) THEN ! Control print … … 269 269 WRITE(numout,*) ' choose the albedo parameterization nn_ice_alb = ', nn_ice_alb 270 270 WRITE(numout,*) ' albedo of bare puddled ice rn_albice = ', rn_albice 271 IF(lflush) CALL flush(numout) 271 272 ENDIF 272 273 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90
r10269 r11101 183 183 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 184 184 IF(lwp) WRITE(numout,*) 185 IF(lflush) CALL flush(numout) 185 186 186 187 ncplmodel = kcplmodel … … 235 236 WRITE(numout,*) ' multiexchg: nldi, nlei, nimpp =', nldi, nlei, nimpp 236 237 WRITE(numout,*) ' multiexchg: nldj, nlej, njmpp =', nldj, nlej, njmpp 238 IF(lflush) CALL flush(numout) 237 239 ENDIF 238 240 … … 313 315 CALL oasis_abort ( ssnd(ji)%nid(jc,jm), 'cpl_define', 'Failure in oasis_def_var' ) 314 316 ENDIF 315 IF( ln_ctl .AND. ssnd(ji)%nid(jc,jm) /= -1 ) WRITE(numout,*) "variable defined in the namcouple" 316 IF( ln_ctl .AND. ssnd(ji)%nid(jc,jm) == -1 ) WRITE(numout,*) "variable NOT defined in the namcouple" 317 IF( ln_ctl) THEN 318 IF(ssnd(ji)%nid(jc,jm) /= -1 ) WRITE(numout,*) "variable defined in the namcouple" 319 IF(ssnd(ji)%nid(jc,jm) == -1 ) WRITE(numout,*) "variable NOT defined in the namcouple" 320 IF(lflush) CALL flush(numout) 321 ENDIF 317 322 END DO 318 323 END DO … … 351 356 END IF 352 357 #endif 353 IF( ln_ctl ) WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_In 354 flush(numout) 358 IF( ln_ctl ) THEN 359 WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_In 360 IF(lflush) CALL flush(numout) 361 ENDIF 355 362 356 363 ! Define 0D (Greenland or Antarctic ice mass) or 1D (river outflow) coupling fields … … 394 401 CALL oasis_abort ( srcv(ji)%nid(jc,jm), 'cpl_define', 'Failure in oasis_def_var' ) 395 402 ENDIF 396 IF( ln_ctl .AND. srcv(ji)%nid(jc,jm) /= -1 ) WRITE(numout,*) "variable defined in the namcouple" 397 IF( ln_ctl .AND. srcv(ji)%nid(jc,jm) == -1 ) WRITE(numout,*) "variable NOT defined in the namcouple" 398 403 IF(ln_ctl) THEN 404 IF(srcv(ji)%nid(jc,jm) /= -1 ) WRITE(numout,*) "variable defined in the namcouple" 405 IF(srcv(ji)%nid(jc,jm) == -1 ) WRITE(numout,*) "variable NOT defined in the namcouple" 406 IF(lflush) CALL flush(numout) 407 ENDIF 399 408 END DO 400 409 END DO … … 447 456 WRITE(numout,*) ' - Sum value is ', SUM(pdata(:,:,jc)) 448 457 WRITE(numout,*) '****************' 458 IF(lflush) CALL flush(numout) 449 459 ENDIF 450 460 ENDIF … … 491 501 & kinfo == OASIS_RecvOut .OR. kinfo == OASIS_FromRestOut 492 502 493 IF ( ln_ctl ) WRITE(numout,*) "llaction, kinfo, kstep, ivarid: " , llaction, kinfo, kstep, srcv(kid)%nid(jc,jm) 503 IF ( ln_ctl ) THEN 504 WRITE(numout,*) "llaction, kinfo, kstep, ivarid: " , llaction, kinfo, kstep, srcv(kid)%nid(jc,jm) 505 IF(lflush) CALL flush(numout) 506 ENDIF 494 507 495 508 IF ( llaction ) THEN … … 513 526 WRITE(numout,*) ' - Sum value is ', SUM(pdata(:,:,jc)) 514 527 WRITE(numout,*) '****************' 528 IF(lflush) CALL flush(numout) 515 529 ENDIF 516 530 … … 603 617 WRITE(numout,*) ' - Start of data is ', pdata(1:number_to_print) 604 618 WRITE(numout,*) '****************' 619 IF(lflush) CALL flush(numout) 605 620 ENDIF 606 621 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_rnf_1d.F90
r10269 r11101 70 70 IF(lwp) WRITE(numout,*) 'cpl_rnf_1d_init : initialization of river runoff coupling' 71 71 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 72 IF(lwp .AND. lflush) CALL flush(numout) 72 73 73 74 REWIND(numnam_cfg) … … 78 79 READ ( numnam_cfg, nam_cpl_rnf_1d, IOSTAT = ios, ERR = 902 ) 79 80 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_cpl_rnf_1d in configuration namelist', lwp ) 80 IF(lwm ) WRITE ( numond, nam_cpl_rnf_1d )81 IF(lwm .AND. nprint > 2) WRITE ( numond, nam_cpl_rnf_1d ) 81 82 82 83 ! ! Parameter control and print … … 87 88 IF(lwp) WRITE(numout,*) ' Print river information = ',ln_print_river_info 88 89 IF(lwp) WRITE(numout,*) ' ' 90 IF(lwp .AND. lflush) CALL flush(numout) 89 91 90 92 ! Assign space for river numbers … … 109 111 110 112 ! Print out the largest river number 111 IF ( ln_print_river_info ) THEN113 IF ( ln_print_river_info .AND. lwp) THEN 112 114 WRITE(numout,*) 'Maximum river number in input file = ',max_river 115 IF(lflush) CALL flush(numout) 113 116 END IF 114 117 … … 137 140 rivers%river_area(:) = REAL(rivers%river_area_c(:),wp) 138 141 139 IF ( ln_print_river_info ) THEN142 IF ( ln_print_river_info .AND. lwp) THEN 140 143 WRITE(numout,*) 'Area of rivers 1 to 10 are ',rivers%river_area(1:10) 144 IF(lflush) CALL flush(numout) 141 145 END IF 142 146 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90
r6487 r11101 275 275 ! 276 276 IF( sd(jf)%ln_tint ) THEN ! temporal interpolation 277 IF(lwp .AND. kt - nit000 <= 100) THEN277 IF(lwp .AND. nprint > 1 ) THEN 278 278 clfmt = "('fld_read: var ', a, ' kt = ', i8, ' (', f9.4,' days), Y/M/D = ', i4.4,'/', i2.2,'/', i2.2," // & 279 279 & "', records b/a: ', i6.4, '/', i6.4, ' (days ', f9.4,'/', f9.4, ')')" … … 281 281 & sd(jf)%nrec_b(1), sd(jf)%nrec_a(1), REAL(sd(jf)%nrec_b(2),wp)/rday, REAL(sd(jf)%nrec_a(2),wp)/rday 282 282 WRITE(numout, *) 'it_offset is : ',it_offset 283 IF(lflush) CALL flush(numout) 283 284 ENDIF 284 285 ! temporal interpolation weights … … 288 289 sd(jf)%fnow(:,:,:) = ztintb * sd(jf)%fdta(:,:,:,1) + ztinta * sd(jf)%fdta(:,:,:,2) 289 290 ELSE ! nothing to do... 290 IF(lwp .AND. kt - nit000 <= 100) THEN291 IF(lwp .AND. nprint > 1 ) THEN 291 292 clfmt = "('fld_read: var ', a, ' kt = ', i8,' (', f9.4,' days), Y/M/D = ', i4.4,'/', i2.2,'/', i2.2," // & 292 293 & "', record: ', i6.4, ' (days ', f9.4, ' <-> ', f9.4, ')')" 293 294 WRITE(numout, clfmt) TRIM(sd(jf)%clvar), kt, REAL(isecsbc,wp)/rday, nyear, nmonth, nday, & 294 295 & sd(jf)%nrec_a(1), REAL(sd(jf)%nrec_b(2),wp)/rday, REAL(sd(jf)%nrec_a(2),wp)/rday 296 IF(lflush) CALL flush(numout) 295 297 ENDIF 296 298 ENDIF … … 414 416 415 417 clfmt = "('fld_init : time-interpolation for ', a, ' read previous record = ', i6, ' at time = ', f7.2, ' days')" 416 IF(lwp) WRITE(numout, clfmt) TRIM(sdjf%clvar), sdjf%nrec_a(1), REAL(sdjf%nrec_a(2),wp)/rday 418 IF(lwp .AND. nprint > 0) THEN 419 WRITE(numout, clfmt) TRIM(sdjf%clvar), sdjf%nrec_a(1), REAL(sdjf%nrec_a(2),wp)/rday 420 ENDIF 417 421 418 422 ENDIF … … 708 712 IF(lwp) WRITE(numout,*) 'Dim size for ',TRIM(clvar),' is ', ilendta 709 713 IF(lwp) WRITE(numout,*) 'Number of levels for ',TRIM(clvar),' is ', ipk 714 IF(lwp .AND. lflush) CALL flush(numout) 710 715 711 716 SELECT CASE( ipk ) … … 779 784 END DO 780 785 sd(ju)%rotn(jn) = .TRUE. ! vector was rotated 781 IF( lwp .AND. kt == nit000 ) WRITE(numout,*) & 786 IF( lwp .AND. kt == nit000 ) THEN 787 WRITE(numout,*) & 782 788 & 'fld_read: vector pair ('//TRIM(sd(ju)%clvar)//', '//TRIM(sd(iv)%clvar)//') rotated on to model grid' 789 IF(lflush) CALL flush(numout) 790 ENDIF 783 791 ENDIF 784 792 ENDIF … … 941 949 & ' data type: ' , sdf(jf)%cltype , & 942 950 & ' land/sea mask:' , TRIM( sdf(jf)%lsmname ) 943 call flush(numout)944 951 END DO 952 IF(lflush) CALL flush(numout) 945 953 ENDIF 946 954 … … 1000 1008 !!---------------------------------------------------------------------- 1001 1009 ! 1002 DO kw = 1, nxt_wgt-1 1003 WRITE(numout,*) 'weight file: ',TRIM(ref_wgts(kw)%wgtname) 1004 WRITE(numout,*) ' ddims: ',ref_wgts(kw)%ddims(1),ref_wgts(kw)%ddims(2) 1005 WRITE(numout,*) ' numwgt: ',ref_wgts(kw)%numwgt 1006 WRITE(numout,*) ' jpiwgt: ',ref_wgts(kw)%jpiwgt 1007 WRITE(numout,*) ' jpjwgt: ',ref_wgts(kw)%jpjwgt 1008 WRITE(numout,*) ' botleft: ',ref_wgts(kw)%botleft 1009 WRITE(numout,*) ' topright: ',ref_wgts(kw)%topright 1010 IF( ref_wgts(kw)%cyclic ) THEN 1011 WRITE(numout,*) ' cyclical' 1012 IF( ref_wgts(kw)%overlap > 0 ) WRITE(numout,*) ' with overlap of ', ref_wgts(kw)%overlap 1013 ELSE 1014 WRITE(numout,*) ' not cyclical' 1015 ENDIF 1016 IF( ASSOCIATED(ref_wgts(kw)%data_wgt) ) WRITE(numout,*) ' allocated' 1017 END DO 1010 IF(nprint > 0) THEN 1011 DO kw = 1, nxt_wgt-1 1012 WRITE(numout,*) 'weight file: ',TRIM(ref_wgts(kw)%wgtname) 1013 IF(nprint > 1) THEN 1014 WRITE(numout,*) ' ddims: ',ref_wgts(kw)%ddims(1),ref_wgts(kw)%ddims(2) 1015 WRITE(numout,*) ' numwgt: ',ref_wgts(kw)%numwgt 1016 WRITE(numout,*) ' jpiwgt: ',ref_wgts(kw)%jpiwgt 1017 WRITE(numout,*) ' jpjwgt: ',ref_wgts(kw)%jpjwgt 1018 WRITE(numout,*) ' botleft: ',ref_wgts(kw)%botleft 1019 WRITE(numout,*) ' topright: ',ref_wgts(kw)%topright 1020 IF( ref_wgts(kw)%cyclic ) THEN 1021 WRITE(numout,*) ' cyclical' 1022 IF( ref_wgts(kw)%overlap > 0 ) WRITE(numout,*) ' with overlap of ', ref_wgts(kw)%overlap 1023 ELSE 1024 WRITE(numout,*) ' not cyclical' 1025 ENDIF 1026 IF( ASSOCIATED(ref_wgts(kw)%data_wgt) ) WRITE(numout,*) ' allocated' 1027 ENDIF 1028 END DO 1029 1030 IF(lflush) CALL flush(numout) 1031 1032 ENDIF 1018 1033 ! 1019 1034 END SUBROUTINE wgt_print -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90
r6486 r11101 81 81 READ ( numnam_cfg, namsbc_ana, IOSTAT = ios, ERR = 902 ) 82 82 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ana in configuration namelist', lwp ) 83 IF(lwm ) WRITE ( numond, namsbc_ana )83 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_ana ) 84 84 ! 85 85 IF(lwp) WRITE(numout,*)' ' -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90
r9321 r11101 84 84 READ ( numnam_cfg, namsbc_apr, IOSTAT = ios, ERR = 902 ) 85 85 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_apr in configuration namelist', lwp ) 86 IF(lwm ) WRITE ( numond, namsbc_apr )86 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_apr ) 87 87 ! 88 88 ALLOCATE( sf_apr(1), STAT=ierror ) !* allocate and fill sf_sst (forcing structure) with sn_sst … … 99 99 WRITE(numout,*) ' Namelist namsbc_apr : Atmospheric PRessure as extrenal forcing' 100 100 WRITE(numout,*) ' ref. pressure: global mean Patm (T) or a constant (F) ln_ref_apr = ', ln_ref_apr 101 IF(lflush) CALL flush(numout) 101 102 ENDIF 102 103 ! … … 114 115 IF(lwp) WRITE(numout,*) ' Inverse barometer added to OBC ssh data' 115 116 ENDIF 117 118 IF(lwp .AND. lflush) CALL flush(numout) 119 116 120 IF( ( ln_apr_obc ) .AND. .NOT. lk_dynspg_ts ) & 117 121 CALL ctl_stop( 'sbc_apr: use inverse barometer ssh at open boundary ONLY possible with time-splitting' ) … … 152 156 ssh_ibb(:,:) = ssh_ib(:,:) 153 157 ENDIF 158 159 IF(lwp .AND. lflush) CALL flush(numout) 160 154 161 ENDIF 155 162 ! ! ---------------------------------------- ! … … 159 166 IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 160 167 IF(lwp) WRITE(numout,*) '~~~~' 168 IF(lwp .AND. lflush) CALL flush(numout) 161 169 IF(nn_timing == 2) CALL timing_start('iom_rstput') 162 170 CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r6498 r11101 151 151 READ ( numnam_cfg, namsbc_clio, IOSTAT = ios, ERR = 902 ) 152 152 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_clio in configuration namelist', lwp ) 153 IF(lwm ) WRITE ( numond, namsbc_clio )153 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_clio ) 154 154 155 155 ! store namelist information in an array -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r6823 r11101 166 166 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_core in configuration namelist', lwp ) 167 167 168 IF(lwm ) WRITE( numond, namsbc_core )168 IF(lwm .AND. nprint > 2) WRITE( numond, namsbc_core ) 169 169 ! ! check: do we plan to use ln_dm2dc with non-daily forcing? 170 170 IF( ln_dm2dc .AND. sn_qsr%nfreqh /= 24 ) & -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90
r6486 r11101 141 141 READ ( numnam_cfg, namsbc_mfs, IOSTAT = ios, ERR = 902 ) 142 142 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_mfs in configuration namelist', lwp ) 143 IF(lwm ) WRITE ( numond, namsbc_mfs )143 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_mfs ) 144 144 ! 145 145 ! store namelist information in an array -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r10269 r11101 276 276 READ ( numnam_cfg, namsbc_cpl, IOSTAT = ios, ERR = 902 ) 277 277 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cpl in configuration namelist', lwp ) 278 IF(lwm ) WRITE ( numond, namsbc_cpl )278 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_cpl ) 279 279 280 280 IF(lwp) THEN ! control print … … 331 331 ENDIF 332 332 333 IF(lwp .AND. lflush) CALL flush(numout) 333 334 ! ! allocate sbccpl arrays 334 335 !IF( sbc_cpl_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_cpl_alloc : unable to allocate arrays' ) … … 480 481 IF(lwp) WRITE(numout,*) 481 482 IF(lwp) WRITE(numout,*) ' runoffs received from oasis -> force ln_rnf = ', ln_rnf 483 IF(lwp .AND. lflush) CALL flush(numout) 482 484 ENDIF 483 485 ! … … 630 632 WRITE(numout,*)' wind stress module' 631 633 WRITE(numout,*) 634 IF(lflush) CALL flush(numout) 632 635 ENDIF 633 636 ENDIF … … 677 680 WRITE(numout,*)' fraction of solar net radiation absorbed in the first ocean level' 678 681 WRITE(numout,*) 682 IF(lflush) CALL flush(numout) 679 683 ENDIF 680 684 ENDIF … … 907 911 WRITE(numout,*)' fraction of solar net radiation absorbed in the first ocean level' 908 912 WRITE(numout,*) 913 IF(lflush) CALL flush(numout) 909 914 ENDIF 910 915 ENDIF … … 947 952 WRITE(numout,*)' wind stress U,V components' 948 953 WRITE(numout,*)' wind stress module' 954 IF(lflush) CALL flush(numout) 949 955 ENDIF 950 956 ENDIF … … 1104 1110 REAL(wp) :: zzx, zzy ! temporary variables 1105 1111 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty, zmsk, zemp, zqns, zqsr, ztx2, zty2 1112 LOGICAL :: ll_wrtstp !write diagnostics? 1106 1113 !!---------------------------------------------------------------------- 1107 1114 1115 ll_wrtstp = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 1108 1116 ! 1109 1117 IF( nn_timing.gt.0 .and. nn_timing .le. 2 ) CALL timing_start('sbc_cpl_rcv') … … 1467 1475 greenland_icesheet_timelapsed = 0.0_wp 1468 1476 ENDIF 1469 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass (kg) read in is ', zgreenland_icesheet_mass_in 1470 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass (kg) used is ', greenland_icesheet_mass 1471 IF(lwp) WRITE(numout,*) 'Greenland icesheet mass rate of change (kg/s) is ', greenland_icesheet_mass_rate_of_change 1472 IF(lwp) WRITE(numout,*) 'Greenland icesheet seconds lapsed since last change is ', greenland_icesheet_timelapsed 1477 IF(lwp .AND. ll_wrtstp) THEN 1478 WRITE(numout,*) 'Greenland icesheet mass (kg) read in is ', zgreenland_icesheet_mass_in 1479 WRITE(numout,*) 'Greenland icesheet mass (kg) used is ', greenland_icesheet_mass 1480 WRITE(numout,*) 'Greenland icesheet mass rate of change (kg/s) is ', greenland_icesheet_mass_rate_of_change 1481 WRITE(numout,*) 'Greenland icesheet seconds lapsed since last change is ', greenland_icesheet_timelapsed 1482 IF(lflush) CALL flush(numout) 1483 ENDIF 1473 1484 ELSE IF ( nn_coupled_iceshelf_fluxes == 2 ) THEN 1474 1485 greenland_icesheet_mass_rate_of_change = rn_greenland_total_fw_flux … … 1514 1525 antarctica_icesheet_timelapsed = 0.0_wp 1515 1526 ENDIF 1516 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass (kg) read in is ', zantarctica_icesheet_mass_in 1517 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass (kg) used is ', antarctica_icesheet_mass 1518 IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass rate of change (kg/s) is ', antarctica_icesheet_mass_rate_of_change 1519 IF(lwp) WRITE(numout,*) 'Antarctica icesheet seconds lapsed since last change is ', antarctica_icesheet_timelapsed 1527 IF(lwp .AND. ll_wrtstp) THEN 1528 WRITE(numout,*) 'Antarctica icesheet mass (kg) read in is ', zantarctica_icesheet_mass_in 1529 WRITE(numout,*) 'Antarctica icesheet mass (kg) used is ', antarctica_icesheet_mass 1530 WRITE(numout,*) 'Antarctica icesheet mass rate of change (kg/s) is ', antarctica_icesheet_mass_rate_of_change 1531 WRITE(numout,*) 'Antarctica icesheet seconds lapsed since last change is ', antarctica_icesheet_timelapsed 1532 IF(lflush) CALL flush(numout) 1533 ENDIF 1520 1534 ELSE IF ( nn_coupled_iceshelf_fluxes == 2 ) THEN 1521 1535 antarctica_icesheet_mass_rate_of_change = rn_antarctica_total_fw_flux -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90
r6486 r11101 99 99 WRITE(numout,*) '~~~~~~~' 100 100 WRITE(numout,*) 101 IF(lflush) CALL flush(numout) 101 102 ENDIF 102 103 ! allocate sbcdcy arrays -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90
r6486 r11101 98 98 READ ( numnam_cfg, namsbc_flx, IOSTAT = ios, ERR = 902 ) 99 99 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_flx in configuration namelist', lwp ) 100 IF(lwm ) WRITE ( numond, namsbc_flx )100 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_flx ) 101 101 ! 102 102 ! ! check: do we plan to use ln_dm2dc with non-daily forcing? … … 159 159 CALL lbc_lnk( taum(:,:), 'T', 1. ) ; CALL lbc_lnk( wndm(:,:), 'T', 1. ) 160 160 161 IF( n itend-nit000 <= 100.AND. lwp ) THEN ! control print (if less than 100 time-step asked)161 IF( nprint > 2 .AND. lwp ) THEN ! control print (if less than 100 time-step asked) 162 162 WRITE(numout,*) 163 163 WRITE(numout,*) ' read daily momentum, heat and freshwater fluxes OK' … … 170 170 CALL prihre( sf(jf)%fnow, jpi, jpj, 1, jpi, 20, 1, jpj, 10, zfact, numout ) 171 171 END DO 172 CALL FLUSH(numout)172 IF(lflush) CALL flush(numout) 173 173 ENDIF 174 174 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90
r6487 r11101 86 86 IF( kn_fwb == 2 ) WRITE(numout,*) ' adjusted from previous year budget' 87 87 IF( kn_fwb == 3 ) WRITE(numout,*) ' fwf set to zero and spread out over erp area' 88 IF(lflush) CALL flush(numout) 88 89 ENDIF 89 90 ! … … 127 128 IF(lwp)WRITE(numout,*)' year = ',iyear-1, ' freshwater budget read = ', a_fwb 128 129 IF(lwp)WRITE(numout,*)' year = ',iyear-2, ' freshwater budget read = ', a_fwb_b 130 IF(lwp .AND. lflush) CALL flush(numout) 129 131 ENDIF 130 132 ! ! Update fwfold if new year start … … 188 190 erp(:,:) = erp(:,:) + zerp_cor(:,:) 189 191 ! 190 IF( nprint == 1.AND. lwp ) THEN ! control print192 IF( nprint > 3 .AND. lwp ) THEN ! control print 191 193 IF( z_fwf < 0._wp ) THEN 192 194 WRITE(numout,*)' z_fwf < 0' … … 201 203 WRITE(numout,*)' MIN(zerp_cor) = ', MINVAL(zerp_cor) 202 204 WRITE(numout,*)' MAX(zerp_cor) = ', MAXVAL(zerp_cor) 205 IF(lflush) CALL flush(numout) 203 206 ENDIF 204 207 ENDIF -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r9975 r11101 182 182 CALL wrk_alloc( jpi,jpj, ztmp1, ztmp2 ) 183 183 ! 184 IF(lwp) WRITE(numout,*)'cice_sbc_init' 184 IF(lwp .AND. nprint>1) THEN 185 WRITE(numout,*)'cice_sbc_init' 186 IF(lflush) CALL flush(numout) 187 ENDIF 185 188 186 189 ji_off = INT ( (jpiglo - nx_global) / 2 ) … … 342 345 343 346 IF( kt == nit000 ) THEN 344 IF(lwp) WRITE(numout,*)'cice_sbc_in' 347 IF(lwp .AND. nprint>1) THEN 348 WRITE(numout,*)'cice_sbc_in' 349 IF(lflush) CALL flush(numout) 350 ENDIF 345 351 ENDIF 346 352 … … 583 589 584 590 IF( kt == nit000 ) THEN 585 IF(lwp) WRITE(numout,*)'cice_sbc_out' 591 IF(lwp .AND. nprint>1) THEN 592 WRITE(numout,*)'cice_sbc_out' 593 IF(lflush) CALL flush(numout) 594 ENDIF 586 595 ENDIF 587 596 … … 820 829 !!--------------------------------------------------------------------- 821 830 822 IF(lwp) WRITE(numout,*)'cice_sbc_final' 831 IF(lwp .AND. nprint > 1) THEN 832 WRITE(numout,*)'cice_sbc_final' 833 IF(lflush) CALL flush(numout) 834 ENDIF 823 835 824 836 CALL CICE_Finalize … … 895 907 READ ( numnam_cfg, namsbc_cice, IOSTAT = ios, ERR = 902 ) 896 908 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cice in configuration namelist', lwp ) 897 IF(lwm ) WRITE ( numond, namsbc_cice )909 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_cice ) 898 910 899 911 ! store namelist information in an array … … 943 955 WRITE(numout,*) 944 956 WRITE(numout,*) ' read forcing fluxes for CICE OK' 945 CALL FLUSH(numout)957 IF(lflush) CALL flush(numout) 946 958 ENDIF 947 959 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_if.F90
r6498 r11101 82 82 READ ( numnam_cfg, namsbc_iif, IOSTAT = ios, ERR = 902 ) 83 83 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_iif in configuration namelist', lwp ) 84 IF(lwm ) WRITE ( numond, namsbc_iif )84 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_iif ) 85 85 86 86 ALLOCATE( sf_ice(1), STAT=ierror ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90
r8046 r11101 91 91 CHARACTER (LEN=32) :: cvarLeff ! variable name for efficient Length scale 92 92 INTEGER :: ios ! Local integer output status for namelist read 93 LOGICAL :: ll_wrtstp !write diagnostics? 93 94 94 95 REAL(wp), DIMENSION(:,:,:), POINTER :: zfwfisf3d, zqhcisf3d, zqlatisf3d … … 100 101 ! 101 102 ! 103 ll_wrtstp = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 102 104 ! ! ====================== ! 103 105 IF( kt == nit000 ) THEN ! First call kt=nit000 ! … … 110 112 READ ( numnam_cfg, namsbc_isf, IOSTAT = ios, ERR = 902 ) 111 113 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_isf in configuration namelist', lwp ) 112 IF(lwm ) WRITE ( numond, namsbc_isf )114 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_isf ) 113 115 114 116 … … 124 126 IF ( lwp ) WRITE(numout,*) ' nn_gammablk = ', nn_gammablk 125 127 IF ( lwp ) WRITE(numout,*) ' rn_tfri2 = ', rn_tfri2 128 IF(lwp .AND. lflush) CALL flush(numout) 126 129 IF (ln_divisf) THEN ! keep it in the namelist ??? used true anyway as for runoff ? (PM) 127 130 rdivisf = 1._wp … … 271 274 272 275 ! check 273 IF(lwp ) WRITE(numout, *) 'Greenland iceshelf melting climatology (kg/s) : ',zgreenland_fwfisf_sum276 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Greenland iceshelf melting climatology (kg/s) : ',zgreenland_fwfisf_sum 274 277 275 278 zgreenland_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * greenland_icesheet_mask(:,:) ) 276 279 277 IF(lwp ) WRITE(numout, *) 'Greenland iceshelf melting adjusted value (kg/s) : ',zgreenland_fwfisf_sum280 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Greenland iceshelf melting adjusted value (kg/s) : ',zgreenland_fwfisf_sum 278 281 279 282 zantarctica_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * antarctica_icesheet_mask(:,:) ) … … 285 288 286 289 ! check 287 IF(lwp ) WRITE(numout, *) 'Antarctica iceshelf melting climatology (kg/s) : ',zantarctica_fwfisf_sum290 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Antarctica iceshelf melting climatology (kg/s) : ',zantarctica_fwfisf_sum 288 291 289 292 zantarctica_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * antarctica_icesheet_mask(:,:) ) 290 293 291 IF(lwp ) WRITE(numout, *) 'Antarctica iceshelf melting adjusted value (kg/s) : ',zantarctica_fwfisf_sum294 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Antarctica iceshelf melting adjusted value (kg/s) : ',zantarctica_fwfisf_sum 292 295 293 296 ENDIF … … 320 323 321 324 ! check 322 IF(lwp ) WRITE(numout, *) 'Greenland iceshelf melting climatology (kg/s) : ',zgreenland_fwfisf_sum325 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Greenland iceshelf melting climatology (kg/s) : ',zgreenland_fwfisf_sum 323 326 324 327 zgreenland_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * greenland_icesheet_mask(:,:) ) 325 328 326 IF(lwp ) WRITE(numout, *) 'Greenland iceshelf melting adjusted value (kg/s) : ',zgreenland_fwfisf_sum329 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Greenland iceshelf melting adjusted value (kg/s) : ',zgreenland_fwfisf_sum 327 330 328 331 zantarctica_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * antarctica_icesheet_mask(:,:) ) … … 334 337 335 338 ! check 336 IF(lwp ) WRITE(numout, *) 'Antarctica iceshelf melting climatology (kg/s) : ',zantarctica_fwfisf_sum339 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Antarctica iceshelf melting climatology (kg/s) : ',zantarctica_fwfisf_sum 337 340 338 341 zantarctica_fwfisf_sum = glob_sum( fwfisf(:,:) * e1t(:,:) * e2t(:,:) * antarctica_icesheet_mask(:,:) ) 339 342 340 IF(lwp ) WRITE(numout, *) 'Antarctica iceshelf melting adjusted value (kg/s) : ',zantarctica_fwfisf_sum343 IF(lwp .AND. ll_wrtstp) WRITE(numout, *) 'Antarctica iceshelf melting adjusted value (kg/s) : ',zantarctica_fwfisf_sum 341 344 342 345 ENDIF … … 409 412 & iom_varid( numror, 'fwf_isf_b', ldstop = .FALSE. ) > 0 ) THEN 410 413 IF(lwp) WRITE(numout,*) ' nit000-1 isf tracer content forcing fields read in the restart file' 414 IF(lwp .AND. lflush) CALL flush(numout) 411 415 CALL iom_get( numror, jpdom_autoglo, 'fwf_isf_b', fwfisf_b(:,:) ) ! before salt content isf_tsc trend 412 416 CALL iom_get( numror, jpdom_autoglo, 'isf_sc_b', risf_tsc_b(:,:,jp_sal) ) ! before salt content isf_tsc trend -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90
r9321 r11101 99 99 WRITE(numout,*) 'sbc_init : surface boundary condition setting' 100 100 WRITE(numout,*) '~~~~~~~~ ' 101 IF(lflush) CALL flush(numout) 101 102 ENDIF 102 103 … … 108 109 READ ( numnam_cfg, namsbc, IOSTAT = ios, ERR = 902 ) 109 110 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc in configuration namelist', lwp ) 110 IF(lwm ) WRITE ( numond, namsbc )111 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc ) 111 112 112 113 ! ! overwrite namelist parameter using CPP key information … … 160 161 END SELECT 161 162 ! 163 IF(lwp .AND. lflush) CALL flush(numout) 164 ! 162 165 IF ( nn_components /= jp_iam_nemo .AND. .NOT. lk_oasis ) & 163 166 & CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but key_oasis3 disabled' ) … … 264 267 IF( nn_components/= jp_iam_nemo ) & 265 268 & WRITE(numout,*) ' + OASIS coupled SAS' 269 IF(lflush) CALL flush(numout) 266 270 ENDIF 267 271 ! … … 283 287 WRITE(numout,*)" OPA-SAS coupled via OASIS : nn_fsbc re-defined from OASIS namcouple ", nn_fsbc 284 288 WRITE(numout,*) 289 IF(lflush) CALL flush(numout) 285 290 ENDIF 286 291 ENDIF … … 421 426 IF( ln_rstart .AND. & !* Restart: read in restart file 422 427 & iom_varid( numror, 'utau_b', ldstop = .FALSE. ) > 0 ) THEN 423 IF(lwp) WRITE(numout,*) ' nit000-1 surface forcing fields red in the restart file' 428 IF(lwp .AND. nprint > 0) THEN 429 WRITE(numout,*) ' nit000-1 surface forcing fields red in the restart file' 430 IF(lflush) CALL flush(numout) 431 ENDIF 424 432 IF(nn_timing == 2) CALL timing_start('iom_rstget') 425 433 CALL iom_get( numror, jpdom_autoglo, 'utau_b', utau_b ) ! before i-stress (U-point) … … 448 456 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 449 457 ! ! ---------------------------------------- ! 450 IF(lwp) WRITE(numout,*) 451 IF(lwp) WRITE(numout,*) 'sbc : ocean surface forcing fields written in ocean restart file ', & 458 IF(lwp .AND. nprint > 0) THEN 459 WRITE(numout,*) 460 WRITE(numout,*) 'sbc : ocean surface forcing fields written in ocean restart file ', & 452 461 & 'at it= ', kt,' date= ', ndastp 453 IF(lwp) WRITE(numout,*) '~~~~' 462 WRITE(numout,*) '~~~~' 463 IF(lflush) CALL flush(numout) 464 ENDIF 454 465 IF(nn_timing == 2) CALL timing_start('iom_rstput') 455 466 CALL iom_rstput( kt, nitrst, numrow, 'utau_b' , utau ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r9321 r11101 148 148 IF( ln_rstart .AND. & !* Restart: read in restart file 149 149 & iom_varid( numror, 'rnf_b', ldstop = .FALSE. ) > 0 ) THEN 150 IF(lwp ) WRITE(numout,*) ' nit000-1 runoff forcing fields red in the restart file'150 IF(lwp .AND. nprint > 0) WRITE(numout,*) ' nit000-1 runoff forcing fields red in the restart file' 151 151 IF(nn_timing == 2) CALL timing_start('iom_rstget') 152 152 CALL iom_get( numror, jpdom_autoglo, 'rnf_b', rnf_b ) ! before runoff … … 159 159 rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 160 160 ENDIF 161 162 IF(lwp .AND. lflush) CALL flush(numout) 163 161 164 ENDIF 162 165 ! ! ---------------------------------------- ! 163 166 IF( lrst_oce ) THEN ! Write in the ocean restart file ! 164 167 ! ! ---------------------------------------- ! 165 IF(lwp) WRITE(numout,*) 166 IF(lwp) WRITE(numout,*) 'sbcrnf : runoff forcing fields written in ocean restart file ', & 168 IF(lwp .AND. nprint > 0) THEN 169 WRITE(numout,*) 170 WRITE(numout,*) 'sbcrnf : runoff forcing fields written in ocean restart file ', & 167 171 & 'at it= ', kt,' date= ', ndastp 168 IF(lwp) WRITE(numout,*) '~~~~' 172 WRITE(numout,*) '~~~~' 173 IF(lflush) CALL flush(numout) 174 ENDIF 169 175 IF(nn_timing == 2) CALL timing_start('iom_rstput') 170 176 CALL iom_rstput( kt, nitrst, numrow, 'rnf_b' , rnf ) … … 282 288 READ ( numnam_cfg, namsbc_rnf, IOSTAT = ios, ERR = 902 ) 283 289 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_rnf in configuration namelist', lwp ) 284 IF(lwm ) WRITE ( numond, namsbc_rnf )290 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_rnf ) 285 291 ! 286 292 ! ! Control print … … 294 300 WRITE(numout,*) ' depth of river mouth additional mixing rn_hrnf = ', rn_hrnf 295 301 WRITE(numout,*) ' multiplicative factor for runoff rn_rfact = ', rn_rfact 302 IF(lflush) CALL flush(numout) 296 303 ENDIF 297 304 ! ! ================== … … 303 310 IF(lwp) WRITE(numout,*) 304 311 IF(lwp) WRITE(numout,*) ' runoffs inflow read in a file' 312 IF(lwp .AND. lflush) CALL flush(numout) 305 313 IF( ierror > 0 ) THEN 306 314 CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' ) ; RETURN … … 314 322 IF(lwp) WRITE(numout,*) 315 323 IF(lwp) WRITE(numout,*) ' runoffs temperatures read in a file' 324 IF(lwp .AND. lflush) CALL flush(numout) 316 325 ALLOCATE( sf_t_rnf(1), STAT=ierror ) 317 326 IF( ierror > 0 ) THEN … … 326 335 IF(lwp) WRITE(numout,*) 327 336 IF(lwp) WRITE(numout,*) ' runoffs salinities read in a file' 337 IF(lwp .AND. lflush) CALL flush(numout) 328 338 ALLOCATE( sf_s_rnf(1), STAT=ierror ) 329 339 IF( ierror > 0 ) THEN … … 338 348 IF(lwp) WRITE(numout,*) 339 349 IF(lwp) WRITE(numout,*) ' runoffs depth read in a file' 350 IF(lwp .AND. lflush) CALL flush(numout) 340 351 rn_dep_file = TRIM( cn_dir )//TRIM( sn_dep_rnf%clname ) 341 352 IF( .NOT. sn_dep_rnf%ln_clim ) THEN ; WRITE(rn_dep_file, '(a,"_y",i4)' ) TRIM( rn_dep_file ), nyear ! add year … … 378 389 IF(lwp) WRITE(numout,*) ' depth over which runoffs is spread rn_dep_max = ', rn_dep_max 379 390 IF(lwp) WRITE(numout,*) ' create (=1) a runoff depth file or not (=0) nn_rnf_depth_file = ', nn_rnf_depth_file 391 IF(lwp .AND. lflush) CALL flush(numout) 380 392 381 393 CALL iom_open( TRIM( sn_rnf%clname ), inum ) ! open runoff file … … 431 443 IF( nn_rnf_depth_file == 1 ) THEN ! save output nb levels for runoff 432 444 IF(lwp) WRITE(numout,*) ' create runoff depht file' 445 IF(lwp .AND. lflush) CALL flush(numout) 433 446 CALL iom_open ( TRIM( sn_dep_rnf%clname ), inum, ldwrt = .TRUE., kiolib = jprstlib ) 434 447 CALL iom_rstput( 0, 0, inum, 'rodepth', h_rnf ) … … 468 481 IF(lwp) WRITE(numout,*) ' - set to zero SSS damping (if ln_ssr=T)' 469 482 IF(lwp) WRITE(numout,*) ' - mixed upstream-centered (if ln_traadv_cen2=T)' 483 IF(lwp .AND. lflush) CALL flush(numout) 470 484 ! 471 485 CALL rnf_mouth ! set river mouth mask … … 474 488 IF(lwp) WRITE(numout,*) 475 489 IF(lwp) WRITE(numout,*) ' No specific treatment at river mouths' 490 IF(lwp .AND. lflush) CALL flush(numout) 476 491 rnfmsk (:,:) = 0._wp 477 492 rnfmsk_z(:) = 0._wp … … 509 524 IF(lwp) WRITE(numout,*) 'rnf_mouth : river mouth mask' 510 525 IF(lwp) WRITE(numout,*) '~~~~~~~~~ ' 526 IF(lwp .AND. lflush) CALL flush(numout) 511 527 ! 512 528 cl_rnfile = TRIM( cn_dir )//TRIM( sn_cnf%clname ) -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r9321 r11101 92 92 IF(lwp) WRITE(numout,*) 93 93 IF(lwp) WRITE(numout,*) '~~~~~~~ mean fields initialised to instantaneous values' 94 IF(lwp .AND. lflush) CALL flush(numout) 94 95 zcoef = REAL( nn_fsbc - 1, wp ) 95 96 ssu_m(:,:) = zcoef * ub(:,:,1) … … 156 157 & 'at it= ', kt,' date= ', ndastp 157 158 IF(lwp) WRITE(numout,*) '~~~~~~~' 159 IF(lwp .AND. lflush) CALL flush(numout) 158 160 zf_sbc = REAL( nn_fsbc, wp ) 159 161 IF(nn_timing == 2) CALL timing_start('iom_rstput') … … 257 259 ENDIF 258 260 ! 261 IF(lwp .AND. lflush) CALL flush(numout) 262 ! 259 263 END SUBROUTINE sbc_ssm_init 260 264 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90
r6486 r11101 175 175 READ ( numnam_cfg, namsbc_ssr, IOSTAT = ios, ERR = 902 ) 176 176 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ssr in configuration namelist', lwp ) 177 IF(lwm ) WRITE ( numond, namsbc_ssr )177 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_ssr ) 178 178 179 179 IF(lwp) THEN !* control print … … 189 189 WRITE(numout,*) ' flag to bound erp term ln_sssr_bnd = ', ln_sssr_bnd 190 190 WRITE(numout,*) ' ABS(Max./Min.) erp threshold rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day' 191 IF(lflush) CALL flush(numout) 191 192 ENDIF 192 193 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90
r6486 r11101 71 71 WRITE(numout,*) Wave(ntide(jk))%cname_tide, utide(jk), ftide(jk), v0tide(jk), omega_tide(jk) 72 72 END DO 73 IF(lflush) CALL flush(numout) 73 74 ENDIF 74 75 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90
r6486 r11101 90 90 READ ( numnam_cfg, namsbc_wave, IOSTAT = ios, ERR = 902 ) 91 91 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_wave in configuration namelist', lwp ) 92 IF(lwm ) WRITE ( numond, namsbc_wave )92 IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_wave ) 93 93 ! 94 94 -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90
r6486 r11101 60 60 WRITE(numout,*) 'tide_init : Initialization of the tidal components' 61 61 WRITE(numout,*) '~~~~~~~~~ ' 62 IF(lflush) CALL flush(numout) 62 63 ENDIF 63 64 ! … … 72 73 READ ( numnam_cfg, nam_tide, IOSTAT = ios, ERR = 902 ) 73 74 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tide in configuration namelist', lwp ) 74 IF(lwm ) WRITE ( numond, nam_tide )75 IF(lwm .AND. nprint > 2) WRITE ( numond, nam_tide ) 75 76 ! 76 77 nb_harmo=0 … … 90 91 WRITE(numout,*) ' ln_tide_ramp = ', ln_tide_ramp 91 92 WRITE(numout,*) ' rdttideramp = ', rdttideramp 93 IF(lflush) CALL flush(numout) 92 94 ENDIF 93 95 IF( ln_tide_ramp.AND.((nitend-nit000+1)*rdt/rday < rdttideramp) ) & -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SOL/solver.F90
r6498 r11101 59 59 ! 60 60 61 IF(lwp ) THEN!* open elliptic solver statistics file (only on the printing processors)61 IF(lwp .AND. (ln_ctl .OR. sn_cfctl%l_runstat)) THEN !* open elliptic solver statistics file (only on the printing processors) 62 62 CALL ctl_opn( numsol, 'solver.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 63 63 ENDIF … … 70 70 READ ( numnam_cfg, namsol, IOSTAT = ios, ERR = 902 ) 71 71 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsol in configuration namelist', lwp ) 72 IF(lwm ) WRITE ( numond, namsol )72 IF(lwm .AND. nprint > 2) WRITE ( numond, namsol ) 73 73 74 74 IF(lwp) THEN !* Control print … … 86 86 WRITE(numout,*) ' optimal coefficient of sor rn_sor = ', rn_sor 87 87 WRITE(numout,*) 88 IF(lflush) CALL flush(numout) 88 89 ENDIF 89 90 eps = rn_eps … … 100 101 CASE ( 1 ) ! preconditioned conjugate gradient solver 101 102 IF(lwp) WRITE(numout,*) ' a preconditioned conjugate gradient solver is used' 103 IF(lwp .AND. lflush) CALL flush(numout) 102 104 IF( jpr2di /= 0 .AND. jpr2dj /= 0 ) CALL ctl_stop( ' jpr2di and jpr2dj should be equal to zero' ) 103 105 ! … … 105 107 IF(lwp) WRITE(numout,*) ' a successive-over-relaxation solver with extra outer halo is used' 106 108 IF(lwp) WRITE(numout,*) ' with jpr2di =', jpr2di, ' and jpr2dj =', jpr2dj 109 IF(lwp .AND. lflush) CALL flush(numout) 107 110 IF( .NOT. lk_mpp .AND. jpr2di /= 0 .AND. jpr2dj /= 0 ) THEN 108 111 CALL ctl_stop( 'jpr2di and jpr2dj are not equal to zero', & -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/STO/stopar.F90
r9321 r11101 268 268 READ ( numnam_cfg, namsto, IOSTAT = ios, ERR = 902 ) 269 269 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsto in configuration namelist', lwp ) 270 IF(lwm ) WRITE ( numond, namsto )270 IF(lwm .AND. nprint > 2) WRITE ( numond, namsto ) 271 271 272 272 !IF(ln_ens_rst_in) cn_storst_in = cn_mem//cn_storst_in -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r8447 r11101 1241 1241 READ ( numnam_cfg, nameos, IOSTAT = ios, ERR = 902 ) 1242 1242 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nameos in configuration namelist', lwp ) 1243 IF(lwm ) WRITE( numond, nameos )1243 IF(lwm .AND. nprint > 2) WRITE( numond, nameos ) 1244 1244 ! 1245 1245 rau0 = 1026._wp !: volumic mass of reference [kg/m3] … … 1258 1258 WRITE(numout,*) ' model does not use Conservative Temperature' 1259 1259 ENDIF 1260 IF(lflush) CALL flush(numout) 1260 1261 ENDIF 1261 1262 ! … … 1272 1273 IF(lwp) WRITE(numout,*) 1273 1274 IF(lwp) WRITE(numout,*) ' use of TEOS-10 equation of state (cons. temp. and abs. salinity)' 1275 IF(lwp .AND. lflush) CALL flush(numout) 1274 1276 ! 1275 1277 rdeltaS = 32._wp … … 1462 1464 IF(lwp) WRITE(numout,*) 1463 1465 IF(lwp) WRITE(numout,*) ' use of EOS-80 equation of state (pot. temp. and pract. salinity)' 1466 IF(lwp .AND. lflush) CALL flush(numout) 1464 1467 ! 1465 1468 rdeltaS = 20._wp … … 1662 1665 WRITE(numout,*) ' 2nd cabbel. coef. rn_nu = ', rn_nu 1663 1666 WRITE(numout,*) ' Caution: rn_beta0=0 incompatible with ddm parameterization ' 1667 IF(lflush) CALL flush(numout) 1664 1668 ENDIF 1665 1669 ! … … 1681 1685 IF(lwp) WRITE(numout,*) ' rau0 * rcp rau0_rcp = ', rau0_rcp 1682 1686 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) r1_rau0_rcp = ', r1_rau0_rcp 1687 IF(lwp .AND. lflush) CALL flush(numout) 1683 1688 ! 1684 1689 END SUBROUTINE eos_init -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90
r7771 r11101 208 208 READ ( numnam_cfg, namtra_adv, IOSTAT = ios, ERR = 902 ) 209 209 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_adv in configuration namelist', lwp ) 210 IF(lwm ) WRITE ( numond, namtra_adv )210 IF(lwm .AND. nprint > 2) WRITE ( numond, namtra_adv ) 211 211 212 212 IF(lwp) THEN ! Namelist print … … 223 223 WRITE(numout,*) ' upstream scheme within muscl ln_traadv_msc_ups = ', ln_traadv_msc_ups 224 224 WRITE(numout,*) ' TVD advection scheme with zts ln_traadv_tvd_zts = ', ln_traadv_tvd_zts 225 IF(lflush) CALL flush(numout) 225 226 ENDIF 226 227 … … 260 261 IF( nadv == 7 ) WRITE(numout,*) ' TVD ZTS scheme is used' 261 262 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection scheme' 263 IF(lflush) CALL flush(numout) 262 264 ENDIF 263 265 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen2.F90
r9321 r11101 139 139 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ ' 140 140 IF(lwp) WRITE(numout,*) 141 IF(lwp .AND. lflush) CALL flush(numout) 141 142 ! 142 143 IF( .NOT. ALLOCATED( upsmsk ) ) THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90
r7179 r11101 97 97 IF(lwp) WRITE(numout,*) 'tra_adv_eiv : eddy induced advection on ', cdtype,' :' 98 98 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ add to velocity fields the eiv component' 99 IF(lwp .AND. lflush) CALL flush(numout) 99 100 # if defined key_diaeiv 100 101 IF( cdtype == 'TRA') THEN -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90
r6486 r11101 290 290 READ ( numnam_cfg, namtra_adv_mle, IOSTAT = ios, ERR = 902 ) 291 291 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_adv_mle in configuration namelist', lwp ) 292 IF(lwm ) WRITE ( numond, namtra_adv_mle )292 IF(lwm .AND. nprint > 2) WRITE ( numond, namtra_adv_mle ) 293 293 294 294 IF(lwp) THEN ! Namelist print … … 306 306 WRITE(numout,*) ' =1 no MLE in case of convection ; =0 always MLE nn_conv = ', nn_conv 307 307 WRITE(numout,*) ' Density difference used to define ML for FK rn_rho_c_mle = ', rn_rho_c_mle 308 IF(lflush) CALL flush(numout) 308 309 ENDIF 309 310 ! … … 317 318 WRITE(numout,*) ' Mixed Layer Eddy parametrisation NOT used' 318 319 ENDIF 320 321 IF(lflush) CALL flush(numout) 322 319 323 ENDIF 320 324 ! … … 325 329 IF(lwp) WRITE(numout,*) ' ML buoyancy criteria = ', rb_c, ' m/s2 ' 326 330 IF(lwp) WRITE(numout,*) ' associated ML density criteria defined in zdfmxl = ', rho_c, 'kg/m3' 331 IF(lwp .AND. lflush) CALL flush(numout) 327 332 ! 328 333 IF( nn_mle == 0 ) THEN ! MLE array allocation & initialisation -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90
r7771 r11101 99 99 IF(lwp) WRITE(numout,*) '~~~~~~~' 100 100 IF(lwp) WRITE(numout,*) 101 IF(lwp .AND. lflush) CALL flush(numout) 101 102 ! 102 103 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90
r7179 r11101 84 84 IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype 85 85 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 86 IF(lwp .AND. lflush) CALL flush(numout) 86 87 ENDIF 87 88 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90
r7179 r11101 101 101 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 102 102 IF(lwp) WRITE(numout,*) 103 IF(lwp .AND. lflush) CALL flush(numout) 103 104 ENDIF 104 105 l_trd = .FALSE. -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90
r10020 r11101 102 102 IF(lwp) WRITE(numout,*) 'tra_adv_tvd : TVD advection scheme on ', cdtype 103 103 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 104 IF(lwp .AND. lflush) CALL flush(numout) 104 105 ! 105 106 ENDIF … … 383 384 IF(lwp) WRITE(numout,*) 'tra_adv_tvd_zts : TVD ZTS advection scheme on ', cdtype 384 385 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 386 IF(lwp .AND. lflush) CALL flush(numout) 385 387 ENDIF 386 388 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90
r7179 r11101 101 101 IF(lwp) WRITE(numout,*) 'tra_adv_ubs : horizontal UBS advection scheme on ', cdtype 102 102 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 103 IF(lwp .AND. lflush) CALL flush(numout) 103 104 ENDIF 104 105 ! -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90
r6486 r11101 147 147 READ ( numnam_cfg, nambbc, IOSTAT = ios, ERR = 902 ) 148 148 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbc in configuration namelist', lwp ) 149 IF(lwm ) WRITE ( numond, nambbc )149 IF(lwm .AND. nprint > 2) WRITE ( numond, nambbc ) 150 150 151 151 IF(lwp) THEN ! Control print … … 158 158 WRITE(numout,*) ' Constant geothermal flux value rn_geoflx_cst = ', rn_geoflx_cst 159 159 WRITE(numout,*) 160 IF(lflush) CALL flush(numout) 160 161 ENDIF 161 162 … … 196 197 IF(lwp) WRITE(numout,*) ' *** no geothermal heat flux' 197 198 ENDIF 199 IF(lwp .AND. lflush) CALL flush(numout) 198 200 ! 199 201 END SUBROUTINE tra_bbc_init -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90
r10681 r11101 358 358 IF(lwp) WRITE(numout,*) 'trabbl:bbl : Compute bbl velocities and diffusive coefficients in ', cdtype 359 359 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 360 IF(lflush) CALL flush(numout) 360 361 ENDIF 361 362 ! !* bottom variables (T, S, alpha, beta, depth, velocity) … … 517 518 READ ( numnam_cfg, nambbl, IOSTAT = ios, ERR = 902 ) 518 519 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbl in configuration namelist', lwp ) 519 IF(lwm ) WRITE ( numond, nambbl )520 IF(lwm .AND. nprint > 2) WRITE ( numond, nambbl ) 520 521 ! 521