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 9971 – NEMO

Changeset 9971


Ignore:
Timestamp:
2018-07-19T14:39:57+02:00 (6 years ago)
Author:
mathiot
Message:

initialization of ztabl, ztabr and ztab to 0 only in case of land suppression as suggested by Sebastien #2118

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/LBC/mpp_nfd_generic.h90

    r9945 r9971  
    9595         ALLOCATE( ztabl(jpimax   ,4,ipk,ipl,ipf) , ztabr(jpimax*jpmaxngh,4,ipk,ipl,ipf) )  
    9696         ! 
    97          ztabr(:,:,:,:,:) = 0._wp 
    98          ztabl(:,:,:,:,:) = 0._wp 
     97         ! when some processors of the north fold are suppressed,  
     98         ! values of ztab* arrays corresponding to these suppressed domain won't be defined  
     99         ! and we need a default definition to 0. 
     100         ! a better test should be: a testing if "suppressed land-processors" belongs to the north-pole folding 
     101         IF ( jpni*jpnj /= jpnij ) THEN 
     102            ztabr(:,:,:,:,:) = 0._wp 
     103            ztabl(:,:,:,:,:) = 0._wp 
     104         END IF 
    99105         ! 
    100106         DO jf = 1, ipf 
     
    183189         ALLOCATE( znorthgloio(jpimax,4,ipk,ipl,ipf,jpni) ) 
    184190         ! 
    185          ztab(:,:,:,:,:)=0._wp 
     191         ! when some processors of the north fold are suppressed, 
     192         ! values of ztab* arrays corresponding to these suppressed domain won't be defined 
     193         ! and we need a default definition to 0. 
     194         ! a better test should be: a testing if "suppressed land-processors" belongs to the north-pole folding 
     195         IF ( jpni*jpnj /= jpnij ) ztab(:,:,:,:,:)=0._wp 
    186196         ! 
    187197         CALL MPI_ALLGATHER( znorthloc  , itaille, MPI_DOUBLE_PRECISION,                & 
Note: See TracChangeset for help on using the changeset viewer.