gdal2tiles.py and ProjNatOriginLongGeoKey

Hello All,

I'm trying to convert GeoTIFF to TMS tiles using gdal2tiles.py script.

listgeo:
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
600449 4934034 0
ModelPixelScaleTag (1,3):
0.5 0.5 0
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,257): "IMAGINE GeoTIFF Support
\nCopyright 1991 - 2006 by Leica Geosystems Geospatial Imaging, LLC.
All Rights Reserved\n@(#)$RCSfile: egtf.c $ IMAGINE 9.1 $Revision:
22.0 $ $Date: 2006/05/24 17:16:00 EST $\nProjection Name = UTM\nUnits
= meters\nGeoTIFF Units = meters"
ProjectedCSTypeGeoKey (Short,1): PCS_WGS84_UTM_zone_36N
PCSCitationGeoKey (Ascii,244): "IMAGINE GeoTIFF Support
\nCopyright 1991 - 2006 by Leica Geosystems Geospatial Imaging, LLC.
All Rights Reserved\n@(#)$RCSfile: egtf.c $ IMAGINE 9.1 $Revision:
22.0 $ $Date: 2006/05/24 17:16:00 EST $\nUTM Zone 36N\nEllipsoid = WGS
84\nDatum = WGS 84"
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.

PCS = 32636 (WGS 84 / UTM zone 36N)
Projection = 16036 (UTM zone 36N)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
ProjNatOriginLongGeoKey: 33.000000 ( 33d 0' 0.00"E)
ProjScaleAtNatOriginGeoKey: 0.999600
ProjFalseEastingGeoKey: 500000.000000 m
ProjFalseNorthingGeoKey: 0.000000 m
GCS: 4326/WGS 84
Datum: 6326/World Geodetic System 1984
Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)

Corner Coordinates:
Upper Left ( 600449.000, 4934034.000) ( 34d15'52.85"E, 44d33'
9.50"N)
Lower Left ( 600449.000, 4931945.500) ( 34d15'51.39"E, 44d32'
1.83"N)
Upper Right ( 602409.500, 4934034.000) ( 34d17'21.69"E, 44d33'
8.51"N)
Lower Right ( 602409.500, 4931945.500) ( 34d17'20.20"E, 44d32'
0.84"N)
Center ( 601429.250, 4932989.750) ( 34d16'36.53"E,
44d32'35.17"N)

I use bash script to convert:
#!/bin/sh

input=1.tif
gdal_translate -of VRT -a_srs EPSG:900913 -gcp 0 0 600449 4934034 -gcp
3921 4177 602409.5 4931945.5 -gcp 3921 0 602409.5 4934034 $input 2.vrt
rm 3.vrt
rm -rf sat
gdalwarp -of VRT -s_srs EPSG:900913 -t_srs EPSG:4326 2.vrt 3.vrt
./gdal2tiles.py -v -z 4-18 --s_srs EPSG:4326 3.vrt sat

The problem is tiles are generated but it looks like they lie far to
the west from place they should be.

I guess ProjNatOriginLongGeoKey tells about this shift but I can't
find out what it means and how use it?

Thanks in advance,
Eduard