pi0 = (ps / b  + a -ps) GF + ps
pinsp = f ps

Solve[pinsp +(pi0 - pinsp) Exp[- γ t] == ps Exp[-h/ hs] (GF/ b - GF +1) +a GF,t]

% /. C[1]->0

alist = {1.1696,1.0,0.8618,0.7562,0.62,0.5043,0.441,0.4,0.375,0.35,0.3295,0.3065,0.2835,0.261,0.248,0.2327}

blist = {0.5578,0.6514,0.7222,0.7825,0.8126,0.8434,0.8693,0.8910,0.9092,0.9222,0.9319,0.9403,0.9477,0.9544,0.9602,0.9653}

γtlist = Table[Log[(alist[[n]] GF /ps+1-f
+GF(1/blist[[n]]-1))/(alist[[n]]/ps GF-f+Exp[-h](1+GF(1/blist[[n]]-1)))], {n,1,Length[alist]}]

T12list = {5.0,8.0,12.5,18.5,27.0,38.3,54.3,77.0,109.0,146.0,187.0,239.0,305.0,390.0,498.0,635.0}

waittimes = Table[γtlist[[n]] T12list[[n]] Log[2]/60,{n,1,Length[γtlist]}]

values = {f-> 0.79, ps-> 1.013, h->2400/8500}

waittimes/. values /. GF->1

waittimes /. values /. GF->.5

waittimes /. values /. GF->.8

nwaittimes = waittimes /. values

Plot[nwaittimes,{GF,0,1},PlotRange->{0,24}, PlotLabels->T12list]

Show[%64,ImageSize->Large]

Plot[Evaluate[waittimes /. h->4500/8500 /. values], {GF,0,1},PlotRange->{0,24}, PlotLabels->T12list, ImageSize->Large]

z

Exp[#]& /@ γtlist /. values /. GF->.8

ps/blist[[16]]  +alist[[16]]/. values

ps Exp[-h] / blist[[16]] + alist[[16]]/. values

Plot[Evaluate[Table[-8500Log[(f - alist[[n]]GF /ps)/(GF/ blist[[n]]-GF +1)],{n,1,Length[alist]}] /. values ],{GF,0,1},PlotRange->{0,10000},PlotLabels-> Placed[T12list,Above]]

?Plot