- Timestamp:
- 2021-03-01T21:25:30+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14122_HPC-08_Mueller_OSMOSIS_streamlining/src/OCE/ZDF/zdfosm.F90
r14564 r14565 1469 1469 ! OSBL is deepening, entrainment > restratification 1470 1470 IF ( zdb_bl(ji,jj) > 1e-15 ) THEN 1471 zgamma_b_nd = MAX( zdbdz_bl_ext(ji,jj), 0.0_wp ) * zdh(ji,jj) / zdb_ml(ji,jj)1471 zgamma_b_nd = MAX( zdbdz_bl_ext(ji,jj), 0.0_wp ) * zdh(ji,jj) / ( zvel_max + MAX( zdb_bl(ji,jj), 1e-15 ) ) 1472 1472 zpsi = ( 1.0_wp - 0.5_wp * zdh(ji,jj) / zhbl(ji,jj) ) * & 1473 1473 & ( zwb0(ji,jj) - MIN( ( zwb_min(ji,jj) + 2.0_wp * zwb_fk_b(ji,jj) ), 0.0_wp ) ) * zdh(ji,jj) / zhbl(ji,jj) … … 1484 1484 zari = MIN( 1.5 * zdb_bl(ji,jj) / ( zhbl(ji,jj) * ( MAX(zdbdz_bl_ext(ji,jj),0._wp) + zdb_bl(ji,jj)**2 / MAX(4.5 * zwstrc(ji,jj)**2 , 1.e-12 )) ), 0.2d0 ) 1485 1485 ENDIF 1486 ! Relaxation to dh_ref = zari * hbl 1487 zddhdt = -a_ddh_2 * ( 1.0 - zdh(ji,jj) / ( zari * zhbl(ji,jj) ) ) * zwb_ent(ji,jj) / zdb_bl(ji,jj) 1486 ! Relaxation to dh_ref = zari * hbl 1487 zddhdt = -1.0_wp * a_ddh_2 * ( 1.0 - zdh(ji,jj) / ( zari * zhbl(ji,jj) ) ) * zwb_ent(ji,jj) / & 1488 & ( zvel_max + MAX( zdb_bl(ji,jj), 1e-15 ) ) 1488 1489 1489 1490 ELSE IF ( j_ddh(ji,jj) == 0 ) THEN 1490 ! Growing shear layer 1491 zddhdt = -1.0_wp * a_ddh * ( 1.0 - zdh(ji,jj) / zhbl(ji,jj) ) * zwb_ent(ji,jj) / zdb_bl(ji,jj) 1491 ! Growing shear layer 1492 zddhdt = -1.0_wp * a_ddh * ( 1.0 - 1.6_wp * zdh(ji,jj) / zhbl(ji,jj) ) * zwb_ent(ji,jj) / & 1493 & ( zvel_max + MAX( zdb_bl(ji,jj), 1e-15 ) ) 1492 1494 zddhdt = EXP( -4.0_wp * ABS( ff_t(ji,jj) ) * zhbl(ji,jj) / MAX(zustar(ji,jj), 1e-8_wp ) ) * zddhdt 1493 1495 ELSE 1494 1496 zddhdt = 0.0_wp 1495 1497 ENDIF ! j_ddh 1496 zdhdt(ji,jj) = zdhdt(ji,jj) + zalpha_b * ( 1.0_wp - 0.5_wp * zdh(ji,jj) / zhbl(ji,jj) ) * zddhdt / ( zvel_max + MAX( zdb_bl(ji,jj), 1e-15 ) ) 1498 zdhdt(ji,jj) = zdhdt(ji,jj) + zalpha_b * ( 1.0_wp - 0.5_wp * zdh(ji,jj) / zhbl(ji,jj) ) * & 1499 & zdb_ml(ji,jj) * MAX( zddhdt, 0.0_wp ) / ( zvel_max + MAX( zdb_bl(ji,jj), 1e-15 ) ) 1497 1500 ELSE ! zdb_bl >0 1498 1501 zdhdt(ji,jj) = -( zwb_ent(ji,jj) + 2.0 * zwb_fk_b(ji,jj) ) / MAX( zvel_max, 1.0e-15) … … 1500 1503 ELSE ! zwb_min + 2*zwb_fk_b < 0 1501 1504 ! OSBL shoaling due to restratification flux. This is the velocity defined in Fox-Kemper et al (2008) 1502 zdhdt(ji,jj) = - zvel_mle(ji,jj)1505 zdhdt(ji,jj) = -1.0_wp * MIN( zvel_mle(ji,jj), hbl(ji,jj) / 10800.0_wp ) 1503 1506 1504 1507 … … 1524 1527 ENDIF 1525 1528 zdhdt(ji,jj) = 2.0 * zdhdt(ji,jj) / MAX(zpert, epsln) 1529 zdhdt(ji,jj) = MAX( zdhdt(ji,jj), -1.0_wp * hbl(ji,jj) / 5400.0_wp ) 1526 1530 ENDIF ! lconv 1527 1531 ELSE ! lshear … … 1547 1551 ELSE 1548 1552 ! OSBL shoaling due to restratification flux. This is the velocity defined in Fox-Kemper et al (2008) 1549 zdhdt(ji,jj) = - zvel_mle(ji,jj)1553 zdhdt(ji,jj) = -1.0_wp * MIN( zvel_mle(ji,jj), hbl(ji,jj) / 10800.0_wp ) 1550 1554 1551 1555 … … 1571 1575 ENDIF 1572 1576 zdhdt(ji,jj) = 2.0 * zdhdt(ji,jj) / MAX(zpert, epsln) 1577 zdhdt(ji,jj) = MAX( zdhdt(ji,jj), -1.0_wp * hbl(ji,jj) / 5400.0_wp ) 1573 1578 ENDIF ! lconv 1574 1579 ENDIF ! lshear
Note: See TracChangeset
for help on using the changeset viewer.