Changeset 14654
- Timestamp:
- 2021-03-29T16:38:08+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/FLO/flodom.F90
r14644 r14654 373 373 !! 374 374 REAL(dp) :: zabt, zbct, zcdt, zdat, zabpt, zbcpt, zcdpt, zdapt 375 !!---------------------------------------------------------------------376 !! Statement function377 REAL(dp) :: fsline378 REAL(dp) :: psax, psay, psbx, psby, psx, psy379 fsline( psax, psay, psbx, psby, psx, psy ) = psy * ( psbx - psax ) &380 & - psx * ( psby - psay ) &381 & + psax * psby - psay * psbx382 !!---------------------------------------------------------------------383 375 384 376 ! 4 semi plane defined by the 4 points and including the T point … … 413 405 ENDIF 414 406 ! 407 !!--------------------------------------------------------------------- 408 !! Statement function 409 CONTAINS 410 FUNCTION fsline(psax, psay, psbx, psby, psx, psy) 411 REAL(dp), INTENT(IN) :: psax, psay, psbx, psby, psx, psy 412 REAL(dp) :: fsline 413 fsline = psy * ( psbx - psax ) & 414 & - psx * ( psby - psay ) & 415 & + psax * psby - psay * psbx 416 END FUNCTION fsline 417 !!--------------------------------------------------------------------- 415 418 END SUBROUTINE flo_findmesh 416 419
Note: See TracChangeset
for help on using the changeset viewer.