- Timestamp:
- 2017-06-13T18:50:51+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/GO6_dyn_vrt_diag/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r7845 r8168 76 76 INTEGER, INTENT(in) :: kt ! ocean time-step index 77 77 ! 78 INTEGER :: id_dia_vrt_ldf_int = 1! TODO remove once flags set properly79 INTEGER :: id_dia_vrt_ldf_mean = 1! TODO remove once flags set properly78 INTEGER :: id_dia_vrt_ldf_int = 31 ! TODO remove once flags set properly 79 INTEGER :: id_dia_vrt_ldf_mean = 32 ! TODO remove once flags set properly 80 80 INTEGER :: ji, jj, jk ! dummy loop indices 81 81 REAL(wp) :: zbt, ze2u, ze2v ! temporary scalar … … 194 194 ze2v = e1v(ji,jj) * fse3v(ji,jj,jk) 195 195 ! horizontal biharmonic diffusive trends 196 zua(ji,jj,jk) = - ( zuf(ji ,jj,jk) - zuf(ji,jj-1,jk) ) / ze2u &197 & + ( zut(ji+1,jj,jk) - zut(ji,jj ,jk) ) / e1u(ji,jj) &196 zua(ji,jj,jk) = ( - ( zuf(ji ,jj,jk) - zuf(ji,jj-1,jk) ) / ze2u & 197 & + ( zut(ji+1,jj,jk) - zut(ji,jj ,jk) ) / e1u(ji,jj) ) & 198 198 & * ( fsahmu(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 199 zva(ji,jj,jk) = + ( zuf(ji,jj ,jk) - zuf(ji-1,jj,jk) ) / ze2v &200 & + ( zut(ji,jj+1,jk) - zut(ji ,jj,jk) ) / e2v(ji,jj) &199 zva(ji,jj,jk) = ( + ( zuf(ji,jj ,jk) - zuf(ji-1,jj,jk) ) / ze2v & 200 & + ( zut(ji,jj+1,jk) - zut(ji ,jj,jk) ) / e2v(ji,jj) ) & 201 201 & * ( fsahmv(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 202 202 ! add it to the general momentum trends … … 209 209 END DO ! End of slab 210 210 ! ! =============== 211 IF ( ( id_dia_vrt_ldf_int == 1 ) .or. ( id_dia_vrt_ldf_mean == 1 ) ) THEN 212 CALL dyn_vrt_dia_3d(zua, zva, id_dia_vrt_ldf_int, id_dia_vrt_ldf_mean) 211 IF ( ( id_dia_vrt_ldf_int == 31 ) .or. ( id_dia_vrt_ldf_mean == 32 ) ) THEN 212 ! TODO - remove kt only used for validation 213 CALL dyn_vrt_dia_3d(zua, zva, id_dia_vrt_ldf_int, id_dia_vrt_ldf_mean, kt) 213 214 END IF 214 215 !
Note: See TracChangeset
for help on using the changeset viewer.