- Timestamp:
- 2020-07-02T16:41:07+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/BDY/bdylib.F90
r12724 r13228 249 249 !!$ zdy_centred = phib(iibm1jp1,ijbm1jp1) - phib(iibm1jm1,ijbm1jm1) 250 250 ! upstream differencing for tangential derivatives 251 zsign_ups = sign( 1. , zdt * zdy_centred )251 zsign_ups = sign( 1.0_wp, zdt * zdy_centred ) 252 252 zsign_ups = 0.5*( zsign_ups + abs(zsign_ups) ) 253 253 zdy = zsign_ups * zdy_1 + (1. - zsign_ups) * zdy_2 … … 257 257 zrx = zdt * zdx / ( zex1 * znor2 ) 258 258 !!$ zrx = min(zrx,2.0_wp) 259 zout = sign( 1. , zrx )259 zout = sign( 1.0_wp, zrx ) 260 260 zout = 0.5*( zout + abs(zout) ) 261 261 zwgt = 2.*rn_Dt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) … … 266 266 & + zwgt * ( phi_ext(jb) - phib(ii,ij) ) ) / ( 1. + zrx ) 267 267 else !! full oblique radiation !! 268 zsign_ups = sign( 1. , zdt * zdy )268 zsign_ups = sign( 1.0_wp, zdt * zdy ) 269 269 zsign_ups = 0.5*( zsign_ups + abs(zsign_ups) ) 270 270 zey = zsign_ups * zey1 + (1.-zsign_ups) * zey2 … … 414 414 !!$ zdy_centred = phib(iibm1jp1,ijbm1jp1,jk) - phib(iibm1jm1,ijbm1jm1,jk) 415 415 ! upstream differencing for tangential derivatives 416 zsign_ups = sign( 1. , zdt * zdy_centred )416 zsign_ups = sign( 1.0_wp, zdt * zdy_centred ) 417 417 zsign_ups = 0.5*( zsign_ups + abs(zsign_ups) ) 418 418 zdy = zsign_ups * zdy_1 + (1. - zsign_ups) * zdy_2 … … 423 423 zrx = zdt * zdx / ( zex1 * znor2 ) 424 424 !!$ zrx = min(zrx,2.0_wp) 425 zout = sign( 1. , zrx )425 zout = sign( 1.0_wp, zrx ) 426 426 zout = 0.5*( zout + abs(zout) ) 427 427 zwgt = 2.*rn_Dt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) … … 432 432 & + zwgt * ( phi_ext(jb,jk) - phib(ii,ij,jk) ) ) / ( 1. + zrx ) 433 433 else !! full oblique radiation !! 434 zsign_ups = sign( 1. , zdt * zdy )434 zsign_ups = sign( 1.0_wp, zdt * zdy ) 435 435 zsign_ups = 0.5*( zsign_ups + abs(zsign_ups) ) 436 436 zey = zsign_ups * zey1 + (1.-zsign_ups) * zey2
Note: See TracChangeset
for help on using the changeset viewer.