Ignore:
Timestamp:
2019-05-23T18:36:06+02:00 (2 years ago)
Author:
girrmann
Message:

dev_r10984_HPC-13 : Step 1, boundary is now detected all over the local domain, this does not change the result. Improve bdy treatment for bdy_rnf in bdytra.F90, this changes the result when keyword runoff is specified in namelist

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdydyn2d.F90

    r11044 r11048  
    9999         ii   = idx%nbi(jb,igrd) 
    100100         ij   = idx%nbj(jb,igrd) 
     101         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove 
    101102         zwgt = idx%nbw(jb,igrd) 
    102103         pua2d(ii,ij) = ( pua2d(ii,ij) + zwgt * ( dta%u2d(jb) - pua2d(ii,ij) ) ) * umask(ii,ij,1) 
     
    107108         ii   = idx%nbi(jb,igrd) 
    108109         ij   = idx%nbj(jb,igrd) 
     110         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove 
    109111         zwgt = idx%nbw(jb,igrd) 
    110112         pva2d(ii,ij) = ( pva2d(ii,ij) + zwgt * ( dta%v2d(jb) - pva2d(ii,ij) ) ) * vmask(ii,ij,1) 
     
    163165         ii = idx%nbi(jb,igrd) 
    164166         ij = idx%nbj(jb,igrd) 
     167         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove 
    165168         IF( ll_wd ) THEN 
    166169            spgu(ii, ij) = dta%ssh(jb)  - ssh_ref  
     
    178181         ii  = idx%nbi(jb,igrd) 
    179182         ij  = idx%nbj(jb,igrd)  
     183         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove 
    180184         flagu => idx%flagu(jb,igrd) 
    181185         iim1 = ii + MAX( 0, INT( flagu ) )   ! T pts i-indice inside the boundary 
     
    196200         ii  = idx%nbi(jb,igrd) 
    197201         ij  = idx%nbj(jb,igrd)  
     202         IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj )  CYCLE   ! to remove 
    198203         flagv => idx%flagv(jb,igrd) 
    199204         ijm1 = ij + MAX( 0, INT( flagv ) )   ! T pts j-indice inside the boundary 
Note: See TracChangeset for help on using the changeset viewer.