Changeset 8464


Ignore:
Timestamp:
2017-08-25T10:31:40+02:00 (3 months ago)
Author:
Tim Graham
Message:

#1934 and #1935 - correct KE calculation in diawri.F90

Location:
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r7494 r8464  
    296296      ! 
    297297      IF ( iom_use("eken") ) THEN 
    298          rke(:,:,jk) = 0._wp                               !      kinetic energy  
     298         rke(:,:,jpk) = 0._wp                               !      kinetic energy  
    299299         DO jk = 1, jpkm1 
    300300            DO jj = 2, jpjm1 
    301301               DO ji = fs_2, fs_jpim1   ! vector opt. 
    302                   zztmp   = 1._wp / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 
    303                   zztmpx  = 0.5 * (  un(ji-1,jj,jk) * un(ji-1,jj,jk) * e2u(ji-1,jj) * fse3u(ji-1,jj,jk)    & 
    304                      &             + un(ji  ,jj,jk) * un(ji  ,jj,jk) * e2u(ji  ,jj) * fse3u(ji  ,jj,jk) )  & 
     302                  zztmp   =  r1_e1e2t(ji,jj) / fse3t(ji,jj,jk)  
     303                  zztmpx  = 0.5 * (  un(ji-1,jj,jk) * un(ji-1,jj,jk) * e1e2u(ji-1,jj) * fse3u(ji-1,jj,jk)    & 
     304                     &             + un(ji  ,jj,jk) * un(ji  ,jj,jk) * e1e2u(ji  ,jj) * fse3u(ji  ,jj,jk) )  & 
    305305                     &          *  zztmp  
    306306                  ! 
    307                   zztmpy  = 0.5 * (  vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1v(ji,jj-1) * fse3v(ji,jj-1,jk)    & 
    308                      &             + vn(ji,jj  ,jk) * vn(ji,jj  ,jk) * e1v(ji,jj  ) * fse3v(ji,jj  ,jk) )  & 
     307                  zztmpy  = 0.5 * (  vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1e2v(ji,jj-1) * fse3v(ji,jj-1,jk)    & 
     308                     &             + vn(ji,jj  ,jk) * vn(ji,jj  ,jk) * e1e2v(ji,jj  ) * fse3v(ji,jj  ,jk) )  & 
    309309                     &          *  zztmp  
    310310                  ! 
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7494 r8464  
    234234      IF( lk_diaar5  )      CALL dia_ar5( kstp )         ! ar5 diag 
    235235      IF( lk_diaharm )      CALL dia_harm( kstp )        ! Tidal harmonic analysis 
    236                             CALL dia_prod( kstp )        ! ocean model: product diagnostics 
     236      IF( lk_diaprod )      CALL dia_prod( kstp )        ! ocean model: product diagnostics 
    237237                            CALL dia_wri( kstp )         ! ocean model: outputs 
    238238      ! 
Note: See TracChangeset for help on using the changeset viewer.