weightfit.py
73 lines
| 3.0 KiB
| text/x-python
|
PythonLexer
Alexander Valdez
|
r1673 | 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 | ||||