Changeset 3562


Ignore:
Timestamp:
2012-11-15T17:17:48+01:00 (9 years ago)
Author:
rblod
Message:

Fix vvl and flux formulae, see ticket #952

Location:
trunk/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r3294 r3562  
    106106         hdivb(:,:,:) = 0._wp   ;   hdivn(:,:,:) = 0._wp 
    107107         ! 
    108          !                                       ! define e3u_b, e3v_b from e3t_b initialized in domzgr 
    109          CALL dom_vvl_2( nit000, fse3u_b(:,:,:), fse3v_b(:,:,:) ) 
    110          ! 
    111108         IF( cp_cfg == 'eel' ) THEN 
    112109            CALL istate_eel                      ! EEL   configuration : start from pre-defined U,V T-S fields 
     
    133130            ENDDO 
    134131         ENDIF 
     132         !                                       ! define e3u_b, e3v_b from e3t_b initialized in domzgr 
     133         CALL dom_vvl_2( nit000, fse3u_b(:,:,:), fse3v_b(:,:,:) ) 
    135134         !  
    136135      ENDIF 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r3294 r3562  
    248248               DO jk = 1, jpkm1                       ! Leap-Frog - Asselin filter and swap:  
    249249                  DO jj = 1, jpj                      !                   applied on thickness weighted velocity 
    250                      DO ji = 1, jpim1                 !                              --------------------------- 
     250                     DO ji = 1, jpi                   !                              --------------------------- 
    251251                        zue3a = ua(ji,jj,jk) * fse3u_a(ji,jj,jk) 
    252252                        zve3a = va(ji,jj,jk) * fse3v_a(ji,jj,jk) 
     
    268268               fse3u_b(:,:,1:jpkm1) = ze3u_f(:,:,1:jpkm1)      ! e3u_b <-- filtered scale factor 
    269269               fse3v_b(:,:,1:jpkm1) = ze3v_f(:,:,1:jpkm1) 
    270                CALL lbc_lnk( ub, 'U', -1. )                    ! lateral boundary conditions 
    271                CALL lbc_lnk( vb, 'V', -1. ) 
    272270            ENDIF 
    273271            ! 
Note: See TracChangeset for help on using the changeset viewer.