Changeset 9075
- Timestamp:
- 2017-12-15T14:50:00+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r9045 r9075 523 523 ! Note that the "unclipped" bottom friction parameter is used even with explicit drag 524 524 IF( ln_wd_il ) THEN 525 zztmp = -1._wp / rdtbt 525 526 DO jj = 2, jpjm1 526 527 DO ji = fs_2, fs_jpim1 ! vector opt. 527 zu_frc(ji,jj) = zu_frc(ji,jj) + MAX(r1_hu_n(ji,jj) * zCdU_u(ji,jj), -1._wp / rdtbt) * zwx(ji,jj) * wdrampu(ji,jj)528 zv_frc(ji,jj) = zv_frc(ji,jj) + MAX(r1_hv_n(ji,jj) * zCdU_v(ji,jj), -1._wp / rdtbt) * zwy(ji,jj) * wdrampv(ji,jj)528 zu_frc(ji,jj) = zu_frc(ji,jj) + MAX(r1_hu_n(ji,jj) * zCdU_u(ji,jj), zztmp ) * zwx(ji,jj) * wdrampu(ji,jj) 529 zv_frc(ji,jj) = zv_frc(ji,jj) + MAX(r1_hv_n(ji,jj) * zCdU_v(ji,jj), zztmp ) * zwy(ji,jj) * wdrampv(ji,jj) 529 530 END DO 530 531 END DO … … 540 541 ! ! Add top stress contribution from baroclinic velocities: 541 542 IF( ln_bt_fw ) THEN 543 zztmp = -1._wp / rdtbt 542 544 DO jj = 2, jpjm1 543 545 DO ji = fs_2, fs_jpim1 ! vector opt.
Note: See TracChangeset
for help on using the changeset viewer.