New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10774 – NEMO

Changeset 10774


Ignore:
Timestamp:
2019-03-19T12:32:19+01:00 (5 years ago)
Author:
andmirek
Message:

GMED 450 add flush after prints

Location:
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM
Files:
132 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref

    r10746 r10774  
    10921092       sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
    10931093       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    nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
    1096    nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
    1097    nn_jctls    =    0      !  start j indice of control               over a subdomain) 
    1098    nn_jctle    =    0      !  end   j indice of control 
    1099    nn_isplt    =    1      !  number of processors in i-direction 
    1100    nn_jsplt    =    1      !  number of processors in j-direction 
    1101    nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    1102                            !     (no physical validity of the results) 
    1103    nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
     1094   nn_print    =    0       !  level of print (0 no extra print) 
     1095   ln_flush    =    .FALSE. !  flust numout after each print 
     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) 
    11041105/ 
    11051106!----------------------------------------------------------------------- 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

    r10477 r10774  
    767767            &                    TRIM( c_asmbal ) // ' at timestep = ', kt 
    768768      ENDIF 
     769 
     770      IF(lwp .AND. lflush) CALL flush(numout) 
    769771                                  
    770772   END SUBROUTINE asm_bgc_bal_wri 
     
    11401142      ENDIF 
    11411143      ! 
     1144      IF(lwp .AND. lflush) CALL flush(numout) 
    11421145   END SUBROUTINE phyto2d_asm_inc 
    11431146 
     
    13271330      ENDIF 
    13281331      ! 
     1332      IF(lwp .AND. lflush) CALL flush(numout) 
    13291333   END SUBROUTINE phyto3d_asm_inc 
    13301334 
     
    15981602      ENDIF 
    15991603      ! 
     1604      IF(lwp .AND. lflush) CALL flush(numout) 
    16001605   END SUBROUTINE pco2_asm_inc 
    16011606 
     
    18211826#endif       
    18221827      ! 
     1828      IF(lwp .AND. lflush) CALL flush(numout) 
     1829 
    18231830   END SUBROUTINE ph_asm_inc 
    18241831 
     
    21862193      ENDIF 
    21872194      ! 
     2195      IF(lwp .AND. lflush) CALL flush(numout) 
    21882196   END SUBROUTINE bgc3d_asm_inc 
    21892197 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90

    r10149 r10774  
    183183      ENDIF 
    184184      !                     
     185      IF(lwp .AND. lflush) CALL flush(numout) 
    185186   END SUBROUTINE asm_bkg_wri 
    186187 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r10759 r10774  
    207207         WRITE(numout,*) '      Choice of MLD for BGC assimilation                  mld_choice_bgc = ', mld_choice_bgc 
    208208         WRITE(numout,*) '      Maximum absolute chlorophyll increment (<=0 = off)    rn_maxchlinc = ', rn_maxchlinc 
     209         IF(lflush) CALL flush(numout) 
    209210      ENDIF 
    210211 
     
    242243         WRITE(numout,*) '       iitiaustr_date = ', iitiaustr_date 
    243244         WRITE(numout,*) '       iitiaufin_date = ', iitiaufin_date 
     245         IF(lflush) CALL flush(numout) 
    244246      ENDIF 
    245247       
     
    369371             WRITE(numout,*) '         Time-integrated weight = ', ztotwgt 
    370372             WRITE(numout,*) '         ===================================' 
     373             IF(lflush) CALL flush(numout) 
    371374          ENDIF 
    372375          
     
    590593      ENDIF 
    591594      ! 
     595      IF(lflush) CALL flush(numout) 
     596 
    592597   END SUBROUTINE asm_inc_init 
    593598 
     
    735740               WRITE(numout,*) 'tra_asm_inc : Tracer IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 
    736741               WRITE(numout,*) '~~~~~~~~~~~~' 
     742               IF(lflush) CALL flush(numout) 
    737743            ENDIF 
    738744 
     
    864870                  &  kt,' with IAU weight = ', wgtiau(it) 
    865871               WRITE(numout,*) '~~~~~~~~~~~~' 
     872               IF(lflush) CALL flush(numout) 
    866873            ENDIF 
    867874 
     
    940947                  &  kt,' with IAU weight = ', wgtiau(it) 
    941948               WRITE(numout,*) '~~~~~~~~~~~~' 
     949               IF(lflush) CALL flush(numout) 
    942950            ENDIF 
    943951 
     
    10301038                  &  kt,' with IAU weight = ', wgtiau(it) 
    10311039               WRITE(numout,*) '~~~~~~~~~~~~' 
     1040               IF(lflush) CALL flush(numout) 
    10321041            ENDIF 
    10331042 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

    r10759 r10774  
    320320          
    321321         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) 
    322323      ENDIF 
    323324      IF( .NOT. ln_bias ) RETURN 
     
    561562         ENDWHERE    
    562563      ENDIF  
    563  
     564       
     565      IF(lwp .AND. lflush) call flush(numout) 
    564566 
    565567   END SUBROUTINE bias_init 
     
    895897            ENDIF 
    896898      ENDIF 
    897  
     899      IF(lwp .AND. lflush) call flush(numout) 
    898900   END SUBROUTINE tra_bias 
    899901 
     
    971973            &          rhd_pc, gru_pc , grv_pc  ) 
    972974      ENDIF 
     975 
     976      IF(lwp .AND. lflush) call flush(numout) 
    973977 
    974978   END SUBROUTINE dyn_bias 
     
    10051009      ENDIF 
    10061010      ! 
     1011      IF(lwp .AND. lflush) call flush(numout) 
     1012 
    10071013   END SUBROUTINE bias_opn 
    10081014 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r10759 r10774  
    894894      ENDDO ! ib_bdy  
    895895 
     896      IF(lwp .AND. lflush) CALL flush(numout) 
     897 
    896898      IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init') 
    897899 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

    r10759 r10774  
    110110      IF( nn_timing == 1 ) CALL timing_start('bdy_init') 
    111111 
    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 
    115118      ! 
    116119 
     
    13481351      ENDIF 
    13491352 
     1353      IF(lwp .AND. lflush) CALL flush(numout) 
     1354 
    13501355      CALL wrk_dealloc(jpi,jpj,zfmask)  
    13511356 
     
    17221727      IF ((itest==0).AND.(lwp)) WRITE(numout,*) 'NO open boundary corner found' 
    17231728      ! 
     1729      IF(lwp. AND. lflush) CALL flush(numout) 
     1730      ! 
    17241731      ! Other tests TBD:  
    17251732      ! segments completly on land 
     
    17571764         IF(lwp) WRITE(numout,*) ' ==========  have different open bdy schemes'                                                   
    17581765         IF(lwp) WRITE(numout,*) 
     1766         IF(lflush) CALL flush(numout) 
    17591767         nstop = nstop + 1 
    17601768      ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90

    r10759 r10774  
    9999         IF(lwp) WRITE(numout,*) 'bdytide_init : initialization of tidal harmonic forcing at open boundaries' 
    100100         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     101         IF(lflush) CALL flush(numout) 
    101102      ENDIF 
    102103 
     
    120121            IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy_tide ) 
    121122            !                                               ! 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_2ddta 
    125             IF(lwp) WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
    126             IF(lwp) WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
    127             IF(lwp) THEN  
    128                     WRITE(numout,*) '             Tidal components: '  
     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: '  
    129130               DO itide = 1, nb_harmo 
    130131                  WRITE(numout,*)  '                 ', Wave(ntide(itide))%cname_tide  
    131132               END DO 
    132             ENDIF  
    133             IF(lwp) WRITE(numout,*) ' ' 
     133               WRITE(numout,*) ' ' 
     134               IF(lflush) CALL flush(numout) 
     135            ENDIF 
    134136 
    135137            ! Allocate space for tidal harmonics data - get size from OBC data arrays 
     
    333335           WRITE(numout,*) 'bdytide_update : (re)Initialization of the tidal bdy forcing at kt=',kt 
    334336           WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
     337           IF(lflush) CALL flush(numout) 
    335338        ENDIF 
    336339        ! 
     
    448451               WRITE(numout,*) 'bdy_tide_dta : Refresh nodal factors for tidal open bdy data at kt=',kt 
    449452               WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
     453               IF(lflush) CALL flush(numout) 
    450454               ENDIF 
    451455               ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/BDY/bdyvol.F90

    r9583 r10774  
    8787         IF(lwp) WRITE(numout,*)'bdy_vol : Correction of velocities along unstructured OBC' 
    8888         IF(lwp) WRITE(numout,*)'~~~~~~~' 
     89         IF(lflush) CALL flush(numout) 
    8990      END IF  
    9091 
     
    164165         IF(lwp) WRITE(numout,*)'          correction velocity zubtpecor =', zubtpecor , '(m/s)' 
    165166         IF(lwp) WRITE(numout,*)'          cumulated transport ztranst   =', ztranst   , '(m3/s)' 
     167         IF(lflush) CALL flush(numout) 
    166168      END IF  
    167169      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/C1D/dyndmp.F90

    r10759 r10774  
    100100         WRITE(numout,*) '      Damping file name               cn_resto  = ', cn_resto 
    101101         WRITE(numout,*) 
     102         IF(lflush) CALL flush(numout) 
    102103      ENDIF 
    103104 
     
    128129         CALL iom_close( imask ) 
    129130      ENDIF 
     131      ! 
     132      IF(lwp .AND. lflush) CALL flush(numout) 
    130133      ! 
    131134   END SUBROUTINE dyn_dmp_init 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r10759 r10774  
    185185     transports_3d(:,:,:,:)=0.0  
    186186     transports_2d(:,:,:)  =0.0  
    187  
     187     ! 
     188     IF(lwp .AND. lflush) CALL flush(numout) 
     189     ! 
    188190     IF( nn_timing == 1 )   CALL timing_stop('dia_dct_init') 
    189191     ! 
     
    301303        CALL wrk_dealloc( nb_sec_max,nb_type_class,nb_class_max , zsum  ) 
    302304     ENDIF     
    303  
     305     ! 
     306     IF(lwp .AND. lflush) CALL flush(numout) 
     307     ! 
    304308     IF( nn_timing == 1 )   CALL timing_stop('dia_dct') 
    305309     ! 
     
    503507  
    504508     nb_sec = jsec-1   !number of section read in the file 
    505  
     509     ! 
     510     IF(lwp .AND. lflush) CALL flush(numout) 
     511  
    506512     CALL wrk_dealloc( nb_point_max, directemp ) 
    507513     ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90

    r10759 r10774  
    160160      ALLOCATE( ana_temp(jpi,jpj,2*nb_ana,3) ) 
    161161      ana_temp(:,:,:,:) = 0._wp 
     162  
     163      IF(lwp .AND. lflush) CALL flush(numout) 
    162164 
    163165   END SUBROUTINE dia_harm_init 
     
    231233      CALL wrk_alloc( jpi , jpj , jpmax_harmo , 2 , ana_amp ) 
    232234 
    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 
    236241 
    237242      ztime_ini = nit000_han*rdt                 ! Initial time in seconds at the beginning of analysis 
     
    400405#endif 
    401406      IF(lwp) WRITE(numout,*) '  ' 
     407      IF(lwp .AND. lflush) CALL flush(numout) 
    402408 
    403409      ! A) Elevation 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90

    r9583 r10774  
    131131         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
    132132         IF(lwp) WRITE(numout,*) 
     133         IF(lwp .AND. lflush) CALL flush(numout) 
    133134      ENDIF 
    134135 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r10759 r10774  
    478478         WRITE(numout,*) '      Poleward heat & salt transport (T) or not (F)      ln_diaptr  = ', ln_diaptr 
    479479         WRITE(numout,*) '      Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins      ln_subbas  = ', ln_subbas 
     480         IF(lflush) CALL flush(numout) 
    480481      ENDIF 
    481482 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r10020 r10774  
    462462      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    463463      !! 
    464       LOGICAL ::   ll_print = .FALSE.                        ! =T print and flush numout 
    465464      CHARACTER (len=40) ::   clhstnam, clop, clmx           ! local names 
    466465      INTEGER  ::   inum = 11                                ! temporary logical unit 
     
    489488      ! ----------------- 
    490489 
    491       ! local variable for debugging 
    492       ll_print = .FALSE. 
    493       ll_print = ll_print .AND. lwp 
    494  
    495490      ! Define frequency of output and means 
    496491      zdt = rdt 
     
    532527         IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma,   & 
    533528                                 ' limit storage in depth = ', ipk 
    534  
     529        
    535530         ! WRITE root name in date.file for use by postpro 
    536531         IF(lwp) THEN 
     
    808803         IF(lwp) WRITE(numout,*) 
    809804         IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization' 
    810          IF(ll_print) CALL FLUSH(numout ) 
     805         IF(lwp .AND. lflush) CALL flush(numout) 
    811806 
    812807      ENDIF 
     
    987982      ! 
    988983      IF( nn_timing == 1 )   CALL timing_stop('dia_wri') 
     984      ! 
     985      IF(lwp .AND. lflush) CALL flush(numout) 
    989986      ! 
    990987   END SUBROUTINE dia_wri 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90

    r6491 r10774  
    7878      IF(lwp) WRITE(numout,*)'dom_clo : closed seas ' 
    7979      IF(lwp) WRITE(numout,*)'~~~~~~~' 
     80      IF(lwp .AND. lflush) CALL flush(numout) 
    8081 
    8182      ! initial values 
     
    403404      CALL lbc_lnk( emp , 'T', 1._wp ) 
    404405      ! 
     406      IF(lwp .AND. lflush) CALL flush(numout) 
     407      ! 
    405408      IF( nn_timing == 1 )  CALL timing_stop('sbc_clo') 
    406409      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90

    r10770 r10774  
    130130      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 = ',   & 
    131131           &                   nyear, '/', nmonth, '/', nday, '  nsec_day:', nsec_day, '  nsec_week:', nsec_week 
     132      IF(lwp .AND. lflush) CALL flush(numout) 
    132133 
    133134      ! Up to now, calendar parameters are related to the end of previous run (nit000-1) 
     
    253254         CALL ymds2ju( nyear, 01, 01, 0.0, fjulstartyear ) 
    254255         ! 
    255          IF(lwp .AND. nprint > 0) 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,   & 
    256258              &   '      New day, DATE Y/M/D = ', nyear, '/', nmonth, '/', nday, '      nday_year = ', nday_year 
    257          IF(lwp .AND. nprint > 0) 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,   & 
    258260              &   '   nsec_month = ', nsec_month, '   nsec_day = ', nsec_day, '   nsec_week = ', nsec_week 
     261             IF(lflush) CALL flush(numout) 
     262         ENDIF 
    259263      ENDIF 
    260264 
     
    321325               CASE ( 2 )   ;   WRITE(numout,*) ' nrstdt = 2 : calendar parameters read in restart' 
    322326               END SELECT 
    323                WRITE(numout,*) 
     327               WRITE(numout,*) ' ' 
     328               IF(lflush) CALL flush(numout) 
    324329            ENDIF 
    325330            ! Control of date 
     
    353358            WRITE(numout,*) '   number of elapsed days since the begining of run : ', adatrj 
    354359            WRITE(numout,*) 
     360            IF(lflush) CALL flush(numout) 
    355361         ENDIF 
    356362         ! 
     
    361367            WRITE(numout,*) 'rst_write : write oce restart file  kt =', kt 
    362368            WRITE(numout,*) '~~~~~~~' 
     369            IF(lflush) CALL flush(numout) 
    363370         ENDIF 
    364371         ! calendar control 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r10759 r10774  
    7979         WRITE(numout,*) 'dom_init : domain initialization' 
    8080         WRITE(numout,*) '~~~~~~~~' 
     81         IF(lflush) CALL flush(numout) 
    8182      ENDIF 
    8283      ! 
     
    193194         WRITE(numout,*) '      overwrite an existing file      ln_clobber = ', ln_clobber 
    194195         WRITE(numout,*) '      NetCDF chunksize (bytes)        nn_chunksz = ', nn_chunksz 
     196         IF(lflush) CALL flush(numout) 
    195197      ENDIF 
    196198 
     
    293295         WRITE(numout,*) '                                      ppkth2            = ', ppkth2 
    294296         WRITE(numout,*) '                                      ppacr2            = ', ppacr2 
     297         IF(lflush) CALL flush(numout) 
    295298      ENDIF 
    296299 
     
    319322         WRITE(numout,*) '   Namelist namcla' 
    320323         WRITE(numout,*) '      cross land advection                 nn_cla    = ', nn_cla 
     324         IF(lflush) CALL flush(numout) 
    321325      ENDIF 
    322326      IF ( nn_cla .EQ. 1 ) THEN 
     
    346350         WRITE(numout,*) '      number of chunks in k-dimension      nn_nchunks_k   = ', nn_nchunks_k 
    347351         WRITE(numout,*) '      apply netcdf4/hdf5 chunking & compression ln_nc4zip = ', ln_nc4zip 
     352         IF(lflush) CALL flush(numout) 
    348353      ENDIF 
    349354 
     
    406411         WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, iima2, ijma2 
    407412         WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2 
     413         IF(lflush) CALL flush(numout) 
    408414      ENDIF 
    409415      ! 
     
    467473         WRITE(numout,*) 'dom_stiff : maximum grid stiffness ratio: ', zrxmax 
    468474         WRITE(numout,*) '~~~~~~~~~' 
     475         IF(lflush) CALL flush(numout) 
    469476      ENDIF 
    470477 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90

    r10770 r10774  
    5656         IF( jperio == 5 )   WRITE(numout,*) '      jperio= 5, north fold with F-point pivot' 
    5757         IF( jperio == 6 )   WRITE(numout,*) '      jperio= 6, cyclic east-west and north fold with F-point pivot' 
     58         IF(lflush) CALL flush(numout) 
    5859      ENDIF 
    5960      ! 
     
    136137         WRITE(numout,*) '             ending index' 
    137138         WRITE(numout,25)              (mj1(jj),jj = 1,jpjdta) 
     139         IF(lflush) CALL flush(numout) 
    138140      ENDIF 
    139141 25   FORMAT( 100(10x,19i4,/) ) 
     
    195197      ENDIF 
    196198      ! 
     199      IF(lwp .AND. lflush) CALL flush(numout) 
     200      ! 
    197201   END SUBROUTINE dom_glo 
    198202 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r10754 r10774  
    600600      ENDIF 
    601601      ! 
     602      IF(lwp .AND. lflush) CALL flush(numout) 
     603      ! 
    602604      IF( nn_timing == 1 )  CALL timing_stop('dom_hgr') 
    603605      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r10759 r10774  
    562562      CALL wrk_dealloc( jpi, jpj, zwf  ) 
    563563      ! 
     564      IF(lwp .AND. lflush) CALL flush(numout) 
     565      ! 
    564566      IF( nn_timing == 1 )  CALL timing_stop('dom_msk') 
    565567      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90

    r6486 r10774  
    105105 
    106106      END SELECT 
     107       
     108      IF(lwp .AND. lflush) CALL flush(numout) 
    107109 
    108110   END SUBROUTINE dom_stp 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r10770 r10774  
    133133      IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
    134134      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     135      IF(lwp .AND. lflush) CALL flush(numout) 
    135136 
    136137      ! choose vertical coordinate (z_star, z_tilde or layer) 
     
    292293         IF(lwp) WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors' 
    293294         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 
     295         IF(lflush) CALL flush(numout) 
    294296      ENDIF 
    295297 
     
    525527         IF( lk_mpp ) CALL mpp_max( z_tmax )                                ! max over the global domain 
    526528         IF( lwp    ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 
     529         IF(lflush) CALL flush(numout) 
    527530      END IF 
    528531 
     
    593596         IF(lwp) WRITE(numout,*) 'dom_vvl_sf_swp : - time filter and swap of scale factors' 
    594597         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   - interpolate scale factors and compute depths for next time step' 
     598         IF(lflush) CALL flush(numout) 
    595599      ENDIF 
    596600 
     
    866870               neuler = 0 
    867871            ENDIF 
     872            ! 
     873            IF(lwp .AND. lflush) CALL flush(numout) 
    868874            !                             ! ----------- ! 
    869875            IF( ln_vvl_zstar ) THEN       ! z_star case ! 
     
    10081014         ! IF(       ln_vvl_kepe ) WRITE(numout,*) '              kinetic to potential energy transfer : option used' 
    10091015         ! IF( .NOT. ln_vvl_kepe ) WRITE(numout,*) '              kinetic to potential energy transfer : option not used' 
     1016         IF(lflush) CALL flush(numout) 
    10101017      ENDIF 
    10111018 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90

    r9830 r10774  
    9090      IF(lwp) WRITE(numout,*) 'dom_wri : create NetCDF mesh and mask information file(s)' 
    9191      IF(lwp) WRITE(numout,*) '~~~~~~~' 
     92      IF(lwp .AND. lflush) CALL flush(numout) 
    9293       
    9394      clnam0 = 'mesh_mask'  ! filename (mesh and mask informations) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r10759 r10774  
    170170            &                   ' w ',   MAXVAL( e3w_0(:,:,:) ) 
    171171      ENDIF 
     172      ! 
     173      IF(lwp .AND. lflush) CALL flush(numout) 
    172174      ! 
    173175      IF( nn_timing == 1 )  CALL timing_stop('dom_zgr') 
     
    265267            ENDIF 
    266268         ENDIF 
     269         IF(lflush) CALL flush(numout) 
    267270      ENDIF 
    268271 
     
    338341         WRITE(numout, "(9x,' level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d  ')" ) 
    339342         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) 
    340344      ENDIF 
    341345      DO jk = 1, jpk                      ! control positivity 
     
    598602      ENDIF 
    599603      ! 
     604      IF(lwp .AND. lflush) CALL flush(numout) 
     605      ! 
    600606      IF( nn_timing == 1 )  CALL timing_stop('zgr_bat') 
    601607      ! 
     
    651657         ! 
    652658      ENDIF 
     659      ! 
     660      IF(lwp .AND. lflush) CALL flush(numout) 
    653661      ! 
    654662   END SUBROUTINE zgr_bat_zoom 
     
    789797      CALL wrk_dealloc( jpi, jpj, zbathy ) 
    790798      ! 
     799      IF(lwp .AND. lflush) CALL flush(numout) 
     800      ! 
    791801      IF( nn_timing == 1 )  CALL timing_stop('zgr_bat_ctl') 
    792802      ! 
     
    818828      IF(lwp) WRITE(numout,*) '    zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ' 
    819829      IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~~~' 
     830      IF(lwp .AND. lflush) CALL flush(numout) 
    820831      ! 
    821832      mbkt(:,:) = MAX( mbathy(:,:) , 1 )    ! bottom k-index of T-level (=1 over land) 
     
    862873      IF(lwp) WRITE(numout,*) '    zgr_top_level : ocean top k-index of T-, U-, V- and W-levels ' 
    863874      IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~~~' 
     875      IF(lwp .AND. lflush) CALL flush(numout) 
    864876      ! 
    865877      mikt(:,:) = MAX( misfdep(:,:) , 1 )    ! top k-index of T-level (=1) 
     
    12461258         WRITE(numout,*) 
    12471259         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) 
    12481261      ENDIF   
    12491262      ! 
     
    12731286      INTEGER  ::   id, jd, nprocd 
    12741287      INTEGER  ::   icompt, ibtest, ibtestim1, ibtestip1, ibtestjm1, ibtestjp1   ! (ISF) 
    1275       LOGICAL  ::   ll_print         ! Allow  control print for debugging 
    12761288      REAL(wp) ::   ze3tp , ze3wp    ! Last ocean level thickness at T- and W-points 
    12771289      REAL(wp) ::   zdepwp, zdepth   ! Ajusted ocean depth to avoid too small e3t 
     
    17531765      ENDIF  
    17541766 
     1767      IF(lwp .AND. lflush) CALL flush(numout) 
     1768 
    17551769      CALL wrk_dealloc( jpi, jpj, zmask, zbathy, zrisfdep ) 
    17561770      CALL wrk_dealloc( jpi, jpj, zmisfdep, zmbathy ) 
     
    18551869         WRITE(numout,*) '        Offset for Zb                                rn_zb_b       = ',rn_zb_b 
    18561870         WRITE(numout,*) '        Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b' 
     1871         IF(lflush) CALL flush(numout) 
    18571872      ENDIF 
    18581873 
     
    22362251      END DO 
    22372252      ! 
     2253      IF(lwp .AND. lflush) CALL flush(numout) 
     2254      ! 
    22382255      CALL wrk_dealloc( jpi, jpj, zenv, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 ) 
    22392256      ! 
     
    25472564      END DO 
    25482565 
     2566      IF(lwp .AND. lflush) CALL flush(numout) 
     2567 
    25492568      CALL wrk_dealloc( jpk, z_gsigw, z_gsigt, z_gsi3w                                      ) 
    25502569      CALL wrk_dealloc( jpk, z_esigt, z_esigw                                               ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90

    r10770 r10774  
    9393            WRITE(numout,*) '   T & S data not used' 
    9494         ENDIF 
     95         IF(lflush) CALL flush(numout) 
    9596      ENDIF 
    9697      ! 
     
    216217            WRITE(numout,*) 
    217218            WRITE(numout,*) 'dta_tsd: interpolates T & S data onto the s- or mixed s-z-coordinate mesh' 
     219            IF(lflush) CALL flush(numout) 
    218220         ENDIF 
    219221         ! 
     
    294296         CALL prihre( ptsd(:,:,jpkm1,jp_sal), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout ) 
    295297         WRITE(numout,*) 
     298         IF(lflush) CALL flush(numout) 
    296299      ENDIF 
    297300      ! 
     
    299302         !                                              (data used only for initialisation) 
    300303         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) 
    301305                                        DEALLOCATE( sf_tsd(jp_tem)%fnow )     ! T arrays in the structure 
    302306         IF( sf_tsd(jp_tem)%ln_tint )   DEALLOCATE( sf_tsd(jp_tem)%fdta ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r6486 r10774  
    7979      IF(lwp) WRITE(numout,*) 'istate_ini : Initialization of the dynamics and tracers' 
    8080      IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 
     81      IF(lwp .AND. lflush) CALL flush(numout) 
    8182 
    8283      CALL dta_tsd_init                       ! Initialisation of T & S input data 
     
    209210      IF(lwp) WRITE(numout,*) 'istate_t_s : Philander s initial temperature profile' 
    210211      IF(lwp) WRITE(numout,*) '~~~~~~~~~~   and constant salinity (',zsal,' psu)' 
     212      IF(lwp .AND. lflush) CALL flush(numout) 
    211213      ! 
    212214      DO jk = 1, jpk 
     
    359361            ! 
    360362      END SELECT 
     363      ! 
     364      IF(lwp .AND. lflush) CALL flush(numout) 
    361365      ! 
    362366   END SUBROUTINE istate_eel 
     
    443447      ENDIF 
    444448      ! 
     449      IF(lwp .AND. lflush) CALL flush(numout) 
     450      ! 
    445451   END SUBROUTINE istate_gyre 
    446452 
     
    566572      CALL wrk_dealloc( jpi, jpj, jpk, zprn) 
    567573      ! 
     574      IF(lwp .AND. lflush) CALL flush(numout) 
     575      ! 
    568576   END SUBROUTINE istate_uvg 
    569577 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

    r6488 r10774  
    207207         WRITE(numout,*) 
    208208         WRITE(numout,*) '          smallest real computer value       rsmall = ', rsmall 
     209         IF(lflush) CALL flush(numout) 
    209210      ENDIF 
    210211 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

    r6487 r10774  
    104104         IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and relative vorticity' 
    105105         IF(lwp) WRITE(numout,*) '~~~~~~~   NOT optimal for auto-tasking case' 
     106         IF(lwp .AND. lflush) CALL flush(numout) 
    106107      ENDIF 
    107108 
     
    285286         IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and' 
    286287         IF(lwp) WRITE(numout,*) '~~~~~~~   relative vorticity' 
     288         IF(lwp .AND. lflush) CALL flush(numout) 
    287289      ENDIF 
    288290 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90

    r10759 r10774  
    152152         IF( nadv ==  3 )   WRITE(numout,*) '         flux form   : UBS       scheme is used' 
    153153         IF( nadv == -1 )   WRITE(numout,*) '         esopa test: use all advection formulation' 
     154         IF(lflush) CALL flush(numout) 
    154155      ENDIF 
    155156      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90

    r6486 r10774  
    6666         WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection' 
    6767         WRITE(numout,*) '~~~~~~~~~~~~' 
     68         IF(lflush) CALL flush(numout) 
    6869      ENDIF 
    6970      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90

    r6486 r10774  
    9090         IF(lwp) WRITE(numout,*) 'dyn_adv_ubs : UBS flux form momentum advection' 
    9191         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     92         IF(lwp .AND. lflush) CALL flush(numout) 
    9293      ENDIF 
    9394      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r10759 r10774  
    136136         IF(lwp) THEN  
    137137         WRITE(numout,*) " ! restore original density" 
     138         IF(lflush) CALL flush(numout) 
    138139         ENDIF 
    139140         rhd(:,:,:) = z_rhd_st(:,:,:)     ! restore original density 
     
    190191         WRITE(numout,*) '      s-coord. (Pressure Jacobian: Cubic polynomial)    ln_hpg_prj    = ', ln_hpg_prj 
    191192         WRITE(numout,*) '      time stepping: centered (F) or semi-implicit (T)  ln_dynhpg_imp = ', ln_dynhpg_imp 
     193         IF(lflush) CALL flush(numout) 
    192194      ENDIF 
    193195      ! 
     
    260262         IF(lwp) WRITE(numout,*) 'dyn:hpg_zco : hydrostatic pressure gradient trend' 
    261263         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   z-coordinate case ' 
     264         IF(lwp .AND. lflush) CALL flush(numout) 
    262265      ENDIF 
    263266 
     
    325328         IF(lwp) WRITE(numout,*) 'dyn:hpg_zps : hydrostatic pressure gradient trend' 
    326329         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   z-coordinate with partial steps - vector optimization' 
     330         IF(lwp .AND. lflush) CALL flush(numout) 
    327331      ENDIF 
    328332 
     
    423427         IF(lwp) WRITE(numout,*) 'dyn:hpg_sco : hydrostatic pressure gradient trend' 
    424428         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, OPA original scheme used' 
     429         IF(lwp .AND. lflush) CALL flush(numout) 
    425430      ENDIF 
    426431 
     
    514519         IF(lwp) WRITE(numout,*) 'dyn:hpg_isf : hydrostatic pressure gradient trend for ice shelf' 
    515520         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, OPA original scheme used' 
     521         IF(lwp .AND. lflush) CALL flush(numout) 
    516522      ENDIF 
    517523 
     
    768774         IF(lwp) WRITE(numout,*) 'dyn:hpg_djc : hydrostatic pressure gradient trend' 
    769775         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, density Jacobian with cubic polynomial scheme' 
     776         IF(lwp .AND. lflush) CALL flush(numout) 
    770777      ENDIF 
    771778 
     
    10011008         IF(lwp) WRITE(numout,*) 'dyn:hpg_prj : hydrostatic pressure gradient trend' 
    10021009         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, cubic spline pressure Jacobian' 
     1010         IF(lwp .AND. lflush) CALL flush(numout) 
    10031011      ENDIF 
    10041012 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

    r6486 r10774  
    8888         IF(lwp) WRITE(numout,*) 'dyn_keg : kinetic energy gradient trend, scheme number=', kscheme 
    8989         IF(lwp) WRITE(numout,*) '~~~~~~~' 
     90         IF(lflush) CALL flush(numout) 
    9091      ENDIF 
    9192 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90

    r6486 r10774  
    100100            IF(lwp) WRITE(numout,*) 'dyn_ldf : no lateral diffusion on momentum setup' 
    101101            IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     102            IF(lflush) CALL flush(numout) 
    102103         ENDIF 
    103104      END SELECT 
     
    139140         WRITE(numout,*) '          horizontal (geopotential)   ln_dynldf_hor   = ', ln_dynldf_hor 
    140141         WRITE(numout,*) '          iso-neutral                 ln_dynldf_iso   = ', ln_dynldf_iso 
     142         IF(lflush) CALL flush(numout) 
    141143      ENDIF 
    142144 
     
    229231         IF( nldf ==  4 )   WRITE(numout,*) '              laplacian and bilaplacian operators' 
    230232         IF( nldf ==  5 )   WRITE(numout,*) '              rotated laplacian and bilaplacian operators' 
     233         IF(lflush) CALL flush(numout) 
    231234      ENDIF 
    232235      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r6486 r10774  
    9090         WRITE(numout,*) 'dyn_ldf_bilap : iso-level bilaplacian operator' 
    9191         WRITE(numout,*) '~~~~~~~~~~~~~' 
     92         IF(lwp .AND. lflush) CALL flush(numout) 
    9293      ENDIF 
    9394 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

    r8280 r10774  
    9898         IF(lwp) WRITE(numout,*) 'dyn_ldf_bilapg : horizontal biharmonic operator in s-coordinate' 
    9999         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 
     100         IF(lwp .AND. lflush) CALL flush(numout) 
    100101         !                                      ! allocate dyn_ldf_bilapg arrays 
    101102         IF( dyn_ldf_bilapg_alloc() /= 0 )   CALL ctl_stop('STOP', 'dyn_ldf_bilapg: failed to allocate arrays') 
     
    123124            IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)),  & 
    124125               &                             ' wi', sqrt(MAXVAL(wslpi))     , ' wj', sqrt(MAXVAL(wslpj)) 
     126            IF(lwp .AND. lflush) CALL flush(numout) 
    125127         endif 
    126128!!end 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_iso.F90

    r6486 r10774  
    127127         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate horizontal diffusive operator' 
    128128         !                                      ! allocate dyn_ldf_bilap arrays 
     129         IF(lflush) CALL flush(numout) 
    129130         IF( dyn_ldf_iso_alloc() /= 0 )   CALL ctl_stop('STOP', 'dyn_ldf_iso: failed to allocate arrays') 
    130131      ENDIF 
     
    151152            IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)),  & 
    152153               &                             ' wi', sqrt(MAXVAL(wslpi))     , ' wj', sqrt(MAXVAL(wslpj)) 
     154            IF(lwp .AND. lflush) CALL flush(numout) 
    153155         endif 
    154156!!end 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90

    r6486 r10774  
    7070         IF(lwp) WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator' 
    7171         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     72         IF(lflush) CALL flush(numout) 
    7273      ENDIF 
    7374      !                                                ! =============== 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90

    r10759 r10774  
    145145            WRITE(numout,*) 
    146146         ENDIF 
     147         IF(lflush) CALL flush(numout) 
    147148      ENDIF 
    148149      ! 
     
    304305      CALL lbc_lnk( zunep, 'U', -1.0_wp ) 
    305306      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 
    309312 
    310313      !!  Compute, once and for all, the horizontal divergence (zhdivnep) 
     
    328331         END DO 
    329332      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 
    333340 
    334341!!    Deallocate temporary workspace arrays, which are all local to 
     
    384391      IF(lwp) WRITE(numout,*) 'horizontal velocity divergence and' 
    385392      IF(lwp) WRITE(numout,*) 'relative vorticity of Neptune flow' 
     393      IF(lwp .AND. lflush) CALL flush(numout) 
    386394#if defined key_noslip_accurate 
    387395   !!---------------------------------------------------------------------- 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r6487 r10774  
    119119         IF(lwp) WRITE(numout,*) 'dyn_nxt : time stepping' 
    120120         IF(lwp) WRITE(numout,*) '~~~~~~~' 
     121         IF(lflush) CALL flush(numout) 
    121122      ENDIF 
    122123 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90

    r6486 r10774  
    231231         WRITE(numout,*) '     Free surface with time splitting       lk_dynspg_ts  = ', lk_dynspg_ts 
    232232         WRITE(numout,*) '     Filtered free surface cst volume       lk_dynspg_flt = ', lk_dynspg_flt 
     233         IF(lflush) CALL flush(numout) 
    233234      ENDIF 
    234235 
     
    266267         IF( nspg ==  1 )   WRITE(numout,*) '     free surface with time splitting scheme' 
    267268         IF( nspg ==  2 )   WRITE(numout,*) '     filtered free surface' 
     269         IF(lflush) CALL flush(numout) 
    268270      ENDIF 
    269271 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_exp.F90

    r9583 r10774  
    7575         ! 
    7676         IF( lk_vvl .AND. lwp ) WRITE(numout,*) '              lk_vvl=T : spg is included in dynhpg' 
     77         IF(lflush) CALL flush(numout) 
    7778      ENDIF 
    7879 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

    r9321 r10774  
    120120         IF(lwp) WRITE(numout,*) 'dyn_spg_flt : surface pressure gradient trend' 
    121121         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   (free surface constant volume case)' 
     122         IF(lwp .AND. lflush) CALL flush(numout) 
    122123        
    123124         ! set to zero free surface specific arrays 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r10759 r10774  
    195195         IF(lwp) WRITE(numout,*) '~~~~~~~~~~   free surface with time splitting' 
    196196         IF(lwp) WRITE(numout,*) 
     197         IF(lflush) CALL flush(numout) 
    197198         ! 
    198199         IF (neuler==0) ll_init=.TRUE. 
     
    11611162      IF(lwp) WRITE(numout,*) '     Barotropic time step [s] is :', rdtbt 
    11621163      IF(lwp) WRITE(numout,*) '     Maximum Courant number is   :', zcmax 
     1164      IF(lwp .AND. lflush) CALL flush(numout) 
    11631165      ! 
    11641166      IF ((.NOT.ln_bt_av).AND.(.NOT.ln_bt_fw)) THEN 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynvor.F90

    r10759 r10774  
    234234         IF(lwp) WRITE(numout,*) 'dyn:vor_ene : vorticity term: energy conserving scheme' 
    235235         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     236         IF(lflush) CALL flush(numout) 
    236237      ENDIF 
    237238 
     
    346347         IF(lwp) WRITE(numout,*) 'dyn:vor_mix : vorticity term: mixed energy/enstrophy conserving scheme' 
    347348         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     349         IF(lflush) CALL flush(numout) 
    348350      ENDIF 
    349351 
     
    463465         IF(lwp) WRITE(numout,*) 'dyn:vor_ens : vorticity term: enstrophy conserving scheme' 
    464466         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     467         IF(lflush) CALL flush(numout) 
    465468      ENDIF 
    466469 
     
    587590         IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
    588591         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     592         IF(lflush) CALL flush(numout) 
    589593#if ! defined key_vvl 
    590594         IF( .NOT.ALLOCATED(ze3f) ) THEN 
     
    743747         WRITE(numout,*) '           enstrophy and energy conserving scheme     ln_dynvor_een = ', ln_dynvor_een 
    744748         WRITE(numout,*) '           enstrophy and energy conserving scheme (old) ln_dynvor_een_old= ', ln_dynvor_een_old 
     749         IF(lflush) CALL flush(numout) 
    745750      ENDIF 
    746751 
     
    798803         IF( nvor ==  3 )   WRITE(numout,*) '         vorticity scheme : energy and enstrophy conserving scheme' 
    799804         IF( nvor == -1 )   WRITE(numout,*) '         esopa test: use all lateral physics options' 
     805         IF(lflush) CALL flush(numout) 
    800806      ENDIF 
    801807      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90

    r6486 r10774  
    7474         IF(lwp)WRITE(numout,*) 
    7575         IF(lwp)WRITE(numout,*) 'dyn_zad : arakawa advection scheme' 
     76         IF(lwp .AND. lflush) CALL flush(numout) 
    7677      ENDIF 
    7778 
     
    188189         IF(lwp)WRITE(numout,*) 
    189190         IF(lwp)WRITE(numout,*) 'dyn_zad_zts : arakawa advection scheme with sub-timesteps' 
     191         IF(lwp .AND. lflush) CALL flush(numout) 
    190192      ENDIF 
    191193 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90

    r6486 r10774  
    136136         IF( nzdf ==  0 )   WRITE(numout,*) '              Explicit time-splitting scheme' 
    137137         IF( nzdf ==  1 )   WRITE(numout,*) '              Implicit (euler backward) scheme' 
     138         IF(lflush) CALL flush(numout) 
    138139      ENDIF 
    139140      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_exp.F90

    r6486 r10774  
    7373         WRITE(numout,*) 'dyn_zdf_exp : vertical momentum diffusion - explicit operator' 
    7474         WRITE(numout,*) '~~~~~~~~~~~ ' 
     75         IF(lflush) CALL flush(numout) 
    7576      ENDIF 
    7677 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

    r6795 r10774  
    8181         IF(lwp) WRITE(numout,*) 'dyn_zdf_imp : vertical momentum diffusion implicit operator' 
    8282         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 
     83         IF(lwp .AND. lflush) CALL flush(numout) 
    8384         ! 
    8485         IF( lk_vvl ) THEN   ;    r_vvl = 1._wp       ! Variable volume indicator 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r8400 r10774  
    8787         IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height' 
    8888         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     89         IF(lwp .AND. lflush) CALL flush(numout) 
    8990         ! 
    9091      ENDIF 
     
    184185         IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity ' 
    185186         IF(lwp) WRITE(numout,*) '~~~~~ ' 
     187         IF(lwp .AND. lflush) CALL flush(numout) 
    186188         ! 
    187189         wn(:,:,jpk) = 0._wp                  ! bottom boundary condition: w=0 (set once for all) 
     
    265267         IF(lwp) WRITE(numout,*) 'ssh_swp : Asselin time filter and swap of sea surface height' 
    266268         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     269         IF(lwp .AND. lflush) CALL flush(numout) 
    267270      ENDIF 
    268271 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90

    r6486 r10774  
    6565      IF(lwp) WRITE(numout,*) '~~~~~~~~' 
    6666      IF(lwp) WRITE(numout,*) '           jpnfl = ',jpnfl 
     67      IF(lwp .AND. lflush) CALL flush(numout) 
    6768       
    6869      !-------------------------! 
     
    7273 
    7374         IF(lwp) WRITE(numout,*) '        float restart file read' 
     75         IF(lwp .AND. lflush) CALL flush(numout) 
    7476          
    7577         ! open the restart file  
     
    9395 
    9496            IF(lwp) WRITE(numout,*) '        add new floats' 
     97            IF(lwp .AND. lflush) CALL flush(numout) 
    9598 
    9699            IF( ln_ariane )THEN  !Add new floats with ariane convention 
     
    148151      DO jfl = kfl_start,kfl_end 
    149152         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)write(numout,*)'read:',jfl,flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash 
     154         IF(lwp .AND. lflush) CALL flush(numout) 
    151155      END DO 
    152156      CLOSE(inum) 
     
    275279               IF(lwp) WRITE(numout,*)'geographical position',flxx(jfl),flyy(jfl),flzz(jfl) 
    276280               IF(lwp) WRITE(numout,*)'index position',tpifl(jfl),tpjfl(jfl),tpkfl(jfl) 
     281               IF(lwp .AND. lflush) CALL flush(numout) 
    277282            ENDIF 
    278283         END DO 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90

    r10745 r10774  
    129129   INTEGER ::   nn_bench     !: benchmark parameter (0/1) 
    130130   INTEGER ::   nn_bit_cmp   =    0    !: bit reproducibility  (0/1) 
     131   LOGICAL ::   ln_flush     !: namelist flush numout after write statement 
    131132   !                                           
     133   LOGICAL ::   lflush       !: flush numout after write statement 
    132134   INTEGER ::   nprint, nictls, nictle, njctls, njctle, isplt, jsplt, nbench    !: OLD namelist names 
    133135 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/IOM/iom_nf90.F90

    r10770 r10774  
    8484         !                 ! ============= 
    8585         IF( ldwrt ) THEN  ! ... in write mode 
    86             IF(lwp .AND. nprint > 0) 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 
    8790            IF( snc4set%luse ) THEN 
    8891               CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_WRITE  , if90id ), clinfo) 
     
    9295            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy                          ), clinfo) 
    9396         ELSE              ! ... in read mode 
    94             IF(lwp .AND. nprint > 0) 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 
    95101            CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_NOWRITE, if90id, chunksize = ichunk ), clinfo) 
    96102         ENDIF 
     
    103109               cdname = TRIM(cltmp) 
    104110            ENDIF 
    105             IF(lwp .AND. nprint > 0) 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 
    106115 
    107116            IF( llclobber ) THEN   ;   imode = IOR( NF90_64BIT_OFFSET, NF90_CLOBBER   ) 
     
    109118            ENDIF 
    110119            IF( snc4set%luse ) THEN 
    111                IF(lwp .AND. nprint > 0) 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 
    112124               CALL GET_NF90_SYMBOL("NF90_HDF5", ihdf5) 
    113125               IF( llclobber ) THEN   ;   imode = IOR(ihdf5, NF90_CLOBBER) 
     
    158170        &                                               len  = iom_file(kiomid)%lenuld ), clinfo ) 
    159171         ENDIF 
    160          IF(lwp .AND. nprint > 0) 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 
    161176      ELSE 
    162177         kiomid = 0               ! return error flag 
     
    431446         ! trick: defined to 0 to say that dimension variables are defined but not yet written 
    432447         iom_file(kiomid)%dimsz(1, 1)  = 0    
    433          IF(lwp .AND. nprint > 1) 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 
    434452      ENDIF 
    435453      ! define the data if it is not already done 
     
    500518            IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' chunked ok. Chunks sizes: ', ichunksz 
    501519         ENDIF 
    502          IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' defined ok' 
     520         IF(lwp .AND. nprint > 1) THEN 
     521            WRITE(numout,*) TRIM(clinfo)//' defined ok' 
     522            IF(lflush) CALL flush(numout) 
     523         ENDIF 
    503524      ELSE 
    504525         idvar = kvid 
     
    543564               CALL iom_nf90_check(NF90_INQUIRE_DIMENSION( if90id, 3, len = iom_file(kiomid)%dimsz(1,3) ), clinfo) 
    544565               iom_file(kiomid)%dimsz(1  , 4) = 1   ! unlimited dimension 
    545                IF(lwp .AND. nprint > 1) 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 
    546570            ENDIF 
    547571         ENDIF 
     
    561585         IF( iom_file(kiomid)%luld(idvar) )   iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar)    & 
    562586               &                            = iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar) + 1 
    563          IF(lwp .AND. nprint > 0) 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 
    564591      ENDIF 
    565592      !      
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/IOM/prtctl.F90

    r6486 r10774  
    369369         WRITE(j_id,*) 
    370370         WRITE(j_id,*) 
     371         IF(lflush) CALL flush(j_id) 
    371372 
    3723739000     FORMAT(a41,i4.4,a14) 
     
    515516      IF(lwp) WRITE(numout,*)' sum ilcitl(1,j) = ', zjdom, ' jpjglo = ', jpjglo 
    516517      IF(lwp) WRITE(numout,*) 
     518      IF(lwp .AND. lflush) CALL flush(numout) 
    517519       
    518520 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

    r10770 r10774  
    116116                  ENDIF 
    117117               ENDIF 
     118               IF(lflush) CALL flush(numout) 
    118119            ENDIF 
    119120            ! 
     
    211212            IF ( snc4set%luse )      WRITE(numout,*) 'rst_read : configured with NetCDF4 support' 
    212213            WRITE(numout,*) '~~~~~~~~' 
     214            IF(lflush) CALL flush(numout) 
    213215         ENDIF 
    214216 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LBC/cla.F90

    r6486 r10774  
    8383         IF(lwp) WRITE(numout,*) 'div_cla : cross land advection on hdiv ' 
    8484         IF(lwp) WRITE(numout,*) '~~~~~~~~' 
     85         IF(lflush) CALL flush(numout)  
    8586         ! 
    8687         IF( nbab == 1 )   CALL cla_bab_el_mandeb('ini')    ! Bab el Mandeb ( Red Sea - Indian ocean ) 
     
    120121         IF(lwp) WRITE(numout,*) 'tra_cla : cross land advection on tracers ' 
    121122         IF(lwp) WRITE(numout,*) '~~~~~~~~' 
     123         IF(lflush) CALL flush(numout) 
    122124      ENDIF 
    123125      ! 
     
    147149         IF(lwp) WRITE(numout,*) 'cla_dynspg : cross land advection on (ua,va) ' 
    148150         IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 
     151         IF(lflush) CALL flush(numout) 
    149152      ENDIF 
    150153      ! 
     
    181184      IF(lwp) WRITE(numout,*) 'cla_init : cross land advection initialisation ' 
    182185      IF(lwp) WRITE(numout,*) '~~~~~~~~~' 
     186      IF(lwp .AND. lflush) CALL flush(numout) 
    183187      ! 
    184188      !                           ! Allocate arrays for this module 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

    r10770 r10774  
    7373         WRITE(numout,*) '         nimpp  = ', nimpp 
    7474         WRITE(numout,*) '         njmpp  = ', njmpp 
     75         IF(lflush) CALL flush(numout) 
    7576      ENDIF 
    7677 
     
    137138      IF(lwp) WRITE(numout,*) 'mpp_init : Message Passing MPI' 
    138139      IF(lwp) WRITE(numout,*) '~~~~~~~~' 
     140      IF(lflush) CALL flush(numout) 
    139141 
    140142 
     
    287289      IF(lwp) WRITE(numout,*)' sum ilcit(1,j)=', zjdom, ' jpjglo=', jpjglo 
    288290      IF(lwp) WRITE(numout,*) 
     291      IF(lwp .AND. lflush) CALL flush(numout) 
    289292 
    290293      IF(lwp .AND. nprint>1) THEN 
     
    305308            il1 = il1+ifreq 
    306309         END DO 
     310         IF(lflush) CALL flush(numout) 
     311      ENDIF 
    307312 9200     FORMAT('     ***',20('*************',a3)) 
    308313 9203     FORMAT('     *     ',20('         *   ',a3)) 
     
    310315 9202     FORMAT(' ',i3,' *  ',20(i3,'  x',i3,'   *   ')) 
    311316 9204     FORMAT('     *  ',20('      ',i3,'   *   ')) 
    312       ENDIF 
    313317 
    314318      ! 5. From global to local 
     
    447451         WRITE(numout,*) ' jprecj = ', jprecj , ' npnw   = ', npnw 
    448452         WRITE(numout,*) 
     453         IF(lflush) CALL flush(numout) 
    449454      ENDIF 
    450455 
     
    456461         CALL mpp_ini_north 
    457462         IF(lwp) WRITE(numout,*) ' mpp_init : North fold boundary prepared for jpni >1' 
     463         IF(lwp .AND. lflush) CALL flush(numout) 
    458464      ENDIF 
    459465 
     
    511517          WRITE(numout,*) '                    ihals = ', ihals(1), ihals(2) 
    512518          WRITE(numout,*) '                    ihale = ', ihale(1), ihale(2) 
     519          IF(lflush) CALL flush(numout) 
    513520      ENDIF 
    514521      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90

    r10759 r10774  
    9393      IF(lwp)WRITE(numout,*) '~~~~~~~~' 
    9494      IF(lwp)WRITE(numout,*) ' ' 
     95      IF(lwp .AND. lflush) CALL flush(numout) 
    9596 
    9697      IF( jpni*jpnj < jpnij ) CALL ctl_stop( ' jpnij > jpni x jpnj impossible' ) 
     
    184185      IF(lwp) WRITE(numout,*) ' sum ilcj(1,j)=',zjdom,' jpjglo=',jpjglo 
    185186      IF(lwp) WRITE(numout,*) 
     187      IF(lwp .AND. lflush) CALL flush(numout) 
    186188 
    187189 
     
    364366            il1 = il1+ifreq 
    365367         END DO 
     368         IF(lflush) CALL flush(numout) 
    366369 9400     FORMAT('     ***',20('*************',a3)) 
    367370 9403     FORMAT('     *     ',20('         *   ',a3)) 
     
    572575         WRITE(numout,*) ' jprecj = ', jprecj , ' npnw   = ', npnw 
    573576         WRITE(numout,*) 
     577         IF(lflush) CALL flush(numout) 
    574578      ENDIF 
    575579 
     
    581585         CALL mpp_ini_north 
    582586         IF(lwp) WRITE(numout,*) ' mpp_init2 : North fold boundary prepared for jpni >1' 
     587         IF(lwp .AND. lflush) CALL flush(numout) 
    583588      ENDIF 
    584589 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r10759 r10774  
    159159     nkahm_smag = 1 
    160160#endif 
    161  
     161      ! 
     162      IF(lwp .AND. lflush) CALL flush(numout) 
    162163      ! 
    163164   END SUBROUTINE ldf_dyn_init 
     
    210211      ENDIF 
    211212      ! 
     213      IF(lwp .AND. lflush) CALL flush(numout) 
     214      ! 
    212215   END SUBROUTINE ldf_zpf_1d 
    213216 
     
    252255      ENDIF 
    253256      ! 
     257      IF(lwp .AND. lflush) CALL flush(numout) 
     258      ! 
    254259   END SUBROUTINE ldf_zpf_1d_3d 
    255260 
     
    293298      ENDIF 
    294299      ! 
     300      IF(lwp .AND. lflush) CALL flush(numout) 
     301      ! 
    295302   END SUBROUTINE ldf_zpf_3d 
    296303 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

    r6486 r10774  
    4141      IF(lwp) WRITE(numout,*) 'ldf_dyn_c2d : 2d lateral eddy viscosity coefficient' 
    4242      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     43      IF(lwp .AND. lflush) CALL flush(numout) 
    4344 
    4445      ! harmonic operator (ahm1, ahm2) : ( T- and F- points) (used for laplacian operators 
     
    8586            CALL prihre(ahm2,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 
    8687         ENDIF 
     88         IF(lwp .AND. lflush) CALL flush(numout) 
    8789      ENDIF 
    8890 
     
    119121            CALL prihre(ahm4,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 
    120122         ENDIF 
     123         IF(lwp .AND. lflush) CALL flush(numout) 
    121124      ENDIF 
    122125      ! 
     
    160163      IF(lwp) WRITE(numout,*) '~~~~~~  --' 
    161164      IF(lwp) WRITE(numout,*) '        orca ocean configuration' 
     165      IF(lwp .AND. lflush) CALL flush(numout) 
    162166 
    163167      IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 
     
    316320      ENDIF 
    317321      ! 
     322      IF(lwp .AND. lflush) CALL flush(numout) 
     323      ! 
    318324      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    319325      ! 
     
    358364      IF(lwp) WRITE(numout,*) '~~~~~~  --' 
    359365      IF(lwp) WRITE(numout,*) '        orca_r1 configuration' 
     366      IF(lwp .AND. lflush) CALL flush(numout) 
    360367 
    361368      IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 
     
    519526      ENDIF 
    520527      ! 
     528      IF(lwp .AND. lflush) CALL flush(numout) 
     529      ! 
    521530      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    522531      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

    r9975 r10774  
    179179         ENDIF 
    180180      ENDIF 
     181      ! 
     182      IF(lwp .AND. lflush) CALL flush(numout) 
    181183      ! 
    182184      CALL wrk_dealloc( jpk, zcoef ) 
     
    432434      ENDIF 
    433435      ! 
     436      IF(lwp .AND. lflush) CALL flush(numout) 
     437      ! 
    434438      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    435439      CALL wrk_dealloc( jpk   ,         zcoef ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_smag.F90

    r9583 r10774  
    106106 
    107107 
    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) WRITE(numout,*) 
     109         IF(lwp) WRITE(numout,*) 'ldf_dyn_smag : 3D lateral eddy viscosity coefficient' 
     110         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     111         IF(lwp .AND. lflush) CALL flush(numout) 
     112 
    112113      ENDIF 
    113114       
     
    174175            ahm1(:,:,jpk) = ahm1(:,:,jpkm1) 
    175176            ahm2(:,:,jpk) = ahm2(:,:,jpkm1) 
    176             IF(lwp.and.kt==nit000) WRITE(numout,'(36x," ahm ", 7x)') 
     177            IF(lwp.and.kt==nit000) THEN 
     178               WRITE(numout,'(36x," ahm ", 7x)')  
     179               IF(lflush) CALL flush(numout) 
     180            ENDIF 
     181         IF(lwp.and.kt==nit000) THEN 
    177182            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 
     183               WRITE(numout,'(30x,E10.2,8x,i3)') ahm1(jpi/2,jpj/2,jk), jk 
    180184            END DO 
     185            IF(lflush) CALL flush(numout) 
     186         ENDIF 
    181187      CALL lbc_lnk( ahm1, 'T', 1. )   ! Lateral boundary conditions on ( ahtt ) 
    182188      CALL lbc_lnk( ahm2, 'F', 1. )   ! Lateral boundary conditions on ( ahtt ) 
     
    245251            ahm4(:,:,jpk) = ahm4(:,:,jpkm1) 
    246252 
    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 
     253      IF(  kt == nit000 .AND. lwp) THEN 
     254         DO jk = 1, jpk 
     255               WRITE(numout,'(30x,E10.2,8x,i3)') ahm3(jpi/2,jpj/2,jk), jk 
     256               IF(lflush) CALL flush(numout) 
     257         END DO 
    251258      ENDIF    
    252       END DO 
     259 
    253260      CALL lbc_lnk( ahm3, 'U', 1. )   ! Lateral boundary conditions 
    254261      CALL lbc_lnk( ahm4, 'V', 1. ) 
     
    269276      IF( lk_mpp )   CALL mpp_max( zdeltat )                 ! max over the global domain 
    270277      ! 
    271       IF( MOD( kt, nwrite ) == 1 .AND. lwp )   WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap:  abs(ahm) max: ', zdeltat 
     278      IF( MOD( kt, nwrite ) == 1 .AND. lwp )  THEN 
     279         WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap:  abs(ahm) max: ', zdeltat 
     280         IF(lflush) CALL flush(numout) 
     281      ENDIF 
    272282    ENDIF 
    273283    If(ln_dynldf_bilap)THEN 
     
    282292      IF( lk_mpp )   CALL mpp_max( zdeltat )                 ! max over the global domain 
    283293      ! 
    284       IF( MOD( kt, nwrite ) == 1 .AND. lwp )   WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 
     294      IF( MOD( kt, nwrite ) == 1 .AND. lwp )  THEN 
     295         WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 
     296         IF(lflush) CALL flush(numout) 
     297      ENDIF 
    285298      ! 
    286299   ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

    r6498 r10774  
    414414         IF(lwp) THEN  
    415415            WRITE(numout,*) '          Horizontal mixing in s-coordinate: slope = slope of s-surfaces'  
     416            IF(lflush) CALL flush(numout) 
    416417         ENDIF  
    417418 
     
    451452            IF(lwp) WRITE(numout,*) ' max slop: u',SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)),  &  
    452453               &                             ' wi', sqrt(MAXVAL(wslpi)), ' wj', sqrt(MAXVAL(wslpj))  
     454            IF(lwp .AND. lflush) CALL flush(numout) 
    453455         endif  
    454456   
     
    842844         WRITE(numout,*) 'ldf_slp_init : direction of lateral mixing' 
    843845         WRITE(numout,*) '~~~~~~~~~~~~' 
     846         IF(lflush) CALL flush(numout) 
    844847      ENDIF 
    845848 
     
    865868         IF(ln_sco .AND.  (ln_traldf_hor .OR. ln_dynldf_hor )) THEN 
    866869            IF(lwp)   WRITE(numout,*) '          Horizontal mixing in s-coordinate: slope = slope of s-surfaces' 
     870            IF(lwp .AND. lflush) CALL flush(numout) 
    867871 
    868872            ! geopotential diffusion in s-coordinates on tracers and/or momentum 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r10759 r10774  
    175175        IF(lwp.and.ln_traldf_bilap) CALL ctl_stop(' SMAGORINSKY + BILAPLACIAN - UNSTABLE OR NON_CONSERVATIVE' ) 
    176176#endif 
    177  
     177      ! 
     178      IF(lwp .AND. lflush) CALL flush(numout) 
    178179      ! 
    179180   END SUBROUTINE ldf_tra_init 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c2d.h90

    r6486 r10774  
    3636         IF(lwp) WRITE(numout,*) ' ldf_tra_c2d : 2D eddy diffusivity and eddy' 
    3737         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~   --  induced velocity coefficients' 
     38         IF(lwp .AND. lflush) CALL flush(numout) 
    3839      ELSE 
    3940         IF(lwp) WRITE(numout,*) 
    4041         IF(lwp) WRITE(numout,*) ' ldf_tra2d : 2D eddy diffusivity coefficient' 
    4142         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~   --' 
     43         IF(lwp .AND. lflush) CALL flush(numout) 
    4244      ENDIF 
    4345 
     
    7577            IF(lwp) WRITE(numout,*) '               Constant values used for eddy diffusivity coefficients' 
    7678            IF(lwp) WRITE(numout,*) '               Variation lat/lon only for eddy induced velocity coefficients' 
     79            IF(lwp .AND. lflush) CALL flush(numout) 
    7780         ENDIF 
    7881 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c3d.h90

    r6486 r10774  
    4646         IF(lwp) WRITE(numout,*) '               Coefficients set to constant' 
    4747      ENDIF 
    48  
     48    
    4949 
    5050      ! biharmonic operator   (T-point) 
     
    121121 
    122122# endif 
     123 
     124   IF(lwp .AND. lflush) CALL flush(numout) 
     125 
    123126   END SUBROUTINE ldf_tra_c3d 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_smag.F90

    r9583 r10774  
    9999         IF(lwp) WRITE(numout,*) 
    100100         IF(lwp) WRITE(numout,*) 
     101         IF(lwp .AND. lflush) CALL flush(numout) 
    101102      ENDIF 
    102103 
     
    108109      ahtt(:,:,:) = rn_aht_0 
    109110       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 
    112116       ENDIF 
    113117 
     
    190194         CALL prihre( ahtw(:,:,1), jpi, jpj, 1, jpi, 1,   & 
    191195            &                                1, jpj, 1, 1.e-1, numout ) 
     196         IF(lflush) CALL flush(numout) 
    192197      ENDIF 
    193198ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90

    r10759 r10774  
    269269         WRITE(numout,*) '      choose the albedo parameterization                  nn_ice_alb = ', nn_ice_alb 
    270270         WRITE(numout,*) '      albedo of bare puddled ice                          rn_albice  = ', rn_albice 
     271         IF(lflush) CALL flush(numout) 
    271272      ENDIF 
    272273      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90

    r10269 r10774  
    183183      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 
    184184      IF(lwp) WRITE(numout,*) 
     185      IF(lflush) CALL flush(numout) 
    185186 
    186187      ncplmodel = kcplmodel 
     
    235236         WRITE(numout,*) ' multiexchg: nldi, nlei, nimpp =', nldi, nlei, nimpp 
    236237         WRITE(numout,*) ' multiexchg: nldj, nlej, njmpp =', nldj, nlej, njmpp 
     238         IF(lflush) CALL flush(numout) 
    237239      ENDIF 
    238240       
     
    313315                     CALL oasis_abort ( ssnd(ji)%nid(jc,jm), 'cpl_define', 'Failure in oasis_def_var' ) 
    314316                  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 
    317322               END DO 
    318323            END DO 
     
    351356                  END IF 
    352357#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 
    355362 
    356363                  ! Define 0D (Greenland or Antarctic ice mass) or 1D (river outflow) coupling fields 
     
    394401                     CALL oasis_abort ( srcv(ji)%nid(jc,jm), 'cpl_define', 'Failure in oasis_def_var' ) 
    395402                  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 
    399408               END DO 
    400409            END DO 
     
    447456                     WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
    448457                     WRITE(numout,*) '****************' 
     458                     IF(lflush) CALL flush(numout) 
    449459                  ENDIF 
    450460               ENDIF 
     
    491501                &        kinfo == OASIS_RecvOut .OR. kinfo == OASIS_FromRestOut 
    492502                
    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 
    494507                
    495508               IF ( llaction ) THEN 
     
    513526                     WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
    514527                     WRITE(numout,*) '****************' 
     528                     IF(lflush) CALL flush(numout) 
    515529                  ENDIF 
    516530 
     
    603617                     WRITE(numout,*) '     - Start of data is ', pdata(1:number_to_print) 
    604618                     WRITE(numout,*) '****************' 
     619                     IF(lflush) CALL flush(numout) 
    605620                  ENDIF 
    606621                   
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_rnf_1d.F90

    r10759 r10774  
    7070      IF(lwp) WRITE(numout,*) 'cpl_rnf_1d_init : initialization of river runoff coupling' 
    7171      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     72      IF(lwp .AND. lflush) CALL flush(numout) 
    7273       
    7374      REWIND(numnam_cfg) 
     
    8788      IF(lwp) WRITE(numout,*) '             Print river information = ',ln_print_river_info 
    8889      IF(lwp) WRITE(numout,*) ' ' 
     90      IF(lwp .AND. lflush) CALL flush(numout) 
    8991       
    9092      ! Assign space for river numbers 
     
    109111       
    110112      ! Print out the largest river number 
    111       IF ( ln_print_river_info ) THEN 
     113      IF ( ln_print_river_info .AND. lwp) THEN 
    112114         WRITE(numout,*) 'Maximum river number in input file = ',max_river 
     115         IF(lflush) CALL flush(numout) 
    113116      END IF 
    114117       
     
    137140      rivers%river_area(:) = REAL(rivers%river_area_c(:),wp) 
    138141       
    139       IF ( ln_print_river_info ) THEN 
     142      IF ( ln_print_river_info .AND. lwp) THEN 
    140143        WRITE(numout,*) 'Area of rivers 1 to 10 are ',rivers%river_area(1:10) 
     144        IF(lflush) CALL flush(numout) 
    141145      END IF 
    142146       
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90

    r10770 r10774  
    281281                     & 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 
    282282                  WRITE(numout, *) 'it_offset is : ',it_offset 
     283                  IF(lflush) CALL flush(numout) 
    283284               ENDIF 
    284285               ! temporal interpolation weights 
     
    293294                  WRITE(numout, clfmt) TRIM(sd(jf)%clvar), kt, REAL(isecsbc,wp)/rday, nyear, nmonth, nday,    & 
    294295                     &                 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) 
    295297               ENDIF 
    296298            ENDIF 
     
    414416 
    415417         clfmt = "('fld_init : time-interpolation for ', a, ' read previous record = ', i6, ' at time = ', f7.2, ' days')" 
    416          IF(lwp .AND. nprint > 0) 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 
    417421 
    418422      ENDIF 
     
    708712      IF(lwp) WRITE(numout,*) 'Dim size for ',TRIM(clvar),' is ', ilendta 
    709713      IF(lwp) WRITE(numout,*) 'Number of levels for ',TRIM(clvar),' is ', ipk 
     714      IF(lwp .AND. lflush) CALL flush(numout) 
    710715 
    711716      SELECT CASE( ipk ) 
     
    779784                     END DO 
    780785                     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,*)   & 
    782788                        &   '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 
    783791                  ENDIF 
    784792               ENDIF 
     
    941949               &                          ' data type: '      ,       sdf(jf)%cltype      ,   & 
    942950               &                          ' land/sea mask:'   , TRIM( sdf(jf)%lsmname    ) 
    943             IF(nprint > 4) call flush(numout) 
    944951         END DO 
     952         IF(lflush) CALL flush(numout) 
    945953      ENDIF 
    946954       
     
    10191027            ENDIF 
    10201028         END DO 
     1029  
     1030        IF(lflush) CALL flush(numout) 
     1031 
    10211032      ENDIF 
    10221033      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90

    r10759 r10774  
    9999            WRITE(numout,*) '   Namelist namsbc_apr : Atmospheric PRessure as extrenal forcing' 
    100100            WRITE(numout,*) '      ref. pressure: global mean Patm (T) or a constant (F)  ln_ref_apr = ', ln_ref_apr 
     101            IF(lflush) CALL flush(numout) 
    101102         ENDIF 
    102103         ! 
     
    114115            IF(lwp) WRITE(numout,*) '         Inverse barometer added to OBC ssh data' 
    115116         ENDIF 
     117          
     118         IF(lwp .AND. lflush) CALL flush(numout) 
     119 
    116120         IF( ( ln_apr_obc ) .AND. .NOT. lk_dynspg_ts )   & 
    117121            CALL ctl_stop( 'sbc_apr: use inverse barometer ssh at open boundary ONLY possible with time-splitting' ) 
     
    152156            ssh_ibb(:,:) = ssh_ib(:,:) 
    153157         ENDIF 
     158   
     159         IF(lwp .AND. lflush) CALL flush(numout) 
     160 
    154161      ENDIF 
    155162      !                                         ! ---------------------------------------- ! 
     
    159166         IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 
    160167         IF(lwp) WRITE(numout,*) '~~~~' 
     168         IF(lwp .AND. lflush) CALL flush(numout) 
    161169         IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
    162170         CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r10759 r10774  
    331331      ENDIF 
    332332 
     333      IF(lwp .AND. lflush) CALL flush(numout) 
    333334      !                                   ! allocate sbccpl arrays 
    334335      !IF( sbc_cpl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_cpl_alloc : unable to allocate arrays' ) 
     
    480481         IF(lwp) WRITE(numout,*) 
    481482         IF(lwp) WRITE(numout,*) '   runoffs received from oasis -> force ln_rnf = ', ln_rnf 
     483         IF(lwp .AND. lflush) CALL flush(numout) 
    482484      ENDIF 
    483485      ! 
     
    630632            WRITE(numout,*)'                  wind stress module' 
    631633            WRITE(numout,*) 
     634            IF(lflush) CALL flush(numout) 
    632635         ENDIF 
    633636      ENDIF 
     
    677680            WRITE(numout,*)'               fraction of solar net radiation absorbed in the first ocean level' 
    678681            WRITE(numout,*) 
     682            IF(lflush) CALL flush(numout) 
    679683         ENDIF 
    680684      ENDIF 
     
    907911            WRITE(numout,*)'               fraction of solar net radiation absorbed in the first ocean level' 
    908912            WRITE(numout,*) 
     913            IF(lflush) CALL flush(numout) 
    909914         ENDIF 
    910915      ENDIF 
     
    947952            WRITE(numout,*)'                  wind stress U,V components' 
    948953            WRITE(numout,*)'                  wind stress module' 
     954            IF(lflush) CALL flush(numout) 
    949955         ENDIF 
    950956      ENDIF 
     
    14711477         IF(lwp) WRITE(numout,*) 'Greenland icesheet mass rate of change (kg/s) is ', greenland_icesheet_mass_rate_of_change 
    14721478         IF(lwp) WRITE(numout,*) 'Greenland icesheet seconds lapsed since last change is ', greenland_icesheet_timelapsed 
     1479         IF(lwp .AND. lflush) CALL flush(numout) 
    14731480      ELSE IF ( nn_coupled_iceshelf_fluxes == 2 ) THEN 
    14741481         greenland_icesheet_mass_rate_of_change = rn_greenland_total_fw_flux 
     
    15181525         IF(lwp) WRITE(numout,*) 'Antarctica icesheet mass rate of change (kg/s) is ', antarctica_icesheet_mass_rate_of_change 
    15191526         IF(lwp) WRITE(numout,*) 'Antarctica icesheet seconds lapsed since last change is ', antarctica_icesheet_timelapsed 
     1527         IF(lwp .AND. lflush) CALL flush(numout) 
    15201528      ELSE IF ( nn_coupled_iceshelf_fluxes == 2 ) THEN 
    15211529         antarctica_icesheet_mass_rate_of_change = rn_antarctica_total_fw_flux 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90

    r6486 r10774  
    9999            WRITE(numout,*) '~~~~~~~' 
    100100            WRITE(numout,*) 
     101            IF(lflush) CALL flush(numout) 
    101102         ENDIF 
    102103         ! allocate sbcdcy arrays 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

    r10759 r10774  
    159159         CALL lbc_lnk( taum(:,:), 'T', 1. )   ;   CALL lbc_lnk( wndm(:,:), 'T', 1. ) 
    160160 
    161          IF( nitend-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) 
    162162            WRITE(numout,*)  
    163163            WRITE(numout,*) '        read daily momentum, heat and freshwater fluxes OK' 
     
    170170               CALL prihre( sf(jf)%fnow, jpi, jpj, 1, jpi, 20, 1, jpj, 10, zfact, numout ) 
    171171            END DO 
    172             CALL FLUSH(numout) 
     172            IF(lflush) CALL flush(numout) 
    173173         ENDIF 
    174174         ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r10754 r10774  
    8686            IF( kn_fwb == 2 )   WRITE(numout,*) '          adjusted from previous year budget' 
    8787            IF( kn_fwb == 3 )   WRITE(numout,*) '          fwf set to zero and spread out over erp area' 
     88            IF(lflush) CALL flush(numout) 
    8889         ENDIF 
    8990         ! 
     
    127128            IF(lwp)WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
    128129            IF(lwp)WRITE(numout,*)'          year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b 
     130            IF(lwp .AND. lflush) CALL flush(numout) 
    129131         ENDIF    
    130132         !                                         ! Update fwfold if new year start 
     
    201203               WRITE(numout,*)'   MIN(zerp_cor) = ', MINVAL(zerp_cor)  
    202204               WRITE(numout,*)'   MAX(zerp_cor) = ', MAXVAL(zerp_cor)  
     205               IF(lflush) CALL flush(numout) 
    203206            ENDIF 
    204207         ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r10770 r10774  
    182182      CALL wrk_alloc( jpi,jpj, ztmp1, ztmp2 ) 
    183183      ! 
    184       IF(lwp .AND. nprint>1) 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 
    185188 
    186189      ji_off = INT ( (jpiglo - nx_global) / 2 ) 
     
    342345 
    343346      IF( kt == nit000 )  THEN 
    344          IF(lwp .AND. nprint>1) 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 
    345351      ENDIF 
    346352 
     
    583589       
    584590      IF( kt == nit000 )  THEN 
    585          IF(lwp .AND. nprint>1) 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 
    586595      ENDIF 
    587596       
     
    820829      !!--------------------------------------------------------------------- 
    821830 
    822       IF(lwp .AND. nprint > 1) 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 
    823835 
    824836      CALL CICE_Finalize 
     
    943955         WRITE(numout,*)  
    944956         WRITE(numout,*) '        read forcing fluxes for CICE OK' 
    945          CALL FLUSH(numout) 
     957         IF(lflush) CALL flush(numout) 
    946958      ENDIF 
    947959 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90

    r10759 r10774  
    124124         IF ( lwp ) WRITE(numout,*) '        nn_gammablk = ', nn_gammablk  
    125125         IF ( lwp ) WRITE(numout,*) '        rn_tfri2    = ', rn_tfri2  
     126         IF(lwp .AND. lflush) CALL flush(numout) 
    126127         IF (ln_divisf) THEN       ! keep it in the namelist ??? used true anyway as for runoff ? (PM) 
    127128            rdivisf = 1._wp 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r10770 r10774  
    9999         WRITE(numout,*) 'sbc_init : surface boundary condition setting' 
    100100         WRITE(numout,*) '~~~~~~~~ ' 
     101         IF(lflush) CALL flush(numout) 
    101102      ENDIF 
    102103 
     
    160161      END SELECT 
    161162      ! 
     163      IF(lwp .AND. lflush) CALL flush(numout) 
     164      ! 
    162165      IF ( nn_components /= jp_iam_nemo .AND. .NOT. lk_oasis )   & 
    163166         &      CALL ctl_stop( 'STOP', 'sbc_init : OPA-SAS coupled via OASIS, but key_oasis3 disabled' ) 
     
    264267         IF( nn_components/= jp_iam_nemo )  & 
    265268            &                       WRITE(numout,*) '              + OASIS coupled SAS' 
     269         IF(lflush) CALL flush(numout) 
    266270      ENDIF 
    267271      ! 
     
    283287            WRITE(numout,*)"   OPA-SAS coupled via OASIS : nn_fsbc re-defined from OASIS namcouple ", nn_fsbc 
    284288            WRITE(numout,*) 
     289            IF(lflush) CALL flush(numout) 
    285290         ENDIF 
    286291      ENDIF 
     
    421426         IF( ln_rstart .AND.    &                               !* Restart: read in restart file 
    422427            & iom_varid( numror, 'utau_b', ldstop = .FALSE. ) > 0 ) THEN  
    423             IF(lwp .AND. nprint > 0) 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 
    424432            IF(nn_timing == 2)  CALL timing_start('iom_rstget') 
    425433            CALL iom_get( numror, jpdom_autoglo, 'utau_b', utau_b )   ! before i-stress  (U-point) 
     
    453461            &                    'at it= ', kt,' date= ', ndastp 
    454462            WRITE(numout,*) '~~~~' 
     463            IF(lflush) CALL flush(numout) 
    455464         ENDIF 
    456465         IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90

    r10770 r10774  
    159159            rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 
    160160         ENDIF 
     161  
     162         IF(lwp .AND. lflush) CALL flush(numout) 
     163 
    161164      ENDIF 
    162165      !                                                ! ---------------------------------------- ! 
     
    168171            &                    'at it= ', kt,' date= ', ndastp 
    169172            WRITE(numout,*) '~~~~' 
     173            IF(lflush) CALL flush(numout) 
    170174         ENDIF 
    171175         IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
     
    296300         WRITE(numout,*) '      depth of river mouth additional mixing     rn_hrnf      = ', rn_hrnf 
    297301         WRITE(numout,*) '      multiplicative factor for runoff           rn_rfact     = ', rn_rfact 
     302         IF(lflush) CALL flush(numout) 
    298303      ENDIF 
    299304      !                                   ! ================== 
     
    305310         IF(lwp) WRITE(numout,*) 
    306311         IF(lwp) WRITE(numout,*) '          runoffs inflow read in a file' 
     312         IF(lwp .AND. lflush) CALL flush(numout) 
    307313         IF( ierror > 0 ) THEN 
    308314            CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' )   ;   RETURN 
     
    316322         IF(lwp) WRITE(numout,*) 
    317323         IF(lwp) WRITE(numout,*) '          runoffs temperatures read in a file' 
     324         IF(lwp .AND. lflush) CALL flush(numout) 
    318325         ALLOCATE( sf_t_rnf(1), STAT=ierror  ) 
    319326         IF( ierror > 0 ) THEN 
     
    328335         IF(lwp) WRITE(numout,*) 
    329336         IF(lwp) WRITE(numout,*) '          runoffs salinities read in a file' 
     337         IF(lwp .AND. lflush) CALL flush(numout) 
    330338         ALLOCATE( sf_s_rnf(1), STAT=ierror  ) 
    331339         IF( ierror > 0 ) THEN 
     
    340348         IF(lwp) WRITE(numout,*) 
    341349         IF(lwp) WRITE(numout,*) '          runoffs depth read in a file' 
     350         IF(lwp .AND. lflush) CALL flush(numout) 
    342351         rn_dep_file = TRIM( cn_dir )//TRIM( sn_dep_rnf%clname ) 
    343352         IF( .NOT. sn_dep_rnf%ln_clim ) THEN   ;   WRITE(rn_dep_file, '(a,"_y",i4)' ) TRIM( rn_dep_file ), nyear    ! add year  
     
    380389         IF(lwp) WRITE(numout,*) '    depth over which runoffs is spread                        rn_dep_max = ', rn_dep_max 
    381390         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) 
    382392 
    383393         CALL iom_open( TRIM( sn_rnf%clname ), inum )    !  open runoff file 
     
    433443         IF( nn_rnf_depth_file == 1 ) THEN      !  save  output nb levels for runoff 
    434444            IF(lwp) WRITE(numout,*) '              create runoff depht file' 
     445            IF(lwp .AND. lflush) CALL flush(numout) 
    435446            CALL iom_open  ( TRIM( sn_dep_rnf%clname ), inum, ldwrt = .TRUE., kiolib = jprstlib ) 
    436447            CALL iom_rstput( 0, 0, inum, 'rodepth', h_rnf ) 
     
    470481         IF(lwp) WRITE(numout,*) '             - set to zero SSS damping       (if ln_ssr=T)' 
    471482         IF(lwp) WRITE(numout,*) '             - mixed upstream-centered       (if ln_traadv_cen2=T)' 
     483         IF(lwp .AND. lflush) CALL flush(numout) 
    472484         ! 
    473485         CALL rnf_mouth                             ! set river mouth mask 
     
    476488         IF(lwp) WRITE(numout,*) 
    477489         IF(lwp) WRITE(numout,*) '          No specific treatment at river mouths' 
     490         IF(lwp .AND. lflush) CALL flush(numout) 
    478491         rnfmsk  (:,:) = 0._wp 
    479492         rnfmsk_z(:)   = 0._wp 
     
    511524      IF(lwp) WRITE(numout,*) 'rnf_mouth : river mouth mask' 
    512525      IF(lwp) WRITE(numout,*) '~~~~~~~~~ ' 
     526      IF(lwp .AND. lflush) CALL flush(numout) 
    513527      ! 
    514528      cl_rnfile = TRIM( cn_dir )//TRIM( sn_cnf%clname ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90

    r9321 r10774  
    9292            IF(lwp) WRITE(numout,*) 
    9393            IF(lwp) WRITE(numout,*) '~~~~~~~   mean fields initialised to instantaneous values' 
     94            IF(lwp .AND. lflush) CALL flush(numout) 
    9495            zcoef = REAL( nn_fsbc - 1, wp ) 
    9596            ssu_m(:,:) = zcoef * ub(:,:,1) 
     
    156157               &                    'at it= ', kt,' date= ', ndastp 
    157158            IF(lwp) WRITE(numout,*) '~~~~~~~' 
     159            IF(lwp .AND. lflush) CALL flush(numout) 
    158160            zf_sbc = REAL( nn_fsbc, wp ) 
    159161            IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
     
    257259      ENDIF 
    258260      ! 
     261      IF(lwp .AND. lflush) CALL flush(numout) 
     262      ! 
    259263   END SUBROUTINE sbc_ssm_init 
    260264 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

    r10759 r10774  
    189189         WRITE(numout,*) '      flag to bound erp term                 ln_sssr_bnd = ', ln_sssr_bnd 
    190190         WRITE(numout,*) '      ABS(Max./Min.) erp threshold           rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day' 
     191         IF(lflush) CALL flush(numout) 
    191192      ENDIF 
    192193      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90

    r6486 r10774  
    7171               WRITE(numout,*) Wave(ntide(jk))%cname_tide, utide(jk), ftide(jk), v0tide(jk), omega_tide(jk) 
    7272            END DO 
     73            IF(lflush) CALL flush(numout) 
    7374         ENDIF 
    7475         ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90

    r10759 r10774  
    6060          WRITE(numout,*) 'tide_init : Initialization of the tidal components' 
    6161          WRITE(numout,*) '~~~~~~~~~ ' 
     62          IF(lflush) CALL flush(numout) 
    6263       ENDIF 
    6364       ! 
     
    9091          WRITE(numout,*) '                                     ln_tide_ramp = ', ln_tide_ramp  
    9192          WRITE(numout,*) '                                     rdttideramp  = ', rdttideramp 
     93          IF(lflush) CALL flush(numout) 
    9294       ENDIF 
    9395       IF( ln_tide_ramp.AND.((nitend-nit000+1)*rdt/rday < rdttideramp) )   & 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/SOL/solver.F90

    r10759 r10774  
    8686         WRITE(numout,*) '      optimal coefficient of sor         rn_sor     = ', rn_sor 
    8787         WRITE(numout,*) 
     88         IF(lflush) CALL flush(numout) 
    8889      ENDIF 
    8990      eps = rn_eps 
     
    100101      CASE ( 1 )                          ! preconditioned conjugate gradient solver 
    101102         IF(lwp) WRITE(numout,*) '   a preconditioned conjugate gradient solver is used' 
     103         IF(lwp .AND. lflush) CALL flush(numout) 
    102104         IF( jpr2di /= 0 .AND. jpr2dj /= 0 )   CALL ctl_stop( ' jpr2di and jpr2dj should be equal to zero' ) 
    103105         ! 
     
    105107         IF(lwp) WRITE(numout,*) '   a successive-over-relaxation solver with extra outer halo is used' 
    106108         IF(lwp) WRITE(numout,*) '   with jpr2di =', jpr2di, ' and  jpr2dj =', jpr2dj 
     109         IF(lwp .AND. lflush) CALL flush(numout) 
    107110         IF( .NOT. lk_mpp .AND. jpr2di /= 0 .AND. jpr2dj /= 0 ) THEN 
    108111             CALL ctl_stop( 'jpr2di and jpr2dj are not equal to zero',   & 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r10759 r10774  
    12581258            WRITE(numout,*) '             model does not use Conservative Temperature' 
    12591259         ENDIF 
     1260         IF(lflush) CALL flush(numout) 
    12601261      ENDIF 
    12611262      ! 
     
    12721273         IF(lwp) WRITE(numout,*) 
    12731274         IF(lwp) WRITE(numout,*) '          use of TEOS-10 equation of state (cons. temp. and abs. salinity)' 
     1275         IF(lwp .AND. lflush) CALL flush(numout) 
    12741276         ! 
    12751277         rdeltaS = 32._wp 
     
    14621464         IF(lwp) WRITE(numout,*) 
    14631465         IF(lwp) WRITE(numout,*) '          use of EOS-80 equation of state (pot. temp. and pract. salinity)' 
     1466         IF(lwp .AND. lflush) CALL flush(numout) 
    14641467         ! 
    14651468         rdeltaS = 20._wp 
     
    16621665            WRITE(numout,*) '             2nd cabbel. coef.     rn_nu      = ', rn_nu 
    16631666            WRITE(numout,*) '               Caution: rn_beta0=0 incompatible with ddm parameterization ' 
     1667            IF(lflush) CALL flush(numout) 
    16641668         ENDIF 
    16651669         ! 
     
    16811685      IF(lwp) WRITE(numout,*) '          rau0 * rcp                       rau0_rcp = ', rau0_rcp 
    16821686      IF(lwp) WRITE(numout,*) '          1. / ( rau0 * rcp )           r1_rau0_rcp = ', r1_rau0_rcp 
     1687      IF(lwp .AND. lflush) CALL flush(numout) 
    16831688      ! 
    16841689   END SUBROUTINE eos_init 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r10759 r10774  
    223223         WRITE(numout,*) '      upstream scheme within muscl   ln_traadv_msc_ups = ', ln_traadv_msc_ups 
    224224         WRITE(numout,*) '      TVD advection scheme with zts  ln_traadv_tvd_zts = ', ln_traadv_tvd_zts 
     225         IF(lflush) CALL flush(numout) 
    225226      ENDIF 
    226227 
     
    260261         IF( nadv ==  7 )   WRITE(numout,*) '         TVD ZTS   scheme is used' 
    261262         IF( nadv == -1 )   WRITE(numout,*) '         esopa test: use all advection scheme' 
     263         IF(lflush) CALL flush(numout) 
    262264      ENDIF 
    263265      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r9321 r10774  
    139139         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ ' 
    140140         IF(lwp) WRITE(numout,*) 
     141         IF(lwp .AND. lflush) CALL flush(numout) 
    141142         ! 
    142143         IF( .NOT. ALLOCATED( upsmsk ) )  THEN 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90

    r7179 r10774  
    9797         IF(lwp) WRITE(numout,*) 'tra_adv_eiv : eddy induced advection on ', cdtype,' :' 
    9898         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   add to velocity fields the eiv component' 
     99         IF(lwp .AND. lflush) CALL flush(numout) 
    99100# if defined key_diaeiv  
    100101         IF( cdtype == 'TRA') THEN 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90

    r10759 r10774  
    306306         WRITE(numout,*) '      =1 no MLE in case of convection ; =0 always MLE              nn_conv   = ', nn_conv 
    307307         WRITE(numout,*) '      Density difference used to define ML for FK              rn_rho_c_mle  = ', rn_rho_c_mle 
     308         IF(lflush) CALL flush(numout) 
    308309      ENDIF 
    309310      ! 
     
    317318            WRITE(numout,*) '   Mixed Layer Eddy parametrisation NOT used' 
    318319         ENDIF 
     320 
     321         IF(lflush) CALL flush(numout) 
     322 
    319323      ENDIF 
    320324      ! 
     
    325329         IF(lwp) WRITE(numout,*) '      ML buoyancy criteria = ', rb_c, ' m/s2 ' 
    326330         IF(lwp) WRITE(numout,*) '      associated ML density criteria defined in zdfmxl = ', rho_c, 'kg/m3' 
     331         IF(lwp .AND. lflush) CALL flush(numout) 
    327332         ! 
    328333         IF( nn_mle == 0 ) THEN           ! MLE array allocation & initialisation 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r7771 r10774  
    9999         IF(lwp) WRITE(numout,*) '~~~~~~~' 
    100100         IF(lwp) WRITE(numout,*) 
     101         IF(lwp .AND. lflush) CALL flush(numout) 
    101102         ! 
    102103         ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r7179 r10774  
    8484         IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype 
    8585         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 
     86         IF(lwp .AND. lflush) CALL flush(numout) 
    8687      ENDIF 
    8788      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r7179 r10774  
    101101         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    102102         IF(lwp) WRITE(numout,*) 
     103         IF(lwp .AND. lflush) CALL flush(numout) 
    103104      ENDIF 
    104105      l_trd = .FALSE. 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r10020 r10774  
    102102         IF(lwp) WRITE(numout,*) 'tra_adv_tvd : TVD advection scheme on ', cdtype 
    103103         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     104         IF(lwp .AND. lflush) CALL flush(numout) 
    104105         ! 
    105106      ENDIF 
     
    383384         IF(lwp) WRITE(numout,*) 'tra_adv_tvd_zts : TVD ZTS advection scheme on ', cdtype 
    384385         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     386         IF(lwp .AND. lflush) CALL flush(numout) 
    385387      ENDIF 
    386388      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r7179 r10774  
    101101         IF(lwp) WRITE(numout,*) 'tra_adv_ubs :  horizontal UBS advection scheme on ', cdtype 
    102102         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     103         IF(lwp .AND. lflush) CALL flush(numout) 
    103104      ENDIF 
    104105      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

    r10759 r10774  
    158158         WRITE(numout,*) '      Constant geothermal flux value               rn_geoflx_cst = ', rn_geoflx_cst 
    159159         WRITE(numout,*) 
     160         IF(lflush) CALL flush(numout) 
    160161      ENDIF 
    161162 
     
    196197         IF(lwp) WRITE(numout,*) '      *** no geothermal heat flux' 
    197198      ENDIF 
     199      IF(lwp .AND. lflush) CALL flush(numout) 
    198200      ! 
    199201   END SUBROUTINE tra_bbc_init 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r10759 r10774  
    358358         IF(lwp)  WRITE(numout,*) 'trabbl:bbl : Compute bbl velocities and diffusive coefficients in ', cdtype 
    359359         IF(lwp)  WRITE(numout,*) '~~~~~~~~~~' 
     360         IF(lflush) CALL flush(numout) 
    360361      ENDIF 
    361362      !                                        !* bottom variables (T, S, alpha, beta, depth, velocity) 
     
    530531         WRITE(numout,*) '          diffusive bbl coefficient           rn_ahtbbl  = ', rn_ahtbbl, ' m2/s' 
    531532         WRITE(numout,*) '          advective bbl coefficient           rn_gambbl  = ', rn_gambbl, ' s' 
     533         IF(lflush) CALL flush(numout) 
    532534      ENDIF 
    533535 
     
    538540         IF( nn_bbl_adv == 1 )    WRITE(numout,*) '       * Advective BBL using upper velocity' 
    539541         IF( nn_bbl_adv == 2 )    WRITE(numout,*) '       * Advective BBL using velocity = F( delta rho)' 
     542         IF(lflush) CALL flush(numout) 
    540543      ENDIF 
    541544 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

    r10759 r10774  
    214214         WRITE(numout,*) '      Damping file name               cn_resto  = ', cn_resto 
    215215         WRITE(numout,*) 
     216         IF(lflush) CALL flush(numout) 
    216217      ENDIF 
    217218 
     
    227228         CASE ( 2 )  ; IF(lwp) WRITE(numout,*) '   no tracer damping in the mixed layer' 
    228229         END SELECT 
     230         IF(lwp .AND. lflush) CALL flush(numout) 
    229231 
    230232         !TG: Initialisation of dtatsd - Would it be better to have dmpdta routine 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf.F90

    r6498 r10774  
    159159         WRITE(numout,*) '   see ldf_tra_init report for lateral mixing parameters' 
    160160         WRITE(numout,*) 
     161         IF(lflush) CALL flush(numout) 
    161162      ENDIF 
    162163 
     
    227228      IF( lk_esopa ) THEN 
    228229         IF(lwp) WRITE(numout,*) '          esopa control: use all lateral physics options' 
     230         IF(lwp .AND. lflush) CALL flush(numout) 
    229231         nldf = -1 
    230232      ENDIF 
     
    238240         IF( nldf ==  2 )   WRITE(numout,*) '          bilaplacian operator' 
    239241         IF( nldf ==  3 )   WRITE(numout,*) '          Rotated bilaplacian' 
     242         IF(lflush) CALL flush(numout) 
    240243      ENDIF 
    241244 
     
    278281         WRITE(numout,*) 'tra:ldf_ano : lateral diffusion acting on anomalies' 
    279282         WRITE(numout,*) '~~~~~~~~~~~' 
     283         IF(lflush) CALL flush(numout) 
    280284      ENDIF 
    281285 
     
    292296      zs_ref(:,:,:) = 35.0 * tmask(:,:,:) 
    293297      IF(lwp) WRITE(numout,*) '              homogeneous ocean T = ', zt0, ' S = ',zs0 
     298      IF(lwp .AND. lflush) CALL flush(numout) 
    294299 
    295300      ! Initialisation of gtui/gtvi in case of no cavity 
     
    359364         WRITE(numout,*) 'tra:ldf_ano : lateral diffusion acting on the full fields' 
    360365         WRITE(numout,*) '~~~~~~~~~~~' 
     366         IF(lflush) CALL flush(numout) 
    361367      ENDIF 
    362368   END SUBROUTINE ldf_ano 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90

    r7179 r10774  
    102102         IF(lwp) WRITE(numout,*) 'tra_ldf_bilap : iso-level biharmonic operator on ', cdtype 
    103103         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~' 
     104         IF(lwp .AND. lflush) CALL flush(numout) 
    104105      ENDIF 
    105106      !                                                          ! =========== 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

    r8280 r10774  
    8888         IF(lwp) WRITE(numout,*) 'tra_ldf_bilapg : horizontal biharmonic operator in s-coordinate on ', cdtype 
    8989         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 
     90         IF(lwp .AND. lflush) CALL flush(numout) 
    9091      ENDIF 
    9192 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90

    r7771 r10774  
    137137         IF(lwp) WRITE(numout,*) 'tra_ldf_iso : rotated laplacian diffusion operator on ', cdtype 
    138138         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     139         IF(lwp .AND. lflush) CALL flush(numout) 
    139140      ENDIF 
    140141      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso_grif.F90

    r7179 r10774  
    128128         IF(lwp) WRITE(numout,*) 'tra_ldf_iso_grif : rotated laplacian diffusion operator on ', cdtype 
    129129         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     130         IF(lwp .AND. lflush) CALL flush(numout) 
    130131         ALLOCATE( ah_wslp2(jpi,jpj,jpk) , zdkt3d(jpi,jpj,0:1), STAT=ierr ) 
    131132         IF( lk_mpp   )   CALL mpp_sum ( ierr ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90

    r7179 r10774  
    8585         IF(lwp) WRITE(numout,*) 'tra_ldf_lap : iso-level laplacian diffusion on ', cdtype 
    8686         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 
     87         IF(lwp .AND. lflush) CALL flush(numout) 
    8788      ENDIF 
    8889 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tranpc.F90

    r6486 r10774  
    164164                              END DO 
    165165                              WRITE(numout,*) 
     166                              IF(lflush) CALL flush(numout) 
    166167                           ENDIF 
    167168                            
     
    169170                           IF( jiter == 1 )   inpcc = inpcc + 1  
    170171 
    171                            IF( lp_monitor_point )   WRITE(numout, *) 'Negative N2 at ikp =',ikp,' for layer #', ilayer 
     172                           IF( lp_monitor_point )   THEN 
     173                              WRITE(numout, *) 'Negative N2 at ikp =',ikp,' for layer #', ilayer 
     174                              IF(lflush) CALL flush(numout) 
     175                           ENDIF 
    172176 
    173177                           !! ikup is the uppermost point where mixing will start: 
     
    223227                              WRITE(numout,*) '  => Mean Alfa  in that portion =', zalfa 
    224228                              WRITE(numout,*) '  => Mean Beta  in that portion =', zbeta 
     229                              IF(lflush) CALL flush(numout) 
    225230                           ENDIF 
    226231 
     
    285290                           END DO 
    286291                           WRITE(numout,*) 
     292                           IF(lflush) CALL flush(numout) 
    287293                        ENDIF 
    288294                        ! 
     
    323329            WRITE(numout,*) 'Exiting tra_npc , kt = ',kt,', => numb. of statically instable water-columns: ', inpcc 
    324330            WRITE(numout,*) 
     331            IF(lflush) CALL flush(numout) 
    325332         ENDIF 
    326333         ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r9163 r10774  
    105105         IF(lwp) WRITE(numout,*) 'tra_nxt : achieve the time stepping by Asselin filter and array swap' 
    106106         IF(lwp) WRITE(numout,*) '~~~~~~~' 
     107         IF(lwp .AND. lflush) CALL flush(numout) 
    107108         ! 
    108109         rbcp = 0.25_wp * (1._wp + atfp) * (1._wp + atfp) * ( 1._wp - atfp)      ! Brown & Campana parameter for semi-implicit hpg 
     
    243244         IF(lwp) WRITE(numout,*) 'tra_nxt_fix : time stepping', cdtype 
    244245         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     246         IF(lwp .AND. lflush) CALL flush(numout) 
    245247      ENDIF 
    246248      ! 
     
    317319         IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping', cdtype 
    318320         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     321         IF(lwp .AND. lflush) CALL flush(numout) 
    319322      ENDIF 
    320323      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

    r10770 r10774  
    119119         IF(lwp) WRITE(numout,*) 'tra_qsr : penetration of the surface solar radiation' 
    120120         IF(lwp) WRITE(numout,*) '~~~~~~~' 
     121         IF(lwp .AND. lflush) CALL flush(numout) 
    121122         IF( .NOT.ln_traqsr )   RETURN 
    122123      ENDIF 
     
    135136         IF( ln_rstart .AND.    &                    ! Restart: read in restart file 
    136137              & iom_varid( numror, 'qsr_hc_b', ldstop = .FALSE. ) > 0 ) THEN 
    137             IF(lwp .AND. nprint >0) WRITE(numout,*) '          nit000-1 qsr tracer content forcing field red in the restart file' 
     138            IF(lwp .AND. nprint >0) THEN 
     139               WRITE(numout,*) '          nit000-1 qsr tracer content forcing field red in the restart file' 
     140               IF(lflush) CALL flush(numout) 
     141            ENDIF 
    138142            zfact = 0.5e0 
    139143            CALL iom_get( numror, jpdom_autoglo, 'qsr_hc_b', qsr_hc_b )   ! before heat content trend due to Qsr flux 
     
    369373            &                    'at it= ', kt,' date= ', ndastp 
    370374            WRITE(numout,*) '~~~~' 
     375            IF(lflush) CALL flush(numout) 
    371376         ENDIF 
    372377         IF(nn_timing == 2)  CALL timing_start('iom_rstput')  
     
    456461         WRITE(numout,*) '      RGB & 2 bands: shortess depth of extinction  rn_si0 = ', rn_si0 
    457462         WRITE(numout,*) '      2 bands: longest depth of extinction         rn_si1 = ', rn_si1 
     463         IF(lflush) CALL flush(numout) 
    458464      ENDIF 
    459465 
     
    487493            IF( nqsr ==  4 )   WRITE(numout,*) '         2 bands light penetration' 
    488494            IF( nqsr ==  5 )   WRITE(numout,*) '         bio-model light penetration' 
     495            IF(lflush) CALL flush(numout) 
    489496         ENDIF 
    490497         ! 
     
    569576            ENDIF 
    570577            ! 
     578            IF(lwp .AND. lflush) CALL flush(numout) 
     579            ! 
    571580         ENDIF 
    572581            !                             ! ---------------------------------- ! 
     
    578587            IF(lwp) THEN 
    579588               WRITE(numout,*) 
    580             IF(lwp) WRITE(numout,*) '        level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m' 
     589               WRITE(numout,*) '        level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m' 
     590               IF(lflush) CALL flush(numout) 
    581591            ENDIF 
    582592            ! 
    583593            IF( lk_vvl ) THEN                   ! variable volume 
    584594               IF(lwp) WRITE(numout,*) '        key_vvl: light distribution will be computed at each time step' 
     595               IF(lwp .AND. lflush) CALL flush(numout) 
    585596            ELSE                                ! constant volume: computes one for all 
    586597               zz0 =        rn_abs   * r1_rau0_rcp 
     
    606617            WRITE(numout,*) 'tra_qsr_init : NO solar flux penetration' 
    607618            WRITE(numout,*) '~~~~~~~~~~~~' 
     619            IF(lflush) CALL flush(numout) 
    608620         ENDIF 
    609621      ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r10770 r10774  
    133133         IF(lwp) WRITE(numout,*) 'tra_sbc : TRAcer Surface Boundary Condition' 
    134134         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     135         IF(lwp .AND. lflush) CALL flush(numout) 
    135136      ENDIF 
    136137 
     
    157158              & iom_varid( numror, 'sbc_hc_b', ldstop = .FALSE. ) > 0 ) THEN 
    158159            IF(lwp) WRITE(numout,*) '          nit000-1 surface tracer content forcing fields red in the restart file' 
     160            IF(lwp .AND. lflush) CALL flush(numout) 
    159161            zfact = 0.5_wp 
    160162            CALL iom_get( numror, jpdom_autoglo, 'sbc_hc_b', sbc_tsc_b(:,:,jp_tem) )   ! before heat content sbc trend 
     
    214216            &                    'at it= ', kt,' date= ', ndastp 
    215217            WRITE(numout,*) '~~~~' 
     218            IF(lflush) CALL flush(numout) 
    216219         ENDIF 
    217220         IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
     
    262265               &                    'at it= ', kt,' date= ', ndastp 
    263266                WRITE(numout,*) '~~~~' 
     267                IF(lflush) CALL flush(numout) 
    264268            ENDIF 
    265269            IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf.F90

    r8104 r10774  
    163163         IF( nzdf ==  0 )   WRITE(numout,*) '              Explicit time-splitting scheme' 
    164164         IF( nzdf ==  1 )   WRITE(numout,*) '              Implicit (euler backward) scheme' 
     165         IF(lflush) CALL flush(numout) 
    165166      ENDIF 
    166167      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90

    r6486 r10774  
    100100         IF(lwp) WRITE(numout,*) 'tra_zdf_exp : explicit vertical mixing on ', cdtype 
    101101         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     102         IF(lwp .AND. lflush) CALL flush(numout) 
    102103      ENDIF 
    103104 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r6486 r10774  
    9999         IF(lwp)WRITE(numout,*) 'tra_zdf_imp : implicit vertical mixing on ', cdtype 
    100100         IF(lwp)WRITE(numout,*) '~~~~~~~~~~~ ' 
     101         IF(lwp .AND. lflush) CALL flush(numout) 
    101102         ! 
    102103         IF( lk_vvl ) THEN   ;    r_vvl = 1._wp       ! Variable volume indicator 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdglo.F90

    r6486 r10774  
    285285            WRITE (numout,9513) umo(jpdyn_tau) / tvolu, vmo(jpdyn_tau) / tvolv 
    286286            IF( ln_bfrimp )   WRITE (numout,9514) umo(jpdyn_bfri) / tvolu, vmo(jpdyn_bfri) / tvolv 
     287            IF(lflush) CALL flush(numout) 
    287288         ENDIF 
    288289 
     
    324325            WRITE (numout,9533) hke(jpdyn_tau) / tvolt 
    325326            IF( ln_bfrimp )   WRITE (numout,9534) hke(jpdyn_bfri) / tvolt 
     327            IF(lflush) CALL flush(numout) 
    326328         ENDIF 
    327329 
     
    356358            WRITE (numout,*) 
    357359            WRITE (numout,*) 
     360            IF(lflush) CALL flush(numout) 
    358361         ENDIF 
    359362 
     
    424427            &                   (  smo(jptra_xad) + smo(jptra_yad) + smo(jptra_zad) + smo(jptra_ldf) + smo(jptra_zdf)   & 
    425428            &                    + smo(jptra_npc) + smo(jptra_dmp)                   + smo(jptra_nsr) ) / tvolt 
     429            IF(lflush) CALL flush(numout) 
    426430         ENDIF 
    427431 
     
    459463            &                   (  s2(jptra_xad) + s2(jptra_yad) + s2(jptra_zad) + s2(jptra_ldf) + s2(jptra_zdf)   & 
    460464            &                    + s2(jptra_npc) + s2(jptra_dmp)                  + s2(jptra_nsr) ) / tvolt 
     465            IF(lflush) CALL flush(numout) 
    461466         ENDIF 
    462467 
     
    491496            WRITE (numout,9448)   t2(jptra_zdf)/tvolt,   s2(jptra_zdf)/tvolt 
    492497            WRITE (numout,9449)   t2(jptra_npc)/tvolt,   s2(jptra_npc)/tvolt 
     498            IF(lflush) CALL flush(numout) 
    493499         ENDIF 
    494500 
     
    524530         WRITE(numout,*) 'trd_glo_init : integral constraints properties trends' 
    525531         WRITE(numout,*) '~~~~~~~~~~~~~' 
     532         IF(lflush) CALL flush(numout) 
    526533      ENDIF 
    527534 
     
    534541 
    535542      IF(lwp) WRITE(numout,*) '                total ocean volume at T-point   tvolt = ',tvolt 
     543      IF(lwp .AND. lflush) CALL flush(numout) 
    536544 
    537545      ! Initialization of potential to kinetic energy conversion 
     
    557565         WRITE(numout,*) '                total ocean volume at U-point   tvolu = ',tvolu 
    558566         WRITE(numout,*) '                total ocean volume at V-point   tvolv = ',tvolv 
     567         IF(lflush) CALL flush(numout) 
    559568      ENDIF 
    560569      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdini.F90

    r10759 r10774  
    7171         ! 
    7272         WRITE(numout,*) '      frequency of trends diagnostics (glo)     nn_trd      = ', nn_trd 
     73         IF(lflush) CALL flush(numout) 
    7374      ENDIF 
    7475      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90

    r10020 r10774  
    287287         WRITE(numout,*) 'trd_ken_init : 3D Kinetic Energy trends' 
    288288         WRITE(numout,*) '~~~~~~~~~~~~~' 
     289         IF(lflush) CALL flush(numout) 
    289290      ENDIF 
    290291      !                           ! allocate box volume arrays 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl.F90

    r10759 r10774  
    393393         IF( ln_ctl ) THEN 
    394394            WRITE(numout,*) '             we reach kt == nit000 + 1 = ', nit000+1 
     395            IF(lflush) CALL flush(numout) 
    395396            CALL prt_ctl(tab2d_1=tmlbb   , clinfo1=' tmlbb   -   : ', mask1=tmask, ovlap=1) 
    396397            CALL prt_ctl(tab2d_1=tmlbn   , clinfo1=' tmlbn   -   : ', mask1=tmask, ovlap=1) 
     
    403404         IF( ln_trdmxl_instant ) THEN 
    404405            WRITE(numout,*) '             restart from kt == nit000 = ', nit000 
     406            IF(lflush) CALL flush(numout) 
    405407            CALL prt_ctl(tab2d_1=tmlbb   , clinfo1=' tmlbb   -   : ', mask1=tmask, ovlap=1) 
    406408            CALL prt_ctl(tab2d_1=tmlbn   , clinfo1=' tmlbn   -   : ', mask1=tmask, ovlap=1) 
     
    408410         ELSE 
    409411            WRITE(numout,*) '             restart from kt == nit000 = ', nit000 
     412            IF(lflush) CALL flush(numout) 
    410413            CALL prt_ctl(tab2d_1=tmlbn          , clinfo1=' tmlbn           -  : ', mask1=tmask, ovlap=1) 
    411414            CALL prt_ctl(tab2d_1=hmxlbn         , clinfo1=' hmxlbn          -  : ', mask1=tmask, ovlap=1) 
     
    633636            WRITE(numout,*) '          TRA zsmlres (jpi/2,jpj/2) : ', zsmlres (jpi/2,jpj/2) 
    634637            WRITE(numout,*) '          TRA zsmlres2(jpi/2,jpj/2) : ', zsmlres2(jpi/2,jpj/2) 
     638            IF(lflush) CALL flush(numout) 
    635639            ! 
    636640         END IF 
     
    778782         WRITE(numout,*) '      unit conversion factor                     rn_ucf             = ', rn_ucf 
    779783         WRITE(numout,*) '      criteria to compute the MLD                rn_rho_c           = ', rn_rho_c 
     784         IF(lflush) CALL flush(numout) 
    780785      ENDIF 
    781786 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl_rst.F90

    r9321 r10774  
    6969            ELSE                          ;   WRITE(numout,*) '             kt = '             , kt,' date= ', ndastp 
    7070            ENDIF 
     71            IF(lflush) CALL flush(numout) 
    7172         ENDIF 
    7273 
     
    7980         WRITE(numout,*) '~~~~~~~~~~' 
    8081         WRITE(numout,*) 
     82         IF(lflush) CALL flush(numout) 
    8183      ENDIF 
    8284      IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
     
    145147         WRITE(numout,*) ' trd_mxl_rst_read : read the NetCDF mixed layer trend restart file' 
    146148         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~' 
     149         IF(lflush) CALL flush(numout) 
    147150      ENDIF 
    148151 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdpen.F90

    r7179 r10774  
    144144         WRITE(numout,*) 'trd_pen_init : 3D Potential ENergy trends' 
    145145         WRITE(numout,*) '~~~~~~~~~~~~~' 
     146         IF(lflush) CALL flush(numout) 
    146147      ENDIF 
    147148      !                           ! allocate box volume arrays 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90

    r6486 r10774  
    209209      IF( ndebug /= 0 ) THEN 
    210210         IF(lwp) WRITE(numout,*) ' debuging trd_vor_zint: I done' 
    211          CALL FLUSH(numout) 
     211         IF(lwp .AND. lflush) CALL flush(numout) 
    212212      ENDIF 
    213213      ! 
     
    305305      IF( ndebug /= 0 ) THEN 
    306306         IF(lwp) WRITE(numout,*) ' debuging trd_vor_zint: I done' 
    307          CALL FLUSH(numout) 
     307         IF(lwp .AND. lflush) CALL flush(numout) 
    308308      ENDIF 
    309309      ! 
     
    434434            WRITE(numout,*) 'trd_vor : write trends in the NetCDF file at kt = ', kt 
    435435            WRITE(numout,*) '~~~~~~~  ' 
     436            IF(lflush) CALL flush(numout) 
    436437         ENDIF 
    437438  
     
    452453         IF( ndebug /= 0 ) THEN 
    453454            WRITE(numout,*) ' debuging trd_vor: III.4 done' 
    454             CALL FLUSH(numout) 
     455            IF(lwp .AND. lflush) CALL flush(numout) 
    455456         ENDIF 
    456457         ! 
     
    497498         WRITE(numout,*) '               ##########################################################################' 
    498499         WRITE(numout,*) ' ' 
     500         IF(lflush) CALL flush(numout) 
    499501      ENDIF 
    500502 
     
    510512      IF( ndebug /= 0 ) THEN 
    511513         WRITE(numout,*) ' debuging trd_vor_init: I. done' 
    512          CALL FLUSH(numout) 
     514         IF(lflush) CALL flush(numout) 
    513515      ENDIF 
    514516 
     
    533535 
    534536      IF(lwp) WRITE(numout,*) '               netCDF initialization' 
     537      IF(lwp .AND. lflush) CALL flush(numout) 
    535538 
    536539      ! II.2 Compute julian date from starting date of the run 
     
    542545         &                   ' YEAR ', nyear,' MONTH '      , nmonth,   & 
    543546         &                   ' DAY ' , nday, 'Julian day : ', zjulian 
     547      IF(lwp .AND. lflush) CALL flush(numout) 
    544548 
    545549      ! II.3 Define the T grid trend file (nidvor) 
     
    547551      CALL dia_nam( clhstnam, nn_trd, 'vort' )                  ! filename 
    548552      IF(lwp) WRITE(numout,*) ' Name of NETCDF file ', clhstnam 
     553      IF(lwp .AND. lflush) CALL flush(numout) 
    549554      CALL histbeg( clhstnam, jpi, glamf, jpj, gphif,1, jpi,   &  ! Horizontal grid : glamt and gphit 
    550555         &          1, jpj, nit000-1, zjulian, rdt, nh_t, nidvor, domain_id=nidom, snc4chunks=snc4set ) 
     
    582587      IF( ndebug /= 0 ) THEN 
    583588         WRITE(numout,*) ' debuging trd_vor_init: II. done' 
    584          CALL FLUSH(numout) 
     589         IF(lflush) CALL flush(numout) 
    585590      ENDIF 
    586591      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

    r10759 r10774  
    105105         WRITE(numout,*) 'zdf_bfr : Set bottom friction coefficient (non-linear case)' 
    106106         WRITE(numout,*) '~~~~~~~~' 
     107         IF(lflush) CALL flush(numout) 
    107108      ENDIF 
    108109      ! 
     
    262263      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~' 
    263264      IF(lwp) WRITE(numout,*) '   Namelist nam_bfr : set bottom friction parameters' 
     265      IF(lwp .AND. lflush) CALL flush(numout) 
    264266      ! 
    265267      SELECT CASE (nn_bfr) 
     
    399401      ! 
    400402      IF(lwp) WRITE(numout,*) '      implicit bottom friction switch                ln_bfrimp  = ', ln_bfrimp 
     403      IF(lwp .AND. lflush) CALL flush(numout) 
    401404      ! 
    402405      !                              ! Make sure ln_zdfexp=.false. when use implicit bfr 
     
    430433                   WRITE(numout,*) 'BFR ', narea, nimpp+ji, njmpp+jj, ikbu 
    431434                   WRITE(numout,*) 'BFR ', ABS( bfrcoef2d(ji,jj) ), zfru 
     435                   IF(lflush) CALL flush(numout) 
    432436                ENDIF 
    433437                ictu = ictu + 1 
     
    437441                     WRITE(numout,*) 'BFR ', narea, nimpp+ji, njmpp+jj, ikbv 
    438442                     WRITE(numout,*) 'BFR ', bfrcoef2d(ji,jj), zfrv 
     443                     IF(lflush) CALL flush(numout) 
    439444                 ENDIF 
    440445                 ictv = ictv + 1 
     
    452457                      WRITE(numout,*) 'TFR ', narea, nimpp+ji, njmpp+jj, ikbu 
    453458                      WRITE(numout,*) 'TFR ', ABS( tfrcoef2d(ji,jj) ), zfru 
     459                      IF(lflush) CALL flush(numout) 
    454460                   ENDIF 
    455461                   ictu = ictu + 1 
     
    459465                      WRITE(numout,*) 'TFR ', narea, nimpp+ji, njmpp+jj, ikbv 
    460466                      WRITE(numout,*) 'TFR ', tfrcoef2d(ji,jj), zfrv 
     467                      IF(lflush) CALL flush(numout) 
    461468                   ENDIF 
    462469                   ictv = ictv + 1 
     
    483490         IF ( ln_isfcav ) WRITE(numout,*) ' Top friction coefficient now ranges from: ', zmintfr, ' to ', zmaxtfr 
    484491         WRITE(numout,*) ' Bottom/Top friction coefficient will be reduced where necessary' 
     492         IF(lflush) CALL flush(numout) 
    485493      ENDIF 
    486494      ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfddm.F90

    r10759 r10774  
    255255         WRITE(numout,*) '      maximum avs for dd mixing      rn_avts   = ', rn_avts 
    256256         WRITE(numout,*) '      heat/salt buoyancy flux ratio  rn_hsbfr  = ', rn_hsbfr 
     257         IF(lflush) CALL flush(numout) 
    257258      ENDIF 
    258259      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90

    r7061 r10774  
    6969         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
    7070         IF(lwp) WRITE(numout,*) 
     71         IF(lwp .AND. lflush) CALL flush(numout) 
    7172      ENDIF 
    7273 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r10759 r10774  
    903903         WRITE(numout,*) '      Surface roughness (m)                         rn_hsro        = ', rn_hsro 
    904904         WRITE(numout,*) '      Bottom roughness (m) (nambfr namelist)        rn_bfrz0       = ', rn_bfrz0 
     905         IF(lflush) CALL flush(numout) 
    905906      ENDIF 
    906907 
     
    990991         ! 
    991992      END SELECT 
    992  
     993      ! 
     994      IF(lwp .AND. lflush) CALL flush(numout) 
    993995      ! 
    994996      SELECT CASE ( nn_stab_func )     !* set the parameters of the stability functions 
     
    10661068         ! 
    10671069      END SELECT 
    1068      
     1070      ! 
     1071      IF(lwp .AND. lflush) CALL flush(numout)  
    10691072      !                                !* Set Schmidt number for psi diffusion in the wave breaking case 
    10701073      !                                     ! See Eq. (13) of Carniel et al, OM, 30, 225-239, 2009 
     
    11181121         WRITE(numout,*) 'rl_sf   = ',rl_sf 
    11191122         WRITE(numout,*) 
     1123         IF(lflush) CALL flush(numout) 
    11201124      ENDIF 
    11211125 
     
    11941198            ELSE                         
    11951199               IF(lwp) WRITE(numout,*) ' ===>>>> : previous run without gls scheme, en and mxln computed by iterative loop' 
     1200               IF(lwp .AND. lflush) CALL flush(numout) 
    11961201               en  (:,:,:) = rn_emin 
    11971202               mxln(:,:,:) = 0.05         
     
    12041209         ELSE                                   !* Start from rest 
    12051210            IF(lwp) WRITE(numout,*) ' ===>>>> : Initialisation of en and mxln by background values' 
     1211            IF(lwp .AND. lflush) CALL flush(numout) 
    12061212            en  (:,:,:) = rn_emin 
    12071213            mxln(:,:,:) = 0.05        
     
    12111217         !                                   ! ------------------- 
    12121218         IF(lwp) WRITE(numout,*) '---- gls-rst ----' 
     1219         IF(lwp .AND. lflush) CALL flush(numout) 
    12131220         IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
    12141221         CALL iom_rstput( kt, nitrst, numrow, 'en'   , en     )  
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfini.F90

    r10759 r10774  
    8383         WRITE(numout,*) '      npc call  frequency                 nn_npc    = ', nn_npc 
    8484         WRITE(numout,*) '      npc print frequency                 nn_npcp   = ', nn_npcp 
     85         IF(lflush) CALL flush(numout) 
    8586      ENDIF 
    8687 
     
    154155         CALL ctl_stop( ' except for TKE, GLS or KPP physics, a convection scheme is',   & 
    155156         &              ' required: ln_zdfevd or ln_zdfnpc logicals' ) 
    156  
     157      ! 
     158      IF(lwp .AND. lflush) CALL flush(numout) 
    157159      !                               !* Background eddy viscosity and diffusivity profil 
    158160      IF( nn_avb == 0 ) THEN                ! Define avmb, avtb from namelist parameter 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfkpp.F90

    r10759 r10774  
    12141214         IF(lwp) WRITE(numout,*) 'tra_kpp : KPP non-local tracer fluxes' 
    12151215         IF(lwp) WRITE(numout,*) '~~~~~~~   ' 
     1216         IF(lflush) CALL flush(numout) 
    12161217      ENDIF 
    12171218 
     
    12861287         IF(lwp) WRITE(numout,*) 'trc_kpp : KPP non-local tracer fluxes' 
    12871288         IF(lwp) WRITE(numout,*) '~~~~~~~   ' 
     1289         IF(lflush) CALL flush(numout) 
    12881290      ENDIF 
    12891291 
     
    13841386         WRITE(numout,*) '     max. mix. in interior convec.                 rn_difcon   = ', rn_difcon 
    13851387         WRITE(numout,*) '     horizontal average flag                       nn_ave      = ', nn_ave 
     1388         IF(lflush) CALL flush(numout) 
    13861389      ENDIF 
    13871390 
     
    14101413         ENDIF 
    14111414      ENDIF 
     1415     
     1416      IF(lwp .AND. lflush) CALL flush(numout) 
    14121417       
    14131418 
     
    14211426         WRITE(numout,*) '     Constant value for unreso. turbul. velocity shear Vtc = ', Vtc 
    14221427         WRITE(numout,*) '     Non-dimensional coef. for nonlocal transport      rcg = ', rcg 
     1428         IF(lflush) CALL flush(numout) 
    14231429       ENDIF 
    14241430 
     
    14381444         IF(lwp) WRITE(numout,*) '          no horizontal average on avt, avmu, avmv' 
    14391445         IF(lwp) WRITE(numout,*) '          only in very high horizontal resolution !' 
     1446         IF(lwp .AND. lflush) CALL flush(numout) 
    14401447         ! weighting mean arrays etmean, eumean and evmean 
    14411448         !           ( 1  1 )                                          ( 1 ) 
     
    14641471      CASE ( 1 )                ! horizontal average  
    14651472         IF(lwp) WRITE(numout,*) '          horizontal average on avt, avmu, avmv' 
     1473         IF(lwp .AND. lflush) CALL flush(numout) 
    14661474         ! weighting mean arrays etmean, eumean and evmean 
    14671475         !           ( 1/2  1  1/2 )              ( 1/2  1/2 )             ( 1/2  1  1/2 ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r10759 r10774  
    108108         IF(lwp) WRITE(numout,*) 'zdf_mxl : mixed layer depth' 
    109109         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     110         IF(lflush) CALL flush(numout) 
    110111         !                             ! allocate zdfmxl arrays 
    111112         IF( zdf_mxl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'zdf_mxl : unable to allocate arrays' ) 
     
    396397            END DO 
    397398         END DO 
    398          WRITE(numout,*) 'zthick_0(jk =',jk,') =',zthick_0(2,2) 
    399          WRITE(numout,*) 'fsdepw(jk+1 =',jk+1,') =',fsdepw(2,2,jk+1) 
     399         IF(lwp) THEN 
     400            WRITE(numout,*) 'zthick_0(jk =',jk,') =',zthick_0(2,2) 
     401            WRITE(numout,*) 'fsdepw(jk+1 =',jk+1,') =',fsdepw(2,2,jk+1) 
     402         ENDIF 
    400403      END DO 
     404       
     405      IF(lflush) CALL flush(numout) 
    401406 
    402407      ! Surface boundary condition 
     
    429434      END DO 
    430435 
    431       WRITE(numout,*) 'htc_mld(after) =',htc_mld(2,2) 
     436      IF(lwp) WRITE(numout,*) 'htc_mld(after) =',htc_mld(2,2) 
     437      IF(lwp .AND. lflush) CALL flush(numout) 
    432438 
    433439      ! Convert to heat content 
     
    491497            END DO 
    492498            WRITE(numout,*) '====================================================================' 
     499            IF(lflush) CALL flush(numout) 
    493500         ENDIF 
    494501      ENDIF 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

    r10759 r10774  
    276276         WRITE(numout,*) '      Vertical eddy Visc. in the ML  rn_wvmix  = ', rn_wvmix 
    277277         WRITE(numout,*) '      Use the MLD parameterization   ln_mldw   = ', ln_mldw 
     278         IF(lflush) CALL flush(numout) 
    278279      ENDIF 
    279280      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r10770 r10774  
    797797         WRITE(numout,*) 'zdf_tke_init :  New tidal mixing case: force rn_emin = 1.e-10 and rmxl_min = 1.e-3 ' 
    798798         WRITE(numout,*) '~~~~~~~~~~~~' 
     799         IF(lflush) CALL flush(numout) 
    799800      ENDIF 
    800801# else 
     
    825826         WRITE(numout,*) 
    826827         WRITE(numout,*) '      critical Richardson nb with your parameters  ri_cri = ', ri_cri 
     828         IF(lflush) CALL flush(numout) 
    827829      ENDIF 
    828830      ! 
     
    838840      IF( ln_mxl0 ) THEN 
    839841         IF(lwp) WRITE(numout,*) '   use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min' 
     842         IF(lwp .AND. lflush) CALL flush(numout) 
    840843         rn_mxl0 = rmxl_min 
    841844      ENDIF 
     
    949952           ELSE                                     ! No TKE array found: initialisation 
    950953              IF(lwp) WRITE(numout,*) ' ===>>>> : previous run without tke scheme, en computed by iterative loop' 
     954              IF(lwp .AND. lflush) CALL flush(numout) 
    951955              en (:,:,:) = rn_emin * tmask(:,:,:) 
    952956              CALL tke_avn                               ! recompute avt, avm, avmu, avmv and dissl (approximation) 
     
    971975     ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    972976        !                                   ! ------------------- 
    973         IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- tke-rst ----' 
     977        IF(lwp .AND. nprint > 0) THEN 
     978           WRITE(numout,*) '---- tke-rst ----' 
     979           IF(lflush) CALL flush(numout) 
     980        ENDIF 
    974981        IF(nn_timing == 2)  CALL timing_start('iom_rstput') 
    975982        CALL iom_rstput( kt, nitrst, numrow, 'en'   , en     ) 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftmx.F90

    r10759 r10774  
    156156         IF(lwp) WRITE(numout,*)  
    157157         IF(lwp) WRITE(numout,*) '          N Total power consumption by av_tide    : ztpc = ', ztpc * 1.e-12 ,'TW' 
     158         IF(lwp .AND. lflush) CALL flush(numout) 
    158159      ENDIF 
    159160        
     
    310311         ztpc= rau0 * ztpc / ( rn_me * rn_tfe_itf ) 
    311312         IF(lwp) WRITE(numout,*) '          N Total power consumption by zavt_itf: ztpc = ', ztpc * 1.e-12 ,'TW' 
     313         IF(lwp .AND. lflush) CALL flush(numout) 
    312314      ENDIF 
    313315 
     
    401403         WRITE(numout,*) '      ITF specific parameterisation         = ', ln_tmx_itf 
    402404         WRITE(numout,*) '      ITF tidal dissipation efficiency      = ', rn_tfe_itf 
     405         IF(lflush) CALL flush(numout) 
    403406      ENDIF 
    404407 
     
    468471         WRITE(numout,*)  
    469472         WRITE(numout,*) '          Total power consumption of the tidally driven part of Kz : ztpc = ', ztpc * 1.e-12 ,'TW' 
     473         IF(lflush) CALL flush(numout) 
    470474 
    471475 
     
    497501         END DO 
    498502         WRITE(numout,*) '          Min de zkz ', ztpc, ' Max = ', maxval(zkz(:,:) ) 
     503         IF(lflush) CALL flush(numout) 
    499504 
    500505         DO jk = 2, jpkm1 
     
    516521         ztpc= rau0 * 1/(rn_tfe * rn_me) * ztpc 
    517522         WRITE(numout,*) '          2 Total power consumption of the tidally driven part of Kz : ztpc = ', ztpc * 1.e-12 ,'TW' 
     523         IF(lflush) CALL flush(numout) 
    518524 
    519525         DO jk = 1, jpk 
     
    528534            WRITE(numout,*) '            N2 min - jk= ', jk,'   ', ze_z * 1.e4,' cm2/s min= ',ztpc*1.e4,   & 
    529535               &       'max= ', MAXVAL(zav_tide(:,:,jk) )*1.e4, ' cm2/s' 
     536            IF(lflush) CALL flush(numout) 
    530537         END DO 
    531538 
     
    533540         WRITE(numout,*)  
    534541         WRITE(numout,*) '          Initial profile of tidal vertical mixing' 
     542         IF(lflush) CALL flush(numout) 
     543 
    535544         DO jk = 1, jpk 
    536545            DO jj = 1,jpj 
     
    543552            WRITE(numout,*) '                jk= ', jk,'   ', ze_z * 1.e4,' cm2/s' 
    544553         END DO 
     554         IF(lflush) CALL flush(numout) 
    545555         DO jk = 1, jpk 
    546556            zkz(:,:) = az_tmx(:,:,jk) /rn_n2min 
     
    551561               &       'max= ', MAXVAL(zkz)*1.e4, ' cm2/s' 
    552562         END DO 
     563         IF(lflush) CALL flush(numout) 
    553564         ! 
    554565      ENDIF 
     
    862873            WRITE(numout,*) 
    863874            WRITE(numout,*) '      Total power consumption by av_wave: ztpc =  ', ztpc * 1.e-12_wp, 'TW' 
     875            IF(lflush) CALL flush(numout) 
    864876         ENDIF 
    865877      ENDIF 
     
    983995         WRITE(numout,*) '      Variable (T) or constant (F) mixing efficiency            = ', ln_mevar 
    984996         WRITE(numout,*) '      Differential internal wave-driven mixing (T) or not (F)   = ', ln_tsdiff 
     997         IF(lflush) CALL flush(numout) 
    985998      ENDIF 
    986999       
     
    9951008         WRITE(numout,*) '   Force the background value applied to avm & avt in TKE to be everywhere ',   & 
    9961009            &               'the viscous molecular value & a very small diffusive value, resp.' 
     1010         IF(lflush) CALL flush(numout) 
    9971011      ENDIF 
    9981012       
     
    10421056         WRITE(numout,*) '      Pycnocline-intensifed wave-breaking energy: ', zpyc * 1.e-12_wp, 'TW' 
    10431057         WRITE(numout,*) '      Critical slope wave-breaking energy:        ', zcri * 1.e-12_wp, 'TW' 
     1058         IF(lflush) CALL flush(numout) 
    10441059      ENDIF 
    10451060      ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r10759 r10774  
    250250      NAMELIST/namctl/ ln_ctl  ,sn_cfctl, nn_print, nn_ictls, nn_ictle,   & 
    251251         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    252          &             nn_bench, nn_timing 
     252         &             nn_bench, nn_timing, ln_flush 
    253253      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    254254         &             jpizoom, jpjzoom, jperio, ln_use_jattr 
     
    296296#endif 
    297297      nprint    = nn_print          ! convert DOCTOR namelist names into OLD names 
     298      lflush    = ln_flush          ! convert namelist variable to model variable 
    298299      !do it here because output.namelist.dyn is opened before nprint was set earlier 
    299300      ! 
     
    405406         END DO 
    406407         WRITE(numout,cform_aaa)                                         ! Flag AAAAAAA 
     408         IF(lflush) CALL flush(numout) 
    407409         ! 
    408410      ENDIF 
     
    510512      IF( lk_asminc     )   CALL asm_inc_init   ! Initialize assimilation increments 
    511513      IF(lwp) WRITE(numout,*) 'Euler time step switch is ', neuler 
     514      IF(lwp .AND. lflush) CALL flush(numout) 
    512515       
    513516      IF (nstop > 0) THEN 
     
    546549         WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr  
    547550         WRITE(numout,*) '      level of print                  nn_print   = ', nn_print 
     551         WRITE(numout,*) '      flush numout/stat               ln_flush   = ', ln_flush 
    548552         WRITE(numout,*) '      Start i indice for SUM control  nn_ictls   = ', nn_ictls 
    549553         WRITE(numout,*) '      End i indice for SUM control    nn_ictle   = ', nn_ictle 
     
    554558         WRITE(numout,*) '      benchmark parameter (0/1)       nn_bench   = ', nn_bench 
    555559         WRITE(numout,*) '      timing activated    (0/1)       nn_timing  = ', nn_timing 
     560         IF(lflush) CALL flush(numout) 
    556561      ENDIF 
    557562      ! 
     
    581586         WRITE(numout,*) '      lateral cond. type (between 0 and 6) jperio = ', jperio    
    582587         WRITE(numout,*) '      use file attribute if exists as i/p j-start ln_use_jattr = ', ln_use_jattr 
     588         IF(lflush) CALL flush(numout) 
    583589      ENDIF 
    584590      !                             ! Parameter control 
     
    596602         IF(lwp) WRITE(numout,*)'          - The total number of processors over which the' 
    597603         IF(lwp) WRITE(numout,*)'            print control will be done is ijsplt : ', ijsplt 
     604         IF(lwp .AND. lflush) CALL flush(numout) 
    598605         ! 
    599606         !                              ! indices used for the SUM control 
     
    767774         WRITE (numout, *) 'WARNING: factorisation of number of PEs failed' 
    768775         WRITE (numout, *) '       : using grid of ',num_pes,' x 1' 
     776         IF(lflush) CALL flush(numout) 
    769777         jpnj = 1 
    770778         jpni = num_pes 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

    r10748 r10774  
    7171         WRITE(numout,*) 'stp_ctl : time-stepping control' 
    7272         WRITE(numout,*) '~~~~~~~' 
     73         IF(lflush) CALL flush(numout) 
    7374         ! open time.step file with special treatment for SAS 
    7475         IF ( nn_components == jp_iam_sas ) THEN 
     
    99100      IF( lk_mpp )   CALL mpp_max( zumax )                 ! max over the global domain 
    100101      ! 
    101       IF( ll_colruns )   WRITE(numout,*) ' ==>> time-step= ',kt,' abs(U) max: ', zumax 
     102      IF( ll_colruns )   THEN 
     103         WRITE(numout,*) ' ==>> time-step= ',kt,' abs(U) max: ', zumax 
     104         IF(lflush) CALL flush(numout) 
     105      ENDIF 
    102106      ! 
    103107      IF( zumax > 20.e0 ) THEN 
     
    117121            WRITE(numout,*) 
    118122            WRITE(numout,*) '          output of last fields in numwso' 
     123            IF(lflush) CALL flush(numout) 
    119124         ENDIF 
    120125         kindic = -3 
     
    133138      IF( lk_mpp )   CALL mpp_min( zsmin )                ! min over the global domain 
    134139      ! 
    135       IF( MOD( kt, nwrite ) == 1 .AND. lwp )   WRITE(numout,*) ' ==>> time-step= ',kt,' SSS min:', zsmin 
     140      IF( MOD( kt, nwrite ) == 1 .AND. lwp )   THEN 
     141         WRITE(numout,*) ' ==>> time-step= ',kt,' SSS min:', zsmin 
     142         IF(lflush) CALL flush(numout) 
     143      ENDIF 
    136144      ! 
    137145      IF( zsmin < 0.) THEN  
     
    188196            WRITE(numout,9600) kt, ztmax, ii, ij 
    189197            WRITE(numout,*) 'stp_ctl:tracer anomaly: ***** END OF WARNING *****' 
     198            IF(lflush) CALL flush(numout) 
    190199         ENDIF 
    191200      ENDIF 
     
    205214            WRITE(numout,9700) kt, ztmin, ii, ij 
    206215            WRITE(numout,*) 'stp_ctl:tracer anomaly: ***** END OF WARNING *****' 
     216            IF(lflush) CALL flush(numout) 
    207217         ENDIF 
    208218      ENDIF 
     
    218228      IF( lk_dynspg_flt ) THEN      ! elliptic solver statistics (if required) 
    219229         ! 
    220          IF(ll_wrtruns) WRITE(numsol,9200) kt, niter, res, SQRT(epsr)/eps       ! Solver 
     230         IF(ll_wrtruns) THEN 
     231            WRITE(numsol,9200) kt, niter, res, SQRT(epsr)/eps       ! Solver 
     232            IF(lflush) CALL flush(numsol) 
     233         ENDIF 
    221234         ! 
    222235         IF( kindic < 0 .AND. zsmin > 0.e0 .AND. zumax <= 20.e0 ) THEN   ! create a abort file if problem found  
     
    228241               WRITE(numout,*) ' stpctl: output of last fields' 
    229242               WRITE(numout,*) ' ======  ' 
     243               IF(lflush) CALL flush(numout) 
    230244            ENDIF 
    231245         ENDIF 
     
    241255         ! 
    242256            WRITE(numsol,9300) kt, zssh2, zumax, zsmin      ! ssh statistics 
     257            IF(lflush) CALL flush(numsol) 
    243258         ENDIF 
    244259         ! 
  • branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

    r10754 r10774  
    115115         WRITE(numout,*) 'trc_oce_rgb : Initialisation of the optical look-up table' 
    116116         WRITE(numout,*) '~~~~~~~~~~~ ' 
     117         IF(lflush) CALL flush(numout) 
    117118      ENDIF 
    118119      ! 
     
    184185      r_si2 = 1.e0 / zrgb(2, 1)        ! blue with the smallest chlorophyll concentration) 
    185186      IF(lwp) WRITE(numout,*) '      RGB longest depth of extinction    r_si2 = ', r_si2 
     187      IF(lwp .AND. lflush) CALL flush(numout) 
    186188      ! 
    187189      DO jc = 1, 61                         ! check 
     
    222224         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~' 
    223225         WRITE(numout,*)  
     226         IF(lflush) CALL flush(numout) 
    224227      ENDIF 
    225228      ! 
     
    229232         irgb = NINT( 41 + 20.* LOG10( zchl ) + 1.e-15 )    
    230233         IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  irgb = ', irgb   
     234         IF(lwp .AND. lflush) CALL flush(numout) 
    231235         IF( irgb /= jc ) THEN   
    232236            IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  Chl class = ', irgb 
     
    238242      r_si2 = 1.e0 / prgb(1, 1)      ! blue with the smallest chlorophyll concentration) 
    239243      IF(lwp) WRITE(numout,*) '      RGB longest depth of extinction    r_si2 = ', r_si2 
     244      IF(lwp .AND. lflush) CALL flush(numout) 
    240245      ! 
    241246   END SUBROUTINE trc_oce_rgb_read 
Note: See TracChangeset for help on using the changeset viewer.