Changeset 14551 for NEMO/branches/2021
- Timestamp:
- 2021-02-26T14:42:53+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
r14550 r14551 2390 2390 zsc_ws_1(:,:) = pws0(A2D(0)) / ( 1.0_wp - 0.56_wp * EXP( phol(A2D(0)) ) ) 2391 2391 WHERE ( ldpyc(A2D(0)) ) ! Pycnocline scales 2392 zsc_wth_pyc(:,:) = -0. 2_wp * pwb0(A2D(0)) * pdt_bl(A2D(0)) / pdb_bl(A2D(0))2393 zsc_ws_pyc(:,:) = -0. 2_wp * pwb0(A2D(0)) * pds_bl(A2D(0)) / pdb_bl(A2D(0))2392 zsc_wth_pyc(:,:) = -0.003_wp * pwstrc(A2D(0)) * ( 1.0_wp - pdh(A2D(0)) / phbl(A2D(0)) ) * pdt_ml(A2D(0)) 2393 zsc_ws_pyc(:,:) = -0.003_wp * pwstrc(A2D(0)) * ( 1.0_wp - pdh(A2D(0)) / phbl(A2D(0)) ) * pds_ml(A2D(0)) 2394 2394 END WHERE 2395 2395 ELSEWHERE … … 2406 2406 & ( 1.0_wp - EXP( -15.0_wp * ( 1.0_wp - zznd_ml ) ) ) 2407 2407 END IF 2408 ! 2409 ! may need to comment out lpyc block 2408 2410 IF ( ldpyc(ji,jj) .AND. ( jk >= kmld(ji,jj) ) .AND. ( jk <= kbld(ji,jj) ) ) THEN ! Pycnocline 2409 2411 zznd_pyc = -1.0_wp * ( gdepw(ji,jj,jk,Kmm) - phbl(ji,jj) ) / pdh(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.