my_get_gamma_H_constant.m

 function [gamma_H,model] = my_get_gamma_H_constant(model)
          
         if ~isfield(model,'gamma_H_data')
              
              if model.time_const
                  time_type='time_const';
              else
                  time_type ='time_var';
              end
              
              if model.cone_const
                  model_type = 'cone_const';
              elseif model.vx_const
                  model_type='vx_const';
              else
                  model_type='';
              end
              
              if isfield(model,'base_model')
                  model.k = model.base_model.k;
                  model.name_output_functional = model.base_model.name_output_functional;
              end
              
              
              
              if strcmp(model.name_output_functional, 'box_mean_time')
                  func_type = 'timeInt';
              else
                  func_type = '';
              end
              
              load(['gamma_data_H',func_type,num2str(model.k),'_',model_type,'_',time_type,'_coarse_',num2str(model.coarse_factor),'.mat']);
              
              model.gamma_H_data = gamma_H_data;
          end
          
          mu = get_mu(model);
          
          el_indx = coord2leaf_element(model.gamma_H_data.gamma_H_grid,mu);
          
          gamma_H = model.gamma_H_data.gamma_H_const(el_indx);