anchor.diff:      elem_field(ne_Vdot,1:num_elems) = &
            elem_field(ne_Vdot,1:num_elems)/elem_field(ne_Vdot,1)
   
  +    close(10)
  +    
       call enter_exit(sub_name,2)
   
     end subroutine evaluate_vent
  
positive.code:                sin(2.0_dp*pi*(0.5_dp+(ttime-Tinsp)/(2.0_dp*Texpn)))/ &
                 (2.0_dp*Texpn)*dt
         endif
         
      case('passive')
         if(ttime.le.Tinsp+0.5_dp*dt)then
            dP_muscle = ventilation_values%P_muscle_estimate*ventilation_values%factor_P_muscle_insp*PI*dt* &
                 sin(pi*ttime/Tinsp)/(2.0_dp*Tinsp)
            sum_dP_muscle = sum_dP_muscle+dP_muscle
            sum_dP_muscle_ei = sum_dP_muscle
         else
            Tpass = 0.1_dp
            dP_muscle = MIN(-sum_dP_muscle_ei/(Tpass*Texpn)*dt,-sum_dP_muscle)
            sum_dP_muscle = sum_dP_muscle+dP_muscle
         endif
         
      end select
      
      P_muscle = P_muscle + dP_muscle !current value for muscle pressure
  
      call enter_exit(sub_name,2)
  
    end subroutine set_driving_pressures
  
negative.code:          Tinsp,ttime
      real(dp) :: dP_muscle,P_muscle
      ! Local variables
      real(dp) :: mu
      character(len=60) :: sub_name
      
      ! --------------------------------------------------------------------------
  