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/DOM – 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/DOM
Files:
13 edited

Legend:

Unmodified
Added
Removed

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