Changeset 10986


Ignore:
Timestamp:
2019-05-16T15:23:56+02:00 (17 months ago)
Author:
andmirek
Message:

GMED 462 add flush

Location:
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE
Files:
128 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/ASM/asmbkg.F90

    r10888 r10986  
    8282         ! 
    8383         IF( .NOT. llok ) THEN 
    84             IF(lwp) WRITE(numout,*) ' Setting up assimilation background file '// TRIM( c_asmbkg ) 
     84            IF(lwp) THEN 
     85               WRITE(numout,*) ' Setting up assimilation background file '// TRIM( c_asmbkg ) 
     86               IF(lflush) CALL FLUSH(numout) 
     87            ENDIF 
    8588            ! 
    8689            !                                      ! Define the output file         
     
    9093               zdate = REAL( ndastp ) 
    9194               IF( ln_zdftke ) THEN                   ! read turbulent kinetic energy ( en ) 
    92                   IF(lwp) WRITE(numout,*) ' Reading TKE (en) from restart...' 
     95                  IF(lwp) THEN 
     96                     WRITE(numout,*) ' Reading TKE (en) from restart...' 
     97                     IF(lflush) CALL FLUSH(numout) 
     98                  ENDIF 
    9399                  CALL tke_rst( nit000, 'READ' ) 
    94100               ENDIF 
     
    120126         ! 
    121127         IF( .NOT. llok ) THEN 
    122             IF(lwp) WRITE(numout,*) ' Setting up assimilation background file '// TRIM( c_asmdin ) 
     128            IF(lwp) THEN 
     129               WRITE(numout,*) ' Setting up assimilation background file '// TRIM( c_asmdin ) 
     130               IF(lflush) CALL FLUSH(numout) 
     131            ENDIF 
    123132            ! 
    124133            !                                      ! Define the output file         
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/ASM/asminc.F90

    r10968 r10986  
    173173         WRITE(numout,*) '      Logical switch for ensuring that the sa > salfixmin      ln_salfix = ', ln_salfix 
    174174         WRITE(numout,*) '      Minimum salinity after applying the increments           salfixmin = ', salfixmin 
     175         IF(lflush) CALL FLUSH(numout) 
    175176      ENDIF 
    176177 
     
    209210         WRITE(numout,*) '       ditiaustr_date = ', ditiaustr_date 
    210211         WRITE(numout,*) '       ditiaufin_date = ', ditiaufin_date 
     212         IF(lflush) CALL FLUSH(numout) 
    211213      ENDIF 
    212214 
     
    311313             WRITE(numout,*) '         Time-integrated weight = ', ztotwgt 
    312314             WRITE(numout,*) '         ===================================' 
     315             IF(lflush) CALL FLUSH(numout) 
    313316          ENDIF 
    314317          
     
    347350            WRITE(numout,*) 'asm_inc_init : Assimilation increments valid between dates ', z_inc_dateb,' and ', z_inc_datef 
    348351            WRITE(numout,*) '~~~~~~~~~~~~' 
     352            IF(lflush) CALL FLUSH(numout) 
    349353         ENDIF 
    350354         ! 
     
    462466            WRITE(numout,*) '   ==>>>  Assimilation background state valid at : ', zdate_bkg 
    463467            WRITE(numout,*) 
     468            IF(lflush) CALL FLUSH(numout) 
    464469         ENDIF 
    465470         ! 
     
    491496      ENDIF 
    492497      ! 
    493       IF(lwp) WRITE(numout,*) '   ==>>>   Euler time step switch is ', neuler 
     498      IF(lwp) THEN 
     499         WRITE(numout,*) '   ==>>>   Euler time step switch is ', neuler 
     500         IF(lflush) CALL FLUSH(numout) 
     501      ENDIF 
    494502      ! 
    495503      IF( lk_asminc ) THEN                            !==  data assimilation  ==! 
     
    542550               WRITE(numout,*) 'tra_asm_inc : Tracer IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 
    543551               WRITE(numout,*) '~~~~~~~~~~~~' 
     552               IF(lflush) CALL FLUSH(numout) 
    544553            ENDIF 
    545554            ! 
     
    657666               WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 
    658667               WRITE(numout,*) '~~~~~~~~~~~~' 
     668               IF(lflush) CALL FLUSH(numout) 
    659669            ENDIF 
    660670            ! 
     
    728738                  &  kt,' with IAU weight = ', wgtiau(it) 
    729739               WRITE(numout,*) '~~~~~~~~~~~~' 
     740               IF(lflush) CALL FLUSH(numout) 
    730741            ENDIF 
    731742            ! 
     
    845856               WRITE(numout,*) 'seaice_asm_inc : sea ice conc IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 
    846857               WRITE(numout,*) '~~~~~~~~~~~~' 
     858               IF(lflush) CALL FLUSH(numout) 
    847859            ENDIF 
    848860            ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/BDY/bdydta.F90

    r10968 r10986  
    449449      !!--------------------------------------------------------------------------- 
    450450      ! 
    451       IF(lwp) WRITE(numout,*) 
    452       IF(lwp) WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries' 
    453       IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 
    454       IF(lwp) WRITE(numout,*) '' 
     451      IF(lwp) THEN 
     452         WRITE(numout,*) 
     453         WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries' 
     454         WRITE(numout,*) '~~~~~~~~~~' 
     455         WRITE(numout,*) '' 
     456         IF(lflush) CALL FLUSH(numout) 
     457      ENDIF 
    455458 
    456459      ! Set nn_dta 
     
    485488         ENDIF 
    486489#endif                
    487          IF(lwp) WRITE(numout,*) 'Maximum number of files to open =', nb_bdy_fld(jbdy) 
     490         IF(lwp) THEN 
     491            WRITE(numout,*) 'Maximum number of files to open =', nb_bdy_fld(jbdy) 
     492            IF(lflush) CALL FLUSH(numout) 
     493         ENDIF 
    488494      END DO             
    489495 
     
    532538 
    533539               IF( dta%ll_ssh ) THEN  
    534                   if(lwp) write(numout,*) '++++++ reading in ssh field' 
     540                  IF(lwp) THEN 
     541                     write(numout,*) '++++++ reading in ssh field' 
     542                     IF(lflush) CALL FLUSH(numout) 
     543                  ENDIF 
    535544                  jfld = jfld + 1 
    536545                  blf_i(jfld) = bn_ssh 
     
    543552 
    544553               IF( dta%ll_u2d .and. .not. ln_full_vel_array(jbdy) ) THEN 
    545                   if(lwp) write(numout,*) '++++++ reading in u2d field' 
     554                  IF(lwp) THEN 
     555                     write(numout,*) '++++++ reading in u2d field' 
     556                     IF(lflush) CALL FLUSH(numout) 
     557                  ENDIF 
    546558                  jfld = jfld + 1 
    547559                  blf_i(jfld) = bn_u2d 
     
    554566 
    555567               IF( dta%ll_v2d .and. .not. ln_full_vel_array(jbdy) ) THEN 
    556                   if(lwp) write(numout,*) '++++++ reading in v2d field' 
     568                  IF(lwp) THEN 
     569                     write(numout,*) '++++++ reading in v2d field' 
     570                     IF(lflush) CALL FLUSH(numout) 
     571                  ENDIF 
    557572                  jfld = jfld + 1 
    558573                  blf_i(jfld) = bn_v2d 
     
    572587 
    573588               IF( dta%ll_u3d .OR. ( ln_full_vel_array(jbdy) .and. dta%ll_u2d ) ) THEN 
    574                   if(lwp) write(numout,*) '++++++ reading in u3d field' 
     589                  IF(lwp) THEN 
     590                     write(numout,*) '++++++ reading in u3d field' 
     591                     IF(lflush) CALL FLUSH(numout) 
     592                  ENDIF 
    575593                  jfld = jfld + 1 
    576594                  blf_i(jfld) = bn_u3d 
     
    583601 
    584602               IF( dta%ll_v3d .OR. ( ln_full_vel_array(jbdy) .and. dta%ll_v2d ) ) THEN 
    585                   if(lwp) write(numout,*) '++++++ reading in v3d field' 
     603                  IF(lwp) THEN 
     604                     write(numout,*) '++++++ reading in v3d field' 
     605                     IF(lflush) CALL FLUSH(numout) 
     606                  ENDIF 
    586607                  jfld = jfld + 1 
    587608                  blf_i(jfld) = bn_v3d 
     
    599620 
    600621               IF( dta%ll_tem ) THEN 
    601                   if(lwp) write(numout,*) '++++++ reading in tem field' 
     622                  IF(lwp) THEN 
     623                     write(numout,*) '++++++ reading in tem field' 
     624                     IF(lflush) CALL FLUSH(numout) 
     625                  ENDIF 
    602626                  jfld = jfld + 1 
    603627                  blf_i(jfld) = bn_tem 
     
    609633 
    610634               IF( dta%ll_sal ) THEN 
    611                   if(lwp) write(numout,*) '++++++ reading in sal field' 
     635                  IF(lwp) THEN 
     636                     write(numout,*) '++++++ reading in sal field' 
     637                     IF(lflush) CALL FLUSH(numout) 
     638                  ENDIF 
    612639                  jfld = jfld + 1 
    613640                  blf_i(jfld) = bn_sal 
     
    732759            write(numout,*) '++++++ dta%ll_tem = ',dta%ll_tem 
    733760            write(numout,*) '++++++ dta%ll_sal = ',dta%ll_sal 
     761            IF(lflush) CALL FLUSH(numout) 
    734762         endif 
    735763 
    736764         IF ( nn_dyn2d_dta(jbdy) == 0 .or. nn_dyn2d_dta(jbdy) == 2 ) THEN 
    737             if(lwp) write(numout,*) '++++++ dta%ssh/u2d/u3d allocated space' 
     765            IF(lwp) THEN 
     766               write(numout,*) '++++++ dta%ssh/u2d/u3d allocated space' 
     767               IF(lflush) CALL FLUSH(numout) 
     768            ENDIF 
    738769            IF( dta%ll_ssh ) ALLOCATE( dta%ssh(nblen(1)) ) 
    739770            IF( dta%ll_u2d ) ALLOCATE( dta%u2d(nblen(2)) ) 
     
    742773         IF ( nn_dyn2d_dta(jbdy) == 1 .or. nn_dyn2d_dta(jbdy) == 3 ) THEN 
    743774            IF( dta%ll_ssh ) THEN 
    744                if(lwp) write(numout,*) '++++++ dta%ssh pointing to fnow' 
     775               IF(lwp) THEN 
     776                  write(numout,*) '++++++ dta%ssh pointing to fnow' 
     777                  IF(lflush) CALL FLUSH(numout) 
     778               ENDIF 
    745779               jfld = jfld + 1 
    746780               dta%ssh => bf(jfld)%fnow(:,1,1) 
     
    748782            IF ( dta%ll_u2d ) THEN 
    749783               IF ( ln_full_vel_array(jbdy) ) THEN 
    750                   if(lwp) write(numout,*) '++++++ dta%u2d allocated space' 
     784                  IF(lwp) THEN 
     785                     write(numout,*) '++++++ dta%u2d allocated space' 
     786                     IF(lflush) CALL FLUSH(numout) 
     787                  ENDIF 
    751788                  ALLOCATE( dta%u2d(nblen(2)) ) 
    752789               ELSE 
    753                   if(lwp) write(numout,*) '++++++ dta%u2d pointing to fnow' 
     790                  IF(lwp) THEN 
     791                     write(numout,*) '++++++ dta%u2d pointing to fnow' 
     792                     IF(lflush) CALL FLUSH(numout) 
     793                  ENDIF 
    754794                  jfld = jfld + 1 
    755795                  dta%u2d => bf(jfld)%fnow(:,1,1) 
     
    758798            IF ( dta%ll_v2d ) THEN 
    759799               IF ( ln_full_vel_array(jbdy) ) THEN 
    760                   if(lwp) write(numout,*) '++++++ dta%v2d allocated space' 
     800                  IF(lwp) THEN 
     801                     write(numout,*) '++++++ dta%v2d allocated space' 
     802                     IF(lflush) CALL FLUSH(numout) 
     803                  ENDIF 
    761804                  ALLOCATE( dta%v2d(nblen(3)) ) 
    762805               ELSE 
    763                   if(lwp) write(numout,*) '++++++ dta%v2d pointing to fnow' 
     806                  IF(lwp) THEN 
     807                     write(numout,*) '++++++ dta%v2d pointing to fnow' 
     808                     IF(lflush) CALL FLUSH(numout) 
     809                  ENDIF 
    764810                  jfld = jfld + 1 
    765811                  dta%v2d => bf(jfld)%fnow(:,1,1) 
     
    769815 
    770816         IF ( nn_dyn3d_dta(jbdy) == 0 ) THEN 
    771             if(lwp) write(numout,*) '++++++ dta%u3d/v3d allocated space' 
     817            IF(lwp) THEN 
     818               write(numout,*) '++++++ dta%u3d/v3d allocated space' 
     819               IF(lflush) CALL FLUSH(numout) 
     820            ENDIF 
    772821            IF( dta%ll_u3d ) ALLOCATE( dta_bdy(jbdy)%u3d(nblen(2),jpk) ) 
    773822            IF( dta%ll_v3d ) ALLOCATE( dta_bdy(jbdy)%v3d(nblen(3),jpk) ) 
     
    776825           &  ( ln_full_vel_array(jbdy) .and. ( nn_dyn2d_dta(jbdy) == 1 .or. nn_dyn2d_dta(jbdy) == 3 ) ) ) THEN 
    777826            IF ( dta%ll_u3d .or. ( ln_full_vel_array(jbdy) .and. dta%ll_u2d ) ) THEN 
    778                if(lwp) write(numout,*) '++++++ dta%u3d pointing to fnow' 
     827               IF(lwp) THEN 
     828                  write(numout,*) '++++++ dta%u3d pointing to fnow' 
     829                  IF(lflush) CALL FLUSH(numout) 
     830               ENDIF 
    779831               jfld = jfld + 1 
    780832               dta_bdy(jbdy)%u3d => bf(jfld)%fnow(:,1,:) 
    781833            ENDIF 
    782834            IF ( dta%ll_v3d .or. ( ln_full_vel_array(jbdy) .and. dta%ll_v2d ) ) THEN 
    783                if(lwp) write(numout,*) '++++++ dta%v3d pointing to fnow' 
     835               IF(lwp) THEN 
     836                  write(numout,*) '++++++ dta%v3d pointing to fnow' 
     837                  IF(lflush) CALL FLUSH(numout) 
     838               ENDIF 
    784839               jfld = jfld + 1 
    785840               dta_bdy(jbdy)%v3d => bf(jfld)%fnow(:,1,:) 
     
    788843 
    789844         IF( nn_tra_dta(jbdy) == 0 ) THEN 
    790             if(lwp) write(numout,*) '++++++ dta%tem/sal allocated space' 
     845            IF(lwp) THEN 
     846               write(numout,*) '++++++ dta%tem/sal allocated space' 
     847               IF(lflush) CALL FLUSH(numout) 
     848            ENDIF 
    791849            IF( dta%ll_tem ) ALLOCATE( dta_bdy(jbdy)%tem(nblen(1),jpk) ) 
    792850            IF( dta%ll_sal ) ALLOCATE( dta_bdy(jbdy)%sal(nblen(1),jpk) ) 
    793851         ELSE 
    794852            IF( dta%ll_tem ) THEN 
    795                if(lwp) write(numout,*) '++++++ dta%tem pointing to fnow' 
     853               IF(lwp) THEN 
     854                  write(numout,*) '++++++ dta%tem pointing to fnow' 
     855                  IF(lflush) CALL FLUSH(numout) 
     856               ENDIF 
    796857               jfld = jfld + 1 
    797858               dta_bdy(jbdy)%tem => bf(jfld)%fnow(:,1,:) 
    798859            ENDIF 
    799860            IF( dta%ll_sal ) THEN  
    800                if(lwp) write(numout,*) '++++++ dta%sal pointing to fnow' 
     861               IF(lwp) THEN 
     862                  write(numout,*) '++++++ dta%sal pointing to fnow' 
     863                  IF(lflush) CALL FLUSH(numout) 
     864               ENDIF 
    801865               jfld = jfld + 1 
    802866               dta_bdy(jbdy)%sal => bf(jfld)%fnow(:,1,:) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/BDY/bdyini.F90

    r10968 r10986  
    9191      ! ----------------------------------------- 
    9292      IF ( ln_bdy ) THEN 
    93          IF(lwp) WRITE(numout,*) 
    94          IF(lwp) WRITE(numout,*) 'bdy_init : initialization of open boundaries' 
    95          IF(lwp) WRITE(numout,*) '~~~~~~~~' 
     93         IF(lwp) THEN 
     94            WRITE(numout,*) 
     95            WRITE(numout,*) 'bdy_init : initialization of open boundaries' 
     96            WRITE(numout,*) '~~~~~~~~' 
     97            IF(lflush) CALL FLUSH(numout) 
     98         ENDIF 
    9699         ! 
    97100         ! Open boundaries definition (arrays and masks) 
     
    105108         ! 
    106109      ELSE 
    107          IF(lwp) WRITE(numout,*) 
    108          IF(lwp) WRITE(numout,*) 'bdy_init : open boundaries not used (ln_bdy = F)' 
    109          IF(lwp) WRITE(numout,*) '~~~~~~~~' 
     110         IF(lwp) THEN 
     111            WRITE(numout,*) 
     112            WRITE(numout,*) 'bdy_init : open boundaries not used (ln_bdy = F)' 
     113            WRITE(numout,*) '~~~~~~~~' 
     114            IF(lflush) CALL FLUSH(numout) 
     115         ENDIF 
    110116         ! 
    111117      ENDIF 
     
    403409        ENDIF 
    404410     ENDIF 
     411 
     412     IF(lflush .AND. lwp) CALL FLUSH(numout) 
    405413 
    406414      ! ------------------------------------------------- 
     
    568576               ! 
    569577               ibr_max = MAXVAL( nbrdta(:,igrd,ib_bdy) ) 
    570                IF(lwp) WRITE(numout,*) 
    571                IF(lwp) WRITE(numout,*) ' Maximum rimwidth in file is ', ibr_max 
    572                IF(lwp) WRITE(numout,*) ' nn_rimwidth from namelist is ', nn_rimwidth(ib_bdy) 
     578               IF(lwp) THEN 
     579                  WRITE(numout,*) 
     580                  WRITE(numout,*) ' Maximum rimwidth in file is ', ibr_max 
     581                  WRITE(numout,*) ' nn_rimwidth from namelist is ', nn_rimwidth(ib_bdy) 
     582                  IF(lflush) CALL FLUSH(numout) 
     583               ENDIF 
    573584               IF (ibr_max < nn_rimwidth(ib_bdy))   & 
    574585                     CALL ctl_stop( 'nn_rimwidth is larger than maximum rimwidth in file',cn_coords_file(ib_bdy) ) 
     
    12051216               IF( i_offset == 1 .and. zefl + zwfl == 2 ) THEN 
    12061217                  icount = icount + 1 
    1207                   IF(lwp) WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 
     1218                  IF(lwp) THEN 
     1219                     WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 
     1220                     IF(lflush) CALL FLUSH(numout) 
     1221                  ENDIF 
    12081222               ELSE 
    12091223                  idx_bdy(ib_bdy)%flagu(ib,igrd) = -zefl + zwfl 
     
    12371251               ! This error check only works if you are using the bdyXmask arrays 
    12381252               IF( j_offset == 1 .and. znfl + zsfl == 2 ) THEN 
    1239                   IF(lwp) WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 
     1253                  IF(lwp) THEN 
     1254                     WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 
     1255                     IF(lflush) CALL FLUSH(numout) 
     1256                  ENDIF 
    12401257                  icount = icount + 1 
    12411258               ELSE 
     
    12731290      !!---------------------------------------------------------------------- 
    12741291      ! 
    1275       IF (lwp) WRITE(numout,*) ' ' 
    1276       IF (lwp) WRITE(numout,*) 'bdy_ctl_seg: Check analytical segments' 
    1277       IF (lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    1278       ! 
    1279       IF(lwp) WRITE(numout,*) 'Number of east  segments     : ', nbdysege 
    1280       IF(lwp) WRITE(numout,*) 'Number of west  segments     : ', nbdysegw 
    1281       IF(lwp) WRITE(numout,*) 'Number of north segments     : ', nbdysegn 
    1282       IF(lwp) WRITE(numout,*) 'Number of south segments     : ', nbdysegs 
     1292      IF (lwp) THEN 
     1293         WRITE(numout,*) ' ' 
     1294         WRITE(numout,*) 'bdy_ctl_seg: Check analytical segments' 
     1295         WRITE(numout,*) '~~~~~~~~~~~~' 
     1296      ! 
     1297         WRITE(numout,*) 'Number of east  segments     : ', nbdysege 
     1298         WRITE(numout,*) 'Number of west  segments     : ', nbdysegw 
     1299         WRITE(numout,*) 'Number of north segments     : ', nbdysegn 
     1300         WRITE(numout,*) 'Number of south segments     : ', nbdysegs 
     1301         IF(lflush) CALL FLUSH(numout) 
     1302      ENDIF 
    12831303      ! 1. Check bounds 
    12841304      !---------------- 
     
    13221342      ! 2. Look for segment crossings 
    13231343      !------------------------------  
    1324       IF (lwp) WRITE(numout,*) '**Look for segments corners  :' 
     1344      IF (lwp) THEN 
     1345         WRITE(numout,*) '**Look for segments corners  :' 
     1346         IF(lflush) CALL FLUSH(numout) 
     1347      ENDIF 
    13251348      ! 
    13261349      itest = 0 ! corner number 
     
    14821505            ELSE 
    14831506               ! This is a corner 
    1484                IF(lwp) WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 
     1507               IF(lwp) THEN 
     1508                  WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 
     1509                  IF(lflush) CALL FLUSH(numout) 
     1510               ENDIF 
    14851511               CALL bdy_ctl_corn(npckgw(ib), icornw(ib,1)) 
    14861512               itest=itest+1 
     
    14941520            ELSE 
    14951521               ! This is a corner 
    1496                IF(lwp) WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 
     1522               IF(lwp) THEN 
     1523                  WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 
     1524                  IF(lflush) CALL FLUSH(numout) 
     1525               ENDIF 
    14971526               CALL bdy_ctl_corn(npckgw(ib), icornw(ib,2)) 
    14981527               itest=itest+1 
     
    15221551            ELSE 
    15231552               ! This is a corner 
    1524                IF(lwp) WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 
     1553               IF(lwp) THEN 
     1554                  WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 
     1555                  IF(lflush) CALL FLUSH(numout) 
     1556               ENDIF 
    15251557               CALL bdy_ctl_corn(npckge(ib), icorne(ib,1)) 
    15261558               itest=itest+1 
     
    15341566            ELSE 
    15351567               ! This is a corner 
    1536                IF(lwp) WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 
     1568               IF(lwp) THEN 
     1569                  WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 
     1570                  IF(lflush) CALL FLUSH(numout) 
     1571               ENDIF 
    15371572               CALL bdy_ctl_corn(npckge(ib), icorne(ib,2)) 
    15381573               itest=itest+1 
     
    15931628      END DO 
    15941629      ! 
    1595       IF ((itest==0).AND.(lwp)) WRITE(numout,*) 'NO open boundary corner found' 
     1630      IF ((itest==0).AND.(lwp)) THEN 
     1631         WRITE(numout,*) 'NO open boundary corner found' 
     1632         IF(lflush) CALL FLUSH(numout) 
     1633      ENDIF 
    15961634      ! 
    15971635      ! Other tests TBD:  
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/BDY/bdytides.F90

    r10968 r10986  
    8282      !!---------------------------------------------------------------------- 
    8383      ! 
    84       IF (nb_bdy>0) THEN 
    85          IF(lwp) WRITE(numout,*) 
    86          IF(lwp) WRITE(numout,*) 'bdytide_init : initialization of tidal harmonic forcing at open boundaries' 
    87          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     84      IF (nb_bdy>0 .AND. lwp) THEN 
     85         WRITE(numout,*) 
     86         WRITE(numout,*) 'bdytide_init : initialization of tidal harmonic forcing at open boundaries' 
     87         WRITE(numout,*) '~~~~~~~~~~~~' 
     88         IF(lflush) CALL FLUSH(numout) 
    8889      ENDIF 
    8990 
     
    107108            IF(lwm .AND. nprint > 2) WRITE ( numond, nambdy_tide ) 
    108109            !                                               ! Parameter control and print 
    109             IF(lwp) WRITE(numout,*) '  ' 
    110             IF(lwp) WRITE(numout,*) '          Namelist nambdy_tide : tidal harmonic forcing at open boundaries' 
    111             IF(lwp) WRITE(numout,*) '             read tidal data in 2d files: ', ln_bdytide_2ddta 
    112             IF(lwp) WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
    113             IF(lwp) WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
    114             IF(lwp) THEN  
    115                     WRITE(numout,*) '             Tidal components: '  
     110            IF(lwp) THEN 
     111               WRITE(numout,*) '  ' 
     112               WRITE(numout,*) '          Namelist nambdy_tide : tidal harmonic forcing at open boundaries' 
     113               WRITE(numout,*) '             read tidal data in 2d files: ', ln_bdytide_2ddta 
     114               WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
     115               WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
     116               WRITE(numout,*) '             Tidal components: '  
    116117               DO itide = 1, nb_harmo 
    117118                  WRITE(numout,*)  '                 ', Wave(ntide(itide))%cname_tide  
    118119               END DO 
     120               WRITE(numout,*) ' ' 
     121               IF(lflush) CALL FLUSH(numout) 
    119122            ENDIF  
    120             IF(lwp) WRITE(numout,*) ' ' 
    121123 
    122124            ! Allocate space for tidal harmonics data - get size from OBC data arrays 
     
    312314           WRITE(numout,*) 'bdytide_update : (re)Initialization of the tidal bdy forcing at kt=',kt 
    313315           WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
     316           IF(lflush) CALL FLUSH(numout) 
    314317        ENDIF 
    315318        ! 
     
    419422               ! 
    420423               IF(lwp) THEN 
    421                WRITE(numout,*) 
    422                WRITE(numout,*) 'bdy_tide_dta : Refresh nodal factors for tidal open bdy data at kt=',kt 
    423                WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
     424                  WRITE(numout,*) 
     425                  WRITE(numout,*) 'bdy_tide_dta : Refresh nodal factors for tidal open bdy data at kt=',kt 
     426                  WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
     427                  IF(lflush) CALL FLUSH(numout) 
    424428               ENDIF 
    425429               ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/BDY/bdyvol.F90

    r10888 r10986  
    163163 
    164164 
    165          IF(lwp) WRITE(numout,*) 
    166          IF(lwp) WRITE(numout,*)'bdy_vol : time step :', kt 
    167          IF(lwp) WRITE(numout,*)'~~~~~~~ ' 
    168          IF(lwp) WRITE(numout,*)'          cumulate flux EMP             =', z_cflxemp  , ' (m3/s)' 
    169          IF(lwp) WRITE(numout,*)'          total lateral surface of OBC  =', bdysurftot, '(m2)' 
    170          IF(lwp) WRITE(numout,*)'          correction velocity zubtpecor =', zubtpecor , '(m/s)' 
    171          IF(lwp) WRITE(numout,*)'          cumulated transport ztranst   =', ztranst   , '(m3/s)' 
     165         IF(lwp) THEN 
     166            WRITE(numout,*) 
     167            WRITE(numout,*)'bdy_vol : time step :', kt 
     168            WRITE(numout,*)'~~~~~~~ ' 
     169            WRITE(numout,*)'          cumulate flux EMP             =', z_cflxemp  , ' (m3/s)' 
     170            WRITE(numout,*)'          total lateral surface of OBC  =', bdysurftot, '(m2)' 
     171            WRITE(numout,*)'          correction velocity zubtpecor =', zubtpecor , '(m/s)' 
     172            WRITE(numout,*)'          cumulated transport ztranst   =', ztranst   , '(m3/s)' 
     173            IF(lflush) CALL FLUSH(numout) 
     174         ENDIF 
    172175      END IF  
    173176      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/C1D/c1d.F90

    r10968 r10986  
    6767         WRITE(numout,*) '      column longitude                rn_lon1d     = ', rn_lon1d 
    6868         WRITE(numout,*) '      column localization in a grid   ln_c1d_locpt = ', ln_c1d_locpt 
     69         IF(lflush) CALL FLUSH(numout) 
    6970      ENDIF 
    7071      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/C1D/dtauvd.F90

    r10968 r10986  
    8484            WRITE(numout,*) '   U & V current data not used' 
    8585         ENDIF 
     86         IF(lflush) CALL FLUSH(numout) 
    8687      ENDIF 
    8788      !                             ! no initialization when restarting 
     
    155156            WRITE(numout,*) 
    156157            WRITE(numout,*) 'dta_uvd: interpolate U & V current data onto the s- or mixed s-z-coordinate mesh' 
     158            IF(lflush) CALL FLUSH(numout) 
    157159         ENDIF 
    158160         ! 
     
    210212      IF( .NOT. ln_uvd_dyndmp    ) THEN   !==   deallocate U & V current structure   ==!  
    211213         !                                !==   (data used only for initialization)  ==! 
    212          IF(lwp) WRITE(numout,*) 'dta_uvd: deallocate U & V current arrays as they are only used to initialize the run' 
     214         IF(lwp) THEN 
     215            WRITE(numout,*) 'dta_uvd: deallocate U & V current arrays as they are only used to initialize the run' 
     216            IF(lflush) CALL FLUSH(numout) 
     217         ENDIF 
    213218                                   DEALLOCATE( sf_uvd(1)%fnow )     ! U current arrays in the structure 
    214219         IF( sf_uvd(1)%ln_tint )   DEALLOCATE( sf_uvd(1)%fdta ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/C1D/dyncor_c1d.F90

    r10888 r10986  
    4848      !!---------------------------------------------------------------------- 
    4949 
    50       IF(lwp) WRITE(numout,*) 
    51       IF(lwp) WRITE(numout,*) 'cor_c1d : Coriolis factor at T-point' 
    52       IF(lwp) WRITE(numout,*) '~~~~~~~' 
     50      IF(lwp) THEN 
     51         WRITE(numout,*) 
     52         WRITE(numout,*) 'cor_c1d : Coriolis factor at T-point' 
     53         WRITE(numout,*) '~~~~~~~' 
     54         IF(lflush) CALL FLUSH(numout) 
     55      ENDIF 
    5356 
    5457      ! 
     
    6871      !!---------------------------------------------------------------------- 
    6972      ! 
    70       IF( kt == nit000 ) THEN 
    71          IF(lwp) WRITE(numout,*) 
    72          IF(lwp) WRITE(numout,*) 'dyn_cor_c1d : total vorticity trend in 1D' 
    73          IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 
     73      IF( kt == nit000 .AND. lwp) THEN 
     74         WRITE(numout,*) 
     75         WRITE(numout,*) 'dyn_cor_c1d : total vorticity trend in 1D' 
     76         WRITE(numout,*) '~~~~~~~~~~' 
     77         IF(lflush) CALL FLUSH(numout) 
    7478      ENDIF 
    7579      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/C1D/dyndmp.F90

    r10968 r10986  
    9898         WRITE(numout,*) '      Damping file name               cn_resto  = ', cn_resto 
    9999         WRITE(numout,*) 
     100         IF(lflush) CALL FLUSH(numout) 
    100101      ENDIF 
    101102      ! 
     
    112113            CALL ctl_stop(ctmp1) 
    113114         END SELECT 
     115         IF(lwp .AND. lflush) CALL FLUSH(numout) 
    114116         ! 
    115117         IF( .NOT. ln_uvd_dyndmp ) THEN      ! force the initialization of U & V current data for damping 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/CRS/crsdom.F90

    r10888 r10986  
    20852085         WRITE(numout,*) ' njmpp_full  = ', njmpp_full 
    20862086         WRITE(numout,*) 
     2087         IF(lflush) CALL FLUSH(numout) 
    20872088      ENDIF 
    20882089       
     
    21472148 
    21482149         CASE ( 2 )  
    2149             WRITE(numout,*)  'crs_init, jperio=2 not supported'  
     2150            IF(lwp) THEN 
     2151               WRITE(numout,*)  'crs_init, jperio=2 not supported'  
     2152               IF(lflush) CALL FLUSH(numout) 
     2153            ENDIF 
    21502154         
    21512155         CASE ( 5, 6 )    ! F-pivot at North Fold 
     
    21782182 
    21792183         CASE DEFAULT 
    2180             WRITE(numout,*) 'crs_init. Only jperio = 0, 1, 3, 4, 5, 6 supported'  
     2184            IF(lwp) THEN 
     2185               WRITE(numout,*) 'crs_init. Only jperio = 0, 1, 3, 4, 5, 6 supported'  
     2186               IF(lflush) CALL FLUSH(numout) 
     2187            ENDIF 
    21812188  
    21822189         END SELECT 
    21832190 
    21842191      CASE (1 ) 
    2185          WRITE(numout,*) 'crs_init.  Equator-centered bins option not yet available'  
     2192         IF(lwp) THEN 
     2193            WRITE(numout,*) 'crs_init.  Equator-centered bins option not yet available'  
     2194            IF(lflush) CALL FLUSH(numout) 
     2195         ENDIF 
    21862196 
    21872197      END SELECT 
     
    22502260 
    22512261      ! 
    2252       IF(lwp) WRITE(numout,*) 
    2253       IF(lwp) WRITE(numout,*) '    crsini : mbkt is ocean bottom k-index of T-, U-, V- and W-levels ' 
    2254       IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~~~' 
     2262      IF(lwp) THEN 
     2263         WRITE(numout,*) 
     2264         WRITE(numout,*) '    crsini : mbkt is ocean bottom k-index of T-, U-, V- and W-levels ' 
     2265         WRITE(numout,*) '    ~~~~~~~~~~~~~' 
     2266         IF(lflush) CALL FLUSH(numout) 
     2267      ENDIF 
    22552268      ! 
    22562269      mbkt_crs(:,:) = MAX( mbathy_crs(:,:) , 1 )    ! bottom k-index of T-level (=1 over land) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/CRS/crsdomwri.F90

    r10888 r10986  
    5858      ! 
    5959      ! 
    60       IF(lwp) WRITE(numout,*) 
    61       IF(lwp) WRITE(numout,*) 'crs_dom_wri : create NetCDF mesh and mask file' 
    62       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     60      IF(lwp) THEN 
     61         WRITE(numout,*) 
     62         WRITE(numout,*) 'crs_dom_wri : create NetCDF mesh and mask file' 
     63         WRITE(numout,*) '~~~~~~~~~~~' 
     64         IF(lflush) CALL FLUSH(numout) 
     65      ENDIF 
    6366       
    6467      clnam = 'mesh_mask_crs'  ! filename (mesh and mask informations) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/CRS/crsini.F90

    r10968 r10986  
    9999        WRITE(numout,*) '      type of Kz coarsening (0,1,2)         nn_crs_kz  = ', nn_crs_kz 
    100100        WRITE(numout,*) '      wn coarsened or computed using hdivn  ln_crs_wn  = ', ln_crs_wn 
     101        IF(lflush) CALL FLUSH(numout) 
    101102     ENDIF 
    102103               
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/dia25h.F90

    r10968 r10986  
    6767         WRITE(numout,*) '   Namelist nam_dia25h : set 25h outputs ' 
    6868         WRITE(numout,*) '      Switch for 25h diagnostics (T) or not (F)  ln_dia25h  = ', ln_dia25h 
     69         IF(lflush) CALL FLUSH(numout) 
    6970      ENDIF 
    7071      IF( .NOT. ln_dia25h )   RETURN 
     
    159160              WRITE(numout,*) 'dia_wri_tide : Summing instantaneous hourly diagnostics at timestep ',kt 
    160161              WRITE(numout,*) '~~~~~~~~~~~~ ' 
     162              IF(lflush) CALL FLUSH(numout) 
    161163         ENDIF 
    162164 
     
    180182         IF (lwp) THEN 
    181183            WRITE(numout,*) 'dia_tide : Summed the following number of hourly values so far',cnt_25h 
     184            IF(lflush) CALL FLUSH(numout) 
    182185         ENDIF 
    183186         ! 
     
    190193            WRITE(numout,*) 'dia_wri_tide : Writing 25 hour mean tide diagnostics at timestep', kt 
    191194            WRITE(numout,*) '~~~~~~~~~~~~ ' 
     195            IF(lflush) CALL FLUSH(numout) 
    192196         ENDIF 
    193197         ! 
     
    208212         ENDIF 
    209213         ! 
    210          IF(lwp)  WRITE(numout,*) 'dia_wri_tide : Mean calculated by dividing 25 hour sums and writing output' 
     214         IF(lwp)  THEN 
     215            WRITE(numout,*) 'dia_wri_tide : Mean calculated by dividing 25 hour sums and writing output' 
     216            IF(lflush) CALL FLUSH(numout) 
     217         ENDIF 
    211218         zmdi=1.e+20 !missing data indicator for masking 
    212219         ! write tracers (instantaneous) 
     
    261268         ENDIF 
    262269         cnt_25h = 1 
    263          IF(lwp)  WRITE(numout,*) 'dia_wri_tide :   & 
     270         IF(lwp)  THEN 
     271            WRITE(numout,*) 'dia_wri_tide :   & 
    264272            &    After 25hr mean write, reset sum to current value and cnt_25h to one for overlapping average', cnt_25h 
     273            IF(lflush) CALL FLUSH(numout) 
     274         ENDIF 
    265275      ENDIF !  cnt_25h .EQ. 25 .AND.  MOD( kt, i_steps * 24) == 0 .AND. kt .NE. nn_it000 
    266276      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diacfl.F90

    r10888 r10986  
    142142         WRITE(numout,*) '   Max Cv = ', rCv_max, ' at (i,j,k) = (',nCv_loc(1),nCv_loc(2),nCv_loc(3),') => dt/C = ', z2dt/rCv_max 
    143143         WRITE(numout,*) '   Max Cw = ', rCw_max, ' at (i,j,k) = (',nCw_loc(1),nCw_loc(2),nCw_loc(3),') => dt/C = ', z2dt/rCw_max 
     144         IF(lflush) CALL FLUSH(numout) 
    144145      ENDIF 
    145146      ! 
     
    166167         WRITE(numcfl,*) 'Timestep  Direction  Max C     i    j    k' 
    167168         WRITE(numcfl,*) '******************************************' 
     169         IF(lflush) CALL FLUSH(numout) 
    168170      ENDIF 
    169171      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diadct.F90

    r10968 r10986  
    158158        IF(nn_dct .GE. nn_dctwri .AND. MOD(nn_dct,nn_dctwri) .NE. 0)  & 
    159159          &  CALL ctl_stop( 'diadct: nn_dct should be smaller and a multiple of nn_dctwri' ) 
    160  
     160        IF(lflush) CALL FLUSH(numout) 
    161161     ENDIF 
    162162 
     
    222222         WRITE(numout,*) "~~~~~~~~~~~~~~~~~~~~~~~~~" 
    223223         WRITE(numout,*) "nb_sec = ",nb_sec 
     224         IF(lflush) CALL FLUSH(numout) 
    224225     ENDIF 
    225226 
     
    241242        IF( MOD(kt,nn_dctwri)==0 )THEN 
    242243 
    243            IF( kt==nit000+nn_dctwri-1 )WRITE(numout,*)"      diadct: average transports and write at kt = ",kt          
     244           IF( kt==nit000+nn_dctwri-1 ) THEN 
     245              WRITE(numout,*)"      diadct: average transports and write at kt = ",kt          
     246              IF(lflush) CALL FLUSH(numout) 
     247           ENDIF 
    244248   
    245249           !! divide arrays by nn_dctwri/nn_dct to obtain average  
     
    321325     DO jsec=1,nb_sec_max      !loop on the nb_sec sections 
    322326 
    323         IF (  jsec==nn_secdebug .OR. nn_secdebug==-1  ) & 
    324            & WRITE(numout,*)'debuging for section number: ',jsec  
     327        IF (  jsec==nn_secdebug .OR. nn_secdebug==-1  ) THEN 
     328           WRITE(numout,*)'debuging for section number: ',jsec  
     329           IF(lflush) CALL FLUSH(numout) 
     330        ENDIF 
    325331 
    326332        !initialization 
     
    341347        IF( jsec .NE. isec )  CALL ctl_stop( cltmp ) 
    342348 
    343         IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )WRITE(numout,*)"isec ",isec  
     349        IF( jsec==nn_secdebug .OR. nn_secdebug==-1 ) THEN 
     350           WRITE(numout,*)"isec ",isec  
     351           IF(lflush) CALL FLUSH(numout) 
     352        ENDIF 
    344353 
    345354        READ(numdct_in)secs(jsec)%name 
     
    376385            WRITE(numout,clformat)"      Temperature classes :             ",secs(jsec)%ztem 
    377386            WRITE(numout,clformat)"      Depth classes :                   ",secs(jsec)%zlay 
     387            IF(lflush) CALL FLUSH(numout) 
    378388        ENDIF                
    379389 
     
    398408                 WRITE(numout,*)'        # I J ',jpt,coordtemp(jpt),directemp(jpt) 
    399409              ENDDO                   
     410              IF(lflush) CALL FLUSH(numout) 
    400411           ENDIF 
    401412  
     
    434445                 WRITE(numout,*)'         # I J : ',iiglo,ijglo 
    435446              ENDDO 
     447              IF(lflush) CALL FLUSH(numout) 
    436448           ENDIF 
    437449 
     
    468480                 ijglo = secs(jsec)%listPoint(jpt)%J + njmpp - 1 
    469481                 WRITE(numout,*)'         # I J : ',iiglo,ijglo 
    470                  CALL FLUSH(numout) 
     482                 IF(lflush) CALL FLUSH(numout) 
    471483              ENDDO 
    472484           ENDIF 
    473485 
    474486        ELSE  ! iptglo = 0 
    475            IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )& 
     487           IF( jsec==nn_secdebug .OR. nn_secdebug==-1 ) THEN 
    476488              WRITE(numout,*)'   No points for this section.' 
     489              IF(lflush) CALL FLUSH(numout) 
     490           ENDIF 
    477491        ENDIF 
    478492 
     
    510524     !---------------------------------------------------------------------------- 
    511525     ! 
    512      IF( ld_debug )WRITE(numout,*)'      -------------------------' 
    513      IF( ld_debug )WRITE(numout,*)'      removepoints in listpoint' 
     526     IF( ld_debug ) THEN 
     527        WRITE(numout,*)'      -------------------------' 
     528        WRITE(numout,*)'      removepoints in listpoint' 
     529        IF(lflush) CALL FLUSH(numout) 
     530     ENDIF 
    514531 
    515532     !iextr=extremity of list_point that we verify 
     
    529546        WRITE(numout,*)'      ', cdind,' ',cdextr,' ',itest 
    530547        WRITE(numout,*)'      Actual number of points: ',sec%nb_point 
     548        IF(lflush) CALL FLUSH(numout) 
    531549     ENDIF 
    532550 
     
    556574        WRITE(numout,*)'      Number of points after removepoints :',sec%nb_point 
    557575        WRITE(numout,*)'      sec%direction after removepoints :',sec%direction(1:sec%nb_point) 
     576        IF(lflush) CALL FLUSH(numout) 
    558577     ENDIF 
    559578      ! 
     
    592611      !!-------------------------------------------------------- 
    593612      ! 
    594       IF( ld_debug )WRITE(numout,*)'      Compute transport' 
     613      IF( ld_debug ) THEN 
     614         WRITE(numout,*)'      Compute transport' 
     615         IF(lflush) CALL FLUSH(numout) 
     616      ENDIF 
    595617 
    596618      !---------------------------! 
     
    626648         IF( sec%slopeSection .GE. 9999. )     isgnv =  1 
    627649 
    628          IF( ld_debug )write(numout,*)"sec%slopeSection isgnu isgnv ",sec%slopeSection,isgnu,isgnv 
     650         IF( ld_debug ) THEN 
     651            write(numout,*)"sec%slopeSection isgnu isgnv ",sec%slopeSection,isgnu,isgnv 
     652            IF(lflush) CALL FLUSH(numout) 
     653         ENDIF 
    629654 
    630655         !--------------------------------------! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diaharm.F90

    r10968 r10986  
    8181         WRITE(numout,*) 'dia_harm_init: Tidal harmonic analysis initialization' 
    8282         WRITE(numout,*) '~~~~~~~ ' 
     83         IF(lflush) CALL FLUSH(numout) 
    8384      ENDIF 
    8485      ! 
     
    99100         WRITE(numout,*) 'Last  time step used for analysis:  nitend_han= ', nitend_han 
    100101         WRITE(numout,*) 'Time step frequency for harmonic analysis:  nstep_han= ', nstep_han 
     102         IF(lflush) CALL FLUSH(numout) 
    101103      ENDIF 
    102104 
     
    123125         WRITE(numout,*) '        Namelist nam_diaharm' 
    124126         WRITE(numout,*) '        nb_ana    = ', nb_ana 
    125          CALL flush(numout) 
     127         IF(lflush) CALL FLUSH(numout) 
    126128      ENDIF 
    127129      ! 
     
    148150      CALL tide_harmo( ana_freq, vt, ut, ft, name, nb_ana ) 
    149151 
    150       IF(lwp) WRITE(numout,*) 'Analysed frequency  : ',nb_ana ,'Frequency ' 
    151  
    152       DO jh = 1, nb_ana 
    153         IF(lwp) WRITE(numout,*) '                    : ',tname(jh),' ',ana_freq(jh) 
    154       END DO 
     152      IF(lwp) THEN 
     153         WRITE(numout,*) 'Analysed frequency  : ',nb_ana ,'Frequency ' 
     154         DO jh = 1, nb_ana 
     155            WRITE(numout,*) '                    : ',tname(jh),' ',ana_freq(jh) 
     156         END DO 
     157         IF(lflush) CALL FLUSH(numout) 
     158      ENDIF 
    155159 
    156160      ! Initialize temporary arrays: 
     
    228232      !!-------------------------------------------------------------------- 
    229233      ! 
    230       IF(lwp) WRITE(numout,*) 
    231       IF(lwp) WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 
    232       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     234      IF(lwp) THEN 
     235         WRITE(numout,*) 
     236         WRITE(numout,*) 'anharmo_end: kt=nitend_han: Perform harmonic analysis' 
     237         WRITE(numout,*) '~~~~~~~~~~~~' 
     238         IF(lflush) CALL FLUSH(numout) 
     239      ENDIF 
    233240 
    234241      ztime_ini = nit000_han*rdt                 ! Initial time in seconds at the beginning of analysis 
     
    382389      !!---------------------------------------------------------------------- 
    383390 
    384       IF(lwp) WRITE(numout,*) '  ' 
    385       IF(lwp) WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 
    386       IF(lwp) WRITE(numout,*) '  ' 
     391      IF(lwp) THEN 
     392         WRITE(numout,*) '  ' 
     393         WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 
     394         WRITE(numout,*) '  ' 
     395         IF(lflush) CALL FLUSH(numout) 
     396      ENDIF 
    387397 
    388398      ! A) Elevation 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diahsb.F90

    r10968 r10986  
    211211         CALL iom_put( 'bgvole3t' , zdiff_v2 * 1.e-9    )              ! volume e3t drift      (km3)   
    212212         ! 
    213          IF( kt == nitend .AND. lwp) THEN 
     213         IF( kt == nitend .AND. lwp .AND. nprint > 0) THEN 
    214214            WRITE(numout,*) 
    215215            WRITE(numout,*) 'dia_hsb : last time step hsb diagnostics: at it= ', kt,' date= ', ndastp 
     
    219219            WRITE(numout,*) '   volume ssh  drift = ', zdiff_v1 * 1.e-9   , ' km^3' 
    220220            WRITE(numout,*) '   volume e3t  drift = ', zdiff_v2 * 1.e-9   , ' km^3' 
     221            IF(lflush) CALL FLUSH(numout) 
    221222         ENDIF 
    222223         ! 
     
    257258         IF( ln_rstart ) THEN                   !* Read the restart file 
    258259            ! 
    259             IF(lwp) WRITE(numout,*) 
    260             IF(lwp) WRITE(numout,*) '   dia_hsb_rst : read hsb restart at it= ', kt,' date= ', ndastp 
    261             IF(lwp) WRITE(numout,*) 
     260            IF(lwp) THEN 
     261               WRITE(numout,*) 
     262               WRITE(numout,*) '   dia_hsb_rst : read hsb restart at it= ', kt,' date= ', ndastp 
     263               WRITE(numout,*) 
     264               IF(lflush) CALL FLUSH(numout) 
     265            ENDIF 
    262266            CALL iom_get( numror, 'frc_v', frc_v, ldxios = lrxios ) 
    263267            CALL iom_get( numror, 'frc_t', frc_t, ldxios = lrxios ) 
     
    277281            ENDIF 
    278282         ELSE 
    279             IF(lwp) WRITE(numout,*) 
    280             IF(lwp) WRITE(numout,*) '   dia_hsb_rst : initialise hsb at initial state ' 
    281             IF(lwp) WRITE(numout,*) 
     283            IF(lwp) THEN 
     284               WRITE(numout,*) 
     285               WRITE(numout,*) '   dia_hsb_rst : initialise hsb at initial state ' 
     286               WRITE(numout,*) 
     287            ENDIF 
    282288            surf_ini(:,:) = e1e2t(:,:) * tmask_i(:,:)         ! initial ocean surface 
    283289            ssh_ini(:,:) = sshn(:,:)                          ! initial ssh 
     
    310316      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    311317         !                                   ! ------------------- 
    312          IF(lwp .AND. nprint > 0) WRITE(numout,*) 
    313          IF(lwp .AND. nprint > 0) WRITE(numout,*) '   dia_hsb_rst : write restart at it= ', kt,' date= ', ndastp 
    314          IF(lwp .AND. nprint > 0) WRITE(numout,*) 
     318         IF(lwp .AND. nprint > 0) THEN 
     319            WRITE(numout,*) 
     320            WRITE(numout,*) '   dia_hsb_rst : write restart at it= ', kt,' date= ', ndastp 
     321            WRITE(numout,*) 
     322            IF(lflush) CALL FLUSH(numout) 
     323         ENDIF 
    315324         ! 
    316325         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
     
    359368         WRITE(numout,*) 'dia_hsb_init : heat and salt budgets diagnostics' 
    360369         WRITE(numout,*) '~~~~~~~~~~~~ ' 
     370         IF(lflush) CALL FLUSH(numout) 
    361371      ENDIF 
    362372      REWIND( numnam_ref )              ! Namelist namhsb in reference namelist 
     
    371381         WRITE(numout,*) '   Namelist  namhsb :' 
    372382         WRITE(numout,*) '      check the heat and salt budgets (T) or not (F)       ln_diahsb = ', ln_diahsb 
     383         IF(lflush) CALL FLUSH(numout) 
    373384      ENDIF 
    374385      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diahth.F90

    r10888 r10986  
    127127         END IF 
    128128 
    129          IF(lwp) WRITE(numout,*) 
    130          IF(lwp) WRITE(numout,*) 'dia_hth : diagnostics of the thermocline depth' 
    131          IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
    132          IF(lwp) WRITE(numout,*) 
     129         IF(lwp) THEN 
     130            WRITE(numout,*) 
     131            WRITE(numout,*) 'dia_hth : diagnostics of the thermocline depth' 
     132            WRITE(numout,*) '~~~~~~~ ' 
     133            WRITE(numout,*) 
     134            IF(lflush) CALL FLUSH(numout) 
     135         ENDIF 
    133136      ENDIF 
    134137 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diaptr.F90

    r10968 r10986  
    407407         WRITE(numout,*) '      Poleward heat & salt transport (T) or not (F)      ln_diaptr  = ', ln_diaptr 
    408408         WRITE(numout,*) '      Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins      ln_subbas  = ', ln_subbas 
     409         IF(lflush) CALL FLUSH(numout) 
    409410      ENDIF 
    410411 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diatmb.F90

    r10968 r10986  
    5656         WRITE(numout,*) '   Namelist nam_diatmb : set tmb outputs ' 
    5757         WRITE(numout,*) '      Switch for TMB diagnostics (T) or not (F)  ln_diatmb  = ', ln_diatmb 
     58         IF(lflush) CALL FLUSH(numout) 
    5859      ENDIF 
    5960      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/diawri.F90

    r10888 r10986  
    430430      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
    431431      ! 
    432       LOGICAL ::   ll_print = .FALSE.                        ! =T print and flush numout 
    433432      CHARACTER (len=40) ::   clhstnam, clop, clmx           ! local names 
    434433      INTEGER  ::   inum = 11                                ! temporary logical unit 
     
    488487         CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
    489488         zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    490          IF(lwp)WRITE(numout,*) 
    491          IF(lwp)WRITE(numout,*) 'Date 0 used :', nit000, ' YEAR ', nyear,   & 
     489         IF(lwp) THEN 
     490            WRITE(numout,*) 
     491            WRITE(numout,*) 'Date 0 used :', nit000, ' YEAR ', nyear,   & 
    492492            &                    ' MONTH ', nmonth, ' DAY ', nday, 'Julian day : ', zjulian 
    493          IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma,   & 
     493            WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma,   & 
    494494                                 ' limit storage in depth = ', ipk 
    495  
    496          ! WRITE root name in date.file for use by postpro 
    497          IF(lwp) THEN 
     495            IF(lflush) CALL FLUSH(numout) 
     496            ! WRITE root name in date.file for use by postpro 
    498497            CALL dia_nam( clhstnam, nwrite,' ' ) 
    499498            CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     
    505504 
    506505         CALL dia_nam( clhstnam, nwrite, 'grid_T' ) 
    507          IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename 
     506         IF(lwp) THEN 
     507            WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename 
     508            IF(lflush) CALL FLUSH(numout) 
     509         ENDIF 
    508510         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    509511            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
     
    541543 
    542544         CALL dia_nam( clhstnam, nwrite, 'grid_U' ) 
    543          IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename 
     545         IF(lwp) THEN 
     546            WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename 
     547            IF(lflush) CALL FLUSH(numout) 
     548         ENDIF 
    544549         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    545550            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
     
    554559 
    555560         CALL dia_nam( clhstnam, nwrite, 'grid_V' )                   ! filename 
    556          IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 
     561         IF(lwp) THEN 
     562            WRITE(numout,*) " Name of NETCDF file ", clhstnam 
     563            IF(lflush) CALL FLUSH(numout) 
     564         ENDIF 
    557565         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    558566            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
     
    567575 
    568576         CALL dia_nam( clhstnam, nwrite, 'grid_W' )                   ! filename 
    569          IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 
     577         IF(lwp) THEN 
     578            WRITE(numout,*) " Name of NETCDF file ", clhstnam 
     579            IF(lflush) CALL FLUSH(numout) 
     580         ENDIF 
    570581         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    571582            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
     
    737748         CALL histend( nid_W, snc4chunks=snc4set ) 
    738749 
    739          IF(lwp) WRITE(numout,*) 
    740          IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization' 
    741          IF(ll_print) CALL FLUSH(numout ) 
     750         IF(lwp) THEN 
     751            WRITE(numout,*) 
     752            WRITE(numout,*) 'End of NetCDF Initialization' 
     753            IF(lflush) CALL FLUSH(numout) 
     754         ENDIF 
    742755 
    743756      ENDIF 
     
    753766         WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 
    754767         WRITE(numout,*) '~~~~~~ ' 
     768         IF(lflush) CALL FLUSH(numout) 
    755769      ENDIF 
    756770 
     
    887901      !!---------------------------------------------------------------------- 
    888902      !  
    889       IF(lwp) WRITE(numout,*) 
    890       IF(lwp) WRITE(numout,*) 'dia_wri_state : single instantaneous ocean state' 
    891       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~   and forcing fields file created ' 
    892       IF(lwp) WRITE(numout,*) '                and named :', cdfile_name, '...nc' 
     903      IF(lwp) THEN 
     904         WRITE(numout,*) 
     905         WRITE(numout,*) 'dia_wri_state : single instantaneous ocean state' 
     906         WRITE(numout,*) '~~~~~~~~~~~~~   and forcing fields file created ' 
     907         WRITE(numout,*) '                and named :', cdfile_name, '...nc' 
     908         IF(lflush) CALL FLUSH(numout) 
     909      ENDIF 
    893910 
    894911#if defined key_si3 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/closea.F90

    r10888 r10986  
    9797      !!---------------------------------------------------------------------- 
    9898      ! 
    99       IF(lwp) WRITE(numout,*) 
    100       IF(lwp) WRITE(numout,*)'dom_clo : read in masks to define closed seas ' 
    101       IF(lwp) WRITE(numout,*)'~~~~~~~' 
     99      IF(lwp) THEN 
     100         WRITE(numout,*) 
     101         WRITE(numout,*)'dom_clo : read in masks to define closed seas ' 
     102         WRITE(numout,*)'~~~~~~~' 
     103         IF(lflush) CALL FLUSH(numout) 
     104      ENDIF 
    102105      ! 
    103106      ! read the closed seas masks (if they exist) from domain_cfg file (if it exists) 
     
    120123            CALL mpp_max('closea', jncs) 
    121124            IF( jncs > 0 ) THEN 
    122                IF( lwp ) WRITE(numout,*) 'Number of closed seas : ',jncs 
     125               IF( lwp ) THEN 
     126                  WRITE(numout,*) 'Number of closed seas : ',jncs 
     127                  IF(lflush) CALL FLUSH(numout) 
     128               ENDIF 
    123129            ELSE 
    124130               CALL ctl_stop( 'Problem with closea_mask field in domain_cfg file. Has no values > 0 so no closed seas defined.') 
    125131            ENDIF 
    126132         ELSE  
    127             IF( lwp ) WRITE(numout,*) 
    128             IF( lwp ) WRITE(numout,*) '   ==>>>   closea_mask field not found in domain_cfg file.' 
    129             IF( lwp ) WRITE(numout,*) '           No closed seas defined.' 
    130             IF( lwp ) WRITE(numout,*) 
     133            IF( lwp ) THEN 
     134               WRITE(numout,*) 
     135               WRITE(numout,*) '   ==>>>   closea_mask field not found in domain_cfg file.' 
     136               WRITE(numout,*) '           No closed seas defined.' 
     137               WRITE(numout,*) 
     138               IF(lflush) CALL FLUSH(numout) 
     139            ENDIF 
    131140            l_sbc_clo = .false. 
    132141            jncs = 0  
     
    148157               CALL mpp_max('closea', jncsr) 
    149158               IF( jncsr > 0 ) THEN 
    150                   IF( lwp ) WRITE(numout,*) 'Number of closed seas rnf mappings : ',jncsr 
     159                  IF( lwp ) THEN 
     160                     WRITE(numout,*) 'Number of closed seas rnf mappings : ',jncsr 
     161                     IF(lflush) CALL FLUSH(numout) 
     162                  ENDIF 
    151163               ELSE 
    152164                  CALL ctl_stop( 'Problem with closea_mask_rnf field in domain_cfg file. Has no values > 0 so no closed seas rnf mappings defined.') 
    153165               ENDIF 
    154166            ELSE  
    155                IF( lwp ) WRITE(numout,*) 'closea_mask_rnf field not found in domain_cfg file. No closed seas rnf mappings defined.' 
     167               IF( lwp ) THEN 
     168                  WRITE(numout,*) 'closea_mask_rnf field not found in domain_cfg file. No closed seas rnf mappings defined.' 
     169                  IF(lflush) CALL FLUSH(numout) 
     170               ENDIF 
    156171               jncsr = 0 
    157172            ENDIF 
     
    168183               CALL mpp_max('closea', jncse) 
    169184               IF( jncse > 0 ) THEN  
    170                   IF( lwp ) WRITE(numout,*) 'Number of closed seas empmr mappings : ',jncse 
     185                  IF( lwp ) THEN 
     186                     WRITE(numout,*) 'Number of closed seas empmr mappings : ',jncse 
     187                     IF(lflush) CALL FLUSH(numout) 
     188                  ENDIF 
    171189               ELSE 
    172190                  CALL ctl_stop( 'Problem with closea_mask_empmr field in domain_cfg file. Has no values > 0 so no closed seas empmr mappings defined.') 
    173191               ENDIF 
    174192            ELSE  
    175                IF( lwp ) WRITE(numout,*) 'closea_mask_empmr field not found in domain_cfg file. No closed seas empmr mappings defined.' 
     193               IF( lwp ) THEN 
     194                  WRITE(numout,*) 'closea_mask_empmr field not found in domain_cfg file. No closed seas empmr mappings defined.' 
     195                  IF(lflush) CALL FLUSH(numout) 
     196               ENDIF 
    176197               jncse = 0 
    177198            ENDIF 
     
    182203         ! 
    183204      ELSE ! ln_read_cfg = .false. so no domain_cfg file 
    184          IF( lwp ) WRITE(numout,*) 'No domain_cfg file so no closed seas defined.' 
     205         IF( lwp ) THEN 
     206            WRITE(numout,*) 'No domain_cfg file so no closed seas defined.' 
     207            IF(lflush) CALL FLUSH(numout) 
     208         ENDIF 
    185209         l_sbc_clo = .false. 
    186210         l_clo_rnf = .false. 
     
    219243      IF( kt == nit000 ) THEN                             !  Initialisation  ! 
    220244         !                                                !------------------! 
    221          IF(lwp) WRITE(numout,*) 
    222          IF(lwp) WRITE(numout,*)'sbc_clo : closed seas ' 
    223          IF(lwp) WRITE(numout,*)'~~~~~~~' 
     245         IF(lwp) THEN 
     246            WRITE(numout,*) 
     247            WRITE(numout,*)'sbc_clo : closed seas ' 
     248            WRITE(numout,*)'~~~~~~~' 
     249            IF(lflush) CALL FLUSH(numout) 
     250         ENDIF 
    224251 
    225252         ALLOCATE( surf(jncs+1) , STAT=ierr ) 
     
    267294         ENDIF 
    268295         ! 
    269          IF(lwp) WRITE(numout,*)'     Closed sea surface areas (km2)' 
    270          DO jc = 1, jncs 
    271             IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jc, surf(jc) * 1.0e-6 
    272          END DO 
    273          IF(lwp) WRITE(numout,FMT='(A,ES12.2)') 'Global surface area excluding closed seas (km2): ', surf(jncs+1) * 1.0e-6 
    274          ! 
    275          IF(jncsr > 0) THEN 
    276             IF(lwp) WRITE(numout,*)'     Closed sea target rnf surface areas (km2)' 
    277             DO jcr = 1, jncsr 
    278                IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jcr, surfr(jcr) * 1.0e-6 
     296         IF(lwp) THEN 
     297            WRITE(numout,*)'     Closed sea surface areas (km2)' 
     298            DO jc = 1, jncs 
     299               WRITE(numout,FMT='(1I3,5X,ES12.2)') jc, surf(jc) * 1.0e-6 
    279300            END DO 
    280          ENDIF 
    281          ! 
    282          IF(jncse > 0) THEN 
    283             IF(lwp) WRITE(numout,*)'     Closed sea target empmr surface areas (km2)' 
    284             DO jce = 1, jncse 
    285                IF(lwp) WRITE(numout,FMT='(1I3,5X,ES12.2)') jce, surfe(jce) * 1.0e-6 
    286             END DO 
     301            WRITE(numout,FMT='(A,ES12.2)') 'Global surface area excluding closed seas (km2): ', surf(jncs+1) * 1.0e-6 
     302         ! 
     303            IF(jncsr > 0) THEN 
     304               WRITE(numout,*)'     Closed sea target rnf surface areas (km2)' 
     305               DO jcr = 1, jncsr 
     306                  WRITE(numout,FMT='(1I3,5X,ES12.2)') jcr, surfr(jcr) * 1.0e-6 
     307               END DO 
     308            ENDIF 
     309         ! 
     310            IF(jncse > 0) THEN 
     311               WRITE(numout,*)'     Closed sea target empmr surface areas (km2)' 
     312               DO jce = 1, jncse 
     313                  WRITE(numout,FMT='(1I3,5X,ES12.2)') jce, surfe(jce) * 1.0e-6 
     314               END DO 
     315            ENDIF 
     316            IF(lflush) CALL FLUSH(numout) 
    287317         ENDIF 
    288318      ENDIF 
     
    453483         WRITE(numout,*) 'clo_bat : suppression of closed seas' 
    454484         WRITE(numout,*) '~~~~~~~' 
     485         IF(lflush) CALL FLUSH(numout) 
    455486      ENDIF 
    456487      ! 
     
    461492         id = iom_varid(inum, 'closea_mask', ldstop = .false.)       
    462493         IF( id > 0 ) THEN 
    463             IF( lwp ) WRITE(numout,*) 'Suppressing closed seas in bathymetry based on closea_mask field,' 
     494            IF( lwp ) THEN 
     495               WRITE(numout,*) 'Suppressing closed seas in bathymetry based on closea_mask field,' 
     496               IF(lflush) CALL FLUSH(numout) 
     497            ENDIF 
    464498            CALL iom_get ( inum, jpdom_data, 'closea_mask', zdata_in ) 
    465499            closea_mask(:,:) = NINT(zdata_in(:,:)) 
     
    469503            ENDWHERE 
    470504         ELSE 
    471             IF( lwp ) WRITE(numout,*) 'No closea_mask field found in domain_cfg file. No suppression of closed seas.' 
     505            IF( lwp ) THEN 
     506               WRITE(numout,*) 'No closea_mask field found in domain_cfg file. No suppression of closed seas.' 
     507               IF(lflush) CALL FLUSH(numout) 
     508            ENDIF 
    472509         ENDIF 
    473510         ! 
     
    475512         ! 
    476513      ELSE 
    477          IF( lwp ) WRITE(numout,*) 'No domain_cfg file => no suppression of closed seas.' 
     514         IF( lwp ) THEN 
     515            WRITE(numout,*) 'No domain_cfg file => no suppression of closed seas.' 
     516            IF(lflush) CALL FLUSH(numout) 
     517         ENDIF 
    478518      ENDIF 
    479519      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/daymod.F90

    r10968 r10986  
    134134 
    135135      ! control print 
    136       IF(lwp) WRITE(numout,'(a,i6,a,i2,a,i2,a,i8,a,i8,a,i8,a,i8)')   & 
     136      IF(lwp) THEN 
     137         WRITE(numout,'(a,i6,a,i2,a,i2,a,i8,a,i8,a,i8,a,i8)')   & 
    137138           &                   ' =======>> 1/2 time step before the start of the run DATE Y/M/D = ',   & 
    138139           &                   nyear, '/', nmonth, '/', nday, '  nsec_day:', nsec_day, '  nsec_week:', nsec_week, '  & 
    139140           &                   nsec_month:', nsec_month , '  nsec_year:' , nsec_year 
     141         IF(lflush) CALL FLUSH(numout) 
     142      ENDIF 
    140143 
    141144      ! Up to now, calendar parameters are related to the end of previous run (nit000-1) 
     
    269272         CALL ymds2ju( nyear, 01, 01, 0.0, fjulstartyear ) 
    270273         ! 
    271          IF(lwp .AND. nprint > 0) WRITE(numout,'(a,i8,a,i4.4,a,i2.2,a,i2.2,a,i3.3)') '======>> time-step =', kt,   & 
    272               &   '      New day, DATE Y/M/D = ', nyear, '/', nmonth, '/', nday, '      nday_year = ', nday_year 
    273          IF(lwp .AND. nprint > 0) WRITE(numout,'(a,i8,a,i7,a,i5)') '         nsec_year = ', nsec_year,   & 
     274         IF(lwp .AND. nprint > 0) THEN 
     275            WRITE(numout,'(a,i8,a,i4.4,a,i2.2,a,i2.2,a,i3.3)') '======>> time-step =', kt,   & 
     276           &   '      New day, DATE Y/M/D = ', nyear, '/', nmonth, '/', nday, '      nday_year = ', nday_year 
     277            WRITE(numout,'(a,i8,a,i7,a,i5)') '         nsec_year = ', nsec_year,   & 
    274278              &   '   nsec_month = ', nsec_month, '   nsec_day = ', nsec_day, '   nsec_week = ', nsec_week 
     279            IF(lflush) CALL FLUSH(numout) 
     280         ENDIF 
    275281      ENDIF 
    276282 
     
    337343               END SELECT 
    338344               WRITE(numout,*) 
     345               IF(lflush) CALL FLUSH(numout) 
    339346            ENDIF 
    340347            ! Control of date 
     
    397404       WRITE(numout,*) '   nn_time0                                         : ',nn_time0 
    398405            WRITE(numout,*) 
     406            IF(lflush) CALL FLUSH(numout) 
    399407         ENDIF 
    400408         ! 
     
    406414               WRITE(numout,*) 'rst_write : write oce restart file  kt =', kt 
    407415               WRITE(numout,*) '~~~~~~~' 
     416               IF(lflush) CALL FLUSH(numout) 
    408417            ENDIF 
    409418         ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domain.F90

    r10968 r10986  
    112112         WRITE(numout,*)     '      Ocean model configuration used:' 
    113113         WRITE(numout,*)     '         cn_cfg = ', TRIM( cn_cfg ), '   nn_cfg = ', nn_cfg 
     114         IF(lflush) CALL FLUSH(numout) 
    114115      ENDIF 
    115116      lwxios = .FALSE. 
     
    129130      IF(cdstr == 'SAS') THEN 
    130131         IF(lrxios) THEN 
    131                IF(lwp) write(numout,*) 'Disable reading restart file using XIOS for SAS' 
     132               IF(lwp) THEN 
     133                  write(numout,*) 'Disable reading restart file using XIOS for SAS' 
     134                  IF(lflush) CALL FLUSH(numout) 
     135               ENDIF 
    132136               lrxios = .FALSE. 
    133137         ENDIF 
     
    203207         WRITE(numout,*) '~~~~~~~~' 
    204208         WRITE(numout,*)  
     209         IF(lflush) CALL FLUSH(numout) 
    205210      ENDIF 
    206211      ! 
     
    268273            WRITE(numout,25)              (mj1(jj),jj = 1,jpjglo) 
    269274         ENDIF 
     275         IF(lflush) CALL FLUSH(numout) 
    270276      ENDIF 
    271277 25   FORMAT( 100(10x,19i4,/) ) 
     
    303309         WRITE(numout,*) 'dom_nam : domain initialization through namelist read' 
    304310         WRITE(numout,*) '~~~~~~~ ' 
     311         IF(lflush) CALL FLUSH(numout) 
    305312      ENDIF 
    306313      ! 
     
    349356            WRITE(numout,*) "      AGRIF: ln_xios_read will be ingored. See setting for parent" 
    350357         ENDIF 
     358         IF(lflush) CALL FLUSH(numout) 
    351359      ENDIF 
    352360 
     
    363371      neuler = nn_euler 
    364372      IF( neuler == 1 .AND. .NOT. ln_rstart ) THEN 
    365          IF(lwp) WRITE(numout,*)   
    366          IF(lwp) WRITE(numout,*)'   ==>>>   Start from rest (ln_rstart=F)' 
    367          IF(lwp) WRITE(numout,*)'           an Euler initial time step is used : nn_euler is forced to 0 '    
     373         IF(lwp) THEN 
     374            WRITE(numout,*)   
     375            WRITE(numout,*)'   ==>>>   Start from rest (ln_rstart=F)' 
     376            WRITE(numout,*)'           an Euler initial time step is used : nn_euler is forced to 0 '    
     377            IF(lflush) CALL FLUSH(numout) 
     378         ENDIF 
    368379         neuler = 0 
    369380      ENDIF 
     
    383394      IF( Agrif_Root() ) THEN 
    384395#endif 
    385       IF(lwp) WRITE(numout,*) 
     396      IF(lwp) THEN 
     397         WRITE(numout,*) 
     398         IF(lflush) CALL FLUSH(numout) 
     399      ENDIF 
    386400      SELECT CASE ( nleapy )        ! Choose calendar for IOIPSL 
    387401      CASE (  1 )  
     
    395409         IF(lwp) WRITE(numout,*) '   ==>>>   The IOIPSL calendar is "360d", i.e. 360 days in a year' 
    396410      END SELECT 
     411      IF(lflush .AND. lwp) CALL FLUSH(numout) 
    397412#if defined key_agrif 
    398413      ENDIF 
     
    416431         WRITE(numout,*) '      asselin time filter parameter           rn_atfp     = ', rn_atfp 
    417432         WRITE(numout,*) '      online coarsening of dynamical fields   ln_crs      = ', ln_crs 
     433         IF(lflush) CALL FLUSH(numout) 
    418434      ENDIF 
    419435      ! 
     
    446462         WRITE(numout,*) '      number of chunks in k-dimension             nn_nchunks_k = ', nn_nchunks_k 
    447463         WRITE(numout,*) '      apply netcdf4/hdf5 chunking & compression   ln_nc4zip    = ', ln_nc4zip 
     464         IF(lflush) CALL FLUSH(numout) 
    448465      ENDIF 
    449466 
     
    506523         WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, ima2(1), ima2(2) 
    507524         WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, imi2(1), imi2(2) 
     525         IF(lflush) CALL FLUSH(numout) 
    508526      ENDIF 
    509527      ! 
     
    597615      !!---------------------------------------------------------------------- 
    598616      ! 
    599       IF(lwp) WRITE(numout,*) 
    600       IF(lwp) WRITE(numout,*) 'cfg_write : create the domain configuration file (', TRIM(cn_domcfg_out),'.nc)' 
    601       IF(lwp) WRITE(numout,*) '~~~~~~~~~' 
     617      IF(lwp) THEN 
     618         WRITE(numout,*) 
     619         WRITE(numout,*) 'cfg_write : create the domain configuration file (', TRIM(cn_domcfg_out),'.nc)' 
     620         WRITE(numout,*) '~~~~~~~~~' 
     621         IF(lflush) CALL FLUSH(numout) 
     622      ENDIF 
    602623      ! 
    603624      !                       ! ============================= ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domhgr.F90

    r10888 r10986  
    8686         WRITE(numout,*) '~~~~~~~   ' 
    8787         WRITE(numout,*) '   namcfg : read (=T) or user defined (=F) configuration    ln_read_cfg  = ', ln_read_cfg 
     88         IF(lflush) CALL FLUSH(numout) 
    8889      ENDIF 
    8990      ! 
    9091      ! 
    9192      IF( ln_read_cfg ) THEN        !==  read in mesh_mask.nc file  ==! 
    92          IF(lwp) WRITE(numout,*) 
    93          IF(lwp) WRITE(numout,*) '   ==>>>   read horizontal mesh in ', TRIM( cn_domcfg ), ' file' 
     93         IF(lwp) THEN 
     94            WRITE(numout,*) 
     95            WRITE(numout,*) '   ==>>>   read horizontal mesh in ', TRIM( cn_domcfg ), ' file' 
     96            IF(lflush) CALL FLUSH(numout) 
     97         ENDIF 
    9498         ! 
    9599         CALL hgr_read   ( glamt , glamu , glamv , glamf ,   &    ! geographic position (required) 
     
    101105         ! 
    102106      ELSE                          !==  User defined configuration  ==!  
    103          IF(lwp) WRITE(numout,*) 
    104          IF(lwp) WRITE(numout,*) '          User defined horizontal mesh (usr_def_hgr)' 
     107         IF(lwp) THEN 
     108            WRITE(numout,*) 
     109            WRITE(numout,*) '          User defined horizontal mesh (usr_def_hgr)' 
     110            IF(lflush) CALL FLUSH(numout) 
     111         ENDIF 
    105112         ! 
    106113         CALL usr_def_hgr( glamt , glamu , glamv , glamf ,   &    ! geographic position (required) 
     
    116123      ! 
    117124      IF( iff == 0 ) THEN                 ! Coriolis parameter has not been defined  
    118          IF(lwp) WRITE(numout,*) '          Coriolis parameter calculated on the sphere from gphif & gphit' 
     125         IF(lwp) THEN 
     126            WRITE(numout,*) '          Coriolis parameter calculated on the sphere from gphif & gphit' 
     127            IF(lflush) CALL FLUSH(numout) 
     128         ENDIF 
    119129         ff_f(:,:) = 2. * omega * SIN( rad * gphif(:,:) )     ! compute it on the sphere at f-point 
    120130         ff_t(:,:) = 2. * omega * SIN( rad * gphit(:,:) )     !    -        -       -    at t-point 
    121131      ELSE 
    122132         IF( ln_read_cfg ) THEN 
    123             IF(lwp) WRITE(numout,*) '          Coriolis parameter have been read in ', TRIM( cn_domcfg ), ' file' 
     133            IF(lwp) THEN 
     134               WRITE(numout,*) '          Coriolis parameter have been read in ', TRIM( cn_domcfg ), ' file' 
     135               IF(lflush) CALL FLUSH(numout) 
     136            ENDIF 
    124137         ELSE 
    125             IF(lwp) WRITE(numout,*) '          Coriolis parameter have been set in usr_def_hgr routine' 
     138            IF(lwp) THEN 
     139               WRITE(numout,*) '          Coriolis parameter have been set in usr_def_hgr routine' 
     140               IF(lflush) CALL FLUSH(numout) 
     141            ENDIF 
    126142         ENDIF 
    127143      ENDIF 
     
    138154      e1e2f (:,:) = e1f(:,:) * e2f(:,:)   ;   r1_e1e2f(:,:) = 1._wp / e1e2f(:,:) 
    139155      IF( ie1e2u_v == 0 ) THEN               ! u- & v-surfaces have not been defined 
    140          IF(lwp) WRITE(numout,*) '          u- & v-surfaces calculated as e1 e2 product' 
     156         IF(lwp .AND. nprint > 1) THEN 
     157            WRITE(numout,*) '          u- & v-surfaces calculated as e1 e2 product' 
     158            IF(lflush) CALL FLUSH(numout) 
     159         ENDIF 
    141160         e1e2u (:,:) = e1u(:,:) * e2u(:,:)         ! compute them 
    142161         e1e2v (:,:) = e1v(:,:) * e2v(:,:)  
    143162      ELSE 
    144          IF(lwp) WRITE(numout,*) '          u- & v-surfaces have been read in "mesh_mask" file:' 
    145          IF(lwp) WRITE(numout,*) '                     grid size reduction in strait(s) is used' 
     163         IF(lwp) THEN 
     164            WRITE(numout,*) '          u- & v-surfaces have been read in "mesh_mask" file:' 
     165            WRITE(numout,*) '                     grid size reduction in strait(s) is used' 
     166            IF(lflush) CALL FLUSH(numout) 
     167         ENDIF 
    146168      ENDIF 
    147169      r1_e1e2u(:,:) = 1._wp / e1e2u(:,:)     ! compute their invert in any cases 
     
    185207         WRITE(numout,*) '   hgr_read : read the horizontal coordinates in mesh_mask' 
    186208         WRITE(numout,*) '   ~~~~~~~~      jpiglo = ', jpiglo, ' jpjglo = ', jpjglo, ' jpk = ', jpk 
     209         IF(lflush) CALL FLUSH(numout) 
    187210      ENDIF 
    188211      ! 
     
    211234      IF(  iom_varid( inum, 'ff_f', ldstop = .FALSE. ) > 0  .AND.  & 
    212235         & iom_varid( inum, 'ff_t', ldstop = .FALSE. ) > 0    ) THEN 
    213          IF(lwp) WRITE(numout,*) '           Coriolis factor at f- and t-points read in ', TRIM( cn_domcfg ), ' file' 
     236         IF(lwp) THEN 
     237            WRITE(numout,*) '           Coriolis factor at f- and t-points read in ', TRIM( cn_domcfg ), ' file' 
     238            IF(lflush) CALL FLUSH(numout) 
     239         ENDIF 
    214240         CALL iom_get( inum, jpdom_data, 'ff_f'  , pff_f  , lrowattr=ln_use_jattr ) 
    215241         CALL iom_get( inum, jpdom_data, 'ff_t'  , pff_t  , lrowattr=ln_use_jattr ) 
     
    220246      ! 
    221247      IF( iom_varid( inum, 'e1e2u', ldstop = .FALSE. ) > 0 ) THEN 
    222          IF(lwp) WRITE(numout,*) '           e1e2u & e1e2v read in ', TRIM( cn_domcfg ), ' file' 
     248         IF(lwp) THEN 
     249            WRITE(numout,*) '           e1e2u & e1e2v read in ', TRIM( cn_domcfg ), ' file' 
     250            IF(lflush) CALL FLUSH(numout) 
     251         ENDIF 
    223252         CALL iom_get( inum, jpdom_data, 'e1e2u'  , pe1e2u  , lrowattr=ln_use_jattr ) 
    224253         CALL iom_get( inum, jpdom_data, 'e1e2v'  , pe1e2v  , lrowattr=ln_use_jattr ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/dommsk.F90

    r10968 r10986  
    119119         WRITE(numout,*) '      lateral momentum boundary cond.    rn_shlat  = ',rn_shlat 
    120120         WRITE(numout,*) '      consistency with analytical form   ln_vorlat = ',ln_vorlat  
     121         IF(lflush) CALL FLUSH(numout) 
    121122      ENDIF 
    122123      ! 
     
    129130         CALL ctl_stop( 'dom_msk: wrong value for rn_shlat (i.e. a negalive value). We stop.' ) 
    130131      ENDIF 
    131  
     132      IF(lwp .AND. lflush) CALL FLUSH(numout) 
    132133      !  Ocean/land mask at t-point  (computed from ko_top and ko_bot) 
    133134      ! ---------------------------- 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domvvl.F90

    r10968 r10986  
    120120      !!---------------------------------------------------------------------- 
    121121      ! 
    122       IF(lwp) WRITE(numout,*) 
    123       IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
    124       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     122      IF(lwp) THEN 
     123         WRITE(numout,*) 
     124         WRITE(numout,*) 'dom_vvl_init : Variable volume activated' 
     125         WRITE(numout,*) '~~~~~~~~~~~~' 
     126         IF(lflush) CALL FLUSH(numout) 
     127      ENDIF 
    125128      ! 
    126129      CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer) 
     
    304307      ! 
    305308      IF( kt == nit000 ) THEN 
    306          IF(lwp) WRITE(numout,*) 
    307          IF(lwp) WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors' 
    308          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 
     309         IF(lwp) THEN 
     310            WRITE(numout,*) 
     311            WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors' 
     312            WRITE(numout,*) '~~~~~~~~~~~~~~' 
     313            IF(lflush) CALL FLUSH(numout) 
     314         ENDIF 
    309315      ENDIF 
    310316 
     
    533539         z_tmax = MAXVAL( tmask(:,:,1) *  ABS( ssha(:,:) ) ) 
    534540         CALL mpp_max( 'domvvl', z_tmax )                                ! max over the global domain 
    535          IF( lwp    ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 
     541         IF( lwp    ) THEN 
     542            WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 
     543            IF(lflush) CALL FLUSH(numout) 
     544         ENDIF 
    536545      END IF 
    537546 
     
    601610         IF(lwp) WRITE(numout,*) 'dom_vvl_sf_swp : - time filter and swap of scale factors' 
    602611         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   - interpolate scale factors and compute depths for next time step' 
     612         IF(lflush) CALL FLUSH(numout) 
    603613      ENDIF 
    604614      ! 
     
    820830               CALL iom_get( numror, jpdom_autoglo, 'e3t_n', e3t_n(:,:,:), ldxios = lrxios ) 
    821831               ! needed to restart if land processor not computed  
    822                IF(lwp) write(numout,*) 'dom_vvl_rst : e3t_b and e3t_n found in restart files' 
     832               IF(lwp) THEN 
     833                  write(numout,*) 'dom_vvl_rst : e3t_b and e3t_n found in restart files' 
     834                  IF(lflush) CALL FLUSH(numout) 
     835               ENDIF 
    823836               WHERE ( tmask(:,:,:) == 0.0_wp )  
    824837                  e3t_n(:,:,:) = e3t_0(:,:,:) 
     
    829842               ENDIF 
    830843            ELSE IF( id1 > 0 ) THEN 
    831                IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart files' 
    832                IF(lwp) write(numout,*) 'e3t_n set equal to e3t_b.' 
    833                IF(lwp) write(numout,*) 'neuler is forced to 0' 
     844               IF(lwp) THEN 
     845                  write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart files' 
     846                  write(numout,*) 'e3t_n set equal to e3t_b.' 
     847                  write(numout,*) 'neuler is forced to 0' 
     848                  IF(lflush) CALL FLUSH(numout) 
     849               ENDIF 
    834850               CALL iom_get( numror, jpdom_autoglo, 'e3t_b', e3t_b(:,:,:), ldxios = lrxios ) 
    835851               e3t_n(:,:,:) = e3t_b(:,:,:) 
    836852               neuler = 0 
    837853            ELSE IF( id2 > 0 ) THEN 
    838                IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_b not found in restart files' 
    839                IF(lwp) write(numout,*) 'e3t_b set equal to e3t_n.' 
    840                IF(lwp) write(numout,*) 'neuler is forced to 0' 
     854               IF(lwp) THEN 
     855                  write(numout,*) 'dom_vvl_rst WARNING : e3t_b not found in restart files' 
     856                  write(numout,*) 'e3t_b set equal to e3t_n.' 
     857                  write(numout,*) 'neuler is forced to 0' 
     858                  IF(lflush) CALL FLUSH(numout) 
     859               ENDIF 
    841860               CALL iom_get( numror, jpdom_autoglo, 'e3t_n', e3t_n(:,:,:), ldxios = lrxios ) 
    842861               e3t_b(:,:,:) = e3t_n(:,:,:) 
    843862               neuler = 0 
    844863            ELSE 
    845                IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart file' 
    846                IF(lwp) write(numout,*) 'Compute scale factor from sshn' 
    847                IF(lwp) write(numout,*) 'neuler is forced to 0' 
     864               IF(lwp) THEN 
     865                  write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart file' 
     866                  write(numout,*) 'Compute scale factor from sshn' 
     867                  write(numout,*) 'neuler is forced to 0' 
     868                  IF(lflush) CALL FLUSH(numout) 
     869               ENDIF 
    848870               DO jk = 1, jpk 
    849871                  e3t_n(:,:,jk) =  e3t_0(:,:,jk) * ( ht_0(:,:) + sshn(:,:) ) & 
     
    952974      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    953975         !                                   ! =================== 
    954          IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- dom_vvl_rst ----' 
     976         IF(lwp .AND. nprint > 0) THEN 
     977            WRITE(numout,*) '---- dom_vvl_rst ----' 
     978            IF(lflush) CALL FLUSH(numout) 
     979         ENDIF 
    955980         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    956981         !                                           ! --------- ! 
     
    10241049         ENDIF 
    10251050         WRITE(numout,*) '         debug prints flag                                 ln_vvl_dbg   = ', ln_vvl_dbg 
     1051         IF(lflush) CALL FLUSH(numout) 
    10261052      ENDIF 
    10271053      ! 
     
    10411067         IF( ln_vvl_layer           ) WRITE(numout,*) '      ==>>>   layer vertical coordinate is used' 
    10421068         IF( ln_vvl_ztilde_as_zstar ) WRITE(numout,*) '      ==>>>   to emulate a zstar coordinate' 
     1069         IF(lflush) CALL FLUSH(numout) 
    10431070      ENDIF 
    10441071      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domwri.F90

    r10888 r10986  
    6363      !!---------------------------------------------------------------------- 
    6464      ! 
    65       IF(lwp) WRITE(numout,*) 
    66       IF(lwp) WRITE(numout,*) 'dom_wri : create NetCDF mesh and mask information file(s)' 
    67       IF(lwp) WRITE(numout,*) '~~~~~~~' 
     65      IF(lwp) THEN 
     66         WRITE(numout,*) 
     67         WRITE(numout,*) 'dom_wri : create NetCDF mesh and mask information file(s)' 
     68         WRITE(numout,*) '~~~~~~~' 
     69         IF(lflush) CALL FLUSH(numout) 
     70      ENDIF 
    6871       
    6972      clnam = 'mesh_mask'  ! filename (mesh and mask informations) 
     
    283286         WRITE(numout,*) 'dom_stiff : maximum grid stiffness ratio: ', zrxmax 
    284287         WRITE(numout,*) '~~~~~~~~~' 
     288         IF(lflush) CALL FLUSH(numout) 
    285289      ENDIF 
    286290      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domzgr.F90

    r10968 r10986  
    8080         WRITE(numout,*) 'dom_zgr : vertical coordinate' 
    8181         WRITE(numout,*) '~~~~~~~' 
     82         IF(lflush) CALL FLUSH(numout) 
    8283      ENDIF 
    8384 
     
    8687 
    8788      IF( ln_read_cfg ) THEN        !==  read in mesh_mask.nc file  ==! 
    88          IF(lwp) WRITE(numout,*) 
    89          IF(lwp) WRITE(numout,*) '   ==>>>   Read vertical mesh in ', TRIM( cn_domcfg ), ' file' 
     89         IF(lwp) THEN 
     90            WRITE(numout,*) 
     91            WRITE(numout,*) '   ==>>>   Read vertical mesh in ', TRIM( cn_domcfg ), ' file' 
     92            IF(lflush) CALL FLUSH(numout) 
     93         ENDIF 
    9094         ! 
    9195         CALL zgr_read   ( ln_zco  , ln_zps  , ln_sco, ln_isfcav,   &  
     
    97101            ! 
    98102      ELSE                          !==  User defined configuration  ==! 
    99          IF(lwp) WRITE(numout,*) 
    100          IF(lwp) WRITE(numout,*) '          User defined vertical mesh (usr_def_zgr)' 
     103         IF(lwp) THEN 
     104            WRITE(numout,*) 
     105            WRITE(numout,*) '          User defined vertical mesh (usr_def_zgr)' 
     106            IF(lflush) CALL FLUSH(numout) 
     107         ENDIF 
    101108         ! 
    102109         CALL usr_def_zgr( ln_zco  , ln_zps  , ln_sco, ln_isfcav,   &  
     
    127134         WRITE(numout,*) '      s- or hybrid z-s-coordinate    ln_sco    = ', ln_sco 
    128135         WRITE(numout,*) '      ice shelf cavities             ln_isfcav = ', ln_isfcav 
     136         IF(lflush) CALL FLUSH(numout) 
    129137      ENDIF 
    130138 
     
    163171            &                          ' uw', MAXVAL(  e3uw_0(:,:,:) ), ' vw', MAXVAL(  e3vw_0(:,:,:) ),  & 
    164172            &                          ' w ', MAXVAL(   e3w_0(:,:,:) ) 
     173         IF(lflush) CALL FLUSH(numout) 
    165174      ENDIF 
    166175      ! 
     
    199208         WRITE(numout,*) '   zgr_read : read the vertical coordinates in ', TRIM( cn_domcfg ), ' file' 
    200209         WRITE(numout,*) '   ~~~~~~~~' 
     210         IF(lflush) CALL FLUSH(numout) 
    201211      ENDIF 
    202212      ! 
     
    248258            WRITE(numout, "(9x,' level  gdept_1d  gdepw_1d  e3t_1d   e3w_1d  ')" ) 
    249259            WRITE(numout, "(10x, i4, 4f9.2)" ) ( jk, pdept_1d(jk), pdepw_1d(jk), pe3t_1d(jk), pe3w_1d(jk), jk = 1, jpk ) 
     260            IF(lflush) CALL FLUSH(numout) 
    250261         ENDIF 
    251262      ENDIF 
     
    286297      !!---------------------------------------------------------------------- 
    287298      ! 
    288       IF(lwp) WRITE(numout,*) 
    289       IF(lwp) WRITE(numout,*) '    zgr_top_bot : ocean top and bottom k-index of T-, U-, V- and W-levels ' 
    290       IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~' 
     299      IF(lwp) THEN 
     300         WRITE(numout,*) 
     301         WRITE(numout,*) '    zgr_top_bot : ocean top and bottom k-index of T-, U-, V- and W-levels ' 
     302         WRITE(numout,*) '    ~~~~~~~~~~~' 
     303         IF(lflush) CALL FLUSH(numout) 
     304      ENDIF 
    291305      ! 
    292306      mikt(:,:) = MAX( k_top(:,:) , 1 )    ! top    ocean k-index of T-level (=1 over land) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/dtatsd.F90

    r10968 r10986  
    8787            WRITE(numout,*) '   ===>>   T & S data not used' 
    8888         ENDIF 
     89         IF(lflush) CALL FLUSH(numout) 
    8990      ENDIF 
    9091      ! 
     
    184185            WRITE(numout,*) 
    185186            WRITE(numout,*) 'dta_tsd: interpolates T & S data onto the s- or mixed s-z-coordinate mesh' 
     187            IF(lflush) CALL FLUSH(numout) 
    186188         ENDIF 
    187189         ! 
     
    243245      IF( .NOT.ln_tsd_dmp ) THEN                   !==   deallocate T & S structure   ==!  
    244246         !                                              (data used only for initialisation) 
    245          IF(lwp .AND. nprint > 1) WRITE(numout,*) 'dta_tsd: deallocte T & S arrays as they are only use to initialize the run' 
     247         IF(lwp .AND. nprint > 1) THEN 
     248            WRITE(numout,*) 'dta_tsd: deallocte T & S arrays as they are only use to initialize the run' 
     249            IF(lflush) CALL FLUSH(numout) 
     250         ENDIF 
    246251                                        DEALLOCATE( sf_tsd(jp_tem)%fnow )     ! T arrays in the structure 
    247252         IF( sf_tsd(jp_tem)%ln_tint )   DEALLOCATE( sf_tsd(jp_tem)%fdta ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/iscplhsb.F90

    r10888 r10986  
    170170                           WRITE(numout,*) '                     ',ji,' ',jj,' ',jk,' ',narea 
    171171                           WRITE(numout,*) ' we are now looking for the closest wet cell on the horizontal ' 
     172                           IF(lflush) CALL FLUSH(numout) 
    172173                           ENDIF 
    173174                        ! We deal with these points later. 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/iscplini.F90

    r10968 r10986  
    8383         WRITE(numout,*) ' coupling time step                       = ', rdt_iscpl 
    8484         WRITE(numout,*) ' number of call of the extrapolation loop = ', nn_drown 
     85         IF(lflush) CALL FLUSH(numout) 
    8586      ENDIF 
    8687      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/istate.F90

    r10888 r10986  
    6363      !!---------------------------------------------------------------------- 
    6464      ! 
    65       IF(lwp) WRITE(numout,*) 
    66       IF(lwp) WRITE(numout,*) 'istate_init : Initialization of the dynamics and tracers' 
    67       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     65      IF(lwp) THEN 
     66         WRITE(numout,*) 
     67         WRITE(numout,*) 'istate_init : Initialization of the dynamics and tracers' 
     68         WRITE(numout,*) '~~~~~~~~~~~' 
     69         IF(lflush) CALL FLUSH(numout) 
     70      ENDIF 
    6871 
    6972!!gm  Why not include in the first call of dta_tsd ?   
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/phycst.F90

    r10888 r10986  
    132132         WRITE(numout,*) 
    133133         WRITE(numout,*) '      smallest real computer value       rsmall = ', rsmall 
     134         IF(lflush) CALL FLUSH(numout) 
    134135      ENDIF 
    135136 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/divhor.F90

    r10888 r10986  
    7272         IF(lwp) WRITE(numout,*) 'div_hor : horizontal velocity divergence ' 
    7373         IF(lwp) WRITE(numout,*) '~~~~~~~   ' 
     74         IF(lflush) CALL FLUSH(numout) 
    7475      ENDIF 
    7576      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynadv.F90

    r10968 r10986  
    102102         WRITE(numout,*) 'dyn_adv_init : choice/control of the momentum advection scheme' 
    103103         WRITE(numout,*) '~~~~~~~~~~~~' 
     104         IF(lflush) CALL FLUSH(numout) 
    104105      ENDIF 
    105106      ! 
     
    119120         WRITE(numout,*) '      flux form: 2nd order centred scheme              ln_dynadv_cen2 = ', ln_dynadv_cen2 
    120121         WRITE(numout,*) '                 3rd order UBS scheme                  ln_dynadv_ubs  = ', ln_dynadv_ubs 
     122         IF(lflush) CALL FLUSH(numout) 
    121123      ENDIF 
    122124 
     
    141143         CASE( np_FLX_ubs )   ;   WRITE(numout,*) '   ==>>>   flux form   : UBS       scheme is used' 
    142144         END SELECT 
     145         IF(lflush) CALL FLUSH(numout) 
    143146      ENDIF 
    144147      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynadv_cen2.F90

    r10888 r10986  
    5757         WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection' 
    5858         WRITE(numout,*) '~~~~~~~~~~~~' 
     59         IF(lflush) CALL FLUSH(numout) 
    5960      ENDIF 
    6061      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynadv_ubs.F90

    r10888 r10986  
    8282         IF(lwp) WRITE(numout,*) 'dyn_adv_ubs : UBS flux form momentum advection' 
    8383         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     84         IF(lflush) CALL FLUSH(numout) 
    8485      ENDIF 
    8586      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynhpg.F90

    r10968 r10986  
    170170         WRITE(numout,*) '      s-coord. (Density Jacobian: Cubic polynomial)     ln_hpg_djc    = ', ln_hpg_djc 
    171171         WRITE(numout,*) '      s-coord. (Pressure Jacobian: Cubic polynomial)    ln_hpg_prj    = ', ln_hpg_prj 
     172         IF(lflush) CALL FLUSH(numout) 
    172173      ENDIF 
    173174      ! 
     
    211212         END SELECT 
    212213         WRITE(numout,*) 
     214         IF(lflush) CALL FLUSH(numout) 
    213215      ENDIF 
    214216      !                           
     
    218220      ELSE                            !--- set an ice shelf load 
    219221         ! 
    220          IF(lwp) WRITE(numout,*) 
    221          IF(lwp) WRITE(numout,*) '   ice shelf case: set the ice-shelf load' 
     222         IF(lwp) THEN 
     223            WRITE(numout,*) 
     224            WRITE(numout,*) '   ice shelf case: set the ice-shelf load' 
     225            IF(lflush) CALL FLUSH(numout) 
     226         ENDIF 
    222227         ALLOCATE( zts_top(jpi,jpj,jpts) , zrhd(jpi,jpj,jpk) , zrhdtop_isf(jpi,jpj) , ziceload(jpi,jpj) )  
    223228         ! 
     
    279284      !!---------------------------------------------------------------------- 
    280285      ! 
    281       IF( kt == nit000 ) THEN 
    282          IF(lwp) WRITE(numout,*) 
    283          IF(lwp) WRITE(numout,*) 'dyn:hpg_zco : hydrostatic pressure gradient trend' 
    284          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   z-coordinate case ' 
     286      IF( kt == nit000 .AND. lwp) THEN 
     287         WRITE(numout,*) 
     288         WRITE(numout,*) 'dyn:hpg_zco : hydrostatic pressure gradient trend' 
     289         WRITE(numout,*) '~~~~~~~~~~~   z-coordinate case ' 
     290         IF(lflush) CALL FLUSH(numout) 
    285291      ENDIF 
    286292 
     
    340346      !!---------------------------------------------------------------------- 
    341347      ! 
    342       IF( kt == nit000 ) THEN 
    343          IF(lwp) WRITE(numout,*) 
    344          IF(lwp) WRITE(numout,*) 'dyn:hpg_zps : hydrostatic pressure gradient trend' 
    345          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   z-coordinate with partial steps - vector optimization' 
     348      IF( kt == nit000 .AND. lwp) THEN 
     349         WRITE(numout,*) 
     350         WRITE(numout,*) 'dyn:hpg_zps : hydrostatic pressure gradient trend' 
     351         WRITE(numout,*) '~~~~~~~~~~~   z-coordinate with partial steps - vector optimization' 
     352         IF(lflush) CALL FLUSH(numout) 
    346353      ENDIF 
    347354 
     
    439446      IF( ln_wd_il ) ALLOCATE(zcpx(jpi,jpj), zcpy(jpi,jpj)) 
    440447      ! 
    441       IF( kt == nit000 ) THEN 
    442          IF(lwp) WRITE(numout,*) 
    443          IF(lwp) WRITE(numout,*) 'dyn:hpg_sco : hydrostatic pressure gradient trend' 
    444          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, OPA original scheme used' 
     448      IF( kt == nit000 .AND. lwp) THEN 
     449         WRITE(numout,*) 
     450         WRITE(numout,*) 'dyn:hpg_sco : hydrostatic pressure gradient trend' 
     451         WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, OPA original scheme used' 
     452         IF(lflush) CALL FLUSH(numout) 
    445453      ENDIF 
    446454      ! 
     
    726734      END IF 
    727735 
    728       IF( kt == nit000 ) THEN 
    729          IF(lwp) WRITE(numout,*) 
    730          IF(lwp) WRITE(numout,*) 'dyn:hpg_djc : hydrostatic pressure gradient trend' 
    731          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, density Jacobian with cubic polynomial scheme' 
     736      IF( kt == nit000 .AND. lwp) THEN 
     737         WRITE(numout,*) 
     738         WRITE(numout,*) 'dyn:hpg_djc : hydrostatic pressure gradient trend' 
     739         WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, density Jacobian with cubic polynomial scheme' 
     740         IF(lflush) CALL FLUSH(numout) 
    732741      ENDIF 
    733742 
     
    960969      !!---------------------------------------------------------------------- 
    961970      ! 
    962       IF( kt == nit000 ) THEN 
    963          IF(lwp) WRITE(numout,*) 
    964          IF(lwp) WRITE(numout,*) 'dyn:hpg_prj : hydrostatic pressure gradient trend' 
    965          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, cubic spline pressure Jacobian' 
     971      IF( kt == nit000 .AND. lwp) THEN 
     972         WRITE(numout,*) 
     973         WRITE(numout,*) 'dyn:hpg_prj : hydrostatic pressure gradient trend' 
     974         WRITE(numout,*) '~~~~~~~~~~~   s-coordinate case, cubic spline pressure Jacobian' 
     975         IF(lflush) CALL FLUSH(numout) 
    966976      ENDIF 
    967977 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynkeg.F90

    r10888 r10986  
    8484      IF( ln_timing )   CALL timing_start('dyn_keg') 
    8585      ! 
    86       IF( kt == nit000 ) THEN 
    87          IF(lwp) WRITE(numout,*) 
    88          IF(lwp) WRITE(numout,*) 'dyn_keg : kinetic energy gradient trend, scheme number=', kscheme 
    89          IF(lwp) WRITE(numout,*) '~~~~~~~' 
     86      IF( kt == nit000 .AND. lwp) THEN 
     87         WRITE(numout,*) 
     88         WRITE(numout,*) 'dyn_keg : kinetic energy gradient trend, scheme number=', kscheme 
     89         WRITE(numout,*) '~~~~~~~' 
     90         IF(lflush) CALL FLUSH(numout) 
    9091      ENDIF 
    9192 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynldf.F90

    r10888 r10986  
    106106         CASE( np_blp    )   ;   WRITE(numout,*) '   ==>>>   iso-level bi-laplacian operator' 
    107107         END SELECT 
     108         IF(lflush) CALL FLUSH(numout) 
    108109      ENDIF 
    109110      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynldf_iso.F90

    r10888 r10986  
    114114      ! 
    115115      IF( kt == nit000 ) THEN 
    116          IF(lwp) WRITE(numout,*) 
    117          IF(lwp) WRITE(numout,*) 'dyn_ldf_iso : iso-neutral laplacian diffusive operator or ' 
    118          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   s-coordinate horizontal diffusive operator' 
     116         IF(lwp) THEN 
     117            WRITE(numout,*) 
     118            WRITE(numout,*) 'dyn_ldf_iso : iso-neutral laplacian diffusive operator or ' 
     119            WRITE(numout,*) '~~~~~~~~~~~   s-coordinate horizontal diffusive operator' 
     120            IF(lflush) CALL FLUSH(numout) 
     121         ENDIF 
    119122         !                                      ! allocate dyn_ldf_bilap arrays 
    120123         IF( dyn_ldf_iso_alloc() /= 0 )   CALL ctl_stop('STOP', 'dyn_ldf_iso: failed to allocate arrays') 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynldf_lap_blp.F90

    r10888 r10986  
    6262         WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator, pass=', kpass 
    6363         WRITE(numout,*) '~~~~~~~ ' 
     64         IF(lflush) CALL FLUSH(numout) 
    6465      ENDIF 
    6566      ! 
     
    129130         IF(lwp) WRITE(numout,*) 'dyn_ldf_blp : bilaplacian operator momentum ' 
    130131         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
     132         IF(lflush) CALL FLUSH(numout) 
    131133      ENDIF 
    132134      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynnxt.F90

    r10888 r10986  
    106106      IF( l_trddyn     )   ALLOCATE( zua(jpi,jpj,jpk) , zva(jpi,jpj,jpk) ) 
    107107      ! 
    108       IF( kt == nit000 ) THEN 
    109          IF(lwp) WRITE(numout,*) 
    110          IF(lwp) WRITE(numout,*) 'dyn_nxt : time stepping' 
    111          IF(lwp) WRITE(numout,*) '~~~~~~~' 
     108      IF( kt == nit000 .AND. lwp) THEN 
     109         WRITE(numout,*) 
     110         WRITE(numout,*) 'dyn_nxt : time stepping' 
     111         WRITE(numout,*) '~~~~~~~' 
     112         IF(lflush) CALL FLUSH(numout) 
    112113      ENDIF 
    113114 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynspg.F90

    r10968 r10986  
    198198         WRITE(numout,*) 'dyn_spg_init : choice of the surface pressure gradient scheme' 
    199199         WRITE(numout,*) '~~~~~~~~~~~~' 
     200         IF(lflush) CALL FLUSH(numout) 
    200201      ENDIF 
    201202      ! 
     
    213214         WRITE(numout,*) '      Explicit free surface                  ln_dynspg_exp = ', ln_dynspg_exp 
    214215         WRITE(numout,*) '      Free surface with time splitting       ln_dynspg_ts  = ', ln_dynspg_ts 
     216         IF(lflush) CALL FLUSH(numout) 
    215217      ENDIF 
    216218      !                          ! Control of surface pressure gradient scheme options 
     
    229231         IF( nspg == np_TS  )   WRITE(numout,*) '   ==>>>   free surface with time splitting scheme' 
    230232         IF( nspg == np_NO  )   WRITE(numout,*) '   ==>>>   No surface surface pressure gradient trend in momentum Eqs.' 
     233         IF(lflush) CALL FLUSH(numout) 
    231234      ENDIF 
    232235      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynspg_exp.F90

    r10888 r10986  
    6161      ! 
    6262      IF( kt == nit000 ) THEN 
    63          IF(lwp) WRITE(numout,*) 
    64          IF(lwp) WRITE(numout,*) 'dyn_spg_exp : surface pressure gradient trend' 
    65          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   (explicit free surface)' 
     63         IF(lwp) THEN 
     64            WRITE(numout,*) 
     65            WRITE(numout,*) 'dyn_spg_exp : surface pressure gradient trend' 
     66            WRITE(numout,*) '~~~~~~~~~~~   (explicit free surface)' 
     67            IF(lflush) CALL FLUSH(numout) 
     68         ENDIF 
    6669         ! 
    6770         spgu(:,:) = 0._wp   ;   spgv(:,:) = 0._wp 
    6871         ! 
    69          IF( .NOT.ln_linssh .AND. lwp ) WRITE(numout,*) '      non linear free surface: spg is included in dynhpg' 
     72         IF( .NOT.ln_linssh .AND. lwp ) THEN 
     73            WRITE(numout,*) '      non linear free surface: spg is included in dynhpg' 
     74            IF(lflush) CALL FLUSH(numout) 
     75         ENDIF 
    7076      ENDIF 
    7177 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynspg_ts.F90

    r10968 r10986  
    197197      IF( kt == nit000 ) THEN                   !* initialisation 
    198198         ! 
    199          IF(lwp) WRITE(numout,*) 
    200          IF(lwp) WRITE(numout,*) 'dyn_spg_ts : surface pressure gradient trend' 
    201          IF(lwp) WRITE(numout,*) '~~~~~~~~~~   free surface with time splitting' 
    202          IF(lwp) WRITE(numout,*) 
     199         IF(lwp) THEN 
     200            WRITE(numout,*) 
     201            WRITE(numout,*) 'dyn_spg_ts : surface pressure gradient trend' 
     202            WRITE(numout,*) '~~~~~~~~~~   free surface with time splitting' 
     203            WRITE(numout,*) 
     204            IF(lflush) CALL FLUSH(numout) 
     205         ENDIF 
    203206         ! 
    204207         IF( neuler == 0 )   ll_init=.TRUE. 
     
    14221425#endif 
    14231426         ELSE                                   !* Start from rest 
    1424             IF(lwp) WRITE(numout,*) 
    1425             IF(lwp) WRITE(numout,*) '   ==>>>   start from rest: set barotropic values to 0' 
     1427            IF(lwp) THEN 
     1428               WRITE(numout,*) 
     1429               WRITE(numout,*) '   ==>>>   start from rest: set barotropic values to 0' 
     1430               IF(lflush) CALL FLUSH(numout)  
     1431            ENDIF 
    14261432            ub2_b (:,:) = 0._wp   ;   vb2_b (:,:) = 0._wp   ! used in the 1st interpol of agrif 
    14271433            un_adv(:,:) = 0._wp   ;   vn_adv(:,:) = 0._wp   ! used in the 1st interpol of agrif 
     
    14361442      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    14371443         !                                   ! ------------------- 
    1438          IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- ts_rst ----' 
     1444         IF(lwp .AND. nprint > 0) THEN 
     1445            WRITE(numout,*) '---- ts_rst ----' 
     1446            IF(lflush) CALL FLUSH(numout) 
     1447         ENDIF 
    14391448         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    14401449         CALL iom_rstput( kt, nitrst, numrow, 'ub2_b'   , ub2_b  (:,:), ldxios = lwxios ) 
     
    14941503      zcmax = zcmax * rdtbt 
    14951504      ! Print results 
    1496       IF(lwp) WRITE(numout,*) 
    1497       IF(lwp) WRITE(numout,*) 'dyn_spg_ts_init : split-explicit free surface' 
    1498       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 
    1499       IF( ln_bt_auto ) THEN 
    1500          IF(lwp) WRITE(numout,*) '     ln_ts_auto =.true. Automatically set nn_baro ' 
    1501          IF(lwp) WRITE(numout,*) '     Max. courant number allowed: ', rn_bt_cmax 
    1502       ELSE 
    1503          IF(lwp) WRITE(numout,*) '     ln_ts_auto=.false.: Use nn_baro in namelist   nn_baro = ', nn_baro 
    1504       ENDIF 
    1505  
    1506       IF(ln_bt_av) THEN 
    1507          IF(lwp) WRITE(numout,*) '     ln_bt_av =.true.  ==> Time averaging over nn_baro time steps is on ' 
    1508       ELSE 
    1509          IF(lwp) WRITE(numout,*) '     ln_bt_av =.false. => No time averaging of barotropic variables ' 
    1510       ENDIF 
    1511       ! 
    1512       ! 
    1513       IF(ln_bt_fw) THEN 
    1514          IF(lwp) WRITE(numout,*) '     ln_bt_fw=.true.  => Forward integration of barotropic variables ' 
    1515       ELSE 
    1516          IF(lwp) WRITE(numout,*) '     ln_bt_fw =.false.=> Centred integration of barotropic variables ' 
     1505      IF(lwp) THEN 
     1506         WRITE(numout,*) 
     1507         WRITE(numout,*) 'dyn_spg_ts_init : split-explicit free surface' 
     1508         WRITE(numout,*) '~~~~~~~~~~~~~~~' 
     1509         IF( ln_bt_auto ) THEN 
     1510            WRITE(numout,*) '     ln_ts_auto =.true. Automatically set nn_baro ' 
     1511            WRITE(numout,*) '     Max. courant number allowed: ', rn_bt_cmax 
     1512         ELSE 
     1513            WRITE(numout,*) '     ln_ts_auto=.false.: Use nn_baro in namelist   nn_baro = ', nn_baro 
     1514         ENDIF 
     1515 
     1516         IF(ln_bt_av) THEN 
     1517            WRITE(numout,*) '     ln_bt_av =.true.  ==> Time averaging over nn_baro time steps is on ' 
     1518         ELSE 
     1519            WRITE(numout,*) '     ln_bt_av =.false. => No time averaging of barotropic variables ' 
     1520         ENDIF 
     1521      ! 
     1522      ! 
     1523         IF(ln_bt_fw) THEN 
     1524            WRITE(numout,*) '     ln_bt_fw=.true.  => Forward integration of barotropic variables ' 
     1525         ELSE 
     1526            WRITE(numout,*) '     ln_bt_fw =.false.=> Centred integration of barotropic variables ' 
     1527         ENDIF 
     1528         IF(lflush) CALL FLUSH(numout) 
    15171529      ENDIF 
    15181530      ! 
     
    15301542      END SELECT 
    15311543      ! 
    1532       IF(lwp) WRITE(numout,*) ' ' 
    1533       IF(lwp) WRITE(numout,*) '     nn_baro = ', nn_baro 
    1534       IF(lwp) WRITE(numout,*) '     Barotropic time step [s] is :', rdtbt 
    1535       IF(lwp) WRITE(numout,*) '     Maximum Courant number is   :', zcmax 
    1536       ! 
    1537       IF(lwp) WRITE(numout,*)    '     Time diffusion parameter rn_bt_alpha: ', rn_bt_alpha 
     1544      IF(lwp) THEN 
     1545         WRITE(numout,*) ' ' 
     1546         WRITE(numout,*) '     nn_baro = ', nn_baro 
     1547         WRITE(numout,*) '     Barotropic time step [s] is :', rdtbt 
     1548         WRITE(numout,*) '     Maximum Courant number is   :', zcmax 
     1549      ! 
     1550         WRITE(numout,*)    '     Time diffusion parameter rn_bt_alpha: ', rn_bt_alpha 
     1551         IF(lflush) CALL FLUSH(numout) 
     1552      ENDIF 
    15381553      IF ((ln_bt_av.AND.nn_bt_flt/=0).AND.(rn_bt_alpha>0._wp)) THEN 
    15391554         CALL ctl_stop( 'dynspg_ts ERROR: if rn_bt_alpha > 0, remove temporal averaging' ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynvor.F90

    r10968 r10986  
    220220         IF(lwp) WRITE(numout,*) 'dyn:vor_enT : vorticity term: t-point energy conserving scheme' 
    221221         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     222         IF(lflush) CALL FLUSH(numout) 
    222223      ENDIF 
    223224      ! 
     
    352353      !!---------------------------------------------------------------------- 
    353354      ! 
    354       IF( kt == nit000 ) THEN 
    355          IF(lwp) WRITE(numout,*) 
    356          IF(lwp) WRITE(numout,*) 'dyn:vor_ene : vorticity term: energy conserving scheme' 
    357          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     355      IF( kt == nit000 .AND. lwp) THEN 
     356         WRITE(numout,*) 
     357         WRITE(numout,*) 'dyn:vor_ene : vorticity term: energy conserving scheme' 
     358         WRITE(numout,*) '~~~~~~~~~~~' 
     359         IF(lflush) CALL FLUSH(numout) 
    358360      ENDIF 
    359361      ! 
     
    460462      !!---------------------------------------------------------------------- 
    461463      ! 
    462       IF( kt == nit000 ) THEN 
    463          IF(lwp) WRITE(numout,*) 
    464          IF(lwp) WRITE(numout,*) 'dyn:vor_ens : vorticity term: enstrophy conserving scheme' 
    465          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     464      IF( kt == nit000 .AND. lwp) THEN 
     465         WRITE(numout,*) 
     466         WRITE(numout,*) 'dyn:vor_ens : vorticity term: enstrophy conserving scheme' 
     467         WRITE(numout,*) '~~~~~~~~~~~' 
     468         IF(lflush) CALL FLUSH(numout) 
    466469      ENDIF 
    467470      !                                                ! =============== 
     
    568571      !!---------------------------------------------------------------------- 
    569572      ! 
    570       IF( kt == nit000 ) THEN 
    571          IF(lwp) WRITE(numout,*) 
    572          IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
    573          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     573      IF( kt == nit000 .AND. lwp) THEN 
     574         WRITE(numout,*) 
     575         WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
     576         WRITE(numout,*) '~~~~~~~~~~~' 
     577         IF(lflush) CALL FLUSH(numout) 
    574578      ENDIF 
    575579      ! 
     
    725729      !!---------------------------------------------------------------------- 
    726730      ! 
    727       IF( kt == nit000 ) THEN 
    728          IF(lwp) WRITE(numout,*) 
    729          IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
    730          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     731      IF( kt == nit000 .AND. lwp) THEN 
     732         WRITE(numout,*) 
     733         WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
     734         WRITE(numout,*) '~~~~~~~~~~~' 
     735         IF(lflush) CALL FLUSH(numout) 
    731736      ENDIF 
    732737      ! 
     
    847852         WRITE(numout,*) 'dyn_vor_init : vorticity term : read namelist and control the consistency' 
    848853         WRITE(numout,*) '~~~~~~~~~~~~' 
     854         IF(lflush) CALL FLUSH(numout) 
    849855      ENDIF 
    850856      ! 
     
    867873         WRITE(numout,*) '      mixed enstrophy/energy conserving scheme       ln_dynvor_mix = ', ln_dynvor_mix 
    868874         WRITE(numout,*) '      masked (=T) or unmasked(=F) vorticity          ln_dynvor_msk = ', ln_dynvor_msk 
     875         IF(lflush) CALL FLUSH(numout) 
    869876      ENDIF 
    870877 
     
    874881      ! If energy, enstrophy or mixed advection of momentum in vector form change the value for masks 
    875882      ! at angles with three ocean points and one land point 
    876       IF(lwp) WRITE(numout,*) 
    877       IF(lwp) WRITE(numout,*) '      change fmask value in the angles (T)           ln_vorlat = ', ln_vorlat 
     883      IF(lwp) THEN 
     884         WRITE(numout,*) 
     885         WRITE(numout,*) '      change fmask value in the angles (T)           ln_vorlat = ', ln_vorlat 
     886         IF(lflush) CALL FLUSH(numout) 
     887      ENDIF 
    878888      IF( ln_vorlat .AND. ( ln_dynvor_ene .OR. ln_dynvor_ens .OR. ln_dynvor_mix ) ) THEN 
    879889         DO jk = 1, jpk 
     
    951961         CASE( np_MIX )   ;   WRITE(numout,*) '   ==>>>   mixed enstrophy/energy conserving scheme (MIX)' 
    952962         END SELECT          
     963         IF(lflush) CALL FLUSH(numout) 
    953964      ENDIF 
    954965      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynzad.F90

    r10888 r10986  
    6666         IF(lwp) WRITE(numout,*) 
    6767         IF(lwp) WRITE(numout,*) 'dyn_zad : 2nd order vertical advection scheme' 
     68         IF(lflush) CALL FLUSH(numout) 
    6869      ENDIF 
    6970 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/dynzdf.F90

    r10888 r10986  
    8282      ! 
    8383      IF( kt == nit000 ) THEN       !* initialization 
    84          IF(lwp) WRITE(numout,*) 
    85          IF(lwp) WRITE(numout,*) 'dyn_zdf_imp : vertical momentum diffusion implicit operator' 
    86          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 
     84         IF(lwp) THEN 
     85            WRITE(numout,*) 
     86            WRITE(numout,*) 'dyn_zdf_imp : vertical momentum diffusion implicit operator' 
     87            WRITE(numout,*) '~~~~~~~~~~~ ' 
     88            IF(lflush) CALL FLUSH(numout) 
     89         ENDIF 
    8790         ! 
    8891         If( ln_linssh ) THEN   ;    r_vvl = 0._wp    ! non-linear free surface indicator 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/sshwzv.F90

    r10888 r10986  
    7777      IF( ln_timing )   CALL timing_start('ssh_nxt') 
    7878      ! 
    79       IF( kt == nit000 ) THEN 
    80          IF(lwp) WRITE(numout,*) 
    81          IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height' 
    82          IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     79      IF( kt == nit000 .AND. lwp) THEN 
     80         WRITE(numout,*) 
     81         WRITE(numout,*) 'ssh_nxt : after sea surface height' 
     82         WRITE(numout,*) '~~~~~~~ ' 
     83         IF(lflush) CALL FLUSH(numout) 
    8384      ENDIF 
    8485      ! 
     
    152153      ! 
    153154      IF( kt == nit000 ) THEN 
    154          IF(lwp) WRITE(numout,*) 
    155          IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity ' 
    156          IF(lwp) WRITE(numout,*) '~~~~~ ' 
     155         IF(lwp) THEN 
     156            WRITE(numout,*) 
     157            WRITE(numout,*) 'wzv : now vertical velocity ' 
     158            WRITE(numout,*) '~~~~~ ' 
     159            IF(lflush) CALL FLUSH(numout) 
     160         ENDIF 
    157161         ! 
    158162         wn(:,:,jpk) = 0._wp                  ! bottom boundary condition: w=0 (set once for all) 
     
    240244      IF( ln_timing )   CALL timing_start('ssh_swp') 
    241245      ! 
    242       IF( kt == nit000 ) THEN 
    243          IF(lwp) WRITE(numout,*) 
    244          IF(lwp) WRITE(numout,*) 'ssh_swp : Asselin time filter and swap of sea surface height' 
    245          IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     246      IF( kt == nit000 .AND. lwp) THEN 
     247         WRITE(numout,*) 
     248         WRITE(numout,*) 'ssh_swp : Asselin time filter and swap of sea surface height' 
     249         WRITE(numout,*) '~~~~~~~ ' 
     250         IF(lflush) CALL FLUSH(numout) 
    246251      ENDIF 
    247252      !              !==  Euler time-stepping: no filter, just swap  ==! 
     
    294299      ! 
    295300      IF( kt == nit000 ) THEN 
    296          IF(lwp) WRITE(numout,*) 
    297          IF(lwp) WRITE(numout,*) 'wAimp : Courant number-based partitioning of now vertical velocity ' 
    298          IF(lwp) WRITE(numout,*) '~~~~~ ' 
     301         IF(lwp) THEN 
     302            WRITE(numout,*) 
     303            WRITE(numout,*) 'wAimp : Courant number-based partitioning of now vertical velocity ' 
     304            WRITE(numout,*) '~~~~~ ' 
     305            IF(lflush) CALL FLUSH(numout) 
     306         ENDIF 
    299307         ! 
    300308         Cu_adv(:,:,jpk) = 0._wp              ! bottom value : Cu_adv=0 (set once for all) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DYN/wet_dry.F90

    r10968 r10986  
    104104         WRITE(numout,*) '      cut off depth sbc for wd   rn_wd_sbcdep       = ', rn_wd_sbcdep 
    105105         WRITE(numout,*) '      fraction to start sbc wgt rn_wd_sbcfra        = ', rn_wd_sbcfra 
     106         IF(lflush) CALL FLUSH(numout) 
    106107      ENDIF 
    107108      IF( .NOT. ln_read_cfg ) THEN 
    108          IF(lwp) WRITE(numout,*) '      No configuration file so seting ssh_ref to zero  ' 
     109         IF(lwp) THEN 
     110            WRITE(numout,*) '      No configuration file so seting ssh_ref to zero  ' 
     111            IF(lflush) CALL FLUSH(numout) 
     112         ENDIF 
    109113         ssh_ref=0._wp 
    110114      ENDIF 
     
    266270!!gm 
    267271      ! 
    268       IF(jflag == 1 .AND. lwp)   WRITE(numout,*) 'Need more iterations in wad_lmt!!!' 
     272      IF(jflag == 1 .AND. lwp)   THEN 
     273         WRITE(numout,*) 'Need more iterations in wad_lmt!!!' 
     274         IF(lflush) CALL FLUSH(numout) 
     275      ENDIF 
    269276      ! 
    270277      !IF( ln_rnf      )   CALL sbc_rnf_div( hdivn )          ! runoffs (update hdivn field) 
     
    390397!!gm end 
    391398      ! 
    392       IF( jflag == 1 .AND. lwp )   WRITE(numout,*) 'Need more iterations in wad_lmt_bt!!!' 
     399      IF( jflag == 1 .AND. lwp )   THEN 
     400         WRITE(numout,*) 'Need more iterations in wad_lmt_bt!!!' 
     401         IF(lflush) CALL FLUSH(numout) 
     402      ENDIF 
    393403      ! 
    394404      !IF( ln_rnf      )   CALL sbc_rnf_div( hdivn )          ! runoffs (update hdivn field) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/FLO/flo4rk.F90

    r10888 r10986  
    5858      ! 
    5959      IF( ierror /= 0 ) THEN 
    60          WRITE(numout,*) 'flo_4rk: allocation of workspace arrays failed' 
     60         IF(lwp) THEN 
     61            WRITE(numout,*) 'flo_4rk: allocation of workspace arrays failed' 
     62            IF(lflush) CALL FLUSH(numout) 
     63         ENDIF 
    6164      ENDIF 
    6265 
    6366     
    64       IF( kt == nit000 ) THEN 
    65          IF(lwp) WRITE(numout,*) 
    66          IF(lwp) WRITE(numout,*) 'flo_4rk : compute Runge Kutta trajectories for floats ' 
    67          IF(lwp) WRITE(numout,*) '~~~~~~~' 
     67      IF( kt == nit000 .AND. lwp) THEN 
     68         WRITE(numout,*) 
     69         WRITE(numout,*) 'flo_4rk : compute Runge Kutta trajectories for floats ' 
     70         WRITE(numout,*) '~~~~~~~' 
     71         IF(lflush) CALL FLUSH(numout) 
    6872      ENDIF 
    6973 
     
    114118         ENDIF 
    115119      END DO 
     120      
     121      IF(lflush .AND. lwp) CALL FLUSH(numout) 
    116122       
    117123      ! 4 steps of Runge-Kutta algorithme 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/FLO/floats.F90

    r10968 r10986  
    8484      !!--------------------------------------------------------------------- 
    8585      ! 
    86       IF(lwp) WRITE(numout,*) 
    87       IF(lwp) WRITE(numout,*) 'flo_stp : call floats routine ' 
    88       IF(lwp) WRITE(numout,*) '~~~~~~~' 
     86      IF(lwp) THEN 
     87         WRITE(numout,*) 
     88         WRITE(numout,*) 'flo_stp : call floats routine ' 
     89         WRITE(numout,*) '~~~~~~~' 
     90         IF(lflush) CALL FLUSH(numout) 
     91      ENDIF 
    8992 
    9093      REWIND( numnam_ref )              ! Namelist namflo in reference namelist : Floats 
     
    109112         WRITE(numout,*) '            Use of ariane convention              ln_ariane    = ', ln_ariane 
    110113         WRITE(numout,*) '            ascii output (T) or netcdf output (F) ln_flo_ascii = ', ln_flo_ascii 
    111  
     114         IF(lflush) CALL FLUSH(numout) 
    112115      ENDIF 
    113116      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/FLO/floblk.F90

    r10888 r10986  
    7171      !!--------------------------------------------------------------------- 
    7272 
    73       IF( kt == nit000 ) THEN 
    74          IF(lwp) WRITE(numout,*) 
    75          IF(lwp) WRITE(numout,*) 'flo_blk : compute Blanke trajectories for floats ' 
    76          IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
     73      IF( kt == nit000 .AND. lwp) THEN 
     74         WRITE(numout,*) 
     75         WRITE(numout,*) 'flo_blk : compute Blanke trajectories for floats ' 
     76         WRITE(numout,*) '~~~~~~~ ' 
     77         IF(lflush) CALL FLUSH(numout) 
    7778      ENDIF 
    7879 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/FLO/flodom.F90

    r10888 r10986  
    5959      ! Initialisation with the geographical position or restart 
    6060       
    61       IF(lwp) WRITE(numout,*) 'flo_dom : compute initial position of floats' 
    62       IF(lwp) WRITE(numout,*) '~~~~~~~~' 
    63       IF(lwp) WRITE(numout,*) '           jpnfl = ',jpnfl 
     61      IF(lwp) THEN 
     62         WRITE(numout,*) 'flo_dom : compute initial position of floats' 
     63         WRITE(numout,*) '~~~~~~~~' 
     64         WRITE(numout,*) '           jpnfl = ',jpnfl 
     65         IF(lflush) CALL FLUSH(numout) 
     66      ENDIF 
    6467       
    6568      !-------------------------! 
     
    8992         IF( jpnfl > jpnrstflo )THEN 
    9093 
    91             IF(lwp) WRITE(numout,*) '        add new floats' 
     94            IF(lwp) THEN 
     95               WRITE(numout,*) '        add new floats' 
     96               IF(lflush) CALL FLUSH(numout) 
     97            ENDIF 
    9298 
    9399            IF( ln_ariane )THEN  !Add new floats with ariane convention 
     
    145151      DO jfl = kfl_start,kfl_end 
    146152         READ(inum,*) flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash 
    147          if(lwp)write(numout,*)'read:',jfl,flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash ; call flush(numout) 
     153         IF(lwp) THEN 
     154            write(numout,*)'read:',jfl,flxx(jfl),flyy(jfl),flzz(jfl), nisobfl(jfl),ngrpfl(jfl),itrash 
     155            IF(lflush) CALL FLUSH(numout) 
     156         ENDIF 
    148157      END DO 
    149158      CLOSE(inum) 
     
    266275         DO jfl =  kfl_start,kfl_end 
    267276            IF( idomfl(jfl) == 1 ) THEN 
    268                IF(lwp) WRITE(numout,*)'*****************************' 
    269                IF(lwp) WRITE(numout,*)'!!!!!!!  WARNING   !!!!!!!!!!' 
    270                IF(lwp) WRITE(numout,*)'*****************************' 
    271                IF(lwp) WRITE(numout,*)'The float number',jfl,'is out of the sea.' 
    272                IF(lwp) WRITE(numout,*)'geographical position',flxx(jfl),flyy(jfl),flzz(jfl) 
    273                IF(lwp) WRITE(numout,*)'index position',tpifl(jfl),tpjfl(jfl),tpkfl(jfl) 
     277               IF(lwp) THEN 
     278                  WRITE(numout,*)'*****************************' 
     279                  WRITE(numout,*)'!!!!!!!  WARNING   !!!!!!!!!!' 
     280                  WRITE(numout,*)'*****************************' 
     281                  WRITE(numout,*)'The float number',jfl,'is out of the sea.' 
     282                  WRITE(numout,*)'geographical position',flxx(jfl),flyy(jfl),flzz(jfl) 
     283                  WRITE(numout,*)'index position',tpifl(jfl),tpjfl(jfl),tpkfl(jfl) 
     284                  IF(lflush) CALL FLUSH(numout) 
     285               ENDIF 
    274286            ENDIF 
    275287         END DO 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/FLO/florst.F90

    r10888 r10986  
    6767            WRITE(numout,*) 'flo_rst : write in  restart_float file ' 
    6868            WRITE(numout,*) '~~~~~~~    ' 
     69            IF(lflush) CALL FLUSH(numout) 
    6970         ENDIF 
    7071 
     
    118119                  ENDIF 
    119120               END DO 
     121               IF(lflush) CALL FLUSH(numout) 
    120122            ENDIF 
    121123         ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/IOM/in_out_manager.F90

    r10888 r10986  
    129129   INTEGER ::   nn_isplt         !: number of processors following i 
    130130   INTEGER ::   nn_jsplt         !: number of processors following j 
     131   LOGICAL ::   ln_flush         !: flush numout flag 
    131132   !                                           
    132133   INTEGER ::   nprint, nictls, nictle, njctls, njctle, isplt, jsplt    !: OLD namelist names 
     134   LOGICAL ::   lflush 
    133135 
    134136   INTEGER ::   ijsplt     =    1      !: nb of local domain = nb of processors 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/IOM/iom.F90

    r10968 r10986  
    321321        ENDIF 
    322322!set name of the restart file and enable available fields 
    323         if(lwp .AND. nprint > 1) WRITE(numout,*) 'Setting restart filename (for XIOS) to: ',rst_file 
     323        if(lwp .AND. nprint > 1) THEN 
     324           WRITE(numout,*) 'Setting restart filename (for XIOS) to: ',rst_file 
     325           IF(lflush) CALL FLUSH(numout) 
     326        ENDIF 
    324327        CALL xios_get_handle("file_definition", filegroup_hdl ) 
    325328        CALL xios_add_child(filegroup_hdl, file_hdl, 'rrestart') 
     
    346349                        scalar_ref = "grid_scalar", operation = "instant") 
    347350                END SELECT 
    348                 IF(lwp .AND. nprint > 1) WRITE(numout,*) 'XIOS read: ', TRIM(rst_rfields(i)%vname), ' enabled in ', TRIM(rst_file) 
     351                IF(lwp .AND. nprint > 1) THEN 
     352                   WRITE(numout,*) 'XIOS read: ', TRIM(rst_rfields(i)%vname), ' enabled in ', TRIM(rst_file) 
     353                   IF(lflush) CALL FLUSH(numout) 
     354                ENDIF 
    349355           ENDIF 
    350356         ENDIF 
     
    537543 
    538544!set name of the restart file and enable available fields 
    539         IF(lwp .AND. nprint > 1) WRITE(numout,*) 'Setting restart filename (for XIOS write) to: ',cdrst_file 
     545        IF(lwp .AND. nprint > 1) THEN 
     546           WRITE(numout,*) 'Setting restart filename (for XIOS write) to: ',cdrst_file 
     547           IF(lflush) CALL FLUSH(numout) 
     548        ENDIF 
    540549        CALL xios_get_handle("file_definition", filegroup_hdl ) 
    541550        CALL xios_add_child(filegroup_hdl, file_hdl, 'wrestart') 
     
    543552           CALL xios_set_file_attr( "wrestart", type="one_file", enabled=.TRUE.,&  
    544553                                    mode="write", output_freq=xios_timestep)  
    545            if(lwp .AND. nprint > 1) write(numout,*) 'OPEN ', trim(cdrst_file), ' in one_file mode'  
     554           IF(lwp .AND. nprint > 1) THEN 
     555              write(numout,*) 'OPEN ', trim(cdrst_file), ' in one_file mode'  
     556              IF(lflush) CALL FLUSH(numout) 
     557           ENDIF 
    546558        ELSE   
    547559           CALL xios_set_file_attr( "wrestart", type="multiple_file", enabled=.TRUE.,&  
    548560                                    mode="write", output_freq=xios_timestep)  
    549            if(lwp .AND. nprint > 1) write(numout,*) 'OPEN ', trim(cdrst_file), ' in multiple_file mode'  
     561           IF(lwp .AND. nprint > 1) THEN 
     562              write(numout,*) 'OPEN ', trim(cdrst_file), ' in multiple_file mode'  
     563              IF(lflush) CALL FLUSH(numout) 
     564           ENDIF 
    550565        ENDIF  
    551566        CALL xios_set_file_attr( "wrestart", name=trim(cdrst_file)) 
     
    824839               iom_file(jf)%nfid       = 0          ! free the id  
    825840               IF( PRESENT(kiomid) )   kiomid = 0   ! return 0 as id to specify that the file was closed 
    826                IF(lwp  .AND. nprint > 1) & 
     841               IF(lwp  .AND. nprint > 1) THEN 
    827842                 WRITE(numout,*) TRIM(clinfo)//' close file: '//TRIM(iom_file(jf)%name)//' ok' 
     843                 IF(lflush) CALL FLUSH(numout) 
     844               ENDIF 
    828845            ELSEIF( PRESENT(kiomid) ) THEN 
    829846               WRITE(ctmp1,*) '--->',  kiomid 
     
    942959      ELSE 
    943960#if defined key_iomput 
    944          IF(lwp .AND. nprint > 1) WRITE(numout,*) 'XIOS RST READ (0D): ', trim(cdvar) 
     961         IF(lwp .AND. nprint > 1) THEN 
     962            WRITE(numout,*) 'XIOS RST READ (0D): ', trim(cdvar) 
     963            IF(lflush) CALL FLUSH(numout) 
     964         ENDIF 
    945965         CALL iom_swap( TRIM(crxios_context) ) 
    946966         CALL xios_recv_field( trim(cdvar), pvar) 
     
    12671287 
    12681288            IF( istop == nstop ) THEN   ! no additional errors until this point... 
    1269                IF(lwp .AND. nprint > 1) & 
     1289               IF(lwp .AND. nprint > 1) THEN 
    12701290                  WRITE(numout,"(10x,' read ',a,' (rec: ',i6,') in ',a,' ok')") TRIM(cdvar), itime, TRIM(iom_file(kiomid)%name) 
     1291                  IF(lflush) CALL FLUSH(numout) 
     1292               ENDIF 
    12711293              
    12721294               !--- overlap areas and extra hallows (mpp) 
     
    12981320         IF( PRESENT(pv_r3d) ) THEN 
    12991321            pv_r3d(:, :, :) = 0. 
    1300             if(lwp .AND. nprint > 1) write(numout,*) 'XIOS RST READ (3D): ',trim(cdvar) 
     1322            IF(lwp .AND. nprint > 1) THEN 
     1323               write(numout,*) 'XIOS RST READ (3D): ',trim(cdvar) 
     1324               IF(lflush) CALL FLUSH(numout) 
     1325            ENDIF 
    13011326            CALL xios_recv_field( trim(cdvar), pv_r3d) 
    13021327            IF(idom /= jpdom_unknown ) then 
     
    13051330         ELSEIF( PRESENT(pv_r2d) ) THEN 
    13061331            pv_r2d(:, :) = 0. 
    1307             if(lwp .AND. nprint > 1) write(numout,*) 'XIOS RST READ (2D): ', trim(cdvar) 
     1332            IF(lwp .AND. nprint > 1) THEN 
     1333               write(numout,*) 'XIOS RST READ (2D): ', trim(cdvar) 
     1334               IF(lflush) CALL FLUSH(numout) 
     1335            ENDIF 
    13081336            CALL xios_recv_field( trim(cdvar), pv_r2d) 
    13091337            IF(idom /= jpdom_unknown ) THEN 
     
    13121340         ELSEIF( PRESENT(pv_r1d) ) THEN 
    13131341            pv_r1d(:) = 0. 
    1314             if(lwp .AND. nprint > 1) write(numout,*) 'XIOS RST READ (1D): ', trim(cdvar) 
     1342            IF(lwp .AND. nprint > 1) THEN 
     1343               write(numout,*) 'XIOS RST READ (1D): ', trim(cdvar) 
     1344               IF(lflush) CALL FLUSH(numout) 
     1345            ENDIF 
    13151346            CALL xios_recv_field( trim(cdvar), pv_r1d) 
    13161347         ENDIF 
     
    15141545#ifdef key_iomput 
    15151546      IF( kt == kwrite ) THEN 
    1516           IF(lwp .AND. nprint > 1) write(numout,*) 'RESTART: write (XIOS 0D) ',trim(cdvar) 
     1547          IF(lwp .AND. nprint > 1) THEN 
     1548             write(numout,*) 'RESTART: write (XIOS 0D) ',trim(cdvar) 
     1549             IF(lflush) CALL FLUSH(numout) 
     1550          ENDIF 
    15171551          CALL xios_send_field(trim(cdvar), pvar) 
    15181552      ENDIF 
     
    15441578#ifdef key_iomput 
    15451579      IF( kt == kwrite ) THEN 
    1546          IF(lwp .AND. nprint > 1) write(numout,*) 'RESTART: write (XIOS 1D) ',trim(cdvar) 
     1580         IF(lwp .AND. nprint > 1) THEN 
     1581            write(numout,*) 'RESTART: write (XIOS 1D) ',trim(cdvar) 
     1582            IF(lflush) CALL FLUSH(numout) 
     1583         ENDIF 
    15471584         CALL xios_send_field(trim(cdvar), pvar) 
    15481585      ENDIF 
     
    15741611#ifdef key_iomput 
    15751612      IF( kt == kwrite ) THEN 
    1576          IF(lwp .AND. nprint > 1) write(numout,*) 'RESTART: write (XIOS 2D) ',trim(cdvar) 
     1613         IF(lwp .AND. nprint > 1) THEN 
     1614            write(numout,*) 'RESTART: write (XIOS 2D) ',trim(cdvar) 
     1615            IF(lflush) CALL FLUSH(numout) 
     1616         ENDIF 
    15771617         CALL xios_send_field(trim(cdvar), pvar) 
    15781618      ENDIF 
     
    16041644#ifdef key_iomput 
    16051645      IF( kt == kwrite ) THEN 
    1606          IF(lwp .AND. nprint > 1) write(numout,*) 'RESTART: write (XIOS 3D) ',trim(cdvar) 
     1646         IF(lwp .AND. nprint > 1) THEN 
     1647            write(numout,*) 'RESTART: write (XIOS 3D) ',trim(cdvar) 
     1648            IF(lflush) CALL FLUSH(numout) 
     1649         ENDIF 
    16071650         CALL xios_send_field(trim(cdvar), pvar) 
    16081651      ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/IOM/iom_nf90.F90

    r10968 r10986  
    8888         !                                       !=========================! 
    8989         IF( ldwrt ) THEN  ! ... in write mode 
    90             IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in WRITE mode' 
     90            IF(lwp .AND. nprint > 0) THEN 
     91               WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in WRITE mode' 
     92               IF(lflush) CALL FLUSH(numout) 
     93            ENDIF 
    9194            IF( snc4set%luse ) THEN 
    9295               CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_WRITE  , if90id ), clinfo) 
     
    9699            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy                          ), clinfo) 
    97100         ELSE              ! ... in read mode 
    98             IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in READ mode' 
     101            IF(lwp .AND. nprint > 0) THEN 
     102               WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in READ mode' 
     103               IF(lflush) CALL FLUSH(numout) 
     104            ENDIF 
    99105            CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_NOWRITE, if90id, chunksize = ichunk ), clinfo) 
    100106         ENDIF 
     
    107113               cdname = TRIM(cltmp) 
    108114            ENDIF 
    109             IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' create new file: '//TRIM(cdname)//' in WRITE mode' 
     115            IF(lwp .AND. nprint > 0) THEN 
     116               WRITE(numout,*) TRIM(clinfo)//' create new file: '//TRIM(cdname)//' in WRITE mode' 
     117               IF(lflush) CALL FLUSH(numout) 
     118            ENDIF 
    110119 
    111120            IF( llclobber ) THEN   ;   imode = IOR( NF90_64BIT_OFFSET, NF90_CLOBBER   ) 
     
    113122            ENDIF 
    114123            IF( snc4set%luse ) THEN 
    115                IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' creating file: '//TRIM(cdname)//' in hdf5 (netcdf4) mode' 
     124               IF(lwp .AND. nprint > 0) THEN 
     125                  WRITE(numout,*) TRIM(clinfo)//' creating file: '//TRIM(cdname)//' in hdf5 (netcdf4) mode' 
     126                  IF(lflush) CALL FLUSH(numout) 
     127               ENDIF 
    116128               CALL GET_NF90_SYMBOL("NF90_HDF5", ihdf5) 
    117129               IF( llclobber ) THEN   ;   imode = IOR(ihdf5, NF90_CLOBBER) 
     
    165177               &                                       len  = iom_file(kiomid)%lenuld ), clinfo ) 
    166178         ENDIF 
    167          IF(lwp .AND. nprint > 0) WRITE(numout,*) '                   ---> '//TRIM(cdname)//' OK' 
     179         IF(lwp .AND. nprint > 0) THEN 
     180            WRITE(numout,*) '                   ---> '//TRIM(cdname)//' OK' 
     181            IF(lflush) CALL FLUSH(numout) 
     182         ENDIF 
    168183      ELSE 
    169184         kiomid = 0               ! return error flag 
     
    562577         ! trick: defined to 0 to say that dimension variables are defined but not yet written 
    563578         iom_file(kiomid)%dimsz(1, 1)  = 0    
    564          IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' define dimension variables done' 
     579         IF(lwp .AND. nprint > 1) THEN 
     580            WRITE(numout,*) TRIM(clinfo)//' define dimension variables done' 
     581            IF(lflush) CALL FLUSH(numout) 
     582          ENDIF 
    565583      ENDIF 
    566584      ! define the data if it is not already done 
     
    637655            CALL iom_nf90_check(SET_NF90_DEF_VAR_CHUNKING(if90id, idvar, ichunkalg, ichunksz), clinfo) 
    638656            CALL iom_nf90_check(SET_NF90_DEF_VAR_DEFLATE(if90id, idvar, ishuffle, ideflate, ideflate_level), clinfo) 
    639             IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' chunked ok. Chunks sizes: ', ichunksz 
    640          ENDIF 
    641          IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' defined ok' 
     657            IF(lwp .AND. nprint > 1) THEN 
     658               WRITE(numout,*) TRIM(clinfo)//' chunked ok. Chunks sizes: ', ichunksz 
     659               IF(lflush) CALL FLUSH(numout) 
     660            ENDIF 
     661         ENDIF 
     662         IF(lwp .AND. nprint > 1) THEN 
     663            WRITE(numout,*) TRIM(clinfo)//' defined ok' 
     664            IF(lflush) CALL FLUSH(numout) 
     665          ENDIF 
    642666      ELSE 
    643667         idvar = kvid 
     
    685709               CALL iom_nf90_check( NF90_INQUIRE_DIMENSION( if90id, 3, len = iom_file(kiomid)%dimsz(1,3) ), clinfo ) 
    686710               iom_file(kiomid)%dimsz(1  , 4) = 1   ! unlimited dimension 
    687                IF(lwp .AND. nprint > 1) WRITE(numout,*) TRIM(clinfo)//' write dimension variables done' 
     711               IF(lwp .AND. nprint > 1) THEN 
     712                  WRITE(numout,*) TRIM(clinfo)//' write dimension variables done' 
     713                  IF(lflush) CALL FLUSH(numout) 
     714               ENDIF 
    688715            ENDIF 
    689716         ENDIF 
     
    703730         IF( iom_file(kiomid)%luld(idvar) )   iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar)    & 
    704731               &                            = iom_file(kiomid)%dimsz(iom_file(kiomid)%ndims(idvar), idvar) + 1 
    705          IF(lwp .AND. nprint > 0) WRITE(numout,*) TRIM(clinfo)//' written ok' 
     732         IF(lwp .AND. nprint > 0) THEN 
     733            WRITE(numout,*) TRIM(clinfo)//' written ok' 
     734            IF(lflush) CALL FLUSH(numout) 
     735         ENDIF 
    706736      ENDIF 
    707737      !      
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/IOM/prtctl.F90

    r10888 r10986  
    493493         zidom = zidom + ilcitl(ji,1) - nrecil 
    494494      END DO 
    495       IF(lwp) WRITE(numout,*) 
    496       IF(lwp) WRITE(numout,*)' sum ilcitl(i,1) = ', zidom, ' jpiglo = ', jpiglo 
     495      IF(lwp) THEN  
     496         WRITE(numout,*) 
     497         WRITE(numout,*)' sum ilcitl(i,1) = ', zidom, ' jpiglo = ', jpiglo 
     498         IF(lflush) CALL FLUSH(numout) 
     499      ENDIF 
    497500       
    498501      zjdom = nrecjl 
     
    500503         zjdom = zjdom + ilcjtl(1,jj) - nrecjl 
    501504      END DO 
    502       IF(lwp) WRITE(numout,*)' sum ilcitl(1,j) = ', zjdom, ' jpjglo = ', jpjglo 
    503       IF(lwp) WRITE(numout,*) 
     505      IF(lwp) THEN 
     506         WRITE(numout,*)' sum ilcitl(1,j) = ', zjdom, ' jpjglo = ', jpjglo 
     507         WRITE(numout,*) 
     508         IF(lflush) CALL FLUSH(numout) 
     509      ENDIF 
    504510       
    505511 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/IOM/restart.F90

    r10968 r10986  
    104104                  ENDIF 
    105105               ENDIF 
     106               IF(lflush) CALL FLUSH(numout) 
    106107            ENDIF 
    107108            ! 
     
    211212            IF ( snc4set%luse )      WRITE(numout,*) 'rst_read : configured with NetCDF4 support' 
    212213            WRITE(numout,*) '~~~~~~~~' 
     214            IF(lflush) CALL FLUSH(numout) 
    213215         ENDIF 
    214216         lxios_sini = .FALSE. 
     
    223225             crxios_context = 'nemo_rst' 
    224226             IF( .NOT.lxios_set ) THEN 
    225                  IF(lwp .AND. nprint > 0) WRITE(numout,*) 'Enable restart reading by XIOS' 
     227                 IF(lwp .AND. nprint > 0) THEN 
     228                    WRITE(numout,*) 'Enable restart reading by XIOS' 
     229                    IF(lflush) CALL FLUSH(numout) 
     230                 ENDIF 
    226231                 CALL iom_init( crxios_context, ld_tmppatch = .false. ) 
    227232                 lxios_set = .TRUE. 
     
    230235         IF( TRIM(Agrif_CFixed()) /= '0' .AND. lrxios) THEN 
    231236             CALL iom_init( crxios_context, ld_tmppatch = .false. ) 
    232              IF(lwp .AND. nprint > 0) WRITE(numout,*) 'Enable restart reading by XIOS for AGRIF' 
     237             IF(lwp .AND. nprint > 0) THEN 
     238                WRITE(numout,*) 'Enable restart reading by XIOS for AGRIF' 
     239                IF(lflush) CALL FLUSH(numout) 
     240             ENDIF 
    233241             lxios_set = .TRUE. 
    234242         ENDIF  
     
    264272      IF( ln_diurnal ) CALL iom_get( numror, jpdom_autoglo, 'Dsst' , x_dsst, ldxios = lrxios )  
    265273      IF ( ln_diurnal_only ) THEN  
    266          IF(lwp) WRITE( numout, * ) & 
    267          &   "rst_read:- ln_diurnal_only set, setting rhop=rau0"  
     274         IF(lwp) THEN 
     275            WRITE( numout, * ) "rst_read:- ln_diurnal_only set, setting rhop=rau0"  
     276            IF(lflush) CALL FLUSH(numout) 
     277         ENDIF 
    268278         rhop = rau0 
    269279         CALL iom_get( numror, jpdom_autoglo, 'tn'     , w3d, ldxios = lrxios )  
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LBC/lib_mpp.F90

    r10968 r10986  
    617617         !                                       -------------------------- 
    618618         IF ( SIZE(todelay(idvar)%z1d) /= isz ) THEN                  ! Check dimension coherence 
    619             IF(lwp) WRITE(numout,*) ' WARNING: the nb of delayed variables in restart file is not the model one' 
     619            IF(lwp) THEN 
     620               WRITE(numout,*) ' WARNING: the nb of delayed variables in restart file is not the model one' 
     621               IF(lflush) CALL FLUSH(numout) 
     622            ENDIF 
    620623            DEALLOCATE(todelay(idvar)%z1d) 
    621624            ndelayid(idvar) = -1                                      ! do as if we had no restart 
     
    684687         !                                       -------------------------- 
    685688         IF ( SIZE(todelay(idvar)%z1d) /= isz ) THEN                  ! Check dimension coherence 
    686             IF(lwp) WRITE(numout,*) ' WARNING: the nb of delayed variables in restart file is not the model one' 
     689            IF(lwp) THEN 
     690               WRITE(numout,*) ' WARNING: the nb of delayed variables in restart file is not the model one' 
     691               IF(lflush) CALL FLUSH(numout) 
     692            ENDIF 
    687693            DEALLOCATE(todelay(idvar)%z1d) 
    688694            ndelayid(idvar) = -1                                      ! do as if we had no restart 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LBC/mppini.F90

    r10968 r10986  
    9494         WRITE(numout,*) '   l_Iperio = ', l_Iperio, '    l_Jperio = ', l_Jperio  
    9595         WRITE(numout,*) '     npolj  = ',   npolj , '      njmpp  = ', njmpp 
     96         IF(lflush) CALL FLUSH(numout) 
    9697      ENDIF 
    9798      ! 
     
    189190         WRITE(numout,*) '~~~~~~~~ ' 
    190191         WRITE(numout,*) 
     192         IF(lflush) CALL FLUSH(numout) 
    191193      ENDIF 
    192194      IF( jpni < 1 .OR. jpnj < 1 ) THEN 
     
    211213            WRITE(numout,9000) '      - has mpi subdomains with a maximum size of (jpi = ',  iimax, ', jpj = ',  ijmax,   & 
    212214               &                                                             ', jpi*jpj = ',  iimax* ijmax, ')' 
     215            IF(lflush) CALL FLUSH(numout) 
    213216         ENDIF 
    214217         IF( iimax*ijmax < jpimax*jpjmax ) THEN   ! chosen subdomain size is larger that the best subdomain size 
     
    264267            ENDIF 
    265268            WRITE(numout,*) 
     269            IF(lflush) CALL FLUSH(numout) 
    266270         ENDIF 
    267271         CALL mpp_init_bestpartition( mppsize, ldlist = .TRUE. )   ! must be done by all core 
     
    284288            WRITE(numout,9003) '   Number of suppressed land subdomains = ', jpni*jpnj - inijmin 
    285289            WRITE(numout,*) 
     290            IF(lflush) CALL FLUSH(numout) 
    286291         ENDIF 
    287292      ENDIF 
     
    337342         WRITE(numout,*) '      sum ilci(i,1) = ', sum(ilci(:,1)), ' jpiglo = ', jpiglo 
    338343         WRITE(numout,*) '      sum ilcj(1,j) = ', sum(ilcj(1,:)), ' jpjglo = ', jpjglo 
     344         IF(lflush) CALL FLUSH(numout) 
    339345      ENDIF 
    340346      
     
    507513            il1 = il1+ifreq 
    508514         END DO 
     515         IF(lflush) CALL FLUSH(numout) 
     516      ENDIF 
    509517 9400    FORMAT('           ***'   ,20('*************',a3)    ) 
    510518 9403    FORMAT('           *     ',20('         *   ',a3)    ) 
     
    512520 9402    FORMAT('       ',i3,' *  ',20(i3,'  x',i3,'   *   ') ) 
    513521 9404    FORMAT('           *  '   ,20('      ',i3,'   *   ') ) 
    514       ENDIF 
    515522          
    516523      ! just to save nono etc for all proc 
     
    647654         WRITE(numout,*) '      nrecj  = ', nrecj   
    648655         WRITE(numout,*) '      nn_hls = ', nn_hls  
     656         IF(lflush) CALL FLUSH(numout) 
    649657      ENDIF 
    650658 
     
    655663            WRITE(numout,*) 
    656664            WRITE(numout,*) '   ==>>>   North fold boundary prepared for jpni >1' 
     665            IF(lflush) CALL FLUSH(numout) 
    657666            ! additional prints in layout.dat 
    658667         ENDIF 
     
    948957            WRITE(numout,         *) '  --------------------------------------', '-----', '--------------' 
    949958            WRITE(numout,*) 
     959            IF(lflush) CALL FLUSH(numout) 
    950960         END IF 
    951961         iitarget = MINLOC( inbi0(:)*inbj0(:), mask = inbi0(:)*inbj0(:) >= knbij, dim = 1 ) 
     
    955965            inbij = COUNT(llisoce) 
    956966            DEALLOCATE( llisoce ) 
    957             IF(lwp) WRITE(numout,'(a, i5, a, i5, a, i4, a, i4, a, i9, a, i5, a, i5, a)')    & 
     967            IF(lwp) THEN 
     968               WRITE(numout,'(a, i5, a, i5, a, i4, a, i4, a, i9, a, i5, a, i5, a)')    & 
    958969               &     'nb_cores ' , inbij,' oce + ', inbi0(ji)*inbj0(ji) - inbij             & 
    959970               &                                , ' land ( ', inbi0(ji),' x ', inbj0(ji),   & 
    960971               & ' ), nb_points ', iszi0(ji)*iszj0(ji),' ( ', iszi0(ji),' x ', iszj0(ji),' )' 
     972               IF(lflush) CALL FLUSH(numout) 
     973            ENDIF 
    961974         END DO 
    962975         DEALLOCATE( inbi0, inbj0, iszi0, iszj0 ) 
     
    11711184          WRITE(numout,*) '                    ihals = ', ihals(1), ihals(2) 
    11721185          WRITE(numout,*) '                    ihale = ', ihale(1), ihale(2) 
     1186          IF(lflush) CALL FLUSH(numout) 
    11731187      ENDIF 
    11741188      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LDF/ldfc1d_c2d.F90

    r10888 r10986  
    6161      !!---------------------------------------------------------------------- 
    6262      ! 
    63       IF(lwp) WRITE(numout,*) 
    64       IF(lwp) WRITE(numout,*) '   ldf_c1d : set a given profile to eddy mixing coefficients' 
     63      IF(lwp) THEN 
     64         WRITE(numout,*) 
     65         WRITE(numout,*) '   ldf_c1d : set a given profile to eddy mixing coefficients' 
     66         IF(lflush) CALL FLUSH(numout) 
     67      ENDIF 
    6568      ! 
    6669      ! initialization of the profile 
     
    134137      !!---------------------------------------------------------------------- 
    135138      ! 
    136       IF(lwp) WRITE(numout,*) 
    137       IF(lwp) WRITE(numout,*) '   ldf_c2d :   aht = Ufac * max(e1,e2)   with Ufac = ', pUfac, ' m/s' 
     139      IF(lwp) THEN 
     140         WRITE(numout,*) 
     141         WRITE(numout,*) '   ldf_c2d :   aht = Ufac * max(e1,e2)   with Ufac = ', pUfac, ' m/s' 
     142         IF(lflush) CALL FLUSH(numout) 
     143      ENDIF 
    138144      ! 
    139145      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LDF/ldfdyn.F90

    r10968 r10986  
    151151         WRITE(numout,*) '            lower limit (default 1.0)         rn_minfac    = ', rn_minfac 
    152152         WRITE(numout,*) '            upper limit (default 1.0)         rn_maxfac    = ', rn_maxfac 
     153         IF(lflush) CALL FLUSH(numout) 
    153154      ENDIF 
    154155 
     
    224225         END SELECT 
    225226         WRITE(numout,*) 
     227         IF(lflush) CALL FLUSH(numout) 
    226228      ENDIF 
    227229       
     
    233235      ! 
    234236      IF( ln_dynldf_OFF ) THEN 
    235          IF(lwp) WRITE(numout,*) '   ==>>>   No viscous operator selected. ahmt and ahmf are not allocated' 
     237         IF(lwp) THEN 
     238            WRITE(numout,*) '   ==>>>   No viscous operator selected. ahmt and ahmf are not allocated' 
     239            IF(lflush) CALL FLUSH(numout) 
     240         ENDIF 
    236241         RETURN 
    237242         ! 
     
    255260         ! 
    256261         CASE(   0  )      !==  constant  ==! 
    257             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity. = constant = ', zah0, cl_Units 
     262            IF(lwp) THEN 
     263               WRITE(numout,*) '   ==>>>   eddy viscosity. = constant = ', zah0, cl_Units 
     264               IF(lflush) CALL FLUSH(numout) 
     265            ENDIF 
    258266            ahmt(:,:,1:jpkm1) = zah0 
    259267            ahmf(:,:,1:jpkm1) = zah0 
    260268            ! 
    261269         CASE(  10  )      !==  fixed profile  ==! 
    262             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F( depth )' 
    263             IF(lwp) WRITE(numout,*) '           surface viscous coef. = constant = ', zah0, cl_Units 
     270            IF(lwp) THEN 
     271               WRITE(numout,*) '   ==>>>   eddy viscosity = F( depth )' 
     272               WRITE(numout,*) '           surface viscous coef. = constant = ', zah0, cl_Units 
     273               IF(lflush) CALL FLUSH(numout) 
     274            ENDIF 
    264275            ahmt(:,:,1) = zah0                        ! constant surface value 
    265276            ahmf(:,:,1) = zah0 
     
    267278            ! 
    268279         CASE ( -20 )      !== fixed horizontal shape read in file  ==! 
    269             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F(i,j) read in eddy_viscosity.nc file' 
     280            IF(lwp) THEN 
     281               WRITE(numout,*) '   ==>>>   eddy viscosity = F(i,j) read in eddy_viscosity.nc file' 
     282               IF(lflush) CALL FLUSH(numout) 
     283            ENDIF 
    270284            CALL iom_open( 'eddy_viscosity_2D.nc', inum ) 
    271285            CALL iom_get ( inum, jpdom_data, 'ahmt_2d', ahmt(:,:,1) ) 
     
    278292            ! 
    279293         CASE(  20  )      !== fixed horizontal shape  ==! 
    280             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)' 
    281             IF(lwp) WRITE(numout,*) '           using a fixed viscous velocity = ', rn_Uv  ,' m/s   and   Lv = Max(e1,e2)' 
    282             IF(lwp) WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     294            IF(lwp) THEN 
     295               WRITE(numout,*) '   ==>>>   eddy viscosity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)' 
     296               WRITE(numout,*) '           using a fixed viscous velocity = ', rn_Uv  ,' m/s   and   Lv = Max(e1,e2)' 
     297               WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     298               IF(lflush) CALL FLUSH(numout) 
     299            ENDIF 
    283300            CALL ldf_c2d( 'DYN', zUfac      , inn        , ahmt, ahmf )         ! surface value proportional to scale factor^inn 
    284301            ! 
    285302         CASE( -30  )      !== fixed 3D shape read in file  ==! 
    286             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F(i,j,k) read in eddy_viscosity_3D.nc file' 
     303            IF(lwp) THEN 
     304               WRITE(numout,*) '   ==>>>   eddy viscosity = F(i,j,k) read in eddy_viscosity_3D.nc file' 
     305               IF(lflush) CALL FLUSH(numout) 
     306            ENDIF 
    287307            CALL iom_open( 'eddy_viscosity_3D.nc', inum ) 
    288308            CALL iom_get ( inum, jpdom_data, 'ahmt_3d', ahmt ) 
     
    291311            ! 
    292312         CASE(  30  )       !==  fixed 3D shape  ==! 
    293             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth )' 
    294             IF(lwp) WRITE(numout,*) '           using a fixed viscous velocity = ', rn_Uv  ,' m/s   and   Ld = Max(e1,e2)' 
    295             IF(lwp) WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     313            IF(lwp) THEN 
     314               WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth )' 
     315               WRITE(numout,*) '           using a fixed viscous velocity = ', rn_Uv  ,' m/s   and   Ld = Max(e1,e2)' 
     316               WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     317               IF(lflush) CALL FLUSH(numout) 
     318            ENDIF 
    296319            CALL ldf_c2d( 'DYN', zUfac      , inn        , ahmt, ahmf )         ! surface value proportional to scale factor^inn 
    297320            CALL ldf_c1d( 'DYN', ahmt(:,:,1), ahmf(:,:,1), ahmt, ahmf )  ! reduction with depth 
    298321            ! 
    299322         CASE(  31  )       !==  time varying 3D field  ==! 
    300             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth , time )' 
    301             IF(lwp) WRITE(numout,*) '           proportional to the local velocity : 1/2 |u|e (lap) or 1/12 |u|e^3 (blp)' 
     323            IF(lwp) THEN 
     324               WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth , time )' 
     325               WRITE(numout,*) '           proportional to the local velocity : 1/2 |u|e (lap) or 1/12 |u|e^3 (blp)' 
     326               IF(lflush) CALL FLUSH(numout) 
     327            ENDIF 
    302328            ! 
    303329            l_ldfdyn_time = .TRUE.     ! will be calculated by call to ldf_dyn routine in step.F90 
    304330            ! 
    305331         CASE(  32  )       !==  time varying 3D field  ==! 
    306             IF(lwp) WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth , time )' 
    307             IF(lwp) WRITE(numout,*) '           proportional to the local deformation rate and gridscale (Smagorinsky)' 
     332            IF(lwp) THEN 
     333               WRITE(numout,*) '   ==>>>   eddy viscosity = F( latitude, longitude, depth , time )' 
     334               WRITE(numout,*) '           proportional to the local deformation rate and gridscale (Smagorinsky)' 
     335               IF(lflush) CALL FLUSH(numout) 
     336            ENDIF 
    308337            ! 
    309338            l_ldfdyn_time = .TRUE.     ! will be calculated by call to ldf_dyn routine in step.F90 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LDF/ldfslp.F90

    r10888 r10986  
    746746         WRITE(numout,*) 'ldf_slp_init : direction of lateral mixing' 
    747747         WRITE(numout,*) '~~~~~~~~~~~~' 
     748         IF(lflush) CALL FLUSH(numout) 
    748749      ENDIF 
    749750      ! 
     
    752753      ! 
    753754      IF( ln_traldf_triad ) THEN        ! Griffies operator : triad of slopes 
    754          IF(lwp) WRITE(numout,*) '   ==>>>   triad) operator (Griffies)' 
     755         IF(lwp) THEN 
     756            WRITE(numout,*) '   ==>>>   triad) operator (Griffies)' 
     757            IF(lflush) CALL FLUSH(numout) 
     758         ENDIF 
    755759         ALLOCATE( triadi_g(jpi,jpj,jpk,0:1,0:1) , triadj_g(jpi,jpj,jpk,0:1,0:1) ,     & 
    756760            &      triadi  (jpi,jpj,jpk,0:1,0:1) , triadj  (jpi,jpj,jpk,0:1,0:1) ,     & 
     
    760764         ! 
    761765      ELSE                             ! Madec operator : slopes at u-, v-, and w-points 
    762          IF(lwp) WRITE(numout,*) '   ==>>>   iso operator (Madec)' 
     766         IF(lwp) THEN 
     767            WRITE(numout,*) '   ==>>>   iso operator (Madec)' 
     768            IF(lflush) CALL FLUSH(numout) 
     769         ENDIF 
    763770         ALLOCATE( omlmask(jpi,jpj,jpk) ,                                                                        & 
    764771            &      uslp(jpi,jpj,jpk) , uslpml(jpi,jpj) , wslpi(jpi,jpj,jpk) , wslpiml(jpi,jpj) ,     & 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LDF/ldftra.F90

    r10968 r10986  
    146146         WRITE(numout,*) 'ldf_tra_init : lateral tracer diffusion' 
    147147         WRITE(numout,*) '~~~~~~~~~~~~ ' 
     148         IF(lflush) CALL FLUSH(numout) 
    148149      ENDIF 
    149150       
     
    180181         WRITE(numout,*) '            lateral diffusive velocity (if cst)  rn_Ud           = ', rn_Ud, ' m/s' 
    181182         WRITE(numout,*) '            lateral diffusive length   (if cst)  rn_Ld           = ', rn_Ld, ' m' 
     183         IF(lflush) CALL FLUSH(numout) 
    182184      ENDIF 
    183185      ! 
     
    271273         END SELECT 
    272274         WRITE(numout,*) 
     275         IF(lflush) CALL FLUSH(numout) 
    273276      ENDIF 
    274277 
     
    281284      IF( ln_traldf_OFF ) THEN               !== no explicit diffusive operator  ==! 
    282285         ! 
    283          IF(lwp) WRITE(numout,*) '   ==>>>   No diffusive operator selected. ahtu and ahtv are not allocated' 
     286         IF(lwp) THEN 
     287            WRITE(numout,*) '   ==>>>   No diffusive operator selected. ahtu and ahtv are not allocated' 
     288            IF(lflush) CALL FLUSH(numout) 
     289         ENDIF 
    284290         RETURN 
    285291         ! 
     
    304310         ! 
    305311         CASE(   0  )      !==  constant  ==! 
    306             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = constant = ', aht0, cl_Units 
     312            IF(lwp) THEN 
     313               WRITE(numout,*) '   ==>>>   eddy diffusivity = constant = ', aht0, cl_Units 
     314               IF(lflush) CALL FLUSH(numout) 
     315            ENDIF 
    307316            ahtu(:,:,1:jpkm1) = aht0 
    308317            ahtv(:,:,1:jpkm1) = aht0 
    309318            ! 
    310319         CASE(  10  )      !==  fixed profile  ==! 
    311             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F( depth )' 
    312             IF(lwp) WRITE(numout,*) '           surface eddy diffusivity = constant = ', aht0, cl_Units 
     320            IF(lwp) THEN 
     321               WRITE(numout,*) '   ==>>>   eddy diffusivity = F( depth )' 
     322               WRITE(numout,*) '           surface eddy diffusivity = constant = ', aht0, cl_Units 
     323               IF(lflush) CALL FLUSH(numout) 
     324            ENDIF 
    313325            ahtu(:,:,1) = aht0                        ! constant surface value 
    314326            ahtv(:,:,1) = aht0 
     
    316328            ! 
    317329         CASE ( -20 )      !== fixed horizontal shape and magnitude read in file  ==! 
    318             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F(i,j) read in eddy_diffusivity.nc file' 
     330            IF(lwp) THEN 
     331               WRITE(numout,*) '   ==>>>   eddy diffusivity = F(i,j) read in eddy_diffusivity.nc file' 
     332               IF(lflush) CALL FLUSH(numout) 
     333            ENDIF 
    319334            CALL iom_open( 'eddy_diffusivity_2D.nc', inum ) 
    320335            CALL iom_get ( inum, jpdom_data, 'ahtu_2D', ahtu(:,:,1) ) 
     
    327342            ! 
    328343         CASE(  20  )      !== fixed horizontal shape  ==! 
    329             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap or blp case)' 
    330             IF(lwp) WRITE(numout,*) '           using a fixed diffusive velocity = ', rn_Ud,' m/s   and   Ld = Max(e1,e2)' 
    331             IF(lwp) WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     344            IF(lwp) THEN 
     345               WRITE(numout,*) '   ==>>>   eddy diffusivity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap or blp case)' 
     346               WRITE(numout,*) '           using a fixed diffusive velocity = ', rn_Ud,' m/s   and   Ld = Max(e1,e2)' 
     347               WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     348               IF(lflush) CALL FLUSH(numout) 
     349            ENDIF 
    332350            CALL ldf_c2d( 'TRA', zUfac      , inn        , ahtu, ahtv )    ! value proportional to scale factor^inn 
    333351            ! 
    334352         CASE(  21  )      !==  time varying 2D field  ==! 
    335             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, time )' 
    336             IF(lwp) WRITE(numout,*) '                            = F( growth rate of baroclinic instability )' 
    337             IF(lwp) WRITE(numout,*) '            min value = 0.2 * aht0 (with aht0= 1/2 rn_Ud*rn_Ld)' 
    338             IF(lwp) WRITE(numout,*) '            max value =       aei0 (with aei0=1/2 rn_Ue*Le  increased to aht0 within 20N-20S' 
     353            IF(lwp) THEN 
     354               WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, time )' 
     355               WRITE(numout,*) '                            = F( growth rate of baroclinic instability )' 
     356               WRITE(numout,*) '            min value = 0.2 * aht0 (with aht0= 1/2 rn_Ud*rn_Ld)' 
     357               WRITE(numout,*) '            max value =       aei0 (with aei0=1/2 rn_Ue*Le  increased to aht0 within 20N-20S' 
     358               IF(lflush) CALL FLUSH(numout) 
     359            ENDIF 
    339360            ! 
    340361            l_ldftra_time = .TRUE.     ! will be calculated by call to ldf_tra routine in step.F90 
     
    344365            ! 
    345366         CASE( -30  )      !== fixed 3D shape read in file  ==! 
    346             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F(i,j,k) read in eddy_diffusivity.nc file' 
     367            IF(lwp) THEN 
     368               WRITE(numout,*) '   ==>>>   eddy diffusivity = F(i,j,k) read in eddy_diffusivity.nc file' 
     369               IF(lflush) CALL FLUSH(numout) 
     370            ENDIF 
    347371            CALL iom_open( 'eddy_diffusivity_3D.nc', inum ) 
    348372            CALL iom_get ( inum, jpdom_data, 'ahtu_3D', ahtu ) 
     
    351375            ! 
    352376         CASE(  30  )      !==  fixed 3D shape  ==! 
    353             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, depth )' 
    354             IF(lwp) WRITE(numout,*) '           using a fixed diffusive velocity = ', rn_Ud,' m/s   and   Ld = Max(e1,e2)' 
    355             IF(lwp) WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     377            IF(lwp) THEN 
     378               WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, depth )' 
     379               WRITE(numout,*) '           using a fixed diffusive velocity = ', rn_Ud,' m/s   and   Ld = Max(e1,e2)' 
     380               WRITE(numout,*) '           maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, '  for e1=1°)' 
     381               IF(lflush) CALL FLUSH(numout) 
     382            ENDIF 
    356383            CALL ldf_c2d( 'TRA', zUfac      , inn        , ahtu, ahtv )    ! surface value proportional to scale factor^inn 
    357384            CALL ldf_c1d( 'TRA', ahtu(:,:,1), ahtv(:,:,1), ahtu, ahtv )    ! reduction with depth 
    358385            ! 
    359386         CASE(  31  )      !==  time varying 3D field  ==! 
    360             IF(lwp) WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, depth , time )' 
    361             IF(lwp) WRITE(numout,*) '           proportional to the velocity : 1/2 |u|e or 1/12 |u|e^3' 
     387            IF(lwp) THEN 
     388               WRITE(numout,*) '   ==>>>   eddy diffusivity = F( latitude, longitude, depth , time )' 
     389               WRITE(numout,*) '           proportional to the velocity : 1/2 |u|e or 1/12 |u|e^3' 
     390               IF(lflush) CALL FLUSH(numout) 
     391            ENDIF 
    362392            ! 
    363393            l_ldftra_time = .TRUE.     ! will be calculated by call to ldf_tra routine in step.F90 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/OBS/obs_write.F90

    r10888 r10986  
    186186         WRITE(numout,*)'~~~~~~~~~~~~~' 
    187187         WRITE(numout,*)'Writing '//TRIM(clfiletype)//' feedback file : ',TRIM(clfname) 
     188         IF(lflush) CALL FLUSH(numout) 
    188189      ENDIF 
    189190 
     
    458459         WRITE(numout,*)'~~~~~~~~~~~~~' 
    459460         WRITE(numout,*)'Writing '//TRIM(surfdata%cvars(1))//' feedback file : ',TRIM(clfname) 
     461         IF(lflush) CALL FLUSH(numout) 
    460462      ENDIF 
    461463 
     
    571573         WRITE(numout,*) 'obs_wri_stats :' 
    572574         WRITE(numout,*) '~~~~~~~~~~~~~~~' 
     575         IF(lflush) CALL FLUSH(numout) 
    573576      ENDIF 
    574577 
     
    600603            WRITE(numout,*) 'Overall RMS obs minus model of the good observations: ',sqrt( zsumx2/inumgoodobsmpp ) 
    601604            WRITE(numout,*) '' 
     605            IF(lflush) CALL FLUSH(numout) 
    602606         ENDIF 
    603607 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/OBS/obsinter_h2d.h90

    r10888 r10986  
    281281                  WRITE(numout,*)'Switching to great circle distance' 
    282282                  WRITE(numout,*) 
     283                  IF(lflush) CALL FLUSH(numout) 
    283284               ENDIF 
    284285               CALL obs_int_h2d_ds1( kpk2, ikmax,                        & 
     
    11291130            WRITE(numout,*)' Iterations for i,j exceed max iteration count!' 
    11301131            WRITE(numout,*) 
     1132            IF(lflush) CALL FLUSH(numout) 
    11311133             
    11321134            ldfail = .TRUE. 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/cpl_oasis3.F90

    r10888 r10986  
    157157         IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 
    158158      ENDIF  
    159       IF(lwp) WRITE(numout,*) 
    160       IF(lwp) WRITE(numout,*) 'cpl_define : initialization in coupled ocean/atmosphere case' 
    161       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 
    162       IF(lwp) WRITE(numout,*) 
     159      IF(lwp) THEN 
     160         WRITE(numout,*) 
     161         WRITE(numout,*) 'cpl_define : initialization in coupled ocean/atmosphere case' 
     162         WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 
     163         WRITE(numout,*) 
     164         IF(lflush) CALL FLUSH(numout) 
     165      ENDIF 
    163166 
    164167      ncplmodel = kcplmodel 
     
    208211         WRITE(numout,*) ' multiexchg: nldi, nlei, nimpp =', nldi, nlei, nimpp 
    209212         WRITE(numout,*) ' multiexchg: nldj, nlej, njmpp =', nldj, nlej, njmpp 
     213         IF(lflush) CALL FLUSH(numout) 
    210214      ENDIF 
    211215    
     
    243247                  END IF 
    244248#endif 
    245                   IF( ln_ctl ) WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_Out 
     249                  IF( ln_ctl ) THEN  
     250                     WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_Out 
     251                     IF(lflush) CALL FLUSH(numout) 
     252                  ENDIF 
    246253                  CALL oasis_def_var (ssnd(ji)%nid(jc,jm), zclname, id_part   , (/ 2, 1 /),   & 
    247254                     &                OASIS_Out          , ishape , OASIS_REAL, nerror ) 
     
    288295                  END IF 
    289296#endif 
    290                   IF( ln_ctl ) WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_In 
     297                  IF( ln_ctl ) THEN 
     298                     WRITE(numout,*) "Define", ji, jc, jm, " "//TRIM(zclname), " for ", OASIS_In 
     299                     IF(lflush) CALL FLUSH(numout) 
     300                  ENDIF 
    291301                  CALL oasis_def_var (srcv(ji)%nid(jc,jm), zclname, id_part   , (/ 2, 1 /),   & 
    292302                     &                OASIS_In           , ishape , OASIS_REAL, nerror ) 
     
    361371                     WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
    362372                     WRITE(numout,*) '****************' 
     373                     IF(lflush) CALL FLUSH(numout) 
    363374                  ENDIF 
    364375               ENDIF 
     
    420431                  &        kinfo == OASIS_RecvOut .OR. kinfo == OASIS_FromRestOut 
    421432                
    422                IF ( ln_ctl )   WRITE(numout,*) "llaction, kinfo, kstep, ivarid: " , llaction, kinfo, kstep, srcv(kid)%nid(jc,jm) 
     433               IF ( ln_ctl )   THEN 
     434                  WRITE(numout,*) "llaction, kinfo, kstep, ivarid: " , llaction, kinfo, kstep, srcv(kid)%nid(jc,jm) 
     435                  IF(lflush) CALL FLUSH(numout) 
     436               ENDIF 
    423437                
    424438               IF ( llaction ) THEN 
     
    442456                     WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
    443457                     WRITE(numout,*) '****************' 
     458                     IF(lflush) CALL FLUSH(numout) 
    444459                  ENDIF 
    445460                   
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/fldread.F90

    r10968 r10986  
    304304                     & 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 
    305305                  WRITE(numout, *) '      it_offset is : ',it_offset 
     306                  IF(lflush) CALL FLUSH(numout) 
    306307               ENDIF 
    307308               ! temporal interpolation weights 
     
    315316                  WRITE(numout, clfmt) TRIM(sd(jf)%clvar), kt, REAL(isecsbc,wp)/rday, nyear, nmonth, nday,    & 
    316317                     &                 sd(jf)%nrec_a(1), REAL(sd(jf)%nrec_b(2),wp)/rday, REAL(sd(jf)%nrec_a(2),wp)/rday 
     318                  IF(lflush) CALL FLUSH(numout) 
    317319               ENDIF 
    318320            ENDIF 
     
    440442         ! 
    441443         clfmt = "('   fld_init : time-interpolation for ', a, ' read previous record = ', i6, ' at time = ', f7.2, ' days')" 
    442          IF(lwp .AND. nprint > 1) WRITE(numout, clfmt) TRIM(sdjf%clvar), sdjf%nrec_a(1), REAL(sdjf%nrec_a(2),wp)/rday 
     444         IF(lwp .AND. nprint > 1) THEN 
     445            WRITE(numout, clfmt) TRIM(sdjf%clvar), sdjf%nrec_a(1), REAL(sdjf%nrec_a(2),wp)/rday 
     446            IF(lflush) CALL FLUSH(numout) 
     447         ENDIF 
    443448         ! 
    444449      ENDIF 
     
    762767      ENDIF 
    763768 
    764       IF(lwp.AND. nprint > 1) WRITE(numout,*) 'Dim size for ',        TRIM(clvar),' is ', ilendta 
    765       IF(lwp.AND. nprint > 1) WRITE(numout,*) 'Number of levels for ',TRIM(clvar),' is ', ipk 
     769      IF(lwp.AND. nprint > 1) THEN 
     770         WRITE(numout,*) 'Dim size for ',        TRIM(clvar),' is ', ilendta 
     771         WRITE(numout,*) 'Number of levels for ',TRIM(clvar),' is ', ipk 
     772         IF(lflush) CALL FLUSH(numout) 
     773      ENDIF 
    766774      ! 
    767775      SELECT CASE( ipk ) 
     
    11741182                     END DO 
    11751183                     sd(ju)%rotn(jn) = .TRUE.               ! vector was rotated  
    1176                      IF( lwp .AND. kt == nit000 .AND. nprint > 1)   WRITE(numout,*)   & 
     1184                     IF( lwp .AND. kt == nit000 .AND. nprint > 1)   THEN 
     1185                         WRITE(numout,*)   & 
    11771186                        &   'fld_read: vector pair ('//TRIM(sd(ju)%clvar)//', '//TRIM(sd(iv)%clvar)//') rotated on to model grid' 
     1187                         IF(lflush) CALL FLUSH(numout) 
     1188                     ENDIF 
    11781189                  ENDIF 
    11791190               ENDIF 
     
    13401351            call flush(numout) 
    13411352         END DO 
     1353         IF(lflush) CALL FLUSH(numout) 
    13421354      ENDIF 
    13431355      ! 
     
    14111423         ENDIF 
    14121424         IF( ASSOCIATED(ref_wgts(kw)%data_wgt) )  WRITE(numout,*) '       allocated' 
     1425         IF(lflush) CALL FLUSH(numout) 
    14131426      END DO 
    14141427      ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/geo2ocean.F90

    r10888 r10986  
    6969      ! 
    7070      IF( lmust_init ) THEN      ! at 1st call only: set  gsin. & gcos. 
    71          IF(lwp) WRITE(numout,*) 
    72          IF(lwp) WRITE(numout,*) ' rot_rep: coordinate transformation : geographic <==> model (i,j)-components' 
    73          IF(lwp) WRITE(numout,*) ' ~~~~~~~~    ' 
     71         IF(lwp) THEN 
     72            WRITE(numout,*) 
     73            WRITE(numout,*) ' rot_rep: coordinate transformation : geographic <==> model (i,j)-components' 
     74            WRITE(numout,*) ' ~~~~~~~~    ' 
     75            IF(lflush) CALL FLUSH(numout) 
     76      ENDIF 
    7477         ! 
    7578         CALL angle( glamt, gphit, glamu, gphiu, glamv, gphiv, glamf, gphif )       ! initialization of the transformation 
     
    450453      ! ----------------------------------------------- 
    451454      IF( lmust_init ) THEN 
    452          IF(lwp) WRITE(numout,*) 
    453          IF(lwp) WRITE(numout,*) ' obs_rot : geographic <--> stretched' 
    454          IF(lwp) WRITE(numout,*) ' ~~~~~~~   coordinate transformation' 
     455         IF(lwp) THEN 
     456            WRITE(numout,*) 
     457            WRITE(numout,*) ' obs_rot : geographic <--> stretched' 
     458            WRITE(numout,*) ' ~~~~~~~   coordinate transformation' 
     459            IF(lflush) CALL FLUSH(numout) 
     460      ENDIF 
    455461         CALL angle( glamt, gphit, glamu, gphiu, glamv, gphiv, glamf, gphif )       ! initialization of the transformation 
    456462         lmust_init = .FALSE. 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcapr.F90

    r10968 r10986  
    9191         WRITE(numout,*) '   Namelist namsbc_apr : Atmospheric PRessure as extrenal forcing' 
    9292         WRITE(numout,*) '      ref. pressure: global mean Patm (T) or a constant (F)  ln_ref_apr = ', ln_ref_apr 
     93         IF(lflush) CALL FLUSH(numout) 
    9394      ENDIF 
    9495      ! 
     
    99100         IF(lwp) WRITE(numout,*) '         Reference Patm used : ', rn_pref, ' N/m2' 
    100101      ENDIF 
     102      IF(lwp .AND. lflush) CALL FLUSH(numout) 
    101103      ! 
    102104      r1_grau = 1.e0 / (grav * rau0)               !* constant for optimization 
    103105      ! 
    104106      !                                            !* control check 
    105       IF ( ln_apr_obc  ) THEN 
    106          IF(lwp) WRITE(numout,*) '         Inverse barometer added to OBC ssh data' 
     107      IF ( ln_apr_obc  .AND. lwp) THEN 
     108         WRITE(numout,*) '         Inverse barometer added to OBC ssh data' 
     109         IF(lflush) CALL FLUSH(numout) 
    107110      ENDIF 
    108111!jc: stop below should rather be a warning  
     
    155158         !                                            !* Restart: read in restart file 
    156159         IF( ln_rstart .AND. iom_varid( numror, 'ssh_ibb', ldstop = .FALSE. ) > 0 ) THEN  
    157             IF(lwp) WRITE(numout,*) 'sbc_apr:   ssh_ibb read in the restart file' 
     160            IF(lwp .AND. nprint > 0) THEN 
     161               WRITE(numout,*) 'sbc_apr:   ssh_ibb read in the restart file' 
     162               IF(lflush) CALL FLUSH(numout) 
     163            ENDIF 
    158164            CALL iom_get( numror, jpdom_autoglo, 'ssh_ibb', ssh_ibb, ldxios = lrxios )   ! before inv. barometer ssh 
    159165            ! 
    160166         ELSE                                         !* no restart: set from nit000 values 
    161             IF(lwp) WRITE(numout,*) 'sbc_apr:   ssh_ibb set to nit000 values' 
     167            IF(lwp .AND. nprint > 0) THEN 
     168               WRITE(numout,*) 'sbc_apr:   ssh_ibb set to nit000 values' 
     169               IF(lflush) CALL FLUSH(numout) 
     170            ENDIF 
    162171            ssh_ibb(:,:) = ssh_ib(:,:) 
    163172         ENDIF 
     
    166175      IF( lrst_oce ) THEN                       !      Write in the ocean restart file     ! 
    167176         !                                      ! ---------------------------------------- ! 
    168          IF(lwp) WRITE(numout,*) 
    169          IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 
    170          IF(lwp) WRITE(numout,*) '~~~~' 
     177         IF(lwp .AND. nprint > 0) THEN 
     178            WRITE(numout,*) 
     179            WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp 
     180            WRITE(numout,*) '~~~~' 
     181            IF(lflush) CALL FLUSH(numout) 
     182         ENDIF 
    171183         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    172184         CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib, ldxios = lwxios ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcblk.F90

    r10968 r10986  
    278278         CASE( np_ECMWF     )   ;   WRITE(numout,*) '   ==>>>   "ECMWF" algorithm       (IFS cycle 31)' 
    279279         END SELECT 
     280         IF(lflush) CALL FLUSH(numout) 
    280281         ! 
    281282      ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbccpl.F90

    r10968 r10986  
    330330         WRITE(numout,*)'  ln_usecplmask                       = ', ln_usecplmask 
    331331         WRITE(numout,*)'  nn_cats_cpl                         = ', nn_cats_cpl 
     332         IF(lflush) CALL FLUSH(numout) 
    332333      ENDIF 
    333334 
     
    466467         l_rnfcpl              = .TRUE.                      ! -> no need to read runoffs in sbcrnf 
    467468         ln_rnf                = nn_components /= jp_iam_sas ! -> force to go through sbcrnf if not sas 
    468          IF(lwp) WRITE(numout,*) 
    469          IF(lwp) WRITE(numout,*) '   runoffs received from oasis -> force ln_rnf = ', ln_rnf 
     469         IF(lwp) THEN 
     470            WRITE(numout,*) 
     471            WRITE(numout,*) '   runoffs received from oasis -> force ln_rnf = ', ln_rnf 
     472            IF(lflush) CALL FLUSH(numout) 
     473         ENDIF 
    470474      ENDIF 
    471475      ! 
     
    476480      IF( srcv(jpr_isf)%laction .AND. ln_isf ) THEN 
    477481         l_isfcpl             = .TRUE.                      ! -> no need to read isf in sbcisf 
    478          IF(lwp) WRITE(numout,*) 
    479          IF(lwp) WRITE(numout,*) '   iceshelf received from oasis ' 
     482         IF(lwp) THEN 
     483            WRITE(numout,*) 
     484            WRITE(numout,*) '   iceshelf received from oasis ' 
     485            IF(lflush) CALL FLUSH(numout) 
     486         ENDIF 
    480487      ENDIF 
    481488      ! 
     
    542549         srcv(jpr_co2 )%laction = .TRUE. 
    543550         l_co2cpl = .TRUE. 
    544          IF(lwp) WRITE(numout,*) 
    545          IF(lwp) WRITE(numout,*) '   Atmospheric pco2 received from oasis ' 
    546          IF(lwp) WRITE(numout,*) 
     551         IF(lwp) THEN 
     552            WRITE(numout,*) 
     553            WRITE(numout,*) '   Atmospheric pco2 received from oasis ' 
     554            WRITE(numout,*) 
     555            IF(lflush) CALL FLUSH(numout) 
     556         ENDIF 
    547557      ENDIF 
    548558      ! 
     
    666676            WRITE(numout,*)'                  wind stress module' 
    667677            WRITE(numout,*) 
     678            IF(lflush) CALL FLUSH(numout) 
    668679         ENDIF 
    669680      ENDIF 
     
    713724            WRITE(numout,*)'               fraction of solar net radiation absorbed in the first ocean level' 
    714725            WRITE(numout,*) 
     726            IF(lflush) CALL FLUSH(numout) 
    715727         ENDIF 
    716728      ENDIF 
     
    980992            WRITE(numout,*)'               fraction of solar net radiation absorbed in the first ocean level' 
    981993            WRITE(numout,*) 
     994            IF(lflush) CALL FLUSH(numout) 
    982995         ENDIF 
    983996      ENDIF 
     
    10201033            WRITE(numout,*)'                  wind stress U,V components' 
    10211034            WRITE(numout,*)'                  wind stress module' 
     1035            IF(lflush) CALL FLUSH(numout) 
    10221036         ENDIF 
    10231037      ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcdcy.F90

    r10888 r10986  
    9797            WRITE(numout,*) '~~~~~~~' 
    9898            WRITE(numout,*) 
     99            IF(lflush) CALL FLUSH(numout) 
    99100         ENDIF 
    100101         ! allocate sbcdcy arrays 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcflx.F90

    r10968 r10986  
    169169               WRITE(numout,*) ' day: ', ndastp , TRIM(sf(jf)%clvar), ' * ', zfact 
    170170            END DO 
     171            IF(lflush) CALL FLUSH(numout) 
    171172         ENDIF 
    172173         ! 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcfwb.F90

    r10968 r10986  
    8383            IF( kn_fwb == 2 )   WRITE(numout,*) '          adjusted from previous year budget' 
    8484            IF( kn_fwb == 3 )   WRITE(numout,*) '          fwf set to zero and spread out over erp area' 
     85            IF(lflush) CALL FLUSH(numout) 
    8586         ENDIF 
    8687         ! 
     
    121122            fwfold = a_fwb                            ! current year freshwater budget correction 
    122123            !                                         ! estimate from the previous year budget 
    123             IF(lwp)WRITE(numout,*) 
    124             IF(lwp)WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', fwfold 
    125             IF(lwp)WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
    126             IF(lwp)WRITE(numout,*)'          year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b 
     124            IF(lwp) THEN 
     125               WRITE(numout,*) 
     126               WRITE(numout,*)'sbc_fwb : year = ',iyear  , ' freshwater budget correction = ', fwfold 
     127               WRITE(numout,*)'          year = ',iyear-1, ' freshwater budget read       = ', a_fwb 
     128               WRITE(numout,*)'          year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b 
     129               IF(lflush) CALL FLUSH(numout) 
     130             ENDIF 
    127131         ENDIF    
    128132         !                                         ! Update fwfold if new year start 
     
    200204               WRITE(numout,*)'   MIN(zerp_cor) = ', MINVAL(zerp_cor)  
    201205               WRITE(numout,*)'   MAX(zerp_cor) = ', MAXVAL(zerp_cor)  
     206               IF(lflush) CALL FLUSH(numout) 
    202207            ENDIF 
    203208         ENDIF 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcisf.F90

    r10968 r10986  
    206206         IF( ln_rstart .AND.    &                     ! Restart: read in restart file 
    207207            &   iom_varid( numror, 'fwf_isf_b', ldstop = .FALSE. ) > 0 ) THEN 
    208             IF(lwp .AND. nprint > 0) WRITE(numout,*) '          nit000-1 isf tracer content forcing fields read in the restart file' 
     208            IF(lwp .AND. nprint > 0) THEN 
     209               WRITE(numout,*) '          nit000-1 isf tracer content forcing fields read in the restart file' 
     210               IF(lflush) CALL FLUSH(numout) 
     211            ENDIF 
    209212            CALL iom_get( numror, jpdom_autoglo, 'fwf_isf_b', fwfisf_b(:,:)         , ldxios = lrxios )   ! before salt content isf_tsc trend 
    210213            CALL iom_get( numror, jpdom_autoglo, 'isf_sc_b' , risf_tsc_b(:,:,jp_sal), ldxios = lrxios )   ! before salt content isf_tsc trend 
     
    217220      !  
    218221      IF( lrst_oce ) THEN 
    219          IF(lwp .AND. nprint > 0) WRITE(numout,*) 
    220          IF(lwp .AND. nprint > 0) WRITE(numout,*) 'sbc : isf surface tracer content forcing fields written in ocean restart file ',   & 
     222         IF(lwp .AND. nprint > 0) THEN 
     223            WRITE(numout,*) 
     224            WRITE(numout,*) 'sbc : isf surface tracer content forcing fields written in ocean restart file ',   & 
    221225            &                    'at it= ', kt,' date= ', ndastp 
    222          IF(lwp .AND. nprint > 0) WRITE(numout,*) '~~~~' 
     226            WRITE(numout,*) '~~~~' 
     227            IF(lflush) CALL FLUSH(numout) 
     228         ENDIF 
    223229         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    224230         CALL iom_rstput( kt, nitrst, numrow, 'fwf_isf_b', fwfisf(:,:)         , ldxios = lwxios ) 
     
    285291      IF(lwm .AND. nprint > 2) WRITE ( numond, namsbc_isf ) 
    286292 
    287       IF(lwp) WRITE(numout,*) 
    288       IF(lwp) WRITE(numout,*) 'sbc_isf_init : heat flux of the ice shelf' 
    289       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    290       IF(lwp) WRITE(numout,*) '   Namelist namsbc_isf :' 
    291       IF(lwp) WRITE(numout,*) '      type ice shelf melting/freezing         nn_isf      = ', nn_isf 
    292       IF(lwp) WRITE(numout,*) '      bulk formulation (nn_isf=1 only)        nn_isfblk   = ', nn_isfblk 
    293       IF(lwp) WRITE(numout,*) '      thickness of the top boundary layer     rn_hisf_tbl = ', rn_hisf_tbl 
    294       IF(lwp) WRITE(numout,*) '      gamma formulation                       nn_gammablk = ', nn_gammablk  
    295       IF(lwp) WRITE(numout,*) '      gammat coefficient                      rn_gammat0  = ', rn_gammat0   
    296       IF(lwp) WRITE(numout,*) '      gammas coefficient                      rn_gammas0  = ', rn_gammas0   
    297       IF(lwp) WRITE(numout,*) '      top drag coef. used (from namdrg_top)   rn_Cd0      = ', r_Cdmin_top  
     293      IF(lwp) THEN 
     294         WRITE(numout,*) 
     295         WRITE(numout,*) 'sbc_isf_init : heat flux of the ice shelf' 
     296         WRITE(numout,*) '~~~~~~~~~~~~' 
     297         WRITE(numout,*) '   Namelist namsbc_isf :' 
     298         WRITE(numout,*) '      type ice shelf melting/freezing         nn_isf      = ', nn_isf 
     299         WRITE(numout,*) '      bulk formulation (nn_isf=1 only)        nn_isfblk   = ', nn_isfblk 
     300         WRITE(numout,*) '      thickness of the top boundary layer     rn_hisf_tbl = ', rn_hisf_tbl 
     301         WRITE(numout,*) '      gamma formulation                       nn_gammablk = ', nn_gammablk  
     302         WRITE(numout,*) '      gammat coefficient                      rn_gammat0  = ', rn_gammat0   
     303         WRITE(numout,*) '      gammas coefficient                      rn_gammas0  = ', rn_gammas0   
     304         WRITE(numout,*) '      top drag coef. used (from namdrg_top)   rn_Cd0      = ', r_Cdmin_top  
     305         IF(lflush) CALL FLUSH(numout) 
     306      ENDIF 
    298307 
    299308 
     
    377386         CALL ctl_stop( 'sbc_isf_init: wrong value of nn_isf' ) 
    378387      END SELECT 
     388      IF(lwp .AND. lflush) CALL FLUSH(numout) 
    379389          
    380390      rhisf_tbl_0(:,:) = rhisf_tbl(:,:) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcmod.F90

    r10968 r10986  
    105105         WRITE(numout,*) 'sbc_init : surface boundary condition setting' 
    106106         WRITE(numout,*) '~~~~~~~~ ' 
     107         IF(lflush) CALL FLUSH(numout) 
    107108      ENDIF 
    108109      ! 
     
    162163         WRITE(numout,*) '               Stokes coriolis term                 ln_stcor      = ', ln_stcor 
    163164         WRITE(numout,*) '               neutral drag coefficient (CORE,NCAR) ln_cdgw       = ', ln_cdgw 
     165         IF(lflush) CALL FLUSH(numout) 
    164166      ENDIF 
    165167      ! 
     
    207209         CALL ctl_stop( 'sbc_init : unsupported value for nn_components' ) 
    208210      END SELECT 
     211      IF(lwp .AND. lflush) CALL FLUSH(numout) 
    209212      !                             !* coupled options 
    210213      IF( ln_cpl ) THEN 
     
    287290         END SELECT 
    288291         IF( ll_not_nemo  )       WRITE(numout,*) '               + OASIS coupled SAS' 
     292         IF(lflush) CALL FLUSH(numout) 
    289293      ENDIF 
    290294      ! 
     
    303307            WRITE(numout,*)"   OPA-SAS coupled via OASIS : nn_fsbc re-defined from OASIS namcouple ", nn_fsbc 
    304308            WRITE(numout,*) 
     309            IF(lflush) CALL FLUSH(numout) 
    305310         ENDIF 
    306311      ENDIF 
     
    495500         IF( ln_rstart .AND.    &                               !* Restart: read in restart file 
    496501            & iom_varid( numror, 'utau_b', ldstop = .FALSE. ) > 0 ) THEN 
    497             IF(lwp .AND. nprint > 0) WRITE(numout,*) '          nit000-1 surface forcing fields red in the restart file' 
     502            IF(lwp .AND. nprint > 0) THEN 
     503               WRITE(numout,*) '          nit000-1 surface forcing fields red in the restart file' 
     504               IF(lflush) CALL FLUSH(numout) 
     505            ENDIF 
    498506            CALL iom_get( numror, jpdom_autoglo, 'utau_b', utau_b, ldxios = lrxios )   ! before i-stress  (U-point) 
    499507            CALL iom_get( numror, jpdom_autoglo, 'vtau_b', vtau_b, ldxios = lrxios )   ! before j-stress  (V-point) 
     
    509517            ENDIF 
    510518         ELSE                                                   !* no restart: set from nit000 values 
    511             IF(lwp) WRITE(numout,*) '          nit000-1 surface forcing fields set to nit000' 
     519            IF(lwp .AND. nprint > 0) THEN 
     520               WRITE(numout,*) '          nit000-1 surface forcing fields set to nit000' 
     521               IF(lflush) CALL FLUSH(numout) 
     522            ENDIF 
    512523            utau_b(:,:) = utau(:,:) 
    513524            vtau_b(:,:) = vtau(:,:) 
     
    520531      IF( lrst_oce ) THEN                              !      Write in the ocean restart file     ! 
    521532         !                                             ! ---------------------------------------- ! 
    522          IF(lwp .AND. nprint > 0) WRITE(numout,*) 
    523          IF(lwp .AND. nprint > 0) WRITE(numout,*) 'sbc : ocean surface forcing fields written in ocean restart file ',   & 
     533         IF(lwp .AND. nprint > 0) THEN 
     534            WRITE(numout,*) 
     535            WRITE(numout,*) 'sbc : ocean surface forcing fields written in ocean restart file ',   & 
    524536            &                    'at it= ', kt,' date= ', ndastp 
    525          IF(lwp .AND. nprint > 0) WRITE(numout,*) '~~~~' 
     537            WRITE(numout,*) '~~~~' 
     538            IF(lflush) CALL FLUSH(numout) 
     539         ENDIF 
    526540         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    527541         CALL iom_rstput( kt, nitrst, numrow, 'utau_b' , utau, ldxios = lwxios ) 
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/SBC/sbcrnf.F90

    r10968 r10986  
    146146         IF( ln_rstart .AND.    &                               !* Restart: read in restart file 
    147147            & iom_varid( numror, 'rnf_b', ldstop = .FALSE. ) > 0 ) THEN 
    148             IF(lwp .AND. nprint > 0) WRITE(numout,*) '          nit000-1 runoff forcing fields red in the restart file', lrxios 
     148            IF(lwp .AND. nprint > 0) THEN 
     149               WRITE(numout,*) '          nit000-1 runoff forcing fields red in the restart file', lrxios 
     150               IF(lflush) CALL FLUSH(numout) 
     151            ENDIF 
    149152            CALL iom_get( numror, jpdom_autoglo, 'rnf_b', rnf_b, ldxios = lrxios )     ! before runoff 
    150153            CALL iom_get( numror, jpdom_autoglo, 'rnf_hc_b', rnf_tsc_b(:,:,jp_tem), ldxios = lrxios )   ! before heat content of runoff 
    151154            CALL iom_get( numror, jpdom_autoglo, 'rnf_sc_b', rnf_tsc_b(:,:,jp_sal), ldxios = lrxios )   ! before salinity content of runoff 
    152155         ELSE                                                   !* no restart: set from nit000 values 
    153             IF(lwp .AND. nprint > 0) WRITE(numout,*) '          nit000-1 runoff forcing fields set to nit000' 
     156            IF(lwp .AND. nprint > 0) THEN 
     157               WRITE(numout,*) '          nit000-1 runoff forcing fields set to nit000' 
     158               IF(lflush) CALL FLUSH(numout) 
     159            ENDIF 
    154160            rnf_b    (:,:  ) = rnf    (:,:  ) 
    155161            rnf_tsc_b(:,:,:) = rnf_tsc(:,:,:) 
     
    159165      IF( lrst_oce ) THEN                              !      Write in the ocean restart file     ! 
    160166         !                                             ! ---------------------------------------- ! 
    161          IF(lwp .AND. nprint > 0) WRITE(numout,*) 
    162          IF(lwp .AND. nprint > 0) WRITE(numout,*) 'sbcrnf : runoff forcing fields written in ocean restart file ',   & 
     167         IF(lwp .AND. nprint > 0) THEN 
     168            WRITE(numout,*) 
     169            WRITE(numout,*) 'sbcrnf : runoff forcing fields written in ocean restart file ',   & 
    163170            &                    'at it= ', kt,' date= ', ndastp 
    164          IF(lwp .AND. nprint > 0) WRITE(numout,*) '~~~~' 
     171            WRITE(numout,*) '~~~~' 
     172            IF(lflush) CALL FLUSH(numout) 
     173         ENDIF 
    165174         IF( lwxios ) CALL iom_swap(      cwxios_context          ) 
    166175         CALL iom_rstput( kt, nitrst, numrow, 'rnf_b' , rnf, ldxios = lwxios ) 
     
    284293         WRITE(numout,*) '      depth of river mouth additional mixing     rn_hrnf      = ', rn_hrnf 
    285294         WRITE(numout,*) '      multiplicative factor for runoff           rn_rfact     = ', rn_rfact 
     295         IF(lflush) CALL FLUSH(numout) 
    286296      ENDIF 
    287297      !                                   ! ================== 
     
    291301      IF( .NOT. l_rnfcpl ) THEN                     
    292302         ALLOCATE( sf_rnf(1), STAT=ierror )         ! Create sf_rnf structure (runoff inflow) 
    293          IF(lwp) WRITE(numout,*) 
    294          IF(lwp) WRITE(numout,*) '   ==>>>   runoffs inflow read in a file' 
     303         IF(lwp) THEN 
     304            WRITE(numout,*) 
     305            WRITE(numout,*) '   ==>>>   runoffs inflow read in a file' 
     306            IF(lflush) CALL FLUSH(numout) 
     307         ENDIF 
    295308         IF( ierror > 0 ) THEN 
    296309            CALL ctl_stop( 'sbc_rnf_init: unable to allocate sf_rnf structure' )   ;   RETURN 
     
    302315      ! 
    303316      IF( ln_rnf_tem ) THEN                      ! Create (if required) sf_t_rnf structure 
    304          IF(lwp) WRITE(numout,*) 
    305          IF(lwp) WRITE(numout,*) '   ==>>>   runoffs temperatures read in a file' 
     317         IF(lwp) THEN 
     318            WRITE(numout,*) 
     319            WRITE(numout,*) '   ==>>>   runoffs temperatures read in a file' 
     320            IF(lflush) CALL FLUSH(numout) 
     321         ENDIF 
    306322         ALLOCATE( sf_t_rnf(1), STAT=ierror  ) 
    307323         IF( ierror > 0 ) THEN 
     
    314330      ! 
    315331      IF( ln_rnf_sal  ) THEN                     ! Create (if required) sf_s_rnf and sf_t_rnf structures 
    316          IF(lwp) WRITE(numout,*) 
    317          IF(lwp) WRITE(numout,*) '   ==>>>   runoffs salinities read in a file' 
     332         IF(lwp) THEN 
     333            WRITE(numout,*) 
     334            WRITE(numout,*) '   ==>>>   runoffs salinities read in a file' 
     335            IF(lflush) CALL FLUSH(numout) 
     336         ENDIF 
    318337         ALLOCATE( sf_s_rnf(1), STAT=ierror  ) 
    319338         IF( ierror > 0 ) THEN 
     
    326345      ! 
    327346      IF( ln_rnf_depth ) THEN                    ! depth of runoffs set from a file 
    328          IF(lwp) WRITE(numout,*) 
    329          IF(lwp) WRITE(numout,*) '   ==>>>   runoffs depth read in a file' 
     347         IF(lwp) THEN 
     348            WRITE(numout,*) 
     349            WRITE(numout,*) '   ==>>>   runoffs depth read in a file' 
     350            IF(lflush) CALL FLUSH(numout) 
     351         ENDIF 
    330352         rn_dep_file = TRIM( cn_dir )//TRIM( sn_dep_rnf%clname ) 
    331353         IF( .NOT. sn_dep_rnf%ln_clim ) THEN   ;   WRITE(rn_dep_file, '(a,"_y",i4)' ) TRIM( rn_dep_file ), nyear    ! add year  
     
    363385      ELSE IF( ln_rnf_depth_ini ) THEN           ! runoffs applied at the surface 
    364386         ! 
    365          IF(lwp) WRITE(numout,*) 
    366          IF(lwp) WRITE(numout,*) '   ==>>>   depth of runoff computed once from max value of runoff' 
    367          IF(lwp) WRITE(numout,*) '        max value of the runoff climatologie (over global domain) rn_rnf_max = ', rn_rnf_max 
    368          IF(lwp) WRITE(numout,*) '        depth over which runoffs is spread                        rn_dep_max = ', rn_dep_max 
    369          IF(lwp) WRITE(numout,*) '        create (=1) a runoff depth file or not (=0)      nn_rnf_depth_file  = ', nn_rnf_depth_file 
     387         IF(lwp) THEN 
     388            WRITE(numout,*) 
     389            WRITE(numout,*) '   ==>>>   depth of runoff computed once from max value of runoff' 
     390            WRITE(numout,*) '        max value of the runoff climatologie (over global domain) rn_rnf_max = ', rn_rnf_max 
     391            WRITE(numout,*) '        depth over which runoffs is spread                        rn_dep_max = ', rn_dep_max 
     392            WRITE(numout,*) '        create (=1) a runoff depth file or not (=0)      nn_rnf_depth_file  = ', nn_rnf_depth_file 
     393            IF(lflush) CALL FLUSH(numout) 
     394         ENDIF 
    370395 
    371396         CALL iom_open( TRIM( sn_rnf%clname ), inum )    !  open runoff file 
     
    417442         ! 
    418443         IF( nn_rnf_depth_file == 1 ) THEN      !  save  output nb levels for runoff 
    419             IF(lwp) WRITE(numout,*) '   ==>>>   create runoff depht file' 
     444            IF(lwp) THEN 
     445               WRITE(numout,*) '   ==>>>   create runoff depht file' 
     446               IF(lflush) CALL FLUSH(numout) 
     447            ENDIF 
    420448            CALL iom_open  ( TRIM( sn_dep_rnf%clname ), inum, ldwrt = .TRUE. ) 
    421449            CALL iom_rstput( 0, 0, inum, 'rodepth', h_rnf ) 
     
    449477            IF( ln_sco )   CALL ctl_warn( 'sbc_rnf_init: number of levels over which Kz is increased is computed for zco...' ) 
    450478         ENDIF 
    451          IF(lwp) WRITE(numout,*) 
    452          IF(lwp) WRITE(numout,*) '   ==>>>   Specific treatment used in vicinity of river mouths :' 
    453          IF(lwp) WRITE(numout,*) '             - Increase Kz in surface layers (if rn_hrnf > 0 )' 
    454          IF(lwp) WRITE(numout,*) '               by ', rn_avt_rnf,' m2/s  over ', nkrnf, ' w-levels' 
    455          IF(lwp) WRITE(numout,*) '             - set to zero SSS damping       (if ln_ssr=T)' 
    456          IF(lwp) WRITE(numout,*) '             - mixed upstream-centered       (if ln_traadv_cen2=T)' 
     479         IF(lwp) THEN 
     480            WRITE(numout,*) 
     481            WRITE(numout,*) '   ==>>>   Specific treatment used in vicinity of river mouths :' 
     482            WRITE(numout,*) '             - Increase Kz in surface layers (if rn_hrnf > 0 )' 
     483            WRITE(numout,*) '               by ', rn_avt_rnf,' m2/s  over ', nkrnf, ' w-levels' 
     484            WRITE(numout,*) '             - set to zero SSS damping       (if ln_ssr=T)' 
     485            WRITE(numout,*) '             - mixed upstream-centered       (if ln_traadv_cen2=T)' 
     486            IF(lflush) CALL FLUSH(numout)