- Timestamp:
- 2019-05-16T15:23:56+02:00 (5 years ago)
- Location:
- NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/dia25h.F90
r10968 r10986 67 67 WRITE(numout,*) ' Namelist nam_dia25h : set 25h outputs ' 68 68 WRITE(numout,*) ' Switch for 25h diagnostics (T) or not (F) ln_dia25h = ', ln_dia25h 69 IF(lflush) CALL FLUSH(numout) 69 70 ENDIF 70 71 IF( .NOT. ln_dia25h ) RETURN … … 159 160 WRITE(numout,*) 'dia_wri_tide : Summing instantaneous hourly diagnostics at timestep ',kt 160 161 WRITE(numout,*) '~~~~~~~~~~~~ ' 162 IF(lflush) CALL FLUSH(numout) 161 163 ENDIF 162 164 … … 180 182 IF (lwp) THEN 181 183 WRITE(numout,*) 'dia_tide : Summed the following number of hourly values so far',cnt_25h 184 IF(lflush) CALL FLUSH(numout) 182 185 ENDIF 183 186 ! … … 190 193 WRITE(numout,*) 'dia_wri_tide : Writing 25 hour mean tide diagnostics at timestep', kt 191 194 WRITE(numout,*) '~~~~~~~~~~~~ ' 195 IF(lflush) CALL FLUSH(numout) 192 196 ENDIF 193 197 ! … … 208 212 ENDIF 209 213 ! 210 IF(lwp) WRITE(numout,*) 'dia_wri_tide : Mean calculated by dividing 25 hour sums and writing output' 214 IF(lwp) THEN 215 WRITE(numout,*) 'dia_wri_tide : Mean calculated by dividing 25 hour sums and writing output' 216 IF(lflush) CALL FLUSH(numout) 217 ENDIF 211 218 zmdi=1.e+20 !missing data indicator for masking 212 219 ! write tracers (instantaneous) … … 261 268 ENDIF 262 269 cnt_25h = 1 263 IF(lwp) WRITE(numout,*) 'dia_wri_tide : & 270 IF(lwp) THEN 271 WRITE(numout,*) 'dia_wri_tide : & 264 272 & After 25hr mean write, reset sum to current value and cnt_25h to one for overlapping average', cnt_25h 273 IF(lflush) CALL FLUSH(numout) 274 ENDIF 265 275 ENDIF ! cnt_25h .EQ. 25 .AND. MOD( kt, i_steps * 24) == 0 .AND. kt .NE. nn_it000 266 276 ! -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diacfl.F90
r10888 r10986 142 142 WRITE(numout,*) ' Max Cv = ', rCv_max, ' at (i,j,k) = (',nCv_loc(1),nCv_loc(2),nCv_loc(3),') => dt/C = ', z2dt/rCv_max 143 143 WRITE(numout,*) ' Max Cw = ', rCw_max, ' at (i,j,k) = (',nCw_loc(1),nCw_loc(2),nCw_loc(3),') => dt/C = ', z2dt/rCw_max 144 IF(lflush) CALL FLUSH(numout) 144 145 ENDIF 145 146 ! … … 166 167 WRITE(numcfl,*) 'Timestep Direction Max C i j k' 167 168 WRITE(numcfl,*) '******************************************' 169 IF(lflush) CALL FLUSH(numout) 168 170 ENDIF 169 171 ! -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diadct.F90
r10968 r10986 158 158 IF(nn_dct .GE. nn_dctwri .AND. MOD(nn_dct,nn_dctwri) .NE. 0) & 159 159 & CALL ctl_stop( 'diadct: nn_dct should be smaller and a multiple of nn_dctwri' ) 160 160 IF(lflush) CALL FLUSH(numout) 161 161 ENDIF 162 162 … … 222 222 WRITE(numout,*) "~~~~~~~~~~~~~~~~~~~~~~~~~" 223 223 WRITE(numout,*) "nb_sec = ",nb_sec 224 IF(lflush) CALL FLUSH(numout) 224 225 ENDIF 225 226 … … 241 242 IF( MOD(kt,nn_dctwri)==0 )THEN 242 243 243 IF( kt==nit000+nn_dctwri-1 )WRITE(numout,*)" diadct: average transports and write at kt = ",kt 244 IF( kt==nit000+nn_dctwri-1 ) THEN 245 WRITE(numout,*)" diadct: average transports and write at kt = ",kt 246 IF(lflush) CALL FLUSH(numout) 247 ENDIF 244 248 245 249 !! divide arrays by nn_dctwri/nn_dct to obtain average … … 321 325 DO jsec=1,nb_sec_max !loop on the nb_sec sections 322 326 323 IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) & 324 & WRITE(numout,*)'debuging for section number: ',jsec 327 IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) THEN 328 WRITE(numout,*)'debuging for section number: ',jsec 329 IF(lflush) CALL FLUSH(numout) 330 ENDIF 325 331 326 332 !initialization … … 341 347 IF( jsec .NE. isec ) CALL ctl_stop( cltmp ) 342 348 343 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )WRITE(numout,*)"isec ",isec 349 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 ) THEN 350 WRITE(numout,*)"isec ",isec 351 IF(lflush) CALL FLUSH(numout) 352 ENDIF 344 353 345 354 READ(numdct_in)secs(jsec)%name … … 376 385 WRITE(numout,clformat)" Temperature classes : ",secs(jsec)%ztem 377 386 WRITE(numout,clformat)" Depth classes : ",secs(jsec)%zlay 387 IF(lflush) CALL FLUSH(numout) 378 388 ENDIF 379 389 … … 398 408 WRITE(numout,*)' # I J ',jpt,coordtemp(jpt),directemp(jpt) 399 409 ENDDO 410 IF(lflush) CALL FLUSH(numout) 400 411 ENDIF 401 412 … … 434 445 WRITE(numout,*)' # I J : ',iiglo,ijglo 435 446 ENDDO 447 IF(lflush) CALL FLUSH(numout) 436 448 ENDIF 437 449 … … 468 480 ijglo = secs(jsec)%listPoint(jpt)%J + njmpp - 1 469 481 WRITE(numout,*)' # I J : ',iiglo,ijglo 470 CALL FLUSH(numout)482 IF(lflush) CALL FLUSH(numout) 471 483 ENDDO 472 484 ENDIF 473 485 474 486 ELSE ! iptglo = 0 475 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 ) &487 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 ) THEN 476 488 WRITE(numout,*)' No points for this section.' 489 IF(lflush) CALL FLUSH(numout) 490 ENDIF 477 491 ENDIF 478 492 … … 510 524 !---------------------------------------------------------------------------- 511 525 ! 512 IF( ld_debug )WRITE(numout,*)' -------------------------' 513 IF( ld_debug )WRITE(numout,*)' removepoints in listpoint' 526 IF( ld_debug ) THEN 527 WRITE(numout,*)' -------------------------' 528 WRITE(numout,*)' removepoints in listpoint' 529 IF(lflush) CALL FLUSH(numout) 530 ENDIF 514 531 515 532 !iextr=extremity of list_point that we verify … … 529 546 WRITE(numout,*)' ', cdind,' ',cdextr,' ',itest 530 547 WRITE(numout,*)' Actual number of points: ',sec%nb_point 548 IF(lflush) CALL FLUSH(numout) 531 549 ENDIF 532 550 … … 556 574 WRITE(numout,*)' Number of points after removepoints :',sec%nb_point 557 575 WRITE(numout,*)' sec%direction after removepoints :',sec%direction(1:sec%nb_point) 576 IF(lflush) CALL FLUSH(numout) 558 577 ENDIF 559 578 ! … … 592 611 !!-------------------------------------------------------- 593 612 ! 594 IF( ld_debug )WRITE(numout,*)' Compute transport' 613 IF( ld_debug ) THEN 614 WRITE(numout,*)' Compute transport' 615 IF(lflush) CALL FLUSH(numout) 616 ENDIF 595 617 596 618 !---------------------------! … … 626 648 IF( sec%slopeSection .GE. 9999. ) isgnv = 1 627 649 628 IF( ld_debug )write(numout,*)"sec%slopeSection isgnu isgnv ",sec%slopeSection,isgnu,isgnv 650 IF( ld_debug ) THEN 651 write(numout,*)"sec%slopeSection isgnu isgnv ",sec%slopeSection,isgnu,isgnv 652 IF(lflush) CALL FLUSH(numout) 653 ENDIF 629 654 630 655 !--------------------------------------! -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diaharm.F90
r10968 r10986 81 81 WRITE(numout,*) 'dia_harm_init: Tidal harmonic analysis initialization' 82 82 WRITE(numout,*) '~~~~~~~ ' 83 IF(lflush) CALL FLUSH(numout) 83 84 ENDIF 84 85 ! … … 99 100 WRITE(numout,*) 'Last time step used for analysis: nitend_han= ', nitend_han 100 101 WRITE(numout,*) 'Time step frequency for harmonic analysis: nstep_han= ', nstep_han 102 IF(lflush) CALL FLUSH(numout) 101 103 ENDIF 102 104 … … 123 125 WRITE(numout,*) ' Namelist nam_diaharm' 124 126 WRITE(numout,*) ' nb_ana = ', nb_ana 125 CALL flush(numout)127 IF(lflush) CALL FLUSH(numout) 126 128 ENDIF 127 129 ! … … 148 150 CALL tide_harmo( ana_freq, vt, ut, ft, name, nb_ana ) 149 151 150 IF(lwp) WRITE(numout,*) 'Analysed frequency : ',nb_ana ,'Frequency ' 151 152 DO jh = 1, nb_ana 153 IF(lwp) WRITE(numout,*) ' : ',tname(jh),' ',ana_freq(jh) 154 END DO 152 IF(lwp) THEN 153 WRITE(numout,*) 'Analysed frequency : ',nb_ana ,'Frequency ' 154 DO jh = 1, nb_ana 155 WRITE(numout,*) ' : ',tname(jh),' ',ana_freq(jh) 156 END DO 157 IF(lflush) CALL FLUSH(numout) 158 ENDIF 155 159 156 160 ! Initialize temporary arrays: … … 228 232 !!-------------------------------------------------------------------- 229 233 ! 230 IF(lwp) WRITE(numout,*) 231 IF(lwp) WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 232 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 234 IF(lwp) THEN 235 WRITE(numout,*) 236 WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 237 WRITE(numout,*) '~~~~~~~~~~~~' 238 IF(lflush) CALL FLUSH(numout) 239 ENDIF 233 240 234 241 ztime_ini = nit000_han*rdt ! Initial time in seconds at the beginning of analysis … … 382 389 !!---------------------------------------------------------------------- 383 390 384 IF(lwp) WRITE(numout,*) ' ' 385 IF(lwp) WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 386 IF(lwp) WRITE(numout,*) ' ' 391 IF(lwp) THEN 392 WRITE(numout,*) ' ' 393 WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 394 WRITE(numout,*) ' ' 395 IF(lflush) CALL FLUSH(numout) 396 ENDIF 387 397 388 398 ! A) Elevation -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diahsb.F90
r10968 r10986 211 211 CALL iom_put( 'bgvole3t' , zdiff_v2 * 1.e-9 ) ! volume e3t drift (km3) 212 212 ! 213 IF( kt == nitend .AND. lwp ) THEN213 IF( kt == nitend .AND. lwp .AND. nprint > 0) THEN 214 214 WRITE(numout,*) 215 215 WRITE(numout,*) 'dia_hsb : last time step hsb diagnostics: at it= ', kt,' date= ', ndastp … … 219 219 WRITE(numout,*) ' volume ssh drift = ', zdiff_v1 * 1.e-9 , ' km^3' 220 220 WRITE(numout,*) ' volume e3t drift = ', zdiff_v2 * 1.e-9 , ' km^3' 221 IF(lflush) CALL FLUSH(numout) 221 222 ENDIF 222 223 ! … … 257 258 IF( ln_rstart ) THEN !* Read the restart file 258 259 ! 259 IF(lwp) WRITE(numout,*) 260 IF(lwp) WRITE(numout,*) ' dia_hsb_rst : read hsb restart at it= ', kt,' date= ', ndastp 261 IF(lwp) WRITE(numout,*) 260 IF(lwp) THEN 261 WRITE(numout,*) 262 WRITE(numout,*) ' dia_hsb_rst : read hsb restart at it= ', kt,' date= ', ndastp 263 WRITE(numout,*) 264 IF(lflush) CALL FLUSH(numout) 265 ENDIF 262 266 CALL iom_get( numror, 'frc_v', frc_v, ldxios = lrxios ) 263 267 CALL iom_get( numror, 'frc_t', frc_t, ldxios = lrxios ) … … 277 281 ENDIF 278 282 ELSE 279 IF(lwp) WRITE(numout,*) 280 IF(lwp) WRITE(numout,*) ' dia_hsb_rst : initialise hsb at initial state ' 281 IF(lwp) WRITE(numout,*) 283 IF(lwp) THEN 284 WRITE(numout,*) 285 WRITE(numout,*) ' dia_hsb_rst : initialise hsb at initial state ' 286 WRITE(numout,*) 287 ENDIF 282 288 surf_ini(:,:) = e1e2t(:,:) * tmask_i(:,:) ! initial ocean surface 283 289 ssh_ini(:,:) = sshn(:,:) ! initial ssh … … 310 316 ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file 311 317 ! ! ------------------- 312 IF(lwp .AND. nprint > 0) WRITE(numout,*) 313 IF(lwp .AND. nprint > 0) WRITE(numout,*) ' dia_hsb_rst : write restart at it= ', kt,' date= ', ndastp 314 IF(lwp .AND. nprint > 0) WRITE(numout,*) 318 IF(lwp .AND. nprint > 0) THEN 319 WRITE(numout,*) 320 WRITE(numout,*) ' dia_hsb_rst : write restart at it= ', kt,' date= ', ndastp 321 WRITE(numout,*) 322 IF(lflush) CALL FLUSH(numout) 323 ENDIF 315 324 ! 316 325 IF( lwxios ) CALL iom_swap( cwxios_context ) … … 359 368 WRITE(numout,*) 'dia_hsb_init : heat and salt budgets diagnostics' 360 369 WRITE(numout,*) '~~~~~~~~~~~~ ' 370 IF(lflush) CALL FLUSH(numout) 361 371 ENDIF 362 372 REWIND( numnam_ref ) ! Namelist namhsb in reference namelist … … 371 381 WRITE(numout,*) ' Namelist namhsb :' 372 382 WRITE(numout,*) ' check the heat and salt budgets (T) or not (F) ln_diahsb = ', ln_diahsb 383 IF(lflush) CALL FLUSH(numout) 373 384 ENDIF 374 385 ! -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diahth.F90
r10888 r10986 127 127 END IF 128 128 129 IF(lwp) WRITE(numout,*) 130 IF(lwp) WRITE(numout,*) 'dia_hth : diagnostics of the thermocline depth' 131 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 132 IF(lwp) WRITE(numout,*) 129 IF(lwp) THEN 130 WRITE(numout,*) 131 WRITE(numout,*) 'dia_hth : diagnostics of the thermocline depth' 132 WRITE(numout,*) '~~~~~~~ ' 133 WRITE(numout,*) 134 IF(lflush) CALL FLUSH(numout) 135 ENDIF 133 136 ENDIF 134 137 -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diaptr.F90
r10968 r10986 407 407 WRITE(numout,*) ' Poleward heat & salt transport (T) or not (F) ln_diaptr = ', ln_diaptr 408 408 WRITE(numout,*) ' Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = ', ln_subbas 409 IF(lflush) CALL FLUSH(numout) 409 410 ENDIF 410 411 -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diatmb.F90
r10968 r10986 56 56 WRITE(numout,*) ' Namelist nam_diatmb : set tmb outputs ' 57 57 WRITE(numout,*) ' Switch for TMB diagnostics (T) or not (F) ln_diatmb = ', ln_diatmb 58 IF(lflush) CALL FLUSH(numout) 58 59 ENDIF 59 60 ! -
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diawri.F90
r10888 r10986 430 430 INTEGER, INTENT( in ) :: kt ! ocean time-step index 431 431 ! 432 LOGICAL :: ll_print = .FALSE. ! =T print and flush numout433 432 CHARACTER (len=40) :: clhstnam, clop, clmx ! local names 434 433 INTEGER :: inum = 11 ! temporary logical unit … … 488 487 CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 489 488 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 490 IF(lwp)WRITE(numout,*) 491 IF(lwp)WRITE(numout,*) 'Date 0 used :', nit000, ' YEAR ', nyear, & 489 IF(lwp) THEN 490 WRITE(numout,*) 491 WRITE(numout,*) 'Date 0 used :', nit000, ' YEAR ', nyear, & 492 492 & ' MONTH ', nmonth, ' DAY ', nday, 'Julian day : ', zjulian 493 IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma, &493 WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma, & 494 494 ' limit storage in depth = ', ipk 495 496 ! WRITE root name in date.file for use by postpro 497 IF(lwp) THEN 495 IF(lflush) CALL FLUSH(numout) 496 ! WRITE root name in date.file for use by postpro 498 497 CALL dia_nam( clhstnam, nwrite,' ' ) 499 498 CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) … … 505 504 506 505 CALL dia_nam( clhstnam, nwrite, 'grid_T' ) 507 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 506 IF(lwp) THEN 507 WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 508 IF(lflush) CALL FLUSH(numout) 509 ENDIF 508 510 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 509 511 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & … … 541 543 542 544 CALL dia_nam( clhstnam, nwrite, 'grid_U' ) 543 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 545 IF(lwp) THEN 546 WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 547 IF(lflush) CALL FLUSH(numout) 548 ENDIF 544 549 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu 545 550 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & … … 554 559 555 560 CALL dia_nam( clhstnam, nwrite, 'grid_V' ) ! filename 556 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 561 IF(lwp) THEN 562 WRITE(numout,*) " Name of NETCDF file ", clhstnam 563 IF(lflush) CALL FLUSH(numout) 564 ENDIF 557 565 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv 558 566 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & … … 567 575 568 576 CALL dia_nam( clhstnam, nwrite, 'grid_W' ) ! filename 569 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 577 IF(lwp) THEN 578 WRITE(numout,*) " Name of NETCDF file ", clhstnam 579 IF(lflush) CALL FLUSH(numout) 580 ENDIF 570 581 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit 571 582 & iimi, iima-iimi+1, ijmi, ijma-ijmi+1, & … … 737 748 CALL histend( nid_W, snc4chunks=snc4set ) 738 749 739 IF(lwp) WRITE(numout,*) 740 IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization' 741 IF(ll_print) CALL FLUSH(numout ) 750 IF(lwp) THEN 751 WRITE(numout,*) 752 WRITE(numout,*) 'End of NetCDF Initialization' 753 IF(lflush) CALL FLUSH(numout) 754 ENDIF 742 755 743 756 ENDIF … … 753 766 WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 754 767 WRITE(numout,*) '~~~~~~ ' 768 IF(lflush) CALL FLUSH(numout) 755 769 ENDIF 756 770 … … 887 901 !!---------------------------------------------------------------------- 888 902 ! 889 IF(lwp) WRITE(numout,*) 890 IF(lwp) WRITE(numout,*) 'dia_wri_state : single instantaneous ocean state' 891 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~ and forcing fields file created ' 892 IF(lwp) WRITE(numout,*) ' and named :', cdfile_name, '...nc' 903 IF(lwp) THEN 904 WRITE(numout,*) 905 WRITE(numout,*) 'dia_wri_state : single instantaneous ocean state' 906 WRITE(numout,*) '~~~~~~~~~~~~~ and forcing fields file created ' 907 WRITE(numout,*) ' and named :', cdfile_name, '...nc' 908 IF(lflush) CALL FLUSH(numout) 909 ENDIF 893 910 894 911 #if defined key_si3
Note: See TracChangeset
for help on using the changeset viewer.