@@ -603,8 +603,8 @@ def plotfitting3D(
603603 r11 = re * r_gen1
604604
605605 h = ana .ext_theis2D (
606- re ,
607- time ,
606+ time = time ,
607+ rad = re ,
608608 TG = np .exp (para [0 ]),
609609 sig2 = para [1 ],
610610 corr = para [2 ],
@@ -615,8 +615,8 @@ def plotfitting3D(
615615 ).reshape (- 1 )
616616 h1 = data [:, ri ]
617617 h2 = ana .ext_theis2D (
618- re ,
619- timarr ,
618+ time = timarr ,
619+ rad = re ,
620620 TG = np .exp (para [0 ]),
621621 sig2 = para [1 ],
622622 corr = para [2 ],
@@ -652,8 +652,8 @@ def plotfitting3D(
652652 for te in time :
653653 t11 = te * t_gen
654654 h = ana .ext_theis2D (
655- radarr ,
656- te ,
655+ time = te ,
656+ rad = radarr ,
657657 TG = np .exp (para [0 ]),
658658 sig2 = para [1 ],
659659 corr = para [2 ],
@@ -693,11 +693,11 @@ def plotfitting3Dtheis(data, para, rad, time, radnames, prate, plotname):
693693 r11 = re * r_gen1
694694
695695 h = ana .theis (
696- re , time , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
696+ time = time , rad = re , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
697697 ).reshape (- 1 )
698698 h1 = data [:, ri ]
699699 h2 = ana .theis (
700- re , timarr , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
700+ time = timarr , rad = re , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
701701 ).reshape (- 1 )
702702
703703 zord = 1000 * (len (rad ) - ri )
@@ -726,7 +726,7 @@ def plotfitting3Dtheis(data, para, rad, time, radnames, prate, plotname):
726726 for te in time :
727727 t11 = te * t_gen
728728 h = ana .theis (
729- radarr , te , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
729+ time = te , rad = radarr , T = np .exp (para [0 ]), S = np .exp (para [1 ]), Qw = prate
730730 ).reshape (- 1 )
731731 ax .plot (radarr , t11 , h , color = "k" , alpha = 0.1 , linestyle = "--" )
732732
0 commit comments