	function HOSDustModel(iD,P,C,S)

	    integer	iD
	    integer	P		! Photometer: 1 <= P <=3
	    integer	C		! Color: 1=U, 2=B, 3=V
	    integer	S		! Sector: -15 <= S <=32

	! Table 2 (Perihelion model; 0.3 AU)

	integer	ZP(16,2,3)	/
     &		25800,22350,17440,13280,10040, 7770, 6110, 4920, 3710, 2720, 2150, 1795, 1501, 1316, 1282, 1340,
     &		28500,24700,19280,14600,11000, 8470, 6680, 5370, 4050, 2955, 2315, 1930, 1604, 1392, 1344, 1404,
     &		31800,27550,21500,16280,12260, 9430, 7370, 5890, 4430, 3230, 2520, 2095, 1720, 1485, 1425, 1484,

     &		 5430, 5220, 4870, 4400, 3920, 3440, 3010, 2645, 2210, 1804, 1528, 1342, 1174, 1076, 1052, 1076,
     &		 6020, 5780, 5360, 4840, 4270, 3760, 3280, 2885, 2405, 1952, 1646, 1436, 1262, 1136, 1102, 1130,
     &		 6740, 6490, 6020, 5430, 4780, 4180, 3650, 3200, 2655, 2155, 1804, 1572, 1358, 1230, 1184, 1198/
  
	integer	PP(     3)	/ 818,  880,  938/

	! Table 1 (Aphelion model: 1.0 AU)

	integer	ZA(16,3,2)	/
     &		 1506, 1306, 1022,  775,  596,  459,  367,  296,  228,  173,  138,  116,   99,   90,   91,   97,
     &		 1652, 1424, 1118,  847,  644,  495,  393,  316,  241,  181,  145,  120,  101,   92,   92,   97,
     &		 1882, 1630, 1272,  963,  726,  555,  439,  353,  269,  200,  160,  132,  110,   95,   05,  101,

     &		  317,  304,  284,  257,  227,  200,  180,  157,  135,  113,   98,   90,   81,   79,   80,   82,
     &		  343,  330,  307,  275,  243,  213,  190,  167,  141,  118,  101,   91,   83,   79,   81,   82,
     &		  387,  372,  344,  310,  273,  239,  211,  187,  158,  132,  111,   99,   89,   84,   85,   86/

	integer	PA(     3)	/  52,   53,   57/

	J = S
	if (S .gt. 16) J = S-32		! Modified sector #
	if (J .le. 0) J = 1-J

	if (iD .eq. 1) then

	    if (P .eq. 3) then
		HOSDustModel = PP(C)
	    else
		HOSDustModel = ZP(J,C,P)
	    end if

	else

	    if (P .eq. 3) then
		HOSDustModel = PA(C)
	    else
		HOSDustModel = ZA(J,C,P)
	    end if

	end if

	return
	end