- Timestamp:
- 2020-01-27T15:31:53+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/sbcblk_algo_coare3p6.F90
r12182 r12340 44 44 45 45 PUBLIC :: SBCBLK_ALGO_COARE3P6_INIT, TURB_COARE3P6 46 !! * Substitutions 47 # include "do_loop_substitute.h90" 46 48 47 49 !! COARE own values for given constants: … … 429 431 !!---------------------------------------------------------------------------------- 430 432 ! 431 DO jj = 1, jpj 432 DO ji = 1, jpi 433 ! 434 zta = pzeta(ji,jj) 435 ! 436 zphi_m = ABS(1. - 15.*zta)**.25 !!Kansas unstable 437 ! 438 zpsi_k = 2.*LOG((1. + zphi_m)/2.) + LOG((1. + zphi_m*zphi_m)/2.) & 439 & - 2.*ATAN(zphi_m) + 0.5*rpi 440 ! 441 zphi_c = ABS(1. - 10.15*zta)**.3333 !!Convective 442 ! 443 zpsi_c = 1.5*LOG((1. + zphi_c + zphi_c*zphi_c)/3.) & 444 & - 1.7320508*ATAN((1. + 2.*zphi_c)/1.7320508) + 1.813799447 445 ! 446 zf = zta*zta 447 zf = zf/(1. + zf) 448 zc = MIN(50._wp, 0.35_wp*zta) 449 zstab = 0.5 + SIGN(0.5_wp, zta) 450 ! 451 psi_m_coare(ji,jj) = (1. - zstab) * ( (1. - zf)*zpsi_k + zf*zpsi_c ) & ! (zta < 0) 452 & - zstab * ( 1. + 1.*zta & ! (zta > 0) 453 & + 0.6667*(zta - 14.28)/EXP(zc) + 8.525 ) ! " 454 ! 455 END DO 456 END DO 433 DO_2D_11_11 434 ! 435 zta = pzeta(ji,jj) 436 ! 437 zphi_m = ABS(1. - 15.*zta)**.25 !!Kansas unstable 438 ! 439 zpsi_k = 2.*LOG((1. + zphi_m)/2.) + LOG((1. + zphi_m*zphi_m)/2.) & 440 & - 2.*ATAN(zphi_m) + 0.5*rpi 441 ! 442 zphi_c = ABS(1. - 10.15*zta)**.3333 !!Convective 443 ! 444 zpsi_c = 1.5*LOG((1. + zphi_c + zphi_c*zphi_c)/3.) & 445 & - 1.7320508*ATAN((1. + 2.*zphi_c)/1.7320508) + 1.813799447 446 ! 447 zf = zta*zta 448 zf = zf/(1. + zf) 449 zc = MIN(50._wp, 0.35_wp*zta) 450 zstab = 0.5 + SIGN(0.5_wp, zta) 451 ! 452 psi_m_coare(ji,jj) = (1. - zstab) * ( (1. - zf)*zpsi_k + zf*zpsi_c ) & ! (zta < 0) 453 & - zstab * ( 1. + 1.*zta & ! (zta > 0) 454 & + 0.6667*(zta - 14.28)/EXP(zc) + 8.525 ) ! " 455 ! 456 END_2D 457 457 ! 458 458 END FUNCTION psi_m_coare … … 482 482 REAL(wp) :: zta, zphi_h, zphi_c, zpsi_k, zpsi_c, zf, zc, zstab 483 483 ! 484 DO jj = 1, jpj 485 DO ji = 1, jpi 486 ! 487 zta = pzeta(ji,jj) 488 ! 489 zphi_h = (ABS(1. - 15.*zta))**.5 !! Kansas unstable (zphi_h = zphi_m**2 when unstable, zphi_m when stable) 490 ! 491 zpsi_k = 2.*LOG((1. + zphi_h)/2.) 492 ! 493 zphi_c = (ABS(1. - 34.15*zta))**.3333 !! Convective 494 ! 495 zpsi_c = 1.5*LOG((1. + zphi_c + zphi_c*zphi_c)/3.) & 496 & -1.7320508*ATAN((1. + 2.*zphi_c)/1.7320508) + 1.813799447 497 ! 498 zf = zta*zta 499 zf = zf/(1. + zf) 500 zc = MIN(50._wp,0.35_wp*zta) 501 zstab = 0.5 + SIGN(0.5_wp, zta) 502 ! 503 psi_h_coare(ji,jj) = (1. - zstab) * ( (1. - zf)*zpsi_k + zf*zpsi_c ) & 504 & - zstab * ( (ABS(1. + 2.*zta/3.))**1.5 & 505 & + .6667*(zta - 14.28)/EXP(zc) + 8.525 ) 506 ! 507 END DO 508 END DO 484 DO_2D_11_11 485 ! 486 zta = pzeta(ji,jj) 487 ! 488 zphi_h = (ABS(1. - 15.*zta))**.5 !! Kansas unstable (zphi_h = zphi_m**2 when unstable, zphi_m when stable) 489 ! 490 zpsi_k = 2.*LOG((1. + zphi_h)/2.) 491 ! 492 zphi_c = (ABS(1. - 34.15*zta))**.3333 !! Convective 493 ! 494 zpsi_c = 1.5*LOG((1. + zphi_c + zphi_c*zphi_c)/3.) & 495 & -1.7320508*ATAN((1. + 2.*zphi_c)/1.7320508) + 1.813799447 496 ! 497 zf = zta*zta 498 zf = zf/(1. + zf) 499 zc = MIN(50._wp,0.35_wp*zta) 500 zstab = 0.5 + SIGN(0.5_wp, zta) 501 ! 502 psi_h_coare(ji,jj) = (1. - zstab) * ( (1. - zf)*zpsi_k + zf*zpsi_c ) & 503 & - zstab * ( (ABS(1. + 2.*zta/3.))**1.5 & 504 & + .6667*(zta - 14.28)/EXP(zc) + 8.525 ) 505 ! 506 END_2D 509 507 ! 510 508 END FUNCTION psi_h_coare
Note: See TracChangeset
for help on using the changeset viewer.