- Timestamp:
- 2015-12-02T17:12:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r5930 r5983 41 41 USE sbcapr ! surface boundary condition: atmospheric pressure 42 42 USE dynadv, ONLY: ln_dynadv_vec 43 USE sbcwave, ONLY: usd2d, vsd2d 43 44 #if defined key_agrif 44 45 USE agrif_opa_interp ! agrif … … 422 423 ENDIF 423 424 ! 425 ! Add Stokes Coriolis if defined 426 IF ( ln_stcor ) THEN 427 DO jj = 1, jpjm1 428 DO ji = 1, fs_jpim1 ! vector opt. 429 430 zy1 = ff(ji ,jj-1) * ( vsd2d(ji ,jj-1) + vsd2d(ji+1,jj-1) ) 431 zy2 = ff(ji ,jj ) * ( vsd2d(ji ,jj ) + vsd2d(ji+1,jj ) ) 432 zx1 = ff(ji-1,jj ) * ( usd2d(ji-1,jj ) + usd2d(ji-1,jj+1) ) 433 zx2 = ff(ji ,jj ) * ( usd2d(ji ,jj ) + usd2d(ji ,jj+1) ) 434 435 zu_frc(ji,jj) = zu_frc(ji,jj) + 0.25 * ( zy1 + zy2 ) * hur(ji,jj) 436 zv_frc(ji,jj) = zv_frc(ji,jj) - 0.25 * ( zx1 + zx2 ) * hvr(ji,jj) 437 END DO 438 END DO 439 ENDIF 440 ! 424 441 IF ( ln_apr_dyn ) THEN ! Add atm pressure forcing 425 442 IF (ln_bt_fw) THEN
Note: See TracChangeset
for help on using the changeset viewer.