PRO loadlicmixr_lc,yrls,jday,ndays,mixr_lic,temp_lic ; yrls---yyyy string, year ; jday---start jday, integ ; ndays-- number of days, integ filename='' filename=!statsfile + yrls+'.stats' load,filename,stat print,48L*ndays mixr_lic=fltarr(48L*ndays)-999.0 temp_lic=mixr_lic n1=48L*(jday-1) ;initial number of the array to read n2=n1+len(mixr_lic)-1 mixr_lic(*)=stat(18,n1:n2) temp_lic(*)=stat(14,n1:n2) r=where(temp_lic NE -999,count) IF count GE 1 THEN temp_lic(r)=temp_lic(r)+273.15 return END