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 84 for trunk/NEMO/OPA_SRC/DYN – NEMO

Changeset 84 for trunk/NEMO/OPA_SRC/DYN


Ignore:
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

Location:
trunk/NEMO/OPA_SRC/DYN
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/DYN/dynhpg.F90

    r32 r84  
    162162      END DO 
    163163 
    164       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     164      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    165165         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    166166         zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    … …  
    317317      END DO 
    318318 
    319       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     319      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    320320         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    321321         zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    … …  
    428428      END DO 
    429429 
    430       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     430      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    431431         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    432432         zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynhpg_atsk.F90

    r32 r84  
    150150      !                                                ! =============== 
    151151 
    152       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     152      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    153153         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    154154         zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    … …  
    295295      END DO                                           !   End of slab 
    296296      !                                                ! =============== 
     297 
     298      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     299         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     300         zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     301         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
     302         u_ctl = zuap   ;   v_ctl = zvap 
     303      ENDIF    
     304 
    297305   END SUBROUTINE dyn_hpg_atsk 
    298306 
    … …  
    402410      !                                                ! =============== 
    403411 
    404       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     412      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    405413         WRITE(numout,*) ' hpg  - Ua: ', SUM(ua*umask), ' Va: ', SUM(va*vmask) 
    406414         zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynkeg.F90

    r3 r84  
    119119      !                                                ! =============== 
    120120 
    121       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     121      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    122122         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    123123         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r3 r84  
    219219      END DO                                           !   End of slab 
    220220      !                                                ! =============== 
     221 
     222      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     223         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     224         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     225         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
     226         u_ctl = zua   ;   v_ctl = zva 
     227      ENDIF 
     228 
    221229   END SUBROUTINE dyn_ldf_bilap 
    222230 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

    r32 r84  
    7373      !! * Local declarations 
    7474      INTEGER ::   ji, jj, jk                 ! dummy loop indices 
     75      REAL(wp) ::   zua, zva                  ! temporary scalars 
    7576      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    7677         wk1, wk2,            &  ! work array used for rotated biharmonic 
    … …  
    126127      END DO                                           !   End of slab 
    127128      !                                                ! =============== 
     129 
     130      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     131         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     132         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     133         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
     134         u_ctl = zua   ;   v_ctl = zva 
     135      ENDIF 
    128136 
    129137   END SUBROUTINE dyn_ldf_bilapg 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_lap.F90

    r3 r84  
    118118      !                                                ! =============== 
    119119 
    120       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     120      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    121121         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    122122         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynnxt.F90

    r32 r84  
    133133      !                                                ! =============== 
    134134 
    135       IF( l_ctl .AND. lwp ) THEN         ! print sum fields (used for debugging) 
    136          WRITE(numout,*) ' nxt  - Un: ', SUM(un*umask), ' Vn: ', SUM(vn*vmask) 
    137       ENDIF 
     135      IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un*umask), ' Vn: ', SUM(vn*vmask) 
    138136 
    139137   END SUBROUTINE dyn_nxt 

  • trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90

    r31 r84  
    347347      END DO 
    348348 
    349       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     349      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    350350         WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:jpim1,2:jpjm1,1:jpkm1)*umask(2:jpim1,2:jpjm1,1:jpkm1) ),   & 
    351351            &                   ' Va: ', SUM( va(2:jpim1,2:jpjm1,1:jpkm1)*vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynspg_fsc_atsk.F90

    r31 r84  
    1818   USE oce             ! ocean dynamics and tracers  
    1919   USE dom_oce         ! ocean space and time domain 
    20    USE trdtra_oce     ! ocean active tracer trend  
    21    USE trddyn_oce     ! ocean dynamics trend 
     20   USE zdf_oce         ! ocean vertical physics 
     21   USE trdtra_oce      ! ocean active tracer trend  
     22   USE trddyn_oce      ! ocean dynamics trend 
    2223   USE in_out_manager  ! I/O manager 
    2324   USE phycst          ! physical constant 
    24    USE ocesbc          ! Ocen Surface Boundary condition 
     25   USE ocesbc          ! Ocean Surface Boundary condition 
    2526   USE flxrnf          ! ??? 
    2627   USE sol_oce         ! ocean elliptic solver 
    … …  
    3031   USE obc_oce         ! Lateral open boundary condition 
    3132   USE obcdyn          ! open boudary condition 
    32    USE obcdyn          !    "              " 
    3333   USE obcvol          !    "              " 
    3434   USE lib_mpp         ! ??? 
    … …  
    362362         END DO 
    363363 
     364         IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     365            WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:jpim1,2:jpjm1,1:jpkm1)*umask(2:jpim1,2:jpjm1,1:jpkm1) ),   & 
     366               &                   ' Va: ', SUM( va(2:jpim1,2:jpjm1,1:jpkm1)*vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     367         ENDIF 
     368 
    364369         ! 8. Sea surface elevation time stepping 
    365370         ! -------------------------------------- 

  • trunk/NEMO/OPA_SRC/DYN/dynzad.F90

    r3 r84  
    126126      !                                                ! =============== 
    127127 
    128       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     128      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    129129         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    130130         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    … …  
    227227      END DO 
    228228 
    229       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     229      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    230230         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    231231         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

    r3 r84  
    342342      END DO 
    343343 
    344       IF( l_ctl .AND. lwp ) THEN         ! print sum trends (used for debugging) 
     344      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    345345         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    346346         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_imp_atsk.F90

    r3 r84  
    312312      END DO                                           !   End of slab 
    313313      !                                                ! =============== 
     314 
     315      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     316         zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     317         zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     318         WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
     319         u_ctl = zua   ;   v_ctl = zva 
     320      ENDIF 
     321 
    314322   END SUBROUTINE dyn_zdf_imp_tsk 
    315323 
Note: See TracChangeset for help on using the changeset viewer.