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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10986 for NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA – NEMO

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

GMED 462 add flush

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

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DIA/dia25h.F90

    r10968 r10986  
    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 
Note: See TracChangeset for help on using the changeset viewer.