;$Id:$
;-----------------------------------------------------------------------------------
;
;-----------------------------------------------------------------------------------
pro compute_prior_from_map_2d, Lon, Lat, Bad_Pixel_Mask, Nx, Ny,  Missing, $
                            Lon_Min_Prior, dLon_Prior, Nlon_Prior, $
                            Lat_Min_Prior, dLat_Prior, Nlat_Prior, $
                            Prior, Month, Diurnal,  $
                            Prior_Values, Prior_Values_Ice, Prior_Values_Water

    Prior_Values = make_array(Nx,Ny,/FLOAT, Value=Missing)
    Prior_Values_Ice = make_array(Nx,Ny,/FLOAT, Value=Missing)
    Prior_Values_Water = make_array(Nx,Ny,/FLOAT, Value=Missing)

    idiurnal = Diurnal - 1    ;1(0)= daily averaged
    imonth = Month - 1

    for i = 0, Nx-1 do begin 
       for j = 0, Ny-1 do begin

        ilon = min([Nlon_Prior-1,max([0,fix((Lon[i,j] - Lon_Min_Prior) / (dLon_Prior))])])
        ilat = min([Nlat_Prior-1,max([0,fix((Lat[i,j] - Lat_Min_Prior) / (dLat_Prior))])])

        if (ilon ge 0 and ilat ge 0 and imonth ge 0) then begin
            Prior_Values[i,j] = Prior.Table[ilon,ilat,imonth,idiurnal]
            Prior_Values_Ice[i,j] = Prior.Ice_Table[ilon,ilat,imonth,idiurnal]
            Prior_Values_Water[i,j] = Prior.Water_Table[ilon,ilat,imonth,idiurnal]
        endif

       endfor
    endfor

end
;-----------------------------------------------------------------------------------
; training data (1d)
;-----------------------------------------------------------------------------------
pro compute_prior_from_map_1d, Lon, Lat, Bad_Pixel_Mask, Nx,  Missing, $
                            Lon_Min_Prior, dLon_Prior, Nlon_Prior, $
                            Lat_Min_Prior, dLat_Prior, Nlat_Prior, $
                            Prior, Month, Diurnal, $
                            Prior_Values,Prior_Values_Ice, Prior_Values_Water

    Prior_Values = make_array(Nx,/FLOAT, Value=Missing)
    Prior_Values_Ice = make_array(Nx,/FLOAT, Value=Missing)
    Prior_Values_Water = make_array(Nx,/FLOAT, Value=Missing)

    idiurnal = Diurnal - 1    ;1(0)= daily averaged

    for i = 0, Nx-1 do begin 

        imonth = Month[i] - 1
        ilon = min([Nlon_Prior-1,max([0,fix((Lon[i] - Lon_Min_Prior) / (dLon_Prior))])])
        ilat = min([Nlat_Prior-1,max([0,fix((Lat[i] - Lat_Min_Prior) / (dLat_Prior))])])

        if (ilon ge 0 and ilat ge 0 and imonth ge 0) then begin
           Prior_Values[i] = Prior.Table[ilon,ilat,imonth,idiurnal]
           Prior_Values_Ice[i] = Prior.Ice_Table[ilon,ilat,imonth,idiurnal]
           Prior_Values_Water[i] = Prior.Water_Table[ilon,ilat,imonth,idiurnal]
        endif

    endfor

end