[Proj] NAD27 lat lon to Lambert Conformal Conic

[Larry E. Ives]

I use Arcview to plot point locations which I receive originally in 
NAD27 lat lon. I use Arcview to reproject the value into State Plane 
FIPS 3301 which is North Dakota North NAD27. Arcview lists it as:

Projected Coordinate System: 
NAD_1927_StatePlane_North_Dakota_North_FIPS_3301
Projection: Lambert_Conformal_Conic
False_Easting: 2000000.00000000
False_Northing: 0.00000000
Central_Meridian: -100.50000000
Standard_Parallel_1: 47.43333333
Standard_Parallel_2: 48.73333333
Latitude_Of_Origin: 47.00000000
Linear Unit: Foot_US

I am trying to use cs2cs to duplicate this process. In short my 
latitudes are spot on but my longitudes are way off.

My input values are: -104.0486 +45.7912
In Arcview I get: 1094396.195 -420136.44
Using cs2cs I get: 5656062.96 -420136.44

Here is my script:

#!/bin/bash

cs2cs +proj=latlong +datum=NAD27 \
+to +datum=NAD27 +units=us-ft +proj=lcc \
+lat_0=+47.000000 \
+lat_1=+47.433333 \
+lat_2=+48.733333 \
+lon_0=-100.500000 \
+x_0=2000000 \
+y_0=0.000000 \
<<EOF
-104.0486 +45.7912
EOF

I have tried may things such as changing the input and script to DMS 
values with no change.
I have also tried a completely new set of input values with the same 
results.
Any help would be greatly appreciated.