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 181 for trunk/NEMO/OPA_SRC/stpctl.F90 – NEMO

Ignore:
Timestamp:
2004-11-05T15:08:18+01:00 (19 years ago)
Author:
opalod
Message:

CT : UPDATE126 : improve MPI send possiblities with mpi_bsen and mpi_isend; update the search of extremum of scale factors in mpp

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/stpctl.F90

    r79 r181  
    117117      ENDIF 
    118118      IF( zumax >  20.) THEN 
    119          ilocu = MAXLOC( ABS( un(:,:,:) ) ) 
    120          ii = ilocu(1) + nimpp - 1 
    121          ij = ilocu(2) + njmpp - 1 
    122          ik = ilocu(3) 
    123          IF( lk_mpp ) THEN 
    124             CALL mpp_isl( ii ) 
    125             CALL mpp_isl( ij ) 
    126             CALL mpp_isl( ik ) 
     119         IF (lk_mpp ) THEN 
     120            CALL mpp_maxloc(un,umask,zumax,ii,ij,ik) 
     121         ELSE 
     122           ilocu = MAXLOC( ABS( un(:,:,:) ) ) 
     123           ii = ilocu(1) + nimpp - 1 
     124           ij = ilocu(2) + njmpp - 1 
     125           ik = ilocu(3) 
    127126         ENDIF 
    128127         IF(lwp) THEN 
     
    156155      ENDIF 
    157156      IF( zsmin < 0.) THEN  
    158          ilocs = MINLOC( sn(:,:,1), mask = tmask(:,:,1) == 1.e0 ) 
    159          ii = ilocs(1) + nimpp - 1 
    160          ij = ilocs(2) + njmpp - 1 
    161          IF( lk_mpp )   CALL mpp_isl( ii ) 
    162          IF( lk_mpp )   CALL mpp_isl( ij ) 
     157         IF (lk_mpp) THEN 
     158            CALL mpp_minloc ( sn(:,:,1),tmask(:,:,1), zsmin, ii,ij ) 
     159         ELSE 
     160            ilocs = MINLOC( sn(:,:,1), mask = tmask(:,:,1) == 1.e0 ) 
     161            ii = ilocs(1) + nimpp - 1 
     162            ij = ilocs(2) + njmpp - 1 
     163         END IF 
    163164 
    164165         IF(lwp) THEN 
Note: See TracChangeset for help on using the changeset viewer.