Changeset 4437


Ignore:
Timestamp:
2014-02-04T13:14:54+01:00 (7 years ago)
Author:
trackstand2
Message:

Trim columns in exchs_generic and init mbkmax in mpp_init3

Location:
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC/exchmod.F90

    r4422 r4437  
    88#endif 
    99  USE profile 
    10 !  USE timing, ONLY: timing_start, timing_stop 
     10  USE timing, ONLY: timing_start, timing_stop 
    1111  ! Make some key parameters from mapcomm_mod available to all who  
    1212  ! USE this module 
     
    34863486                           isrcsend, jsrcsend, idesrecv, jdesrecv,     & 
    34873487                           isrcsendp,jsrcsendp,idesrecvp,jdesrecvp,    & 
    3488                            nxrecv,source, iesub, jesub,         & 
     3488                           nxrecv,source, iesub, jesub,                & 
    34893489                           MaxCommDir, MaxComm, idessend, jdessend,    & 
    34903490                           nxsendp, nysendp, nzsendp,                  & 
     
    34953495    USE lib_mpp,     ONLY: mpi_comm_opa 
    34963496#endif 
    3497     USE dom_oce,     ONLY: narea 
     3497    USE dom_oce,     ONLY: narea, mbkmax 
    34983498    USE in_out_manager, ONLY: numout 
    34993499    IMPLICIT none 
     
    35433543 
    35443544    !CALL prof_region_begin(ARPEXCHS_GENERIC, "Exchs_indiv", iprofStat) 
    3545     !CALL timing_start('exchs_generic') 
     3545    CALL timing_start('exchs_generic') 
    35463546 
    35473547    ierr = 0 
     
    37643764                DO j=jstart, jend, 1 
    37653765                   DO i=istart, iend, 1 
    3766                       DO k=1,nzsendp(ipatch,isend,1),1 
     3766                      DO k=1,mbkmax(i,j),1 
     3767                      !DO k=1, nzsendp(ipatch,isend,1),1 
    37673768#else 
    37683769                DO k=1,nzsendp(ipatch,isend,1),1 
     
    38013802                 DO j=jstart, jend, 1 
    38023803                    DO i=istart, iend, 1 
    3803                        DO k=1,nzsendp(ipatch,isend,1),1 
     3804                       DO k=1,mbkmax(i,j),1 ! nzsendp(ipatch,isend,1),1 
    38043805#else 
    38053806                 DO k=1,nzsendp(ipatch,isend,1),1 
     
    39183919                DO j=jstart, jend, 1 
    39193920                   DO i=istart, iend, 1 
    3920                       DO k=1,nzrecvp(ipatch,irecv,1),1 
     3921                      DO k=1,mbkmax(i,j),1 ! nzrecvp(ipatch,irecv,1),1 
     3922                      !DO k=1, nzrecvp(ipatch,irecv,1),1 
    39213923#else 
    39223924                DO k=1,nzrecvp(ipatch,irecv,1),1 
     
    39293931#if defined key_z_first 
    39303932                      ! ARPDBG - wipe anything below the ocean bottom 
    3931                       DO k=nzrecvp(ipatch,irecv,1)+1,jpk,1 
    3932                          b3(i,j,k) = 0.0_wp 
    3933                       END DO 
     3933                      !b3(i,j,nzrecvp(ipatch,irecv,1)+1:jpk) = 0.0_wp 
     3934                      b3(i,j,mbkmax(i,j)+1:jpk) = 0.0_wp 
    39343935#endif 
    39353936                   END DO 
     
    39633964                DO j=jstart, jend, 1 
    39643965                   DO i=istart, iend, 1 
    3965                       DO k=1,nzrecvp(ipatch,irecv,1),1 
     3966                      DO k=1,mbkmax(i,j),1 ! nzrecvp(ipatch,irecv,1),1 
    39663967#else 
    39673968                DO k=1,nzrecvp(ipatch,irecv,1),1 
     
    40824083    IF( ALLOCATED(recvIBuff) )DEALLOCATE(recvIBuff) 
    40834084 
    4084     ! CALL timing_stop('exchs_generic') 
     4085    CALL timing_stop('exchs_generic') 
    40854086    !CALL prof_region_end(ARPEXCHS_GENERIC, iprofStat) 
    40864087 
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/LBC/partition_mod.F90

    r3849 r4437  
    1414                      nperio,      &  ! Local periodicity  
    1515                      nwidthmax,   &  ! Width of widest northern domain 
    16                       narea           ! ID of local area (= rank + 1) 
     16                      narea,       &  ! ID of local area (= rank + 1) 
     17                      mbkmax          ! Deepest level above ocean floor 
    1718#if defined key_mpp_mpi 
    1819   USE lib_mpp,        ONLY: mppsize, mppsync, mpi_comm_opa,                & 
     
    115116      INTEGER :: iproc, jproc                  ! Loop index 
    116117      INTEGER :: ierr                          ! Error flag 
    117       INTEGER :: i,j !,npoints ! ARPDBG for debugging only 
     118      INTEGER :: ji, jj 
    118119      CHARACTER(LEN=4) :: intStr 
    119120 
     
    181182            WRITE (997,FMT='("0 1 0 setrgbcolor")') 
    182183            WRITE (997,FMT='("newpath")') 
    183             DO j = 1,jpjglo,1 
    184                DO i = 1,jpiglo,1 
    185                   IF(imask(i,j) /= 1)THEN 
    186                      WRITE (997,FMT='(I3," u ",I3," u 1 u 0 360 arc closepath")') i, j 
     184            DO jj = 1,jpjglo,1 
     185               DO ji = 1,jpiglo,1 
     186                  IF(imask(ji,jj) /= 1)THEN 
     187                     WRITE (997,FMT='(I3," u ",I3," u 1 u 0 360 arc closepath")') ji, jj 
    187188                     ! Use 'fill' instead of 'stroke' to get filled circles 
    188189                     WRITE (997,FMT='("fill")') 
     
    298299      ! read 
    299300      CALL write_partition() 
     301 
     302      ! Initialise mbkmax because it's needed in any halo swap calls but we didn't have 
     303      ! jpi and jpj until the partitioning was done. 
     304      DO jj=1,jpj,1 
     305         DO ji=1,jpi,1 
     306            mbkmax(ji,jj) = ibotlevel(MIN(jpiglo,ji+nimpp-1), & 
     307                                      MIN(jpjglo,jj+njmpp-1)) 
     308         END DO 
     309      END DO 
    300310 
    301311      ! ARPDBG - test comms setup 
Note: See TracChangeset for help on using the changeset viewer.