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/dynhpg.F90 – NEMO

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

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