comparison MAP/MAP1_14.m @ 9:ecad0ea62b43

May27 mainly better parameters

 IHCs1= IHC_cilia_RPParams.s1;
 IHCGmax= IHC_cilia_RPParams.Gmax;
 IHCGa= IHC_cilia_RPParams.Ga; % (leakage)
 
 IHCGu0 = IHCGa+IHCGmax./(1+exp(IHCu0/IHCs0).*(1+exp(IHCu1/IHCs1)));
-
-
+IHCrestingCiliaCond=IHCGu0;
 
 % Receptor potential
 IHC_Cab= IHC_cilia_RPParams.Cab;
 IHC_Gk= IHC_cilia_RPParams.Gk;
 IHC_Et= IHC_cilia_RPParams.Et;
 IHC_Ek= IHC_cilia_RPParams.Ek;
 IHC_Ekp= IHC_Ek+IHC_Et*IHC_cilia_RPParams.Rpc;
 
-IHCrestingV= -0.06;
-
-IHCrestingV= (IHC_Gk*IHC_Ek+IHCGu0*IHC_Et)/(IHCGu0+IHC_Gk);
+IHCrestingV= (IHC_Gk*IHC_Ekp+IHCGu0*IHC_Et)/(IHCGu0+IHC_Gk);
 
 IHC_Vnow= IHCrestingV*ones(nBFs,1); % initial voltage
 IHC_RP= zeros(nBFs,segmentLength);
 
 % complete record of IHC receptor potential (V)

     % and save it
     IHC_cilia_output(:,segmentStartPTR:segmentStartPTR+segmentLength-1)=...
         IHCciliaDisplacement;
 
     % compute apical conductance
-    G=1./(1+exp(-(IHCciliaDisplacement-IHCu0)/IHCs0).*...
+    G=IHCGmax./(1+exp(-(IHCciliaDisplacement-IHCu0)/IHCs0).*...
         (1+exp(-(IHCciliaDisplacement-IHCu1)/IHCs1)));
-    Gu=IHCGmax*G + IHCGa;
+    Gu=G + IHCGa;
 
     % Compute receptor potential
     for idx=1:segmentLength
         IHC_Vnow=IHC_Vnow+ (-Gu(:, idx).*(IHC_Vnow-IHC_Et)-...
             IHC_Gk*(IHC_Vnow-IHC_Ekp))*  dt/IHC_Cab;