Changeset 13497 for NEMO/trunk/src/ICE/icedyn_adv_umx.F90
- Timestamp:
- 2020-09-21T14:37:46+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icedyn_adv_umx.F90
r13472 r13497 1006 1006 ! !-- Laplacian in j-direction --! 1007 1007 DO jl = 1, jpl 1008 DO_2D( 1, 0, 0, 0 ) 1008 DO_2D( 1, 0, 0, 0 ) ! First derivative (gradient) 1009 1009 ztv1(ji,jj,jl) = ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 1010 1010 END_2D 1011 DO_2D( 0, 0, 0, 0 ) 1011 DO_2D( 0, 0, 0, 0 ) ! Second derivative (Laplacian) 1012 1012 ztv2(ji,jj,jl) = ( ztv1(ji,jj,jl) - ztv1(ji,jj-1,jl) ) * r1_e2t(ji,jj) 1013 1013 END_2D … … 1017 1017 ! !-- BiLaplacian in j-direction --! 1018 1018 DO jl = 1, jpl 1019 DO_2D( 1, 0, 0, 0 ) 1019 DO_2D( 1, 0, 0, 0 ) ! First derivative 1020 1020 ztv3(ji,jj,jl) = ( ztv2(ji,jj+1,jl) - ztv2(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 1021 1021 END_2D 1022 DO_2D( 0, 0, 0, 0 ) 1022 DO_2D( 0, 0, 0, 0 ) ! Second derivative 1023 1023 ztv4(ji,jj,jl) = ( ztv3(ji,jj,jl) - ztv3(ji,jj-1,jl) ) * r1_e2t(ji,jj) 1024 1024 END_2D
Note: See TracChangeset
for help on using the changeset viewer.