Ignore:
Timestamp:
2019-11-27T18:32:01+01:00 (9 months ago)
Author:
mathiot
Message:

ENHANCE-03_closea: rm useless variable, rename some variable and add comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-03_domcfg/src/domutl.F90

    r11604 r11991  
    144144      REAL(wp),                     INTENT(in)    :: rfill          ! filling value 
    145145      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: rdta           ! input data 
    146       REAL(wp), DIMENSION(jpi,jpj) :: rseedmap, rseedmap_b          ! location of new seed (used for processor exchange) 
     146      REAL(wp), DIMENSION(jpi,jpj) :: rseedmap                      ! location of new seed (used for processor exchange) 
    147147  
    148148      INTEGER :: ii  , ij  , jj, kii, kjj     ! working integer 
     
    175175         DO WHILE ( nseed .NE. 0 ) 
    176176            DO WHILE ( ASSOCIATED(seed) ) 
    177                ii=seed%i; ij=seed%j ; rseedmap(ii,ij)=1. 
     177               ii=seed%i; ij=seed%j 
    178178               ! update stack size 
    179179               CALL del_head_idx(seed)  
     
    197197            nseed=SUM(rseedmap); IF( lk_mpp )   CALL mpp_sum('domutil', nseed )  ! this is the sum of all the point check this iteration 
    198198            ! 
    199             rseedmap_b(:,:)=rseedmap(:,:) 
    200199            CALL lbc_lnk('domutil', rseedmap, 'T', 1.) 
    201200            ! 
    202201            ! build new list of seed 
     202            ! new seed only if data > 0 (ie not land), data not already filled and rseedmap > 0 
    203203            DO ii=1,jpi 
    204204               DO jj=1,jpj 
    205                   IF (rseedmap(ii,jj) > 0.0 .AND. rseedmap(ii,jj) /= rseedmap_b(ii,jj) .AND. rdta(ii,jj) > 0 .AND. rdta(ii,jj) /= rfill) THEN 
     205                  IF (rseedmap(ii,jj) > 0.0 .AND. rdta(ii,jj) > 0 .AND. rdta(ii,jj) /= rfill) THEN 
    206206                     CALL add_head_idx(seed, ii, jj, 1) 
    207207                  END IF 
Note: See TracChangeset for help on using the changeset viewer.