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

Changeset 258 for trunk/NEMO/OPA_SRC/DYN


Ignore:
Timestamp:
2005-09-02T17:48:36+02:00 (19 years ago)
Author:
opalod
Message:

nemo_v1_update_004 : CT : Integration of the control print option for debugging work

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

Legend:

Unmodified
Added
Removed

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

    r247 r258  
    2121   USE trdmod          ! ocean dynamics trends  
    2222   USE trdmod_oce      ! ocean variables trends 
     23   USE prtctl          ! Print control 
    2324 
    2425   IMPLICIT NONE 
    … …  
    172173      ENDIF 
    173174 
    174       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    175          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    176          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    177          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    178          u_ctl = zuap   ;   v_ctl = zvap 
     175      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     176         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     177            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    179178      ENDIF 
    180179 
    … …  
    325324      ENDIF 
    326325 
    327       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    328          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    329          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    330          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    331          u_ctl = zuap   ;   v_ctl = zvap 
     326      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     327         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     328            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    332329      ENDIF 
    333330 
    … …  
    443440      ENDIF 
    444441 
    445       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    446          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    447          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    448          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    449          u_ctl = zuap   ;   v_ctl = zvap 
     442      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     443         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     444            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    450445      ENDIF 
    451446 

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

    r247 r258  
    2121   USE trdmod          ! ocean dynamics trends  
    2222   USE trdmod_oce      ! ocean variables trends 
     23   USE prtctl          ! Print control 
    2324 
    2425   IMPLICIT NONE 
    … …  
    162163      ENDIF 
    163164 
    164       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    165          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    166          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    167          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    168          u_ctl = zuap   ;   v_ctl = zvap 
     165      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     166         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     167            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    169168      ENDIF 
    170169 
    … …  
    306305      ENDIF 
    307306 
    308       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    309          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    310          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    311          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    312          u_ctl = zuap   ;   v_ctl = zvap 
     307      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     308         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     309            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    313310      ENDIF    
    314311 
    … …  
    429426      ENDIF 
    430427 
    431       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    432          zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    433          zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    434          WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    435          u_ctl = zuap   ;   v_ctl = zvap 
     428      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     429         CALL prt_ctl(tab3d_1=ua, clinfo1=' hpg  - Ua: ', mask1=umask, & 
     430            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    436431      ENDIF 
    437432 

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

    r247 r258  
    1414   USE trdmod          ! ocean dynamics trends  
    1515   USE trdmod_oce      ! ocean variables trends 
     16   USE prtctl          ! Print control 
    1617 
    1718   IMPLICIT NONE 
    … …  
    119120      ENDIF 
    120121 
    121       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    122          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    123          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    124          WRITE(numout,*) ' keg  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    125          u_ctl = zua   ;   v_ctl = zva 
     122      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     123         CALL prt_ctl(tab3d_1=ua, clinfo1=' keg  - Ua: ', mask1=umask, & 
     124            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    126125      ENDIF 
    127126 

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

    r247 r258  
    1717   USE trdmod_oce      ! ocean variables trends 
    1818   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     19   USE prtctl          ! Print control 
    1920 
    2021   IMPLICIT NONE 
    … …  
    237238      ENDIF 
    238239 
    239       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    240          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    241          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    242          WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    243          u_ctl = zua   ;   v_ctl = zva 
     240      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     241         CALL prt_ctl(tab3d_1=ua, clinfo1=' ldf  - Ua: ', mask1=umask, & 
     242            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    244243      ENDIF 
    245244 

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

    r247 r258  
    2222   USE ldfslp          ! iso-neutral slopes available 
    2323   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     24   USE prtctl          ! Print control 
    2425 
    2526   IMPLICIT NONE 
    … …  
    8182      !! * Local declarations 
    8283      INTEGER ::   ji, jj, jk                 ! dummy loop indices 
    83       REAL(wp) ::   zua, zva                  ! temporary scalars 
    8484      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    8585         zwk1, zwk2                ! work array used for rotated biharmonic 
    … …  
    138138      ENDIF 
    139139 
    140       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    141          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    142          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    143          WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    144          u_ctl = zua   ;   v_ctl = zva 
     140      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     141         CALL prt_ctl(tab3d_1=ua, clinfo1=' ldf  - Ua: ', mask1=umask, & 
     142            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    145143      ENDIF 
    146144 

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

    r247 r258  
    2222   USE ldfslp          ! iso-neutral slopes  
    2323   USE in_out_manager  ! I/O manager 
     24   USE prtctl          ! Print control 
    2425 
    2526   IMPLICIT NONE 
    … …  
    9899         zabe1, zabe2, zcof1, zcof2,   &  ! temporary scalars 
    99100         zmskt, zmskf, zbu, zbv,       & 
    100          zuah, zvah, zua, zva 
     101         zuah, zvah 
    101102      REAL(wp), DIMENSION(jpi,jpj) ::   & 
    102103         ziut, zjuf, zjvt, zivf,        & ! temporary workspace 
    … …  
    259260      ENDIF 
    260261 
    261       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    262          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    263          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    264          WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    265          u_ctl = zua   ;   v_ctl = zva 
     262      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     263         CALL prt_ctl(tab3d_1=ua, clinfo1=' ldf  - Ua: ', mask1=umask, & 
     264            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    266265      ENDIF 
    267266 

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

    r247 r258  
    1818   USE trdmod_oce      ! ocean variables trends 
    1919   USE ldfslp          ! iso-neutral slopes  
     20   USE prtctl          ! Print control 
    2021 
    2122   IMPLICIT NONE 
    … …  
    136137      ENDIF 
    137138 
    138       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    139          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    140          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    141          WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    142          u_ctl = zua   ;   v_ctl = zva 
     139      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     140         CALL prt_ctl(tab3d_1=ua, clinfo1=' ldf  - Ua: ', mask1=umask, & 
     141            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    143142      ENDIF 
    144143 

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

    r247 r258  
    1414   USE obcdyn          ! open boundary condition for momentum (obc_dyn routine) 
    1515   USE lbclnk          ! lateral boundary condition (or mpp link) 
     16   USE prtctl          ! Print control 
    1617 
    1718   IMPLICIT NONE 
    … …  
    135136      !                                                ! =============== 
    136137 
    137       IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 
    138       &                                  ' Vn: ', SUM(vn(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
     138      IF(ln_ctl)   THEN 
     139         CALL prt_ctl(tab3d_1=un, clinfo1=' nxt  - Un: ', mask1=umask, & 
     140            &         tab3d_2=vn, clinfo2=' Vn: ', mask2=vmask) 
     141      ENDIF 
    139142 
    140143   END SUBROUTINE dyn_nxt 

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

    r247 r258  
    3333   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    3434   USE cla_dynspg      ! cross land advection 
     35   USE prtctl          ! Print control 
    3536 
    3637   IMPLICIT NONE 
    … …  
    337338      ENDIF 
    338339 
    339       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    340          WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ),  & 
    341             &                   ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 
     340      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     341         CALL prt_ctl(tab3d_1=ua, clinfo1=' spg  - Ua: ', mask1=umask, & 
     342            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask) 
    342343      ENDIF 
    343344 
    … …  
    370371 
    371372 
    372       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    373          WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 
     373      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     374         CALL prt_ctl(tab2d_1=sshn, clinfo1=' spg  - ssh: ', mask1=tmask) 
    374375      ENDIF 
    375376 

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

    r247 r258  
    3737   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    3838   USE cla_dynspg      ! cross land advection 
     39   USE prtctl          ! Print control 
    3940 
    4041   IMPLICIT NONE 
    … …  
    354355         END DO 
    355356 
    356          IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    357             WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ),  & 
    358                &                   ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 
     357         IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     358            CALL prt_ctl(tab3d_1=ua, clinfo1=' spg  - Ua: ', mask1=umask, & 
     359               &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask) 
    359360         ENDIF 
    360361 
    … …  
    398399      CALL lbc_lnk( sshn, 'T', 1. ) 
    399400 
    400       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    401          WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 
     401      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     402         CALL prt_ctl(tab2d_1=sshn, clinfo1=' spg  - ssh: ', mask1=tmask) 
    402403      ENDIF 
    403404 

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

    r247 r258  
    2020   USE trdmod_oce      ! ocean variables trends 
    2121   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     22   USE prtctl          ! Print control 
    2223 
    2324   IMPLICIT NONE 
    … …  
    188189      ENDIF 
    189190 
    190       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    191          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    192          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    193          WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    194          u_ctl = zua   ;   v_ctl = zva 
     191      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     192         CALL prt_ctl(tab3d_1=ua, clinfo1=' vor  - Ua: ', mask1=umask, & 
     193            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    195194      ENDIF 
    196195 
    … …  
    342341      ENDIF 
    343342 
    344       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    345          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    346          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    347          WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    348          u_ctl = zua   ;   v_ctl = zva 
     343      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     344         CALL prt_ctl(tab3d_1=ua, clinfo1=' vor  - Ua: ', mask1=umask, & 
     345            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    349346      ENDIF 
    350347 
    … …  
    499496      ENDIF 
    500497 
    501       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    502          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    503          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    504          WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    505          u_ctl = zua   ;   v_ctl = zva 
     498      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     499         CALL prt_ctl(tab3d_1=ua, clinfo1=' vor  - Ua: ', mask1=umask, & 
     500            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    506501      ENDIF 
    507502 
    … …  
    662657      ENDIF 
    663658 
    664       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    665          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    666          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    667          WRITE(numout,*) ' vor  een - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    668          u_ctl = zua   ;   v_ctl = zva 
     659      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     660         CALL prt_ctl(tab3d_1=ua, clinfo1=' vor  - Ua: ', mask1=umask, & 
     661            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    669662      ENDIF 
    670663 

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

    r247 r258  
    1515   USE trdmod_oce      ! ocean variables trends 
    1616   USE flxrnf          ! ocean runoffs 
     17   USE prtctl          ! Print control 
    1718 
    1819   IMPLICIT NONE 
    … …  
    142143      ENDIF 
    143144 
    144       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    145          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    146          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    147          WRITE(numout,*) ' zad  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    148          u_ctl = zua   ;   v_ctl = zva 
     145      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     146         CALL prt_ctl(tab3d_1=ua, clinfo1=' zad  - Ua: ', mask1=umask, & 
     147            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    149148      ENDIF 
    150149 
    … …  
    255254      ENDIF 
    256255 
    257       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    258          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    259          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    260          WRITE(numout,*) ' zad  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    261          u_ctl = zua   ;   v_ctl = zva 
     256      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     257         CALL prt_ctl(tab3d_1=ua, clinfo1=' zad  - Ua: ', mask1=umask, & 
     258            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    262259      ENDIF 
    263260 

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

    r247 r258  
    1818   USE trdmod          ! ocean dynamics trends  
    1919   USE trdmod_oce      ! ocean variables trends 
     20   USE prtctl          ! Print control 
    2021 
    2122   IMPLICIT NONE 
    … …  
    197198      ENDIF 
    198199 
    199       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    200          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    201          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    202          WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    203          u_ctl = zua   ;   v_ctl = zva 
     200      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     201         CALL prt_ctl(tab3d_1=ua, clinfo1=' zdf  - Ua: ', mask1=umask, & 
     202            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    204203      ENDIF 
    205204 

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

    r247 r258  
    2222   USE trdmod          ! ocean dynamics trends  
    2323   USE trdmod_oce      ! ocean variables trends 
     24   USE prtctl          ! Print control 
    2425 
    2526   IMPLICIT NONE 
    … …  
    7980         ikbu, ikbum1, ikbv, ikbvm1      ! temporary integers 
    8081      REAL(wp) ::   & 
    81          zrau0r, z2dt, zua, zva,      &  ! temporary scalars 
     82         zrau0r, z2dt,                &  ! temporary scalars 
    8283         z2dtf, zcoef, zzws, zrhs        !    "         " 
    8384      REAL(wp), DIMENSION(jpi,jpj) ::   & 
    … …  
    355356      ENDIF 
    356357 
    357       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    358          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    359          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    360          WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    361          u_ctl = zua   ;   v_ctl = zva 
     358      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     359         CALL prt_ctl(tab3d_1=ua, clinfo1=' zdf  - Ua: ', mask1=umask, & 
     360            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    362361      ENDIF 
    363362 

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

    r247 r258  
    1919   USE trdmod          ! ocean dynamics trends  
    2020   USE trdmod_oce      ! ocean variables trends 
     21   USE prtctl          ! Print control 
    2122 
    2223   IMPLICIT NONE 
    … …  
    7475         ikbu, ikbum1, ikbv, ikbvm1      !    "         "       
    7576      REAL(wp) ::   & 
    76          zrau0r, z2dt, zua, zva,       & !temporary scalars 
     77         zrau0r, z2dt,                 & !temporary scalars 
    7778         z2dtf, zcoef, zzws 
    7879      REAL(wp), DIMENSION(jpi,jpk) ::  & 
    … …  
    338339      ENDIF 
    339340 
    340       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    341          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    342          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    343          WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    344          u_ctl = zua   ;   v_ctl = zva 
     341      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     342         CALL prt_ctl(tab3d_1=ua, clinfo1=' zdf  - Ua: ', mask1=umask, & 
     343            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    345344      ENDIF 
    346345 

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

    r247 r258  
    2121   USE trdmod          ! ocean dynamics trends  
    2222   USE trdmod_oce      ! ocean variables trends 
     23   USE prtctl          ! Print control 
    2324 
    2425   IMPLICIT NONE 
    … …  
    110111      REAL(wp) ::   & 
    111112         zrau0r, z2dt,                       & ! temporary scalars 
    112          z2dtf, zua, zva, zcoef, zzws 
     113         z2dtf, zcoef, zzws 
    113114      REAL(wp) ::   & 
    114115         zcoef0, zcoef3, zcoef4, zbu, zbv, zmkt, zmkf,   & 
    … …  
    406407      ENDIF 
    407408 
    408       IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    409          zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    410          zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    411          WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    412          u_ctl = zua   ;   v_ctl = zva 
     409      IF(ln_ctl) THEN         ! print sum trends (used for debugging) 
     410         CALL prt_ctl(tab3d_1=ua, clinfo1=' zdf  - Ua: ', mask1=umask, & 
     411            &         tab3d_2=va, clinfo2=' Va: ', mask2=vmask, clinfo3='dyn') 
    413412      ENDIF 
    414413 

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

    r247 r258  
    1212   USE dom_oce         ! ocean space and time domain variables  
    1313   USE in_out_manager  ! I/O manager 
     14   USE prtctl          ! Print control 
    1415 
    1516   IMPLICIT NONE 
    … …  
    8081      !                                                ! =============== 
    8182 
    82       IF(l_ctl) THEN                ! print mean trends (used for debugging) 
    83          WRITE(numout,*) ' w**2 -   : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn   (2:nictl,2:njctl,1:jpkm1) ) 
    84       ENDIF 
     83      IF(ln_ctl)   CALL prt_ctl(tab3d_1=wn, clinfo1=' w**2 -   : ', mask1=wn) 
    8584 
    8685   END SUBROUTINE wzv 
    … …  
    131130      END DO 
    132131 
    133       IF(l_ctl) THEN                ! print mean trends (used for debugging) 
    134          WRITE(numout,*) ' w**2 -   : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn   (2:nictl,2:njctl,1:jpkm1) ) 
    135       ENDIF 
     132      IF(ln_ctl)   CALL prt_ctl(tab3d_1=wn, clinfo1=' w**2 -   : ', mask1=wn) 
    136133 
    137134   END SUBROUTINE wzv 
Note: See TracChangeset for help on using the changeset viewer.