- Timestamp:
- 2014-09-29T15:58:28+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r4624 r4794 388 388 ENDIF 389 389 390 #if defined key_surge 391 ! Surface value 392 DO jj = 2, jpjm1 393 DO ji = fs_2, fs_jpim1 ! vector opt. 394 ! hydrostatic pressure gradient along s-surfaces 395 zhpi(ji,jj,1) = zcoef0 / e1u(ji,jj) * ( fse3w(ji+1,jj ,1) * ( znad ) & 396 & - fse3w(ji ,jj ,1) * ( znad ) ) 397 zhpj(ji,jj,1) = zcoef0 / e2v(ji,jj) * ( fse3w(ji ,jj+1,1) * ( znad ) & 398 & - fse3w(ji ,jj ,1) * ( znad ) ) 399 ! s-coordinate pressure gradient correction 400 zuap = -zcoef0 * ( 2._wp * znad ) & 401 & * ( fsde3w(ji+1,jj,1) - fsde3w(ji,jj,1) ) / e1u(ji,jj) 402 zvap = -zcoef0 * ( 2._wp * znad ) & 403 & * ( fsde3w(ji,jj+1,1) - fsde3w(ji,jj,1) ) / e2v(ji,jj) 404 ! add to the general momentum trend 405 ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) + zuap 406 va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) + zvap 407 END DO 408 END DO 409 410 ! interior value (2=<jk=<jpkm1) 411 DO jk = 2, jpkm1 412 DO jj = 2, jpjm1 413 DO ji = fs_2, fs_jpim1 ! vector opt. 414 ! hydrostatic pressure gradient along s-surfaces 415 zhpi(ji,jj,jk) = zhpi(ji,jj,jk-1) + zcoef0 / e1u(ji,jj) & 416 & * ( fse3w(ji+1,jj,jk) * ( 2*znad ) & 417 & - fse3w(ji ,jj,jk) * ( 2*znad ) ) 418 zhpj(ji,jj,jk) = zhpj(ji,jj,jk-1) + zcoef0 / e2v(ji,jj) & 419 & * ( fse3w(ji,jj+1,jk) * ( 2*znad ) & 420 & - fse3w(ji,jj ,jk) * ( 2*znad ) ) 421 ! s-coordinate pressure gradient correction 422 zuap = -zcoef0 * ( 2._wp * znad ) & 423 & * ( fsde3w(ji+1,jj ,jk) - fsde3w(ji,jj,jk) ) / e1u(ji,jj) 424 zvap = -zcoef0 * ( 2._wp * znad ) & 425 & * ( fsde3w(ji ,jj+1,jk) - fsde3w(ji,jj,jk) ) / e2v(ji,jj) 426 ! add to the general momentum trend 427 ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) + zuap 428 va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) + zvap 429 END DO 430 END DO 431 END DO 432 ! 433 #else 390 434 ! Surface value 391 435 DO jj = 2, jpjm1 … … 429 473 END DO 430 474 END DO 475 #endif 431 476 ! 432 477 CALL wrk_dealloc( jpi,jpj,jpk, zhpi, zhpj )
Note: See TracChangeset
for help on using the changeset viewer.