ESM1.5 Reconciliation - esm15_ friction

[penguian]

keyword_maygit owner:jxs599@nci.org.au resolution_fixed type_defect | by yh4968

---

see attachment

---

Issue migrated from trac:347 at 2023-11-27 11:43:15 +1100

[penguian]

@yh4968@nci.org.au uploaded file switches -friciton.docx (49.5 KiB)

[penguian]

@rml599@nci.org.au commented

---

https://trac.nci.org.au/svn/cable/branches/Share/ESM15-CABLE3_draft1/science/canopy/cbl_friction_vel.F90IF( cable_runtime%esm15_friction ) THEN psim_1 = psim( zetar(:,iter), mp, CPI_C ) ELSE
psim_1 = psim( zetar(:,iter) * zref_uv/zref_tq, mp, CPI_C ) ENDIF

rescale = CVONK * MAX( ua, CUMIN )

z_eff = zref_uv / z0m

IF( cable_runtime%esm15_friction ) THEN psim_arg = zetar(:,iter) / z_eff ELSE
psim_arg = zetar(:,iter) * z0m / zref_tq ENDIF

psim_2 = psim( psim_arg, mp, CPI_C )

lower_limit = rescale / ( LOG(z_eff) - psim_1 + psim_2 )

IF( cable_runtime%esm15_friction ) THEN friction_vel = MAX(1.e-6, lower_limit ) ELSE
friction_vel= MIN( MAX(1.e-6, lower_limit), 10.0 ) ENDIF ]}

This change described in #138 - consistency of height between t,q,u,v

Given this is a bug fix, accept the change. In short test, seems to have had some impact on temperature (lower) and GPP/NPP (higher) but may just be different meteorology trajectory.

[{

[penguian]

@jxs599@nci.org.au changed status from new to closed

[penguian]

@jxs599@nci.org.au set resolution to fixed

[penguian]

@jxs599@nci.org.au set milestone to 1. Closed

[penguian]

@ccc561@nci.org.au set keywords to maygit