Changeset 84 for trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

Timestamp:

2004-04-22T15:32:41+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE057 : # General syntax, alignement, comments corrections

# l_ctl alone replace the set (l_ctl .AND. lwp)
# Add of diagnostics which are activated when using l_ctl logical

File:

• trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

@@ -73 +73 @@
       !! * Local declarations
       INTEGER ::   ji, jj, jk                 ! dummy loop indices
+      REAL(wp) ::   zua, zva                  ! temporary scalars
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   &
          wk1, wk2,            &  ! work array used for rotated biharmonic
@@ -126 +127 @@
       END DO                                           !   End of slab
       !                                                ! ===============
+
+      IF(l_ctl) THEN         ! print sum trends (used for debugging)
+         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) )
+         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) )
+         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl
+         u_ctl = zua   ;   v_ctl = zva
+      ENDIF
 
    END SUBROUTINE dyn_ldf_bilapg