Ignore:
Timestamp:
2018-11-15T17:27:18+01:00 (23 months ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 2: add generic glob_min/max/sum and locmin/max, complete timing and report (including bdy and icb), see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/DOM/domain.F90

    r10068 r10314  
    469469      !! ** Method  :   compute and print extrema of masked scale factors 
    470470      !!---------------------------------------------------------------------- 
    471       INTEGER ::   iimi1, ijmi1, iimi2, ijmi2, iima1, ijma1, iima2, ijma2 
     471      INTEGER, DIMENSION(2) ::   imi1, imi2, ima1, ima2 
    472472      INTEGER, DIMENSION(2) ::   iloc   !  
    473473      REAL(wp) ::   ze1min, ze1max, ze2min, ze2max 
     
    475475      ! 
    476476      IF(lk_mpp) THEN 
    477          CALL mpp_minloc( e1t(:,:), tmask_i(:,:), ze1min, iimi1,ijmi1 ) 
    478          CALL mpp_minloc( e2t(:,:), tmask_i(:,:), ze2min, iimi2,ijmi2 ) 
    479          CALL mpp_maxloc( e1t(:,:), tmask_i(:,:), ze1max, iima1,ijma1 ) 
    480          CALL mpp_maxloc( e2t(:,:), tmask_i(:,:), ze2max, iima2,ijma2 ) 
     477         CALL mpp_minloc( 'domain', e1t(:,:), tmask_i(:,:), ze1min, imi1 ) 
     478         CALL mpp_minloc( 'domain', e2t(:,:), tmask_i(:,:), ze2min, imi2 ) 
     479         CALL mpp_maxloc( 'domain', e1t(:,:), tmask_i(:,:), ze1max, ima1 ) 
     480         CALL mpp_maxloc( 'domain', e2t(:,:), tmask_i(:,:), ze2max, ima2 ) 
    481481      ELSE 
    482482         ze1min = MINVAL( e1t(:,:), mask = tmask_i(:,:) == 1._wp )     
     
    486486         ! 
    487487         iloc  = MINLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 
    488          iimi1 = iloc(1) + nimpp - 1 
    489          ijmi1 = iloc(2) + njmpp - 1 
     488         imi1(1) = iloc(1) + nimpp - 1 
     489         imi1(2) = iloc(2) + njmpp - 1 
    490490         iloc  = MINLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 
    491          iimi2 = iloc(1) + nimpp - 1 
    492          ijmi2 = iloc(2) + njmpp - 1 
     491         imi2(1) = iloc(1) + nimpp - 1 
     492         imi2(2) = iloc(2) + njmpp - 1 
    493493         iloc  = MAXLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 
    494          iima1 = iloc(1) + nimpp - 1 
    495          ijma1 = iloc(2) + njmpp - 1 
     494         ima1(1) = iloc(1) + nimpp - 1 
     495         ima1(2) = iloc(2) + njmpp - 1 
    496496         iloc  = MAXLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 
    497          iima2 = iloc(1) + nimpp - 1 
    498          ijma2 = iloc(2) + njmpp - 1 
     497         ima2(1) = iloc(1) + nimpp - 1 
     498         ima2(2) = iloc(2) + njmpp - 1 
    499499      ENDIF 
    500500      IF(lwp) THEN 
     
    502502         WRITE(numout,*) 'dom_ctl : extrema of the masked scale factors' 
    503503         WRITE(numout,*) '~~~~~~~' 
    504          WRITE(numout,"(14x,'e1t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1max, iima1, ijma1 
    505          WRITE(numout,"(14x,'e1t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1min, iimi1, ijmi1 
    506          WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, iima2, ijma2 
    507          WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2 
     504         WRITE(numout,"(14x,'e1t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1max, ima1(1), ima1(2) 
     505         WRITE(numout,"(14x,'e1t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1min, imi1(1), imi1(2) 
     506         WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, ima2(1), ima2(2) 
     507         WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, imi2(1), imi2(2) 
    508508      ENDIF 
    509509      ! 
Note: See TracChangeset for help on using the changeset viewer.