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

Changeset 106 for trunk/NEMO/OPA_SRC/DYN


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/OPA_SRC/DYN
Files:
14 edited

Legend:

Unmodified
Added
Removed

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