Changeset 11991 for NEMO/branches/2019/ENHANCE-03_domcfg/src/domutl.F90

Timestamp:

2019-11-27T18:32:01+01:00 (4 years ago)

Author:

mathiot

Message:

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

File:

• NEMO/branches/2019/ENHANCE-03_domcfg/src/domutl.F90 (modified) (3 diffs)

NEMO/branches/2019/ENHANCE-03_domcfg/src/domutl.F90

@@ -144 +144 @@
       REAL(wp),                     INTENT(in)    :: rfill          ! filling value
       REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: rdta           ! input data
-      REAL(wp), DIMENSION(jpi,jpj) :: rseedmap, rseedmap_b          ! location of new seed (used for processor exchange)
+      REAL(wp), DIMENSION(jpi,jpj) :: rseedmap                      ! location of new seed (used for processor exchange)
  
       INTEGER :: ii  , ij  , jj, kii, kjj     ! working integer
@@ -175 +175 @@
          DO WHILE ( nseed .NE. 0 )
             DO WHILE ( ASSOCIATED(seed) )
-               ii=seed%i; ij=seed%j ; rseedmap(ii,ij)=1.
+               ii=seed%i; ij=seed%j
                ! update stack size
                CALL del_head_idx(seed) 
@@ -197 +197 @@
             nseed=SUM(rseedmap); IF( lk_mpp )   CALL mpp_sum('domutil', nseed )  ! this is the sum of all the point check this iteration
             !
-            rseedmap_b(:,:)=rseedmap(:,:)
             CALL lbc_lnk('domutil', rseedmap, 'T', 1.)
             !
             ! build new list of seed
+            ! new seed only if data > 0 (ie not land), data not already filled and rseedmap > 0
             DO ii=1,jpi
                DO jj=1,jpj
-                  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
+                  IF (rseedmap(ii,jj) > 0.0 .AND. rdta(ii,jj) > 0 .AND. rdta(ii,jj) /= rfill) THEN
                      CALL add_head_idx(seed, ii, jj, 1)
                   END IF