Changeset 6140 for trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90
r5836 r6140 81 81 82 82 !! * Substitutions 83 # include "domzgr_substitute.h90"84 83 # include "vectopt_loop_substitute.h90" 85 84 !!---------------------------------------------------------------------- … … 515 514 ! internal Rossby radius Ro = .5 * sum_jpk(N) / f 516 515 zn2 = MAX( rn2b(ji,jj,jk), 0._wp ) 517 zn(ji,jj) = zn(ji,jj) + SQRT( zn2 ) * fse3w(ji,jj,jk)516 zn(ji,jj) = zn(ji,jj) + SQRT( zn2 ) * e3w_n(ji,jj,jk) 518 517 ! Compute elements required for the inverse time scale of baroclinic 519 518 ! eddies using the isopycnal slopes calculated in ldfslp.F : 520 519 ! T^-1 = sqrt(m_jpk(N^2*(r1^2+r2^2)*e3w)) 521 ze3w = fse3w(ji,jj,jk) * tmask(ji,jj,jk)520 ze3w = e3w_n(ji,jj,jk) * tmask(ji,jj,jk) 522 521 zah(ji,jj) = zah(ji,jj) + zn2 * wslp2(ji,jj,jk) * ze3w 523 522 zhw(ji,jj) = zhw(ji,jj) + ze3w … … 533 532 ! internal Rossby radius Ro = .5 * sum_jpk(N) / f 534 533 zn2 = MAX( rn2b(ji,jj,jk), 0._wp ) 535 zn(ji,jj) = zn(ji,jj) + SQRT( zn2 ) * fse3w(ji,jj,jk)534 zn(ji,jj) = zn(ji,jj) + SQRT( zn2 ) * e3w_n(ji,jj,jk) 536 535 ! Compute elements required for the inverse time scale of baroclinic 537 536 ! eddies using the isopycnal slopes calculated in ldfslp.F : 538 537 ! T^-1 = sqrt(m_jpk(N^2*(r1^2+r2^2)*e3w)) 539 ze3w = fse3w(ji,jj,jk) * tmask(ji,jj,jk)538 ze3w = e3w_n(ji,jj,jk) * tmask(ji,jj,jk) 540 539 zah(ji,jj) = zah(ji,jj) + zn2 * ( wslpi(ji,jj,jk) * wslpi(ji,jj,jk) & 541 540 & + wslpj(ji,jj,jk) * wslpj(ji,jj,jk) ) * ze3w … … 711 710 ! 712 711 DO jk = 1, jpkm1 ! e2u e3u u_eiv = -dk[psi_uw] 713 zw3d(:,:,jk) = ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) ) / ( e2u(:,:) * fse3u(:,:,jk) )712 zw3d(:,:,jk) = ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) ) / ( e2u(:,:) * e3u_n(:,:,jk) ) 714 713 END DO 715 714 CALL iom_put( "uoce_eiv", zw3d ) 716 715 ! 717 716 DO jk = 1, jpkm1 ! e1v e3v v_eiv = -dk[psi_vw] 718 zw3d(:,:,jk) = ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) ) / ( e1v(:,:) * fse3v(:,:,jk) )717 zw3d(:,:,jk) = ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) ) / ( e1v(:,:) * e3v_n(:,:,jk) ) 719 718 END DO 720 719 CALL iom_put( "voce_eiv", zw3d )
Note: See TracChangeset
for help on using the changeset viewer.