- Timestamp:
- 2018-07-29T11:23:51+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/TRA/trasbc.F90
r9939 r10009 125 125 DO ji = fs_2, fs_jpim1 ! vector opt. 126 126 IF ( ll_wd ) THEN ! If near WAD point limit the flux for now 127 IF ( ssh n(ji,jj) + ht_0(ji,jj) > 2._wp * rn_wdmin1 ) THEN127 IF ( ssh(ji,jj,Nnn) + ht_0(ji,jj) > 2._wp * rn_wdmin1 ) THEN 128 128 sbc_tsc(ji,jj,jp_tem) = r1_rho0_rcp * qns(ji,jj) ! non solar heat flux 129 ELSE IF ( ssh n(ji,jj) + ht_0(ji,jj) > rn_wdmin1 ) THEN129 ELSE IF ( ssh(ji,jj,Nnn) + ht_0(ji,jj) > rn_wdmin1 ) THEN 130 130 sbc_tsc(ji,jj,jp_tem) = r1_rho0_rcp * qns(ji,jj) & 131 & * tanh ( 5._wp * ( ( sshn(ji,jj) + ht_0(ji,jj) -rn_wdmin1 ) * r_rn_wdmin1 ) )131 & * tanh ( 5._wp * ( ( ssh(ji,jj,Nnn) + ht_0(ji,jj) - rn_wdmin1 ) * r_rn_wdmin1 ) ) 132 132 ELSE 133 133 sbc_tsc(ji,jj,jp_tem) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.