doc for LST equation observation operator routine

HPSCTerrSys · Jun 20, 2024 · e9f7fb8 · e9f7fb8
1 parent 756f38c
commit e9f7fb8
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/bldsva/intf_DA/pdaf/framework/obs_op_pdaf.F90 b/bldsva/intf_DA/pdaf/framework/obs_op_pdaf.F90
@@ -115,6 +115,8 @@ SUBROUTINE obs_op_pdaf(step, dim_p, dim_obs_p, state_p, m_state_p)
   lpointobs = .false.
 
   DO i = 1, dim_obs_p
+    ! Equation for LST computation from Kustas2009, Eq(7)
+    ! http://dx.doi.org/10.1016/j.agrformet.2009.05.016
      m_state_p(i) &
     = (exp(-0.5*clm_paramarr(obs_index_p(i))) &
                      *state_p(obs_index_p(i))**4 &