Changeset 3301


Ignore:
Timestamp:
2012-02-08T17:56:27+01:00 (9 years ago)
Author:
cbricaud
Message:

undo a bug correction, for ticket 924

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r3300 r3301  
    177177               DO ji = fs_2, fs_jpim1   ! vector opt.    
    178178                  zdir = 0.5 + SIGN( 0.5, pun(ji,jj,jk) )   ! if pun > 0 : zdir = 1 otherwise zdir = 0  
    179                   zdx  = zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) 
     179                  zdx = ( zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) ) * e2u(ji,jj) * fse3u(ji,jj,jk) 
    180180                  zwx(ji,jj,jk)  = ABS( pun(ji,jj,jk) ) * zdt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
    181181                  zfc(ji,jj,jk)  = zdir * ptb(ji  ,jj,jk,jn) + ( 1. - zdir ) * ptb(ji+1,jj,jk,jn)  ! FC in the x-direction for T 
     
    302302               DO ji = fs_2, fs_jpim1   ! vector opt.    
    303303                  zdir = 0.5 + SIGN( 0.5, pvn(ji,jj,jk) )   ! if pun > 0 : zdir = 1 otherwise zdir = 0  
    304                   zdx  = zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) 
     304                  zdx = ( zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) ) * e1v(ji,jj) * fse3v(ji,jj,jk) 
    305305                  zwy(ji,jj,jk)  = ABS( pvn(ji,jj,jk) ) * zdt / zdx    ! (0<zc_cfl<1 : Courant number on x-direction) 
    306306                  zfc(ji,jj,jk)  = zdir * ptb(ji,jj  ,jk,jn) + ( 1. - zdir ) * ptb(ji,jj+1,jk,jn)  ! FC in the x-direction for T 
Note: See TracChangeset for help on using the changeset viewer.