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 106 – NEMO

Changeset 106


Ignore:
Timestamp:
2004-06-28T10:08:06+02:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE067 : Add control indices nictl, njctl used in SUM function output to compare mono versus multi procs runs

Location:
trunk/NEMO
Files:
51 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/LIM_SRC/icestp.F90

    r88 r106  
    184184          
    185185 
    186          IF( l_ctl ) THEN         ! print mean trends (used for debugging) 
     186         IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    187187            WRITE(numout,*) 'Ice Forcings ' 
    188             WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (:,:) ), ' qsr_ice  : ', SUM( qsr_ice (:,:) ) 
    189             WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(:,:) ), ' qnsr_ice : ', SUM( qnsr_ice(:,:) ) 
    190             WRITE(numout,*) ' evap     : ', SUM( evap (:,:) ) 
    191             WRITE(numout,*) ' precip   : ', SUM( tprecip(:,:)  ), ' Snow     : ', SUM( sprecip (:,:) ) 
    192             WRITE(numout,*) ' u-stress : ', SUM( gtaux  (:,:)  ), ' v-stress : ', SUM( gtauy  (:,:) ) 
    193             WRITE(numout,*) ' sst      : ', SUM( sst_io (:,:)  ), ' sss      : ', SUM( sss_io (:,:) ) 
    194             WRITE(numout,*) ' u_io     : ', SUM( u_io   (:,:)  ), ' v_io     : ', SUM( v_io   (:,:) ) 
    195             WRITE(numout,*) ' hsnif  1 : ', SUM( hsnif  (:,:)  ), ' hicnif   : ', SUM( hicif  (:,:) ) 
    196             WRITE(numout,*) ' frld   1 : ', SUM( frld   (:,:)  ), ' sist     : ', SUM( sist   (:,:) ) 
     188            WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (2:nictl,2:njctl) ), ' qsr_ice  : ', SUM( qsr_ice (2:nictl,2:njctl) ) 
     189            WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(2:nictl,2:njctl) ), ' qnsr_ice : ', SUM( qnsr_ice(2:nictl,2:njctl) ) 
     190            WRITE(numout,*) ' evap     : ', SUM( evap    (2:nictl,2:njctl) ) 
     191            WRITE(numout,*) ' precip   : ', SUM( tprecip(2:nictl,2:njctl)  ), ' Snow     : ', SUM( sprecip (2:nictl,2:njctl) ) 
     192            WRITE(numout,*) ' u-stress : ', SUM( gtaux  (2:nictl,2:njctl)  ), ' v-stress : ', SUM( gtauy  (2:nictl,2:njctl) ) 
     193            WRITE(numout,*) ' sst      : ', SUM( sst_io (2:nictl,2:njctl)  ), ' sss      : ', SUM( sss_io (2:nictl,2:njctl) ) 
     194            WRITE(numout,*) ' u_io     : ', SUM( u_io   (2:nictl,2:njctl)  ), ' v_io     : ', SUM( v_io   (2:nictl,2:njctl) ) 
     195            WRITE(numout,*) ' hsnif  1 : ', SUM( hsnif  (2:nictl,2:njctl)  ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
     196            WRITE(numout,*) ' frld   1 : ', SUM( frld   (2:nictl,2:njctl)  ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    197197         ENDIF 
    198198 
    … …  
    203203         CALL lim_dyn                                                ! Ice dynamics !   ( rheology/dynamics ) 
    204204         !                                                           !--------------! 
    205          IF( l_ctl ) THEN 
    206             WRITE(numout,*) ' hsnif  2 : ', SUM( hsnif  (:,:) ), ' hicnif   : ', SUM( hicif  (:,:) ) 
    207             WRITE(numout,*) ' frld   2 : ', SUM( frld   (:,:) ), ' sist     : ', SUM( sist   (:,:) ) 
     205         IF(l_ctl) THEN 
     206            WRITE(numout,*) ' hsnif  2 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
     207            WRITE(numout,*) ' frld   2 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    208208         ENDIF 
    209209 
    … …  
    212212         CALL lim_trp                                                ! Ice transport !  ( Advection/diffusion ) 
    213213         !                                                           !---------------! 
    214          IF( l_ctl ) THEN 
    215             WRITE(numout,*) ' hsnif  3 : ', SUM( hsnif  (:,:) ), ' hicnif   : ', SUM( hicif  (:,:) ) 
    216             WRITE(numout,*) ' frld   3 : ', SUM( frld   (:,:) ), ' sist     : ', SUM( sist   (:,:) ) 
     214         IF(l_ctl) THEN 
     215            WRITE(numout,*) ' hsnif  3 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
     216            WRITE(numout,*) ' frld   3 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    217217         ENDIF 
    218218 
    … …  
    221221         CALL lim_thd                                                ! Ice thermodynamics ! 
    222222         !                                                           !--------------------! 
    223          IF( l_ctl ) THEN 
    224             WRITE(numout,*) ' hsnif  4 : ', SUM( hsnif  (:,:) ), ' hicnif   : ', SUM( hicif  (:,:) ) 
    225             WRITE(numout,*) ' frld   4 : ', SUM( frld   (:,:) ), ' sist     : ', SUM( sist   (:,:) ) 
     223         IF(l_ctl) THEN 
     224            WRITE(numout,*) ' hsnif  4 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
     225            WRITE(numout,*) ' frld   4 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
    226226         ENDIF 
    227227 

  • trunk/NEMO/LIM_SRC/limadv.F90

    r88 r106  
    223223 
    224224      IF(l_ctl) THEN 
    225          WRITE(numout,*) ' lim_adv_x: psm  ', SUM( psm  ), ' ps0  ', SUM( ps0 ) 
    226          WRITE(numout,*) ' lim_adv_x: psx  ', SUM( psx  ), ' psxx ', SUM( psxx ) 
    227          WRITE(numout,*) ' lim_adv_x: psy  ', SUM( psy  ), ' psyy ', SUM( psyy ) 
    228          WRITE(numout,*) ' lim_adv_x: psxy ', SUM( psxy ) 
     225         WRITE(numout,*) ' lim_adv_x: psm  ', SUM( psm (2:nictl,2:njctl) ), ' ps0  ', SUM( ps0 (2:nictl,2:njctl) ) 
     226         WRITE(numout,*) ' lim_adv_x: psx  ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 
     227         WRITE(numout,*) ' lim_adv_x: psy  ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 
     228         WRITE(numout,*) ' lim_adv_x: psxy ', SUM( psxy(2:nictl,2:njctl) ) 
    229229      ENDIF 
    230230 
    … …  
    423423 
    424424      IF(l_ctl) THEN 
    425          WRITE(numout,*) ' lim_adv_y: psm  ', SUM( psm  ), ' ps0  ', SUM( ps0 ) 
    426          WRITE(numout,*) ' lim_adv_y: psx  ', SUM( psx  ), ' psxx ', SUM( psxx ) 
    427          WRITE(numout,*) ' lim_adv_y: psy  ', SUM( psy  ), ' psyy ', SUM( psyy ) 
    428          WRITE(numout,*) ' lim_adv_y: psxy ', SUM( psxy ) 
     425         WRITE(numout,*) ' lim_adv_y: psm  ', SUM( psm (2:nictl,2:njctl) ), ' ps0  ', SUM( ps0 (2:nictl,2:njctl) ) 
     426         WRITE(numout,*) ' lim_adv_y: psx  ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 
     427         WRITE(numout,*) ' lim_adv_y: psy  ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 
     428         WRITE(numout,*) ' lim_adv_y: psxy ', SUM( psxy(2:nictl,2:njctl) ) 
    429429      ENDIF 
    430430 

  • trunk/NEMO/LIM_SRC/limdyn.F90

    r76 r106  
    9090            i_j1 = 1    
    9191            i_jpj = jpj 
    92             IF(l_ctl)   WRITE(numout,*) 'lim_rhg :    i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     92            IF(l_ctl)   WRITE(numout,*) 'lim_dyn :    i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    9393            CALL lim_rhg( i_j1, i_jpj ) 
    9494 
    … …  
    9898               zind(jj) = SUM( frld (:,jj  ) )   ! = FLOAT(jpj) if ocean everywhere on a j-line 
    9999               zmsk(jj) = SUM( tmask(:,jj,1) )   ! = 0          if land  everywhere on a j-line 
    100    !!i         write(numout,*) narea, 'limrhg' , jj, zind(jj), zmsk(jj) 
     100   !!i         write(numout,*) narea, 'limdyn' , jj, zind(jj), zmsk(jj) 
    101101            END DO 
    102102 
    … …  
    111111               END DO 
    112112               i_j1 = MAX( 1, i_j1-1 ) 
    113                IF(l_ctl)   WRITE(numout,*) 'lim_rhg : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     113               IF(l_ctl)   WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    114114     
    115115               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    122122               END DO 
    123123               i_jpj = MIN( jpj, i_jpj+2 ) 
    124                IF(l_ctl)   WRITE(numout,*) 'lim_rhg : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     124               IF(l_ctl)   WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    125125     
    126126               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    141141               i_jpj = MIN( jpj, i_jpj+2) 
    142142     
    143                IF(l_ctl)   WRITE(numout,*) 'lim_rhg : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     143               IF(l_ctl)   WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    144144     
    145145               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    153153 
    154154         IF(l_ctl) THEN  
    155             WRITE(numout,*) ' lim_dyn  : u_oce ', SUM( u_oce ), ' v_oce ', SUM( v_oce ) 
    156             WRITE(numout,*) ' lim_dyn  : u_ice ', SUM( u_ice ), ' v_ice ', SUM( v_ice ) 
     155            WRITE(numout,*) ' lim_dyn  : u_oce ', SUM( u_oce(2:nictl,2:njctl) ), ' v_oce ', SUM( v_oce(2:nictl,2:njctl) ) 
     156            WRITE(numout,*) ' lim_dyn  : u_ice ', SUM( u_ice(2:nictl,2:njctl) ), ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 
    157157         ENDIF 
    158158          
    … …  
    243243 
    244244      IF(l_ctl) THEN  
    245          WRITE(numout,*) ' lim_dyn  : tio_u ', SUM( tio_u ), ' tio_v ', SUM( tio_v ) 
    246          WRITE(numout,*) ' lim_dyn  : ust2s ', SUM( ust2s ) 
     245         WRITE(numout,*) ' lim_dyn  : tio_u ', SUM( tio_u(2:nictl,2:njctl) ), ' tio_v ', SUM( tio_v(2:nictl,2:njctl) ) 
     246         WRITE(numout,*) ' lim_dyn  : ust2s ', SUM( ust2s(2:nictl,2:njctl) ) 
    247247      ENDIF 
    248248 

  • trunk/NEMO/LIM_SRC/limflx.F90

    r88 r106  
    2828 
    2929   !! * Module variables 
    30    REAL(wp)  ::            &  ! constant values 
     30   REAL(wp)  ::           &  ! constant values 
    3131      epsi16 = 1.e-16  ,  & 
    3232      rzero  = 0.e0    ,  & 
    … …  
    194194      IF(l_ctl) THEN 
    195195         WRITE(numout,*) ' lim_flx  ' 
    196          WRITE(numout,*) ' fsolar ', SUM(fsolar), ' fnsolar', SUM( fnsolar ) 
    197          WRITE(numout,*) ' fmass  ', SUM(fmass ), ' fsalt  ', SUM(fsalt) 
    198          WRITE(numout,*) ' ftaux  ', SUM(ftaux ), ' ftauy  ', SUM(ftauy) 
    199          WRITE(numout,*) ' freeze ', SUM(freeze), ' tn_ice ', SUM(tn_ice) 
     196         WRITE(numout,*) ' fsolar ', SUM( fsolar(2:nictl,2:njctl) ), ' fnsolar', SUM( fnsolar(2:nictl,2:njctl) ) 
     197         WRITE(numout,*) ' fmass  ', SUM( fmass (2:nictl,2:njctl) ), ' fsalt  ', SUM( fsalt  (2:nictl,2:njctl) ) 
     198         WRITE(numout,*) ' ftaux  ', SUM( ftaux (2:nictl,2:njctl) ), ' ftauy  ', SUM( ftauy  (2:nictl,2:njctl) ) 
     199         WRITE(numout,*) ' freeze ', SUM( freeze(2:nictl,2:njctl) ), ' tn_ice ', SUM( tn_ice (2:nictl,2:njctl) ) 
    200200      ENDIF  
    201201    

  • trunk/NEMO/LIM_SRC/limhdf.F90

    r88 r106  
    167167      ptab(:,:) = ptab(:,:) 
    168168 
    169       IF(l_ctl)   WRITE(numout,*) ' lim_hdf  : ', SUM( ptab-ptab0 ), ' zconv= ', zconv, ' iter= ', iter 
     169      IF(l_ctl)   WRITE(numout,*) ' lim_hdf  : ', SUM( ptab(2:nictl,2:njctl)-ptab0(2:nictl,2:njctl) ), & 
     170      &                                        ' zconv= ', zconv, ' iter= ', iter 
    170171 
    171172   END SUBROUTINE lim_hdf 

  • trunk/NEMO/LIM_SRC/limrhg.F90

    r77 r106  
    2727 
    2828   !! * Module variables 
    29    REAL(wp)  ::            &  ! constant values 
     29   REAL(wp)  ::           &  ! constant values 
    3030      rzero   = 0.e0   ,  & 
    3131      rone    = 1.e0 
    … …  
    631631 
    632632      IF(l_ctl) THEN 
    633          WRITE(numout,*) ' lim_rhg  : res= ', zresm, 'iter= ', jter,' u_ice ', SUM( u_ice ) , ' v_ice ', SUM( v_ice ) 
     633         WRITE(numout,*) ' lim_rhg  : res= ', zresm, 'iter= ', jter,' u_ice ', SUM( u_ice(2:nictl,2:njctl) ) , & 
     634         &                                                          ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 
    634635      ENDIF 
    635636 

  • trunk/NEMO/LIM_SRC/limthd.F90

    r70 r106  
    116116         END DO 
    117117      END DO 
    118       IF(l_ctl)   WRITE(numout,*) 'lim_thd  : ', SUM( hsnif(:,:) ) , ' hsnif' 
     118      IF(l_ctl)   WRITE(numout,*) 'lim_thd  : ', SUM( hsnif(2:nictl,2:njctl) ) , ' hsnif' 
    119119       
    120120       
    … …  
    149149      END DO 
    150150      IF(l_ctl) THEN 
    151          WRITE(numout,*) 'lim_thd: hicif : ', SUM( hicif ), ' hsnif  ', SUM( hsnif  ) 
    152          WRITE(numout,*) 'lim_thd: dmgwi : ', SUM( dmgwi ), ' qstoif ', SUM( qstoif ) 
    153          WRITE(numout,*) 'lim_thd: frld  : ', SUM( frld  ) 
     151         WRITE(numout,*) 'lim_thd: hicif : ', SUM( hicif (2:nictl,2:njctl) ) 
     152         WRITE(numout,*) 'lim_thd: hsnif : ', SUM( hsnif (2:nictl,2:njctl) ) 
     153         WRITE(numout,*) 'lim_thd: dmgwi : ', SUM( dmgwi (2:nictl,2:njctl) ) 
     154         WRITE(numout,*) 'lim_thd: qstoif: ', SUM( qstoif(2:nictl,2:njctl) ) 
     155         WRITE(numout,*) 'lim_thd: frld  : ', SUM( frld  (2:nictl,2:njctl) ) 
    154156      ENDIF 
    155157 
    … …  
    218220      END DO 
    219221      IF(l_ctl) THEN 
    220          WRITE(numout,*) 'lim_thd: pfrld ', SUM( pfrld  ), ' thcm    ', SUM( thcm    ) 
    221          WRITE(numout,*) 'lim_thd: fdtcn ', SUM( fdtcn  ), ' qdtcn   ', SUM( qdtcn  ) 
    222          WRITE(numout,*) 'lim_thd: qldif ', SUM( qldif  ), ' zqlbsbq ', SUM( zqlbsbq ) 
    223          WRITE(numout,*) 'lim_thd: qcmif ', SUM( qcmif  ), ' fbif    ', SUM( fbif    ) 
    224          WRITE(numout,*) 'lim_thd: qcmif ', SUM( qcmif*tms ) 
    225          WRITE(numout,*) 'lim_thd: zhicifp', SUM( zhicifp ) 
     222         WRITE(numout,*) 'lim_thd: pfrld ' , SUM( pfrld  (2:nictl,2:njctl) ), ' thcm    ', SUM( thcm   (2:nictl,2:njctl) ) 
     223         WRITE(numout,*) 'lim_thd: fdtcn ' , SUM( fdtcn  (2:nictl,2:njctl) ), ' qdtcn   ', SUM( qdtcn  (2:nictl,2:njctl) ) 
     224         WRITE(numout,*) 'lim_thd: qldif ' , SUM( qldif  (2:nictl,2:njctl) ), ' zqlbsbq ', SUM( zqlbsbq(2:nictl,2:njctl) ) 
     225         WRITE(numout,*) 'lim_thd: qcmif ' , SUM( qcmif  (2:nictl,2:njctl) ), ' fbif    ', SUM( fbif   (2:nictl,2:njctl) ) 
     226         WRITE(numout,*) 'lim_thd: qcmif  ', SUM( qcmif  (2:nictl,2:njctl)*tms(2:nictl,2:njctl) ) 
     227         WRITE(numout,*) 'lim_thd: zhicifp', SUM( zhicifp(2:nictl,2:njctl) ) 
    226228         WRITE(numout,*) 'limthd : nbpb = ', nbpb 
    227229      ENDIF 
    … …  
    320322       
    321323      IF(l_ctl) THEN 
    322          WRITE(numout,*) 'lim_thd : phicif ', SUM( phicif ), ' hicif ', SUM( hicif ) 
     324         WRITE(numout,*) 'lim_thd : phicif ', SUM( phicif(2:nictl,2:njctl) ), ' hicif ', SUM( hicif(2:nictl,2:njctl) ) 
    323325         WRITE(numout,*) 'lim_thd : nbpac = ', nbpac 
    324326      ENDIF 
    … …  
    375377      IF(l_ctl) THEN 
    376378         WRITE(numout,*) ' lim_thd  end  ' 
    377          WRITE(numout,*) '  hicif ', SUM( hicif  ), '  hsnif ', SUM( hsnif  ) 
    378          WRITE(numout,*) '  frld  ', SUM( frld   ), '  hicifp', SUM( hicifp ) 
    379          WRITE(numout,*) '  phicif', SUM( phicif ), '  pfrld ', SUM( pfrld  ) 
    380          WRITE(numout,*) '  sist  ', SUM( sist   ), '  tbif 1', SUM( tbif  (:,:,1) ) 
    381          WRITE(numout,*) '  tbif 2', SUM( tbif(:,:,2) ), '  tbif 3', SUM( tbif  (:,:,3) ) 
    382          WRITE(numout,*) '  fdtcn ', SUM( fdtcn  ), ' qdtcn ', SUM( qdtcn ) 
    383          WRITE(numout,*) '  qstoif', SUM( qstoif ), '  fsbbq ', SUM( fsbbq  ) 
     379         WRITE(numout,*) '  hicif ', SUM( hicif (2:nictl,2:njctl)   ), '  hsnif ', SUM( hsnif (2:nictl,2:njctl)   ) 
     380         WRITE(numout,*) '  frld  ', SUM( frld  (2:nictl,2:njctl)   ), '  hicifp', SUM( hicifp(2:nictl,2:njctl)  ) 
     381         WRITE(numout,*) '  phicif', SUM( phicif(2:nictl,2:njctl)   ), '  pfrld ', SUM( pfrld (2:nictl,2:njctl)   ) 
     382         WRITE(numout,*) '  sist  ', SUM( sist  (2:nictl,2:njctl)   ), '  tbif 1', SUM( tbif  (2:nictl,2:njctl,1) ) 
     383         WRITE(numout,*) '  tbif 2', SUM( tbif  (2:nictl,2:njctl,2) ), '  tbif 3', SUM( tbif  (2:nictl,2:njctl,3) ) 
     384         WRITE(numout,*) '  fdtcn ', SUM( fdtcn (2:nictl,2:njctl)   ), '  qdtcn ', SUM( qdtcn (2:nictl,2:njctl)  ) 
     385         WRITE(numout,*) '  qstoif', SUM( qstoif(2:nictl,2:njctl)   ), '  fsbbq ', SUM( fsbbq (2:nictl,2:njctl)   ) 
    384386      ENDIF 
    385387 

  • trunk/NEMO/LIM_SRC/limthd_lac.F90

    r88 r106  
    2323 
    2424   !! * Module variables 
    25    REAL(wp)  ::            &  ! constant values 
     25   REAL(wp)  ::           &  ! constant values 
    2626      epsi20 = 1.e-20  ,  & 
    2727      epsi13 = 1.e-13  ,  & 

  • trunk/NEMO/OPA_SRC/DOM/domain.F90

    r93 r106  
    141141      USE ioipsl 
    142142#endif 
    143       NAMELIST/namrun/ no    , cexper   , ln_rstart , nrstdt , nit000,   & 
    144          &             nitend, ndate0   , nleapy   , ninist , nstock,   & 
    145          &             nprint, nwrite   , nrunoff  , ln_ctl , nbench 
     143      NAMELIST/namrun/ no    , cexper   , ln_rstart , nrstdt , nit000,         & 
     144         &             nitend, ndate0   , nleapy   , ninist , nstock,          & 
     145         &             nprint, nwrite   , nrunoff  , ln_ctl , nictl , njctl,   & 
     146         &             nbench 
    146147      NAMELIST/namdom/ ntopo , e3zps_min, e3zps_rat, ngrid  , nmsh  ,   & 
    147148         &             nacc  , atfp     , rdt      , rdtmin , rdtmax,   & 
    … …  
    176177         WRITE(numout,*) '           runoff option                   nrunoff   = ', nrunoff 
    177178         WRITE(numout,*) '           run control (for debugging)     ln_ctl    = ', ln_ctl 
     179         WRITE(numout,*) '           Max i indice for SUM control    nictl     = ', nictl 
     180         WRITE(numout,*) '           Max j indice for SUM control    njctl     = ', njctl 
    178181         WRITE(numout,*) '           benchmark parameter (0/1)       nbench    = ', nbench 
    179182      ENDIF 
    … …  
    232235         WRITE(numout,*) '           nb of seconds per year   raass = ', raass, ' s' 
    233236         WRITE(numout,*) '           nb of seconds per month  rmoss = ', rmoss, ' s' 
     237      ENDIF 
     238 
     239! ... Control the Max i and j indices used for the SUM control (i.e. when ln_ctl=.true.) 
     240      IF ( nictl > jpim1 ) THEN 
     241          IF(lwp) THEN 
     242             WRITE(numout,cform_war) 
     243             WRITE(numout,*) '           nictl = ', nictl, ' must be <= to jpim1 ' 
     244             WRITE(numout,*) '           nictl forced to be equal to jpim1 ' 
     245          ENDIF 
     246          nwarn = nwarn + 1 
     247          nictl = jpim1 
     248      ENDIF 
     249 
     250      IF ( njctl > jpjm1 ) THEN 
     251          IF(lwp) THEN 
     252             WRITE(numout,cform_war) 
     253             WRITE(numout,*) '           njctl = ', njctl, ' must be <= to jpjm1 ' 
     254             WRITE(numout,*) '           njctl forced to be equal to jpjm1 ' 
     255          ENDIF 
     256          nwarn = nwarn + 1 
     257          njctl = jpjm1 
    234258      ENDIF 
    235259 
    … …  
    257281         WRITE(numout,*) '           depth variation tracer step    rdth      = ', rdth 
    258282         IF( lk_ice_lim ) THEN 
    259             WRITE(numout,*) '           ice model coupling frequency      nfice  = ', nfice 
     283             IF(lwp)WRITE(numout,*) '           ice model coupling frequency      nfice  = ', nfice 
    260284            nfbulk = nfice 
    261285         ENDIF 
    262          WRITE(numout,*) '           bulk computation frequency       nfbulk  = ', nfbulk, ' = nfice if ice model used' 
    263          WRITE(numout,*) '           flag closed sea or not           nclosea = ', nclosea 
     286         IF(lwp)WRITE(numout,*) '           bulk computation frequency       nfbulk  = ', nfbulk, ' = nfice if ice model used' 
     287         IF(lwp)WRITE(numout,*) '           flag closed sea or not           nclosea = ', nclosea 
    264288      ENDIF 
    265289 

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

    r84 r106  
    163163 
    164164      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    165          zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    166          zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     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) ) 
    167167         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    168168         u_ctl = zuap   ;   v_ctl = zvap 
    … …  
    318318 
    319319      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    320          zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    321          zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     320         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     321         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    322322         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    323323         u_ctl = zuap   ;   v_ctl = zvap 
    … …  
    429429 
    430430      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    431          zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    432          zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     431         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     432         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    433433         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    434434         u_ctl = zuap   ;   v_ctl = zvap 

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

    r84 r106  
    151151 
    152152      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    153          zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    154          zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     153         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     154         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    155155         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    156156         u_ctl = zuap   ;   v_ctl = zvap 
    … …  
    297297 
    298298      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) ) 
     299         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     300         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    301301         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    302302         u_ctl = zuap   ;   v_ctl = zvap 
    … …  
    411411 
    412412      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    413          WRITE(numout,*) ' hpg  - Ua: ', SUM(ua*umask), ' Va: ', SUM(va*vmask) 
    414          zuap = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    415          zvap = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     413         WRITE(numout,*) ' hpg  - Ua: ', SUM(ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)) 
     414         WRITE(numout,*) '      - Va: ', SUM(va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
     415         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     416         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    416417         WRITE(numout,*) ' hpg  - Ua: ', zuap-u_ctl, ' Va: ', zvap-v_ctl 
    417418         u_ctl = zuap   ;   v_ctl = zvap 

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

    r84 r106  
    120120 
    121121      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    122          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    123          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     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) ) 
    124124         WRITE(numout,*) ' keg  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    125125         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    221221 
    222222      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) ) 
     223         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     224         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    225225         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    226226         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    129129 
    130130      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) ) 
     131         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     132         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    133133         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    134134         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    119119 
    120120      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    121          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    122          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     121         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     122         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    123123         WRITE(numout,*) ' ldf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    124124         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    133133      !                                                ! =============== 
    134134 
    135       IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un*umask), ' Vn: ', SUM(vn*vmask) 
     135      IF(l_ctl)   WRITE(numout,*) ' nxt  - Un: ', SUM(un(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 
     136      &                                  ' Vn: ', SUM(vn(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
    136137 
    137138   END SUBROUTINE dyn_nxt 

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

    r84 r106  
    266266      ! --------------------- 
    267267      rnorme =0. 
    268       DO jj = 1, jpj 
    269          DO ji = 1, jpi 
    270             zgwgt  = gcdmat(ji,jj) * gcb(ji,jj) 
    271             rnorme = rnorme + gcb(ji,jj) * zgwgt 
    272          END DO 
    273       END DO 
     268      rnorme = SUM( gcb(:,:) * gcdmat(:,:) * gcb(:,:) ) 
    274269      IF( lk_mpp )   CALL mpp_sum( rnorme )   ! sum over the global domain 
    275270 
    … …  
    348343 
    349344      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    350          WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:jpim1,2:jpjm1,1:jpkm1)*umask(2:jpim1,2:jpjm1,1:jpkm1) ),   & 
    351             &                   ' Va: ', SUM( va(2:jpim1,2:jpjm1,1:jpkm1)*vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     345         WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ),   & 
     346            &                   ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 
    352347      ENDIF 
    353348 
    … …  
    382377      CALL lbc_lnk( sshn, 'T', 1. ) 
    383378 
     379      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     380         WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 
     381      ENDIF 
     382 
     383 
    384384   END SUBROUTINE dyn_spg_fsc 
    385385 

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

    r84 r106  
    363363 
    364364         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) ) 
     365            WRITE(numout,*) ' spg  - Ua: ', SUM( ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1) ),   & 
     366               &                   ' Va: ', SUM( va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1) ) 
    367367         ENDIF 
    368368 
    … …  
    395395      CALL lbc_lnk( sshn, 'T', 1. ) 
    396396 
     397      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
     398         WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) ) 
     399      ENDIF 
     400 
     401 
    397402   END SUBROUTINE dyn_spg_fsc_atsk 
    398403 

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

    r52 r106  
    146146 
    147147      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    148          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    149          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     148         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     149         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    150150         WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    151151         u_ctl = zua   ;   v_ctl = zva 
    … …  
    261261 
    262262      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    263          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    264          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     263         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     264         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    265265         WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    266266         u_ctl = zua   ;   v_ctl = zva 
    … …  
    395395 
    396396      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    397          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    398          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     397         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     398         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    399399         WRITE(numout,*) ' vor  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    400400         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    127127 
    128128      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    129          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    130          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     129         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     130         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    131131         WRITE(numout,*) ' zad  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    132132         u_ctl = zua   ;   v_ctl = zva 
    … …  
    228228 
    229229      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    230          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    231          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     230         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     231         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    232232         WRITE(numout,*) ' zad  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    233233         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    343343 
    344344      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    345          zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    346          zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     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) ) 
    347347         WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    348348         u_ctl = zua   ;   v_ctl = zva 

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

    r84 r106  
    314314 
    315315      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) ) 
     316         zua = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
     317         zva = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 
    318318         WRITE(numout,*) ' zdf  - Ua: ', zua-u_ctl, ' Va: ', zva-v_ctl 
    319319         u_ctl = zua   ;   v_ctl = zva 

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

    r3 r106  
    7979 
    8080      IF(l_ctl) THEN                ! print mean trends (used for debugging) 
    81          WRITE(numout,*) ' w**2 -   : ', SUM( wn(:,:,:) * wn   (:,:,:) ) 
     81         WRITE(numout,*) ' w**2 -   : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn   (2:nictl,2:njctl,1:jpkm1) ) 
    8282      ENDIF 
    8383 
    … …  
    130130 
    131131      IF(l_ctl) THEN                ! print mean trends (used for debugging) 
    132          WRITE(numout,*) ' w**2 -   : ', SUM( wn(:,:,:) * wn   (:,:,:) ) 
     132         WRITE(numout,*) ' w**2 -   : ', SUM( wn(2:nictl,2:njctl,1:jpkm1) * wn   (2:nictl,2:njctl,1:jpkm1) ) 
    133133      ENDIF 
    134134 

  • trunk/NEMO/OPA_SRC/LDF/ldfeiv.F90

    r80 r106  
    158158      CALL lbc_lnk( aeiv, 'V', 1. ) 
    159159 
    160       IF(l_ctl)   WRITE(numout,*) ' eiv  - u: ', SUM( aeiu(:,:) ), ' v: ', SUM( aeiv(:,:) ) 
     160      IF(l_ctl)   WRITE(numout,*) ' eiv  - u: ', SUM( aeiu(1:nictl+1,1:njctl+1) ), ' v: ', SUM( aeiv(1:nictl+1,1:njctl+1) ) 
    161161 
    162162 
    … …  
    174174            END DO 
    175175         END DO 
    176          IF(l_ctl)   WRITE(numout,*) ' aht  - u: ', SUM( ahtu(:,:) ), ' v: ', SUM( ahtv(:,:) ), ' w: ', SUM( ahtw(:,:) ) 
     176         IF(l_ctl) THEN 
     177            WRITE(numout,*)' aht  - u: ', SUM( ahtu(1:nictl+1,1:njctl+1) ) 
     178            WRITE(numout,*)'      - v: ', SUM( ahtv(1:nictl+1,1:njctl+1) ) 
     179            WRITE(numout,*)'      - w: ', SUM( ahtw(1:nictl+1,1:njctl+1) ) 
     180         ENDIF 
    177181      ENDIF 
    178182 
    … …  
    319323      CALL lbc_lnk( aeiv, 'V', 1. ) 
    320324 
    321       IF(l_ctl)   WRITE(numout,*) ' eiv  - u: ', SUM( aeiu(:,:) ), ' v: ', SUM( aeiv(:,:) ) 
     325      IF(l_ctl)   WRITE(numout,*) ' eiv  - u: ', SUM( aeiu(1:nictl+1,1:njctl+1) ), ' v: ', SUM( aeiv(1:nictl+1,1:njctl+1) ) 
    322326 
    323327      ! ORCA R05: add a space variation on aht (=aeiv except at the equator and river mouth) 
    … …  
    334338            END DO 
    335339         END DO 
    336          IF(l_ctl)   WRITE(numout,*) ' aht  - u: ', SUM( ahtu(:,:) ), ' v: ', SUM( ahtv(:,:) ), ' w: ', SUM( ahtw(:,:) ) 
     340         IF(l_ctl) THEN 
     341            WRITE(numout,*) ' aht  - u: ', SUM( ahtu(1:nictl+1,1:njctl+1) ) 
     342            WRITE(numout,*) '      - v: ', SUM( ahtv(1:nictl+1,1:njctl+1) ) 
     343            WRITE(numout,*) '      - w: ', SUM( ahtw(1:nictl+1,1:njctl+1) ) 
     344         ENDIF 
    337345      ENDIF 
    338346 

  • trunk/NEMO/OPA_SRC/LDF/ldfslp.F90

    r49 r106  
    479479 
    480480      IF(l_ctl) THEN 
    481          WRITE(numout,*) ' slp  - u : ', SUM( uslp (:,:,:) ), ' v : ', SUM( vslp (:,:,:) ) 
    482          WRITE(numout,*) '      - wi: ', SUM( wslpi(:,:,:) ), ' wj: ', SUM( wslpj(:,:,:) ) 
     481         WRITE(numout,*) ' slp  - u : ', SUM( uslp (1:nictl,1:njctl,:) ), ' v : ', SUM( vslp (1:nictl,1:njctl,:) ) 
     482         WRITE(numout,*) '      - wi: ', SUM( wslpi(1:nictl,1:njctl,:) ), ' wj: ', SUM( wslpj(1:nictl,1:njctl,:) ) 
    483483      ENDIF 
    484484 

  • trunk/NEMO/OPA_SRC/SBC/bulk.F90

    r84 r106  
    9797         IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    9898            WRITE(numout,*) ' Forcings ' 
    99             WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (:,:) * tmask(:,:,1) ) 
    100             WRITE(numout,*) ' qsr_ice  : ', SUM( qsr_ice (:,:) * tmask(:,:,1) ) 
    101             WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(:,:) * tmask(:,:,1) ) 
    102             WRITE(numout,*) ' qnsr_ice : ', SUM( qnsr_ice(:,:) * tmask(:,:,1) ) 
    103             WRITE(numout,*) ' evap     : ', SUM( evap    (:,:) * tmask(:,:,1) ) 
    104             WRITE(numout,*) ' precip   : ', SUM( tprecip (:,:) * tmask(:,:,1) ) / rday 
    105             WRITE(numout,*) ' Snow     : ', SUM( sprecip (:,:) * tmask(:,:,1) ) / rday 
    106             WRITE(numout,*) ' u-stress : ', SUM( taux    (:,:) * umask(:,:,1) ) 
    107             WRITE(numout,*) ' v-stress : ', SUM( tauy    (:,:) * vmask(:,:,1) ) 
    108             WRITE(numout,*) ' sst      : ', SUM( zsst    (:,:) * tmask(:,:,1) ) 
    109             WRITE(numout,*) ' sss      : ', SUM( zsss    (:,:) * tmask(:,:,1) ) 
     99            WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     100            WRITE(numout,*) ' qsr_ice  : ', SUM( qsr_ice (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     101            WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     102            WRITE(numout,*) ' qnsr_ice : ', SUM( qnsr_ice(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     103            WRITE(numout,*) ' evap     : ', SUM( evap    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     104            WRITE(numout,*) ' precip   : ', SUM( tprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday 
     105            WRITE(numout,*) ' Snow     : ', SUM( sprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday 
     106            WRITE(numout,*) ' u-stress : ', SUM( taux    (1:nictl+1,1:njctl+1) * umask(1:nictl+1,1:njctl+1,1) ) 
     107            WRITE(numout,*) ' v-stress : ', SUM( tauy    (1:nictl+1,1:njctl+1) * vmask(1:nictl+1,1:njctl+1,1) ) 
     108            WRITE(numout,*) ' sst      : ', SUM( zsst    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     109            WRITE(numout,*) ' sss      : ', SUM( zsss    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) 
    110110            WRITE(numout,*) 
    111111         ENDIF 

  • trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r74 r106  
    301301 
    302302      IF(l_ctl) THEN 
    303          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    304          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     303         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     304         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    305305         WRITE(numout,*) ' had  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    306306         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    396396 
    397397      IF(l_ctl) THEN 
    398          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    399          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     398         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     399         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    400400         WRITE(numout,*) ' zad  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' centered2' 
    401401         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traadv_cen2_atsk.h90

    r74 r106  
    249249      !                                                ! =============== 
    250250 
    251       IF( l_ctl .AND. lwp ) THEN 
    252          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    253          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     251      IF(l_ctl) THEN 
     252         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     253         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    254254         WRITE(numout,*) ' had  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    255255         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    342342      !                                                ! =============== 
    343343 
    344       IF( l_ctl .AND. lwp ) THEN 
    345          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    346          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     344      IF(l_ctl) THEN 
     345         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     346         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    347347         WRITE(numout,*) ' zad  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' centered2 autotsk' 
    348348         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r97 r106  
    268268 
    269269      IF(l_ctl) THEN 
    270          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    271          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     270         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     271         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    272272         WRITE(numout,*) ' had  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 
    273273         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    394394      END DO 
    395395 
    396       IF( l_ctl .AND. lwp ) THEN 
    397          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    398          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     396      IF(l_ctl) THEN 
     397         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     398         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    399399         WRITE(numout,*) ' zad  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 
    400400         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r97 r106  
    329329 
    330330      IF(l_ctl) THEN 
    331          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    332          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     331         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     332         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    333333         WRITE(numout,*) ' had  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl' 
    334334         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    468468      END DO 
    469469 
    470       IF( l_ctl .AND. lwp ) THEN 
    471          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    472          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     470      IF(l_ctl) THEN 
     471         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     472         zsa = SUM( sa(2:jpim1,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    473473         WRITE(numout,*) ' zad  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' muscl2' 
    474474         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r74 r106  
    260260      END DO 
    261261 
    262       IF( l_ctl .AND. lwp ) THEN 
    263          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    264          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     262      IF(l_ctl) THEN 
     263         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     264         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    265265         WRITE(numout,*) ' zad  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl, ' tvd' 
    266266         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/trabbc.F90

    r32 r106  
    104104         END DO 
    105105#endif 
    106       IF( l_ctl .AND. lwp ) THEN 
    107          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     106 
     107      IF(l_ctl) THEN 
     108         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    108109         WRITE(numout,*) ' bbc  - Ta: ', zta-t_ctl 
    109110         t_ctl = zta 

  • trunk/NEMO/OPA_SRC/TRA/trabbl.F90

    r87 r106  
    340340 
    341341      IF(l_ctl) THEN 
    342          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    343          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     342         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     343         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    344344         WRITE(numout,*) ' bbl  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    345345         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90

    r21 r106  
    407407      END DO 
    408408 
    409       IF( l_ctl .AND. lwp ) THEN 
    410          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    411          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     409      IF(l_ctl) THEN 
     410         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     411         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    412412         WRITE(numout,*) ' bbl  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    413413         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/tradmp.F90

    r61 r106  
    178178 
    179179      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    180          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    181          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     180         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     181         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    182182         WRITE(numout,*) ' dmp  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    183183         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90

    r74 r106  
    232232      !                                             ! =============== 
    233233 
    234       IF( l_ctl ) THEN         ! print mean trends (used for debugging) 
    235          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    236          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     234      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
     235         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     236         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    237237         WRITE(numout,*) ' ldf  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    238238         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traldf_iso_zps.F90

    r63 r106  
    268268      !                                             ! =============== 
    269269 
    270       IF( l_ctl ) THEN         ! print mean trends (used for debugging) 
    271          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    272          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     270      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
     271         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     272         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    273273         WRITE(numout,*) ' ldf  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    274274         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/traldf_lap.F90

    r74 r106  
    142142      !                                                  ! ============= 
    143143 
    144       IF( l_ctl ) THEN         ! print mean trends (used for debugging) 
    145          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    146          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     144      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
     145         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     146         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    147147         WRITE(numout,*) ' ldf  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    148148         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/tranxt.F90

    r87 r106  
    164164 
    165165      IF(l_ctl) THEN         ! print mean field (used for debugging) 
    166          WRITE(numout,*) ' nxt  - Tn: ', SUM(tn*tmask), ' Sn: ', SUM(sn*tmask) 
     166         WRITE(numout,*) ' nxt  - Tn: ', SUM(tn(2:nictl,2:njctl,1:jpkm1)*tmask(2:nictl,2:njctl,1:jpkm1)), & 
     167         &                      ' Sn: ', SUM(sn(2:nictl,2:njctl,1:jpkm1)*tmask(2:nictl,2:njctl,1:jpkm1))  
    167168      ENDIF 
    168169 

  • trunk/NEMO/OPA_SRC/TRA/traqsr.F90

    r87 r106  
    173173 
    174174      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    175          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    176          zta = SUM( ta * tmask )  
     175         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    177176         WRITE(numout,*) ' qsr  - Ta: ', zta-t_ctl 
    178177         t_ctl = zta  

  • trunk/NEMO/OPA_SRC/TRA/trasbc.F90

    r87 r106  
    122122       
    123123      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    124          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    125          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     124         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     125         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    126126         WRITE(numout,*) ' sbc  - Ta: ', zta, ' Sa: ', zsa 
    127127         t_ctl = zta   ;   s_ctl = zsa 

  • trunk/NEMO/OPA_SRC/TRA/trazdf_iso.F90

    r87 r106  
    371371 
    372372      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    373          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    374          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     373         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     374         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    375375         WRITE(numout,*) ' zdf 1- Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    376376         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    534534 
    535535      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    536          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    537          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     536         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     537         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    538538         WRITE(numout,*) ' zdf 2- Ta: ', zta, ' Sa: ', zsa 
    539539      ENDIF 

  • trunk/NEMO/OPA_SRC/TRA/trazdf_iso_vopt.F90

    r67 r106  
    8181 
    8282      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    83          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    84          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
     83         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     84         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
    8585         WRITE(numout,*) ' zdf 1- Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl 
    8686         t_ctl = zta   ;   s_ctl = zsa 
    … …  
    9393 
    9494      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    95          zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    96          zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 
    97          WRITE(numout,*) ' zdf 1- Ta: ', zta, ' Sa: ', zsa 
     95         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     96         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) ) 
     97         WRITE(numout,*) ' zdf 2- Ta: ', zta, ' Sa: ', zsa 
    9898      ENDIF 
    9999 

  • trunk/NEMO/OPA_SRC/ZDF/zdfbfr.F90

    r49 r106  
    168168      CALL lbc_lnk( avmv, 'V', 1. ) 
    169169 
    170       IF( l_ctl )   WRITE(numout,*) ' bfr  u : ', SUM( avmu(:,:,:) ), ' v : ', SUM( avmv(:,:,:) ) 
     170      IF(l_ctl)   WRITE(numout,*) ' bfr  u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 
    171171 
    172172   END SUBROUTINE zdf_bfr 

  • trunk/NEMO/OPA_SRC/ZDF/zdfddm.F90

    r49 r106  
    199199      CALL lbc_lnk( avmv, 'V', 1. ) 
    200200 
    201       IF( l_ctl ) THEN 
    202          WRITE(numout,*) ' ddm  t : ', SUM( avt (:,:,:) ), ' s : ', SUM( avs (:,:,:) ) 
    203          WRITE(numout,*) '      u : ', SUM( avmu(:,:,:) ), ' v : ', SUM( avmv(:,:,:) ) 
     201      IF(l_ctl) THEN 
     202         WRITE(numout,*) ' ddm  t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ), ' s : ', SUM( avs (1:nictl+1,1:njctl+1,:) ) 
     203         WRITE(numout,*) '      u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 
    204204      ENDIF 
    205205       

  • trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r49 r106  
    132132      !                                                ! =============== 
    133133 
    134       IF( l_ctl )   WRITE(numout,*) ' nmln : ', SUM( nmln(:,:) ), ' hmld: ', SUM( hmld(:,:) ) 
     134      IF(l_ctl)   WRITE(numout,*) ' nmln : ', SUM( nmln(1:nictl+1,1:njctl+1) ), ' hmld: ', SUM( hmld(1:nictl+1,1:njctl+1) ) 
    135135 
    136136   END SUBROUTINE zdf_mxl 
    … …  
    175175         IF(lwp) WRITE(numout,*) 'zdf_mxl : mixed layer depth' 
    176176         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
    177          IF(lwp) WRITE(numout,*) 
    178177      ENDIF 
    179178 
    … …  
    234233      END DO 
    235234 
    236       IF( l_ctl )   WRITE(numout,*) ' nmln : ', SUM( nmln(:,:) ), ' hmld: ', SUM( hmld(:,:) ) 
     235      IF(l_ctl)   WRITE(numout,*) ' nmln : ', SUM( nmln(1:nictl+1,1:njctl+1) ), ' hmld: ', SUM( hmld(1:nictl+1,1:njctl+1) ) 
    237236 
    238237   END SUBROUTINE zdf_mxl 

  • trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

    r49 r106  
    578578      CALL lbc_lnk( avt, 'W', 1. ) 
    579579 
    580       IF( l_ctl ) THEN 
    581          WRITE(numout,*) ' tke  e : ', SUM( en  (:,:,:) ), ' t : ', SUM( avt (:,:,:) ) 
    582          WRITE(numout,*) '      u : ', SUM( avmu(:,:,:) ), ' v : ', SUM( avmv(:,:,:) ) 
     580      IF(l_ctl) THEN 
     581         WRITE(numout,*) ' tke  e : ', SUM( en  (1:nictl+1,1:njctl+1,:) ), ' t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ) 
     582         WRITE(numout,*) '      u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 
    583583      ENDIF 
    584584 

  • trunk/NEMO/OPA_SRC/ZDF/zdftke_atsk.h90

    r3 r106  
    476476      CALL lbc_lnk( avt, 'W', 1. ) 
    477477 
     478      IF(l_ctl) THEN 
     479         WRITE(numout,*) ' tke  e : ', SUM( en  (1:nictl+1,1:njctl+1,:) ), ' t : ', SUM( avt (1:nictl+1,1:njctl+1,:) ) 
     480         WRITE(numout,*) '      u : ', SUM( avmu(1:nictl+1,1:njctl+1,:) ), ' v : ', SUM( avmv(1:nictl+1,1:njctl+1,:) ) 
     481      ENDIF 
     482 
     483 
    478484   END SUBROUTINE zdf_tke 

  • trunk/NEMO/OPA_SRC/eosbn2.F90

    r99 r106  
    222222      END SELECT 
    223223 
    224       IF(l_ctl)   WRITE(numout,*) ' eos  : ', SUM( prd(:,:,:) ) 
     224      IF(l_ctl)   WRITE(numout,*) ' eos  : ', SUM( prd(1:nictl+1,1:njctl+1,:) ) 
    225225 
    226226   END SUBROUTINE eos_insitu 
    … …  
    408408      END SELECT 
    409409 
    410       IF(l_ctl)   WRITE(numout,*) ' eos-p: ', SUM( prd(:,:,:) ), ' pot : ', SUM( prhop(:,:,:) ) 
     410      IF(l_ctl)   WRITE(numout,*) ' eos-p: ', SUM( prd(1:nictl+1,1:njctl+1,:) ), ' pot : ', SUM( prhop(1:nictl+1,1:njctl+1,:) ) 
    411411 
    412412   END SUBROUTINE eos_insitu_pot 
    … …  
    584584      END SELECT 
    585585 
    586       IF(l_ctl)   WRITE(numout,*) ' eos2d: ', SUM( prd(1:jpim1,1:jpjm1) ) 
     586      IF(l_ctl)   WRITE(numout,*) ' eos2d: ', SUM( prd(1:nictl+1,1:njctl+1) ) 
    587587 
    588588   END SUBROUTINE eos_insitu_2d 
    … …  
    759759      END SELECT 
    760760 
    761       IF(l_ctl)   WRITE(numout,*) ' bn2  : ', SUM( pn2 (:,:,:) ) 
     761      IF(l_ctl)   WRITE(numout,*) ' bn2  : ', SUM( pn2 (1:nictl+1,1:njctl+1,:) ) 
    762762#if defined key_zdfddm 
    763       IF(l_ctl)   WRITE(numout,*) ' rrau : ', SUM( rrau(:,:,:) ) 
     763      IF(l_ctl)   WRITE(numout,*) ' rrau : ', SUM( rrau(1:nictl+1,1:njctl+1,:) ) 
    764764#endif 
    765765 

  • trunk/NEMO/OPA_SRC/in_out_manager.F90

    r100 r106  
    33   USE lib_print         ! formated print library 
    44   USE par_kind 
     5   USE par_oce 
    56 
    67   PUBLIC 
    … …  
    4950 
    5051   INTEGER ::                &  !: 
    51       nstock = 10 ,          &  !: restart file frequency 
    52       nprint =  0 ,          &  !: level of print (0 no print) 
    53       nwrite = 10               !: restart file frequency 
     52      nstock =   10 ,        &  !: restart file frequency 
     53      nprint =    0 ,        &  !: level of print (0 no print) 
     54      nwrite =   10 ,        &  !: restart file frequency 
     55      nictl  = jpim1,        &  !: max i indice for the SUM control 
     56      njctl  = jpjm1            !: max j indice for the SUM control 
    5457   !!---------------------------------------------------------------------- 
    5558   !! logical units 

  • trunk/NEMO/OPA_SRC/step.F90

    r75 r106  
    214214 
    215215      IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    216          WRITE(numout,*) ' emp  -   : ', SUM( emp       * tmask(:,:,1) ) 
    217          WRITE(numout,*) ' emps -   : ', SUM( emps      * tmask(:,:,1) ) 
    218          WRITE(numout,*) ' qt   -   : ', SUM( qt        * tmask(:,:,1) ) 
    219          WRITE(numout,*) ' qsr  -   : ', SUM( qsr       * tmask(:,:,1) ) 
    220          WRITE(numout,*) ' runoff   : ', SUM( runoff    * tmask(:,:,1) ) 
    221          WRITE(numout,*) ' tmask    : ', SUM( tmask                   ) 
    222          WRITE(numout,*) ' sst  -   : ', SUM( tn(:,:,1) * tmask(:,:,1) ) 
    223          WRITE(numout,*) ' sss  -   : ', SUM( sn(:,:,1) * tmask(:,:,1) ) 
    224          WRITE(numout,*) ' tau  - x : ', SUM( taux ), ' - y : ', SUM( tauy ) 
     216         WRITE(numout,*) ' emp  -   : ', SUM( emp   (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
     217         WRITE(numout,*) ' emps -   : ', SUM( emps  (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
     218         WRITE(numout,*) ' qt   -   : ', SUM( qt    (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
     219         WRITE(numout,*) ' qsr  -   : ', SUM( qsr   (1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
     220         WRITE(numout,*) ' runoff   : ', SUM( runoff(1:nictl+1,1:njctl+1)   * tmask(1:nictl+1,1:njctl+1,1) ) 
     221         WRITE(numout,*) ' tmask    : ', SUM( tmask (1:nictl+1,1:njctl+1,:) ) 
     222         WRITE(numout,*) ' sst  -   : ', SUM( tn    (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     223         WRITE(numout,*) ' sss  -   : ', SUM( sn    (1:nictl+1,1:njctl+1,1) * tmask(1:nictl+1,1:njctl+1,1) ) 
     224         WRITE(numout,*) ' tau  - x : ', SUM( taux  (1:nictl+1,1:njctl+1) ), ' - y : ', SUM( tauy(1:nictl+1,1:njctl+1) ) 
    225225      ENDIF 
    226226 
Note: See TracChangeset for help on using the changeset viewer.