get density FABM standard variable and if needed by FABM allocate memory
subroutine link_density() !! get density FABM standard variable and if needed by FABM !! allocate memory integer rc ! status variable id_density = model%get_interior_variable_id( & fabm_standard_variables%density) if (model%variable_needs_values(id_density)) then allocate(rho_fabm(imt,km,jmt),stat=rc) if (rc /= 0) stop 'link_density(): & Error allocating (rho_fabm)' rho_fabm = 0._rke call model%link_interior_data(id_density,rho_fabm) end if end subroutine link_density