Changeset 11888
 Timestamp:
 20191112T10:09:51+01:00 (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icedyn_rhg_evp.F90
r11467 r11888 401 401 CALL lbc_lnk( zp_delt, 'T', 1. ) 402 402 403 DO jj = 1, jpjm1404 DO j i = 1, jpim1405 406 ! alpha & beta for aEVP 407 IF( ln_aEVP ) THEN403 IF( ln_aEVP ) THEN 404 DO jj = 1, jpjm1 405 DO ji = 1, jpim1 406 407 ! alpha & beta for aEVP 408 408 zalph2 = MAX( 50._wp, rpi * SQRT( 0.5_wp * zp_delt(ji,jj) * r1_e1e2t(ji,jj) * zdt_m(ji,jj) ) ) 409 409 z1_alph2 = 1._wp / ( zalph2 + 1._wp ) 410 410 zbeta(ji,jj) = zalph2 411 ENDIF 412 411 412 ! P/delta at F points 413 zp_delf = 0.25_wp * ( zp_delt(ji,jj) + zp_delt(ji+1,jj) + zp_delt(ji,jj+1) + zp_delt(ji+1,jj+1) ) 414 415 ! stress at F points 416 zs12(ji,jj)= ( zs12(ji,jj) * zalph2 + zp_delf * ( zds(ji,jj) * z1_ecc2 ) * 0.5_wp ) * z1_alph2 417 418 END DO 419 END DO 420 ELSE 421 DO jj = 1, jpjm1 422 DO ji = 1, jpim1 423 413 424 ! P/delta at F points 414 425 zp_delf = 0.25_wp * ( zp_delt(ji,jj) + zp_delt(ji+1,jj) + zp_delt(ji,jj+1) + zp_delt(ji+1,jj+1) ) 415 426 416 427 ! stress at F points 417 428 zs12(ji,jj)= ( zs12(ji,jj) * zalph2 + zp_delf * ( zds(ji,jj) * z1_ecc2 ) * 0.5_wp ) * z1_alph2 418 429 419 END DO 420 END DO 430 END DO 431 END DO 432 ENDIF 421 433 422 434 !  Ice internal stresses (Appendix C of Hunke and Dukowicz, 2002)  !
Note: See TracChangeset
for help on using the changeset viewer.