Changeset 13899 for NEMO/branches/2020/tickets_icb_1900/src/SWE/diawri.F90
- 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/diawri.F90
r12983 r13899 180 180 IF( iom_use("hu") ) THEN ! water column at u-point 181 181 z2d(:,:) = 0._wp 182 DO_2D _10_10182 DO_2D( 1, 0, 1, 0 ) 183 183 z2d(ji,jj) = 0.5_wp * ( e3t(ji ,jj,1,Kmm) * e1e2t(ji ,jj) & 184 184 & + e3t(ji+1,jj,1,Kmm) * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) … … 190 190 IF( iom_use("hv") ) THEN ! water column at v-point 191 191 z2d(:,:) = 0._wp 192 DO_2D _10_10192 DO_2D( 1, 0, 1, 0 ) 193 193 z2d(ji,jj) = 0.5_wp * ( e3t(ji,jj+1,1,Kmm) * e1e2t(ji,jj+1) & 194 194 & + e3t(ji,jj ,1,Kmm) * e1e2t(ji,jj ) ) * r1_e1e2v(ji,jj) … … 200 200 IF( iom_use("hf") ) THEN ! water column at f-point 201 201 z2d(:,:) = 0._wp 202 DO_2D _10_10202 DO_2D( 1, 0, 1, 0 ) 203 203 z2d(ji,jj) = 0.25_wp * ( e3t(ji,jj+1,1,Kmm) * e1e2t(ji,jj+1) + e3t(ji+1,jj+1,1,Kmm) * e1e2t(ji+1,jj+1) & 204 204 & + e3t(ji,jj ,1,Kmm) * e1e2t(ji,jj ) + e3t(ji+1,jj ,1,Kmm) * e1e2t(ji+1,jj ) ) * r1_e1e2f(ji,jj) … … 216 216 zztmp = rho0 * 0.25 217 217 z2d(:,:) = 0._wp 218 DO_2D _00_00218 DO_2D( 0, 0, 0, 0 ) 219 219 zztmp2 = ( ( rCdU_bot(ji+1,jj)+rCdU_bot(ji ,jj) ) * uu(ji ,jj,mbku(ji ,jj),Kmm) )**2 & 220 220 & + ( ( rCdU_bot(ji ,jj)+rCdU_bot(ji-1,jj) ) * uu(ji-1,jj,mbku(ji-1,jj),Kmm) )**2 & … … 235 235 IF ( iom_use("sKE") ) THEN ! surface kinetic energy at T point 236 236 z2d(:,:) = 0._wp 237 DO_2D _00_00237 DO_2D( 0, 0, 0, 0 ) 238 238 z2d(ji,jj) = 0.25_wp * ( uu(ji ,jj,1,Kmm) * uu(ji ,jj,1,Kmm) * e1e2u(ji ,jj) * e3u(ji ,jj,1,Kmm) & 239 239 & + uu(ji-1,jj,1,Kmm) * uu(ji-1,jj,1,Kmm) * e1e2u(ji-1,jj) * e3u(ji-1,jj,1,Kmm) & … … 251 251 IF ( iom_use("sKEf") ) THEN ! surface kinetic energy at F point 252 252 z2d(:,:) = 0._wp ! CAUTION : only valid in SWE, not with bathymetry 253 DO_2D _00_00253 DO_2D( 0, 0, 0, 0 ) 254 254 z2d(ji,jj) = 0.25_wp * ( uu(ji,jj ,1,Kmm) * uu(ji,jj ,1,Kmm) * e1e2u(ji,jj ) * e3u(ji,jj ,1,Kmm) & 255 255 & + uu(ji,jj+1,1,Kmm) * uu(ji,jj+1,1,Kmm) * e1e2u(ji,jj+1) * e3u(ji,jj+1,1,Kmm) & … … 273 273 z2d(:,:) = 0._wp 274 274 ze3 = 0._wp 275 DO_2D _10_10275 DO_2D( 1, 0, 1, 0 ) 276 276 z2d(ji,jj) = ( e2v(ji+1,jj ) * vv(ji+1,jj ,1,Kmm) - e2v(ji,jj) * vv(ji,jj,1,Kmm) & 277 277 & - e1u(ji ,jj+1) * uu(ji ,jj+1,1,Kmm) + e1u(ji,jj) * uu(ji,jj,1,Kmm) ) * r1_e1e2f(ji,jj) … … 282 282 CALL iom_put( "plavor", ff_f ) ! planetary vorticity ( f ) 283 283 ! 284 DO_2D _10_10284 DO_2D( 1, 0, 1, 0 ) 285 285 ze3 = ( e3t(ji,jj+1,1,Kmm) * e1e2t(ji,jj+1) + e3t(ji+1,jj+1,1,Kmm) * e1e2t(ji+1,jj+1) & 286 286 & + e3t(ji,jj ,1,Kmm) * e1e2t(ji,jj ) + e3t(ji+1,jj ,1,Kmm) * e1e2t(ji+1,jj ) ) * r1_e1e2f(ji,jj) … … 293 293 CALL iom_put( "relpotvor", z2d ) ! relative potential vorticity (zeta/h) 294 294 ! 295 DO_2D _10_10295 DO_2D( 1, 0, 1, 0 ) 296 296 ze3 = ( e3t(ji,jj+1,1,Kmm) * e1e2t(ji,jj+1) + e3t(ji+1,jj+1,1,Kmm) * e1e2t(ji+1,jj+1) & 297 297 & + e3t(ji,jj ,1,Kmm) * e1e2t(ji,jj ) + e3t(ji+1,jj ,1,Kmm) * e1e2t(ji+1,jj ) ) * r1_e1e2f(ji,jj) … … 304 304 CALL iom_put( "abspotvor", z2d ) ! absolute potential vorticity ( q ) 305 305 ! 306 DO_2D _10_10306 DO_2D( 1, 0, 1, 0 ) 307 307 z2d(ji,jj) = 0.5_wp * z2d(ji,jj) * z2d(ji,jj) 308 308 END_2D
Note: See TracChangeset
for help on using the changeset viewer.