import numpy def weightfit(w,year,doy,index,ht,beam): if len(w) <= 0 : return w if year == 2020 and doy == 265 : if beam == 0 : if index >= 120 and index <= 183 and ht >=30: w[0,100:120] =0 if beam == 1 : if index >= 132 and index <= 180 and ht ==11: w[0,90:120] = 0 if index >= 120 and index <= 183 and ht >=30: w[0,100:120] = 0 if index >= 176 and index <= 176 and ht ==20: w[0,60:100] = 0.0 if index >= 177 and index <= 177 and ht >=24 and ht <=29: w[0,70:80] = 0.0 if index >= 177 and index <= 177 and ht >=30: w[0,70:80] = 0.0 if index >= 178 and index <= 178 and ht>=18: w[0,100:120] = 0.0 print('entra a la funcion w') # If index eq 96 then stop ;stop 96 return w def Vrfit(p0,year,doy,index,h,beam): if year == 2020 and doy == 265 : if beam == 0 : if index == 177 and h==30 : p0[3]=20 if index == 184 and h>=27 : p0[3]=30 if index == 186 and h>=11 : p0[3]=30 if index == 190 and h==27 : p0[3]=27 if index == 191 and h==26 : p0[3]=28 if index == 191 and h>=35 : p0[3]=24 if index == 131 and h==29 : p0[3]=24 if index == 138 and h==11 : p0[3]=30 if index == 138 and h==14 : p0[3]=30 if index == 139 and h>=34 and h<=36 : p0[3]=20 if index == 142 and h>=34 : p0[3]=35 if index == 155 and h>=27 and h<=40: p0[3]=36 if index == 167 and h>=13 and h<=13: p0[3]=33 if index == 168 and h>=13 and h<=13: p0[3]=37 if index == 172 and h>=13 and h<=13: p0[3]=37 if index == 173 and h>=28 and h<=29: p0[3]=37 print('entra a la funcion vr') if beam == 1 : if index == 177 and h>=31 and h<=32: p0[3]=30 if index == 179 and h==27 : p0[3]=38 if index == 180 and h==23 : p0[3]=33 if index == 180 and h==26 : p0[3]=30 if index == 180 and h==29 : p0[3]=35 if index == 184 and h>=23 : p0[3]=30 if index == 185 and h==11 : p0[3]=20 if index == 185 and h>=29 : p0[3]=30 if index == 186 and h>=11 : p0[3]=30 if index == 190 and h==11 : p0[3]=28 if index == 191 and h>=33 : p0[3]=20 if index == 131 and h==11 : p0[3]=17 if index == 131 and h>=29 and h<=38 : p0[3]=18 if index == 138 and h>=31 : p0[3]=15 if index == 139 and h>=11 and h<=16 : p0[3]=25 if index == 140 and h==33 : p0[3]=26 if index == 142 and h==34 : p0[3]=20 if index == 149 and h>=28 and h<=32: p0[3]=28 if index == 163 and h>=11 and h<=11: p0[3]=30 if index == 167 and h>=11 and h<=16: p0[3]=33 if index == 167 and h>=34 and h<=40: p0[3]=36 if index == 169 and h>=26 and h<=49: p0[3]=33 if index == 172 and h>=28 and h<=28: p0[3]=37 if index == 173 and h>=20 and h<=38: p0[3]=37 if index == 174 and h>=25 and h<=25: p0[3]=34 if index == 175 and h>=27 and h<=35: p0[3]=30 #if h == 11: print('entra ',p0[3]) return p0