- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/SWE/dynldf_lap_blp.F90
r12983 r13899 97 97 DO jk = 1, jpkm1 ! Horizontal slab 98 98 ! 99 DO_2D _01_0199 DO_2D( 0, 1, 0, 1 ) 100 100 ! ! ahm * e3 * curl (computed from 1 to jpim1/jpjm1) 101 101 !!gm open question here : e3f at before or now ? probably now... … … 112 112 END_2D 113 113 ! 114 DO_2D _00_00114 DO_2D( 0, 0, 0, 0 ) 115 115 pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * ( & 116 116 & - ( zcur(ji ,jj) - zcur(ji,jj-1) ) * r1_e2u(ji,jj) / e3u(ji,jj,jk,Kmm) & … … 128 128 DO jk = 1, jpkm1 ! Horizontal slab 129 129 ! 130 DO_2D _01_01130 DO_2D( 0, 1, 0, 1 ) 131 131 ! ! shearing stress component (F-point) NB : ahmf has already been multiplied by fmask 132 132 zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk) & … … 143 143 END_2D 144 144 ! 145 DO_2D _00_00145 DO_2D( 0, 0, 0, 0 ) 146 146 pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm) & 147 147 & * ( ( zten(ji+1,jj ) * e2t(ji+1,jj )*e2t(ji+1,jj ) * e3t(ji+1,jj ,jk,Kmm) & … … 164 164 DO jk = 1, jpkm1 ! Horizontal slab 165 165 ! 166 DO_2D _01_01166 DO_2D( 0, 1, 0, 1 ) 167 167 ! ! shearing stress component (F-point) NB : ahmf has already been multiplied by fmask 168 168 zshe(ji-1,jj-1) = ahmf(ji-1,jj-1,jk) & … … 175 175 END_2D 176 176 ! 177 DO_2D _00_00177 DO_2D( 0, 0, 0, 0 ) 178 178 pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * r1_e1e2u(ji,jj) / e3u(ji,jj,jk,Kmm) & 179 179 & * ( zten(ji+1,jj ) * e2t(ji+1,jj ) * e3t(ji+1,jj ,jk,Kmm) &
Note: See TracChangeset
for help on using the changeset viewer.