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 for trunk/NEMO/LIM_SRC – NEMO

Changeset 106 for trunk/NEMO/LIM_SRC


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/LIM_SRC
Files:
8 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  ,  & 
Note: See TracChangeset for help on using the changeset viewer.