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/LIM_SRC/limdyn.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/LIM_SRC/limdyn.F90

    r247 r258  
    2323   USE lbclnk 
    2424   USE lib_mpp 
     25   USE prtctl          ! Print control 
    2526 
    2627   IMPLICIT NONE 
     
    9293            i_j1 = 1    
    9394            i_jpj = jpj 
    94             IF(l_ctl)   WRITE(numout,*) 'lim_dyn :    i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     95            IF(ln_ctl)    THEN 
     96               CALL prt_ctl_info('lim_dyn  :    i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj) 
     97            ENDIF 
    9598            CALL lim_rhg( i_j1, i_jpj ) 
    9699 
     
    113116               END DO 
    114117               i_j1 = MAX( 1, i_j1-1 ) 
    115                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     118               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    116119     
    117120               CALL lim_rhg( i_j1, i_jpj ) 
     
    124127               END DO 
    125128               i_jpj = MIN( jpj, i_jpj+2 ) 
    126                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     129               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    127130     
    128131               CALL lim_rhg( i_j1, i_jpj ) 
     
    143146               i_jpj = MIN( jpj, i_jpj+2) 
    144147     
    145                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     148               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    146149     
    147150               CALL lim_rhg( i_j1, i_jpj ) 
     
    154157         v_ice(:,1) = 0.e0 
    155158 
    156          IF(l_ctl) THEN  
    157             WRITE(numout,*) ' lim_dyn  : u_oce ', SUM( u_oce(2:nictl,2:njctl) ), ' v_oce ', SUM( v_oce(2:nictl,2:njctl) ) 
    158             WRITE(numout,*) ' lim_dyn  : u_ice ', SUM( u_ice(2:nictl,2:njctl) ), ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 
     159         IF(ln_ctl)  THEN  
     160            CALL prt_ctl(tab2d_1=u_oce , clinfo1=' lim_dyn  : u_oce :', tab2d_2=v_oce , clinfo2=' v_oce :') 
     161            CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_dyn  : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 
    159162         ENDIF 
    160163          
     
    244247      CALL lbc_lnk( tio_v, 'I', -1. )   ! I-point (i.e. ice U-V point) 
    245248 
    246       IF(l_ctl) THEN  
    247          WRITE(numout,*) ' lim_dyn  : tio_u ', SUM( tio_u(2:nictl,2:njctl) ), ' tio_v ', SUM( tio_v(2:nictl,2:njctl) ) 
    248          WRITE(numout,*) ' lim_dyn  : ust2s ', SUM( ust2s(2:nictl,2:njctl) ) 
     249      IF(ln_ctl) THEN  
     250            CALL prt_ctl(tab2d_1=tio_u , clinfo1=' lim_dyn  : tio_u :', tab2d_2=tio_v , clinfo2=' tio_v :') 
     251            CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn  : ust2s :') 
    249252      ENDIF 
    250253 
Note: See TracChangeset for help on using the changeset viewer.