- Timestamp:
- 2017-12-06T09:16:26+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/ROMS_WAD_7832/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r8870 r8907 980 980 ! Add bottom stresses: 981 981 !jth do implicitly instead 982 ! zu_trd(:,:) = zu_trd(:,:) + bfrua(:,:) * un_e(:,:) * hur_e(:,:) 983 ! zv_trd(:,:) = zv_trd(:,:) + bfrva(:,:) * vn_e(:,:) * hvr_e(:,:) 982 IF ( .NOT. ln_wd_il .and. .NOT. ln_wd_dl) THEN ! Revert to explicit for bit comparison tests in non wad runs 983 zu_trd(:,:) = zu_trd(:,:) + bfrua(:,:) * un_e(:,:) * hur_e(:,:) 984 zv_trd(:,:) = zv_trd(:,:) + bfrva(:,:) * vn_e(:,:) * hvr_e(:,:) 985 ENDIF 984 986 ! 985 987 ! Add top stresses: … … 1029 1031 1030 1032 !jth implicit bottom friction: 1031 ua_e(ji,jj) = ua_e(ji,jj) /(1.0 - rdtbt * bfrua(ji,jj) * hur_e(ji,jj)) 1032 va_e(ji,jj) = va_e(ji,jj) /(1.0 - rdtbt * bfrva(ji,jj) * hvr_e(ji,jj)) 1033 IF ( ln_wd_il .or. ln_wd_dl ) THEN ! revert to explicit for bit comparison tests in non wad runs 1034 ua_e(ji,jj) = ua_e(ji,jj) /(1.0 - rdtbt * bfrua(ji,jj) * hur_e(ji,jj)) 1035 va_e(ji,jj) = va_e(ji,jj) /(1.0 - rdtbt * bfrva(ji,jj) * hvr_e(ji,jj)) 1036 ENDIF 1033 1037 1034 1038 END DO
Note: See TracChangeset
for help on using the changeset viewer.