Changeset 11247


Ignore:
Timestamp:
2019-07-10T18:43:31+02:00 (13 months ago)
Author:
mathiot
Message:

correction of tmask/wmask used into eiv code to take into account ice shelf cavity

Location:
branches/UKMO/dev_isf_gm_UKESM_GO6package_r9321/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_isf_gm_UKESM_GO6package_r9321/NEMOGCM/NEMO/OPA_SRC/LDF/ldfeiv.F90

    r6793 r11247  
    9494                  ! eddies using the isopycnal slopes calculated in ldfslp.F :  
    9595                  ! T^-1 = sqrt(m_jpk(N^2*(r1^2+r2^2)*e3w)) 
    96                   ze3w = fse3w(ji,jj,jk) * tmask(ji,jj,jk) 
     96                  ze3w = fse3w(ji,jj,jk) * wmask(ji,jj,jk) 
    9797                  zah(ji,jj) = zah(ji,jj) + zn2 * wslp2(ji,jj,jk) * ze3w 
    9898                  zhw(ji,jj) = zhw(ji,jj) + ze3w 
     
    112112                  ! eddies using the isopycnal slopes calculated in ldfslp.F :  
    113113                  ! T^-1 = sqrt(m_jpk(N^2*(r1^2+r2^2)*e3w)) 
    114                   ze3w = fse3w(ji,jj,jk) * tmask(ji,jj,jk) 
     114                  ze3w = fse3w(ji,jj,jk) * wmask(ji,jj,jk) 
    115115                  zah(ji,jj) = zah(ji,jj) + zn2 * ( wslpi(ji,jj,jk) * wslpi(ji,jj,jk)   & 
    116116                     &                            + wslpj(ji,jj,jk) * wslpj(ji,jj,jk) ) * ze3w 
  • branches/UKMO/dev_isf_gm_UKESM_GO6package_r9321/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r8447 r11247  
    923923               pn2(ji,jj,jk) = grav * (  zaw * ( pts(ji,jj,jk-1,jp_tem) - pts(ji,jj,jk,jp_tem) )     & 
    924924                  &                    - zbw * ( pts(ji,jj,jk-1,jp_sal) - pts(ji,jj,jk,jp_sal) )  )  & 
    925                   &            / fse3w(ji,jj,jk) * tmask(ji,jj,jk) 
     925                  &            / fse3w(ji,jj,jk) * wmask(ji,jj,jk) 
    926926            END DO 
    927927         END DO 
  • branches/UKMO/dev_isf_gm_UKESM_GO6package_r9321/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90

    r7179 r11247  
    113113         DO jj = 1, jpjm1 
    114114            DO ji = 1, fs_jpim1   ! vector opt. 
    115                zuwk = ( wslpi(ji,jj,jk  ) + wslpi(ji+1,jj,jk  ) ) * fsaeiu(ji,jj,jk  ) * umask(ji,jj,jk  ) 
    116                zuwk1= ( wslpi(ji,jj,jk+1) + wslpi(ji+1,jj,jk+1) ) * fsaeiu(ji,jj,jk+1) * umask(ji,jj,jk+1) 
    117                zvwk = ( wslpj(ji,jj,jk  ) + wslpj(ji,jj+1,jk  ) ) * fsaeiv(ji,jj,jk  ) * vmask(ji,jj,jk  ) 
    118                zvwk1= ( wslpj(ji,jj,jk+1) + wslpj(ji,jj+1,jk+1) ) * fsaeiv(ji,jj,jk+1) * vmask(ji,jj,jk+1) 
     115               zuwk = ( wslpi(ji,jj,jk  ) + wslpi(ji+1,jj,jk  ) ) * fsaeiu(ji,jj,jk  ) * wumask(ji,jj,jk  ) 
     116               zuwk1= ( wslpi(ji,jj,jk+1) + wslpi(ji+1,jj,jk+1) ) * fsaeiu(ji,jj,jk+1) * wumask(ji,jj,jk+1) 
     117               zvwk = ( wslpj(ji,jj,jk  ) + wslpj(ji,jj+1,jk  ) ) * fsaeiv(ji,jj,jk  ) * wvmask(ji,jj,jk  ) 
     118               zvwk1= ( wslpj(ji,jj,jk+1) + wslpj(ji,jj+1,jk+1) ) * fsaeiv(ji,jj,jk+1) * wvmask(ji,jj,jk+1) 
    119119 
    120120               zu_eiv(ji,jj) = 0.5 * umask(ji,jj,jk) * ( zuwk - zuwk1 )  
     
    135135               DO ji = fs_2, fs_jpim1   ! vector opt. 
    136136# if defined key_traldf_c2d || defined key_traldf_c3d 
    137                   zuwi  = ( wslpi(ji,jj,jk)+wslpi(ji-1,jj,jk) ) * fsaeiu(ji-1,jj,jk) * e2u(ji-1,jj) * umask(ji-1,jj,jk) 
    138                   zuwi1 = ( wslpi(ji,jj,jk)+wslpi(ji+1,jj,jk) ) * fsaeiu(ji  ,jj,jk) * e2u(ji  ,jj) * umask(ji  ,jj,jk) 
    139                   zvwj  = ( wslpj(ji,jj,jk)+wslpj(ji,jj-1,jk) ) * fsaeiv(ji,jj-1,jk) * e1v(ji,jj-1) * vmask(ji,jj-1,jk) 
    140                   zvwj1 = ( wslpj(ji,jj,jk)+wslpj(ji,jj+1,jk) ) * fsaeiv(ji,jj  ,jk) * e1v(ji  ,jj) * vmask(ji  ,jj,jk) 
     137                  zuwi  = ( wslpi(ji,jj,jk)+wslpi(ji-1,jj,jk) ) * fsaeiu(ji-1,jj,jk) * e2u(ji-1,jj) * wumask(ji-1,jj,jk) 
     138                  zuwi1 = ( wslpi(ji,jj,jk)+wslpi(ji+1,jj,jk) ) * fsaeiu(ji  ,jj,jk) * e2u(ji  ,jj) * wumask(ji  ,jj,jk) 
     139                  zvwj  = ( wslpj(ji,jj,jk)+wslpj(ji,jj-1,jk) ) * fsaeiv(ji,jj-1,jk) * e1v(ji,jj-1) * wvmask(ji,jj-1,jk) 
     140                  zvwj1 = ( wslpj(ji,jj,jk)+wslpj(ji,jj+1,jk) ) * fsaeiv(ji,jj  ,jk) * e1v(ji  ,jj) * wvmask(ji  ,jj,jk) 
    141141   
    142                   zw_eiv(ji,jj) = - 0.5 * tmask(ji,jj,jk) * ( zuwi1 - zuwi + zvwj1 - zvwj )  
     142                  zw_eiv(ji,jj) = - 0.5 * wmask(ji,jj,jk) * ( zuwi1 - zuwi + zvwj1 - zvwj )  
    143143# else 
    144                   zuwi  = ( wslpi(ji,jj,jk) + wslpi(ji-1,jj,jk) ) * e2u(ji-1,jj) * umask(ji-1,jj,jk) 
    145                   zuwi1 = ( wslpi(ji,jj,jk) + wslpi(ji+1,jj,jk) ) * e2u(ji  ,jj) * umask(ji  ,jj,jk) 
    146                   zvwj  = ( wslpj(ji,jj,jk) + wslpj(ji,jj-1,jk) ) * e1v(ji,jj-1) * vmask(ji,jj-1,jk) 
    147                   zvwj1 = ( wslpj(ji,jj,jk) + wslpj(ji,jj+1,jk) ) * e1v(ji  ,jj) * vmask(ji  ,jj,jk) 
    148  
    149                   zw_eiv(ji,jj) = - 0.5 * tmask(ji,jj,jk) * fsaeiw(ji,jj,jk) * ( zuwi1 - zuwi + zvwj1 - zvwj ) 
     144                  zuwi  = ( wslpi(ji,jj,jk) + wslpi(ji-1,jj,jk) ) * e2u(ji-1,jj) * wumask(ji-1,jj,jk) 
     145                  zuwi1 = ( wslpi(ji,jj,jk) + wslpi(ji+1,jj,jk) ) * e2u(ji  ,jj) * wumask(ji  ,jj,jk) 
     146                  zvwj  = ( wslpj(ji,jj,jk) + wslpj(ji,jj-1,jk) ) * e1v(ji,jj-1) * wvmask(ji,jj-1,jk) 
     147                  zvwj1 = ( wslpj(ji,jj,jk) + wslpj(ji,jj+1,jk) ) * e1v(ji  ,jj) * wvmask(ji  ,jj,jk) 
     148 
     149                  zw_eiv(ji,jj) = - 0.5 * wmask(ji,jj,jk) * fsaeiw(ji,jj,jk) * ( zuwi1 - zuwi + zvwj1 - zvwj ) 
    150150# endif 
    151151                  pwn(ji,jj,jk) = pwn(ji,jj,jk) + zw_eiv(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.