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

Changeset 8663


Ignore:
Timestamp:
2017-10-27T12:51:56+02:00 (7 years ago)
Author:
jchanut
Message:

Revert fmask to free slip case at outermost child grid points + slight change in tracer update near bathymetry

Location:
branches/2017/dev_r7963_nemo_v3_6_AGRIF-3_AGRIFVVL/NEMOGCM/NEMO
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7963_nemo_v3_6_AGRIF-3_AGRIFVVL/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90

    r8010 r8663  
    177177      IF (lwp.AND.lk_agrif_debug) Write(*,*) 'Update e3 from grid Number',Agrif_Fixed(), 'Step', Agrif_Nb_Step() 
    178178      ! 
    179       Agrif_UseSpecialValueInUpdate = .FALSE. 
     179      ! JC 
     180!      Agrif_UseSpecialValueInUpdate = .FALSE. 
     181      Agrif_UseSpecialValueInUpdate = .TRUE. 
    180182      Agrif_SpecialValueFineGrid = 0. 
    181183      !  
     
    293295!> jc tmp 
    294296         DO jn = n1,n2 
    295             tabres(i1:i2,j1:j2,k1:k2,jn) =  tabres(i1:i2,j1:j2,k1:k2,jn) * e3t_0(i1:i2,j1:j2,k1:k2) 
     297            tabres(i1:i2,j1:j2,k1:k2,jn) =  tabres(i1:i2,j1:j2,k1:k2,jn) * e3t_0(i1:i2,j1:j2,k1:k2) & 
     298                                         & * tmask(i1:i2,j1:j2,k1:k2) 
    296299         ENDDO 
    297300!< jc tmp 
     
    883886!> jc tmp: 
    884887!         ptab(i1:i2,j1:j2,k1:k2) = fse3t_n(i1:i2,j1:j2,k1:k2) 
    885          ptab(i1:i2,j1:j2,k1:k2) = fse3t_n(i1:i2,j1:j2,k1:k2) / e3t_0(i1:i2,j1:j2,k1:k2) 
     888         ptab(i1:i2,j1:j2,k1:k2) = fse3t_n(i1:i2,j1:j2,k1:k2) / e3t_0(i1:i2,j1:j2,k1:k2) * tmask(i1:i2,j1:j2,k1:k2) 
    886889!< jc tmp: 
    887890      ELSE 
     
    891894         ! 
    892895!> jc tmp: 
    893          ptab(i1:i2,j1:j2,k1:k2) = ptab(i1:i2,j1:j2,k1:k2) * e3t_0(i1:i2,j1:j2,k1:k2) 
     896         DO jk = 1, jpkm1 
     897            DO jj=j1,j2 
     898               DO ji=i1,i2 
     899                  IF (tmask(ji,jj,jk)==1) THEN 
     900                     ptab(ji,jj,jk) = ptab(ji,jj,jk) * e3t_0(ji,jj,jk) 
     901                  ELSE 
     902                     ptab(ji,jj,jk) = e3t_0(ji,jj,jk) 
     903                  ENDIF 
     904               END DO 
     905            END DO 
     906         END DO 
     907!!         ptab(i1:i2,j1:j2,k1:k2) = ptab(i1:i2,j1:j2,k1:k2) * e3t_0(i1:i2,j1:j2,k1:k2) 
    894908!< jc tmp: 
    895909 
  • branches/2017/dev_r7963_nemo_v3_6_AGRIF-3_AGRIFVVL/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r5551 r8663  
    373373            ENDIF 
    374374         END DO 
     375#if defined key_agrif 
     376      IF( .NOT. AGRIF_Root() ) THEN 
     377         IF ((nbondi ==  1).OR.(nbondi == 2)) fmask(nlci-1 , :     ,jk) = 0.e0      ! east 
     378         IF ((nbondi == -1).OR.(nbondi == 2)) fmask(1      , :     ,jk) = 0.e0      ! west 
     379         IF ((nbondj ==  1).OR.(nbondj == 2)) fmask(:      ,nlcj-1 ,jk) = 0.e0      ! north 
     380         IF ((nbondj == -1).OR.(nbondj == 2)) fmask(:      ,1      ,jk) = 0.e0      ! south 
     381      ENDIF 
     382#endif 
    375383      END DO 
    376384      ! 
  • branches/2017/dev_r7963_nemo_v3_6_AGRIF-3_AGRIFVVL/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

    r8010 r8663  
    323323            END DO 
    324324         END DO 
    325  
    326          IF( .NOT. AGRIF_Root() ) THEN 
    327             IF ((nbondi ==  1).OR.(nbondi == 2)) rotn(nlci-1 , :     ,jk) = 0.e0      ! east 
    328             IF ((nbondi == -1).OR.(nbondi == 2)) rotn(1      , :     ,jk) = 0.e0      ! west 
    329             IF ((nbondj ==  1).OR.(nbondj == 2)) rotn(:      ,nlcj-1 ,jk) = 0.e0      ! north 
    330             IF ((nbondj == -1).OR.(nbondj == 2)) rotn(:      ,1      ,jk) = 0.e0      ! south 
    331          ENDIF 
    332325         !                                             ! =============== 
    333326      END DO                                           !   End of slab 
Note: See TracChangeset for help on using the changeset viewer.