Changeset 10986 for NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domvvl.F90

Timestamp:

2019-05-16T15:23:56+02:00 (5 years ago)

Author:

andmirek

Message:

GMED 462 add flush

File:

• NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domvvl.F90 (modified) (9 diffs)

NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/DOM/domvvl.F90

@@ -120 +120 @@
       !!----------------------------------------------------------------------
       !
-      IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) 'dom_vvl_init : Variable volume activated'
-      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~'
+      IF(lwp) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) 'dom_vvl_init : Variable volume activated'
+         WRITE(numout,*) '~~~~~~~~~~~~'
+         IF(lflush) CALL FLUSH(numout)
+      ENDIF
       !
       CALL dom_vvl_ctl     ! choose vertical coordinate (z_star, z_tilde or layer)
@@ -304 +307 @@
       !
       IF( kt == nit000 ) THEN
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors'
-         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~'
+         IF(lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*) 'dom_vvl_sf_nxt : compute after scale factors'
+            WRITE(numout,*) '~~~~~~~~~~~~~~'
+            IF(lflush) CALL FLUSH(numout)
+         ENDIF
       ENDIF
 
@@ -533 +539 @@
          z_tmax = MAXVAL( tmask(:,:,1) *  ABS( ssha(:,:) ) )
          CALL mpp_max( 'domvvl', z_tmax )                                ! max over the global domain
-         IF( lwp    ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax
+         IF( lwp    ) THEN
+            WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax
+            IF(lflush) CALL FLUSH(numout)
+         ENDIF
       END IF
 
@@ -601 +610 @@
          IF(lwp) WRITE(numout,*) 'dom_vvl_sf_swp : - time filter and swap of scale factors'
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   - interpolate scale factors and compute depths for next time step'
+         IF(lflush) CALL FLUSH(numout)
       ENDIF
       !
@@ -820 +830 @@
                CALL iom_get( numror, jpdom_autoglo, 'e3t_n', e3t_n(:,:,:), ldxios = lrxios )
                ! needed to restart if land processor not computed 
-               IF(lwp) write(numout,*) 'dom_vvl_rst : e3t_b and e3t_n found in restart files'
+               IF(lwp) THEN
+                  write(numout,*) 'dom_vvl_rst : e3t_b and e3t_n found in restart files'
+                  IF(lflush) CALL FLUSH(numout)
+               ENDIF
                WHERE ( tmask(:,:,:) == 0.0_wp ) 
                   e3t_n(:,:,:) = e3t_0(:,:,:)
@@ -829 +842 @@
                ENDIF
             ELSE IF( id1 > 0 ) THEN
-               IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart files'
-               IF(lwp) write(numout,*) 'e3t_n set equal to e3t_b.'
-               IF(lwp) write(numout,*) 'neuler is forced to 0'
+               IF(lwp) THEN
+                  write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart files'
+                  write(numout,*) 'e3t_n set equal to e3t_b.'
+                  write(numout,*) 'neuler is forced to 0'
+                  IF(lflush) CALL FLUSH(numout)
+               ENDIF
                CALL iom_get( numror, jpdom_autoglo, 'e3t_b', e3t_b(:,:,:), ldxios = lrxios )
                e3t_n(:,:,:) = e3t_b(:,:,:)
                neuler = 0
             ELSE IF( id2 > 0 ) THEN
-               IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_b not found in restart files'
-               IF(lwp) write(numout,*) 'e3t_b set equal to e3t_n.'
-               IF(lwp) write(numout,*) 'neuler is forced to 0'
+               IF(lwp) THEN
+                  write(numout,*) 'dom_vvl_rst WARNING : e3t_b not found in restart files'
+                  write(numout,*) 'e3t_b set equal to e3t_n.'
+                  write(numout,*) 'neuler is forced to 0'
+                  IF(lflush) CALL FLUSH(numout)
+               ENDIF
                CALL iom_get( numror, jpdom_autoglo, 'e3t_n', e3t_n(:,:,:), ldxios = lrxios )
                e3t_b(:,:,:) = e3t_n(:,:,:)
                neuler = 0
             ELSE
-               IF(lwp) write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart file'
-               IF(lwp) write(numout,*) 'Compute scale factor from sshn'
-               IF(lwp) write(numout,*) 'neuler is forced to 0'
+               IF(lwp) THEN
+                  write(numout,*) 'dom_vvl_rst WARNING : e3t_n not found in restart file'
+                  write(numout,*) 'Compute scale factor from sshn'
+                  write(numout,*) 'neuler is forced to 0'
+                  IF(lflush) CALL FLUSH(numout)
+               ENDIF
                DO jk = 1, jpk
                   e3t_n(:,:,jk) =  e3t_0(:,:,jk) * ( ht_0(:,:) + sshn(:,:) ) &
@@ -952 +974 @@
       ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file
          !                                   ! ===================
-         IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- dom_vvl_rst ----'
+         IF(lwp .AND. nprint > 0) THEN
+            WRITE(numout,*) '---- dom_vvl_rst ----'
+            IF(lflush) CALL FLUSH(numout)
+         ENDIF
          IF( lwxios ) CALL iom_swap(      cwxios_context          )
          !                                           ! --------- !
@@ -1024 +1049 @@
          ENDIF
          WRITE(numout,*) '         debug prints flag                                 ln_vvl_dbg   = ', ln_vvl_dbg
+         IF(lflush) CALL FLUSH(numout)
       ENDIF
       !
@@ -1041 +1067 @@
          IF( ln_vvl_layer           ) WRITE(numout,*) '      ==>>>   layer vertical coordinate is used'
          IF( ln_vvl_ztilde_as_zstar ) WRITE(numout,*) '      ==>>>   to emulate a zstar coordinate'
+         IF(lflush) CALL FLUSH(numout)
       ENDIF
       !