Changeset 4457 for branches/2011
- Timestamp:
- 2014-02-04T13:16:53+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r3211 r4457 405 405 ENDIF 406 406 407 ! Surface value 407 #if ! defined key_z_first 408 ! Surface value 408 409 DO jj = 2, jpjm1 409 410 DO ji = fs_2, fs_jpim1 ! vector opt. … … 423 424 END DO 424 425 END DO 425 426 #endif 427 426 428 ! interior value (2=<jk=<jpkm1) 427 429 #if defined key_z_first 428 430 DO jj = 2, jpjm1 429 431 DO ji = 2, jpim1 430 DO jk = 2, jpkm1 432 433 ! Surface value 434 435 ! hydrostatic pressure gradient along s-surfaces 436 zhpi(ji,jj,1) = zcoef0 / e1u(ji,jj) * ( fse3w(ji+1,jj ,1) * ( znad + rhd(ji+1,jj ,1) ) & 437 & - fse3w(ji ,jj ,1) * ( znad + rhd(ji ,jj ,1) ) ) 438 zhpj(ji,jj,1) = zcoef0 / e2v(ji,jj) * ( fse3w(ji ,jj+1,1) * ( znad + rhd(ji ,jj+1,1) ) & 439 & - fse3w(ji ,jj ,1) * ( znad + rhd(ji ,jj ,1) ) ) 440 ! s-coordinate pressure gradient correction 441 zuap = -zcoef0 * ( rhd (ji+1,jj,1) + rhd (ji,jj,1) + 2._wp * znad ) & 442 & * ( fsde3w(ji+1,jj,1) - fsde3w(ji,jj,1) ) / e1u(ji,jj) 443 zvap = -zcoef0 * ( rhd (ji,jj+1,1) + rhd (ji,jj,1) + 2._wp * znad ) & 444 & * ( fsde3w(ji,jj+1,1) - fsde3w(ji,jj,1) ) / e2v(ji,jj) 445 ! add to the general momentum trend 446 ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) + zuap 447 va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) + zvap 448 449 DO jk = 2, mbkmax(ji,jj)-1 ! jpkm1 431 450 #else 432 451 DO jk = 2, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.