Changeset 5861


Ignore:
Timestamp:
2015-11-05T12:59:33+01:00 (5 years ago)
Author:
gm
Message:

#1612 : correct a bug in dynldf_lap_blp in vvl case: forgot en e3t_b

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90

    r5860 r5861  
    9393                  &        - e1u(ji-1,jj  ) * pub(ji-1,jj  ,jk) + e1u(ji-1,jj-1) * pub(ji-1,jj-1,jk)  )  
    9494               !                                      ! ahm * div        (computed from 2 to jpi/jpj) 
    95                zdiv(ji,jj)     = ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / fse3t(ji,jj,jk) * tmask(ji,jj,jk)                           & 
     95               zdiv(ji,jj)     = ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / fse3t_b(ji,jj,jk) * tmask(ji,jj,jk)                           & 
    9696                  &     * (  e2u(ji,jj)*fse3u_b(ji,jj,jk) * pub(ji,jj,jk) - e2u(ji-1,jj)*fse3u_b(ji-1,jj,jk) * pub(ji-1,jj,jk)  & 
    9797                  &        + e1v(ji,jj)*fse3v_b(ji,jj,jk) * pvb(ji,jj,jk) - e1v(ji,jj-1)*fse3v_b(ji,jj-1,jk) * pvb(ji,jj-1,jk)  ) 
Note: See TracChangeset for help on using the changeset viewer.