wradlib_warp_translate.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Warp Radolan data to GTIFF   \n",
    "https://docs.wradlib.org/en/stable/georef.html#raster-functions  \n",
    "\n",
    "Use this cell to create geotiff in gk3 or ll projection\n",
    "change line 77 and 83 with correct coordinates and crs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#For Files located in single folder\n",
    "#Import Modules\n",
    "import wradlib as wrl\n",
    "import numpy as np\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "import os \n",
    "from osgeo import osr\n",
    "\n",
    "\n",
    "#Working dir paths must be clean otherwise you have to make list comprehension that only bin files in files list\n",
    "wdir = \"/media/elisabeth/0A581C5C581C4939/master/radolan_data/out/bin/2009-2019/\"\n",
    "out = \"out\"\n",
    "#Change to Workingdir\n",
    "os.chdir(wdir)\n",
    "\n",
    "if not os.path.exists(out):\n",
    "    os.mkdir(out)\n",
    "\n",
    "\n",
    "#create the necessary Spatial Reference Objects for the RADOLAN-projection \n",
    "proj_stereo = wrl.georef.create_osr(\"dwd-radolan\") #create osr object\n",
    "proj_wgs = osr.SpatialReference()\n",
    "proj_wgs.ImportFromEPSG(4326)\n",
    "radolan_grid_xy = wrl.georef.get_radolan_grid(900,900)\n",
    "\n",
    "# create Gauss Krueger zone 3 projection osr object\n",
    "proj_gk3 = osr.SpatialReference()\n",
    "proj_gk3.ImportFromEPSG(31467) #dhdn3\n",
    "\n",
    "#Bug handling\n",
    "epsg_31467 = wrl.georef.epsg_to_osr(31467) \n",
    "epsg_4326 = wrl.georef.epsg_to_osr(4326) \n",
    "\n",
    "\n",
    "# transform radolan polar stereographic projection to wgs84 and then to gk3\n",
    "radolan_grid_ll = wrl.georef.reproject(radolan_grid_xy,\n",
    "                                       projection_source=proj_stereo,\n",
    "                                       projection_target=proj_wgs)\n",
    "radolan_grid_gk = wrl.georef.reproject(radolan_grid_ll,\n",
    "                                       projection_source=proj_wgs,\n",
    "                                       projection_target=proj_gk3)\n",
    "\n",
    "#print(radolan_grid_gk)\n",
    "#Get list of files inside the folder\n",
    "files = os.listdir()\n",
    "print(len(files),\" files inside the folder\")\n",
    "\n",
    "#1 because of out is first file later do a list comprehension maybe more beautiful ^^\n",
    "for file in files[1:]:\n",
    "    proj_stereo = wrl.georef.create_osr(\"dwd-radolan\") #create osr object\n",
    "    proj_wgs = osr.SpatialReference()\n",
    "    proj_wgs.ImportFromEPSG(4326)\n",
    "    radolan_grid_xy = wrl.georef.get_radolan_grid(900,900)\n",
    "    #proj_gk3 = osr.SpatialReference()\n",
    "    #proj_gk3.ImportFromEPSG(31467) #dhdn3\n",
    "    radolan_grid_ll = wrl.georef.reproject(radolan_grid_xy,\n",
    "                                       projection_source=proj_stereo,\n",
    "                                       projection_target=proj_wgs)\n",
    "    #radolan_grid_gk = wrl.georef.reproject(radolan_grid_ll,\n",
    "                                      # projection_source=proj_wgs,\n",
    "                                       #projection_target=proj_gk3)\n",
    "    #print(\"Load: \",file ,\" data, georeferencing, reprojection and translate to Gtif...\")\n",
    "    data_raw = None\n",
    "    meta = None\n",
    "    data = None\n",
    "    xy = None\n",
    "    ds_gk3 = None\n",
    "    data_raw, meta  = wrl.io.read_radolan_composite(file)\n",
    "    data_raw = data_raw.round(decimals=3)\n",
    "    data_raw = data_raw.astype('float32')\n",
    "    #print(radolan_grid_ll[1])\n",
    "    #print(data_raw.dtype)\n",
    "    #print(np.median(data_raw))\n",
    "    #flip the data and the coords up-down (from lower left to upper left corner)\n",
    "    #data, xy = wrl.georef.set_raster_origin(data_raw, radolan_grid_xy, 'upper')\n",
    "    data, xy = wrl.georef.set_raster_origin(data_raw, radolan_grid_ll, 'upper')\n",
    "    #print(xy[-1])\n",
    "    #print(xy)\n",
    "    #print(radolan_grid_gk)\n",
    "    #Create raster dataset: Geotiff with radolan projection\n",
    "    #ds = wrl.georef.create_raster_dataset(data, xy, projection=proj_stereo, nodata=-9999)\n",
    "    ds_gk3 = wrl.georef.create_raster_dataset(data, xy, proj_wgs, -9999) # https://docs.wradlib.org/en/stable/generated/wradlib.georef.raster.create_raster_dataset.html\n",
    "    #Reproject dataset to GK3\n",
    "    #ds2 = wrl.georef.reproject_raster_dataset(ds,projection_target=proj_gk3, size=(900,900)) # align=(3052995.645672334, 6067805.270329222) spacing=1000 size=(900,900)\n",
    "    #Bug appears it will fixed in future https://groups.google.com/forum/#!topic/wradlib-users/-dvRhDCjgV0\n",
    "    #This line is needed until fixed\n",
    "    #ds.SetProjection(epsg_31467.ExportToWkt())\n",
    "    #ds.SetProjection(epsg_4326.ExportToWkt())\n",
    "    #write raster dataset\n",
    "    filename = file[15:25] + \".tif\"\n",
    "    wrl.io.write_raster_dataset(wdir + out + \"/\" + filename, ds_gk3, \"GTiff\")\n",
    "    #print(\"Save\", filename, \"to\", wdir + out)\n",
    "\n",
    "print('finish')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "ename": "SyntaxError",
     "evalue": "invalid syntax (<ipython-input-2-5f3b253fd6a5>, line 42)",
     "output_type": "error",
     "traceback": [
      "\u001b[0;36m  File \u001b[0;32m\"<ipython-input-2-5f3b253fd6a5>\"\u001b[0;36m, line \u001b[0;32m42\u001b[0m\n\u001b[0;31m    filename = file1[15:25] + \"_wgs_tr + out2 + \"/\" + filename, rad_wgs, \"GTiff\")\u001b[0m\n\u001b[0m                                                                              ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
     ]
    }
   ],
   "source": [
    "from osgeo import osr\n",
    "import wradlib as wrl\n",
    "\n",
    "#Data Paths\n",
    "wdir = \"/media/elisabeth/0A581C5C581C4939/master/radolan_data/out/bin/2009-2019/\"\n",
    "out2 = \"test/\"\n",
    "file1 = \"raa01-sf_10000-0901012350-dwd---bin\"\n",
    "file2 = \"raa01-sf_10000-1907312350-dwd---bin\"\n",
    "\n",
    "proj_stereo = wrl.georef.create_osr(\"dwd-radolan\") #(osr.SpatialReference) – GDAL/OSR object defining projection\n",
    "proj_wgs = osr.SpatialReference() #Initialize a SpatialReference object\n",
    "proj_wgs.ImportFromEPSG(4326) # Spatial Referenz for EPSG 3426\n",
    "#print(proj_wgs)\n",
    "\n",
    "radolan_grid_xy = wrl.georef.get_radolan_grid(900,900) \n",
    "radolan_grid_ll = wrl.georef.projection.reproject(radolan_grid_xy, projection_source=proj_stereo, projection_target=proj_wgs)\n",
    "#print(\"{0}, ({1:.4f}, {2:.4f})\".format(radolan_grid_ll.shape, *radolan_grid_ll[0,0,:]))\n",
    "\n",
    "# create Gauss Krueger zone 3 projection osr object\n",
    "proj_gk3 = osr.SpatialReference()  #Initialize a SpatialReference object\n",
    "proj_gk3.ImportFromEPSG(31467) #Spatial Referenz for EPSG 31467\n",
    "#print(proj_gk3)\n",
    "\n",
    "# transform radolan polar stereographic projection to wgs84 and then to gk3\n",
    "radolan_grid_gk = wrl.georef.reproject(radolan_grid_ll,projection_source=proj_wgs,projection_target=proj_gk3)\n",
    "#print(\"{0}, ({1:.4f}, {2:.4f})\".format(radolan_grid_gk.shape, *radolan_grid_gk[0,0,:]))\n",
    "\n",
    "data_raw, attributes = wrl.io.radolan.read_radolan_composite(wdir + file1, missing=-9999, loaddata=True) #data : numpy.array of shape (number of rows, number of columns) | xarray.Dataset >// attrs : dictionary of metadata information from the file heade\n",
    "data, xy = wrl.georef.set_raster_origin(data_raw, radolan_grid_xy, 'upper')\n",
    "\n",
    "rad_stereo = wrl.georef.create_raster_dataset(data, xy, projection=proj_stereo, nodata=-9999) #dataset (gdal.Dataset) – In-Memory raster dataset\n",
    "\n",
    "cols = rad_stereo.RasterXSize\n",
    "rows = rad_stereo.RasterYSize\n",
    "\n",
    "filename = file1[15:25] + \"_stereo\" + \".tif\"\n",
    "\n",
    "wrl.io.write_raster_dataset(wdir + out2 + \"/\" + filename, rad_stereo, \"GTiff\")\n",
    "\n",
    "rad_wgs =  wrl.georef.raster.reproject_raster_dataset(rad_stereo, size = (cols,rows), projection_target = proj_wgs) \n",
    "\n",
    "filename = file1[15:25] + \"_wgs_tr + out2 + \"/\" + filename, rad_wgs, \"GTiff\")\n",
    "\n",
    "rad_gk3 =  wrl.georef.raster.reproject_raster_dataset(rad_wgs, size = (cols,rows), projection_target = proj_gk3) \n",
    "\n",
    "filename = file1[15:25] + \"_gk3_1\" + \".tif\"\n",
    "\n",
    "wrl.io.write_raster_dataset(wdir + out2 + \"/\" + filename, rad_gk3, \"GTiff\")\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3862\n",
      "1/3862\n",
      "2/3862\n",
      "3/3862\n",
      "4/3862\n",
      "5/3862\n",
      "6/3862\n",
      "7/3862\n",
      "8/3862\n",
      "9/3862\n",
      "10/3862\n",
      "11/3862\n",
      "12/3862\n",
      "13/3862\n",
      "14/3862\n",
      "15/3862\n",
      "16/3862\n",
      "17/3862\n",
      "18/3862\n",
      "19/3862\n",
      "20/3862\n",
      "21/3862\n",
      "22/3862\n",
      "23/3862\n",
      "24/3862\n",
      "25/3862\n",
      "26/3862\n",
      "27/3862\n",
      "28/3862\n",
      "29/3862\n",
      "30/3862\n",
      "31/3862\n",
      "32/3862\n",
      "33/3862\n",
      "34/3862\n",
      "35/3862\n",
      "36/3862\n",
      "37/3862\n",
      "38/3862\n",
      "39/3862\n",
      "40/3862\n",
      "41/3862\n",
      "42/3862\n",
      "43/3862\n",
      "44/3862\n",
      "45/3862\n",
      "46/3862\n",
      "47/3862\n",
      "48/3862\n",
      "49/3862\n",
      "50/3862\n",
      "51/3862\n",
      "52/3862\n",
      "53/3862\n",
      "54/3862\n",
      "55/3862\n",
      "56/3862\n",
      "57/3862\n",
      "58/3862\n",
      "59/3862\n",
      "60/3862\n",
      "61/3862\n",
      "62/3862\n",
      "63/3862\n",
      "64/3862\n",
      "65/3862\n",
      "66/3862\n",
      "67/3862\n",
      "68/3862\n",
      "69/3862\n",
      "70/3862\n",
      "71/3862\n",
      "72/3862\n",
      "73/3862\n",
      "74/3862\n",
      "75/3862\n",
      "76/3862\n",
      "77/3862\n",
      "78/3862\n",
      "79/3862\n",
      "80/3862\n",
      "81/3862\n",
      "82/3862\n",
      "83/3862\n",
      "84/3862\n",
      "85/3862\n",
      "86/3862\n",
      "87/3862\n",
      "88/3862\n",
      "89/3862\n",
      "90/3862\n",
      "91/3862\n",
      "92/3862\n",
      "93/3862\n",
      "94/3862\n",
      "95/3862\n",
      "96/3862\n",
      "97/3862\n",
      "98/3862\n",
      "99/3862\n",
      "100/3862\n",
      "101/3862\n",
      "102/3862\n",
      "103/3862\n",
      "104/3862\n",
      "105/3862\n",
      "106/3862\n",
      "107/3862\n",
      "108/3862\n",
      "109/3862\n",
      "110/3862\n",
      "111/3862\n",
      "112/3862\n",
      "113/3862\n",
      "114/3862\n",
      "115/3862\n",
      "116/3862\n",
      "117/3862\n",
      "118/3862\n",
      "119/3862\n",
      "120/3862\n",
      "121/3862\n",
      "122/3862\n",
      "123/3862\n",
      "124/3862\n",
      "125/3862\n",
      "126/3862\n",
      "127/3862\n",
      "128/3862\n",
      "129/3862\n",
      "130/3862\n",
      "131/3862\n",
      "132/3862\n",
      "133/3862\n",
      "134/3862\n",
      "135/3862\n",
      "136/3862\n",
      "137/3862\n",
      "138/3862\n",
      "139/3862\n",
      "140/3862\n",
      "141/3862\n",
      "142/3862\n",
      "143/3862\n",
      "144/3862\n",
      "145/3862\n",
      "146/3862\n",
      "147/3862\n",
      "148/3862\n",
      "149/3862\n",
      "150/3862\n",
      "151/3862\n",
      "152/3862\n",
      "153/3862\n",
      "154/3862\n",
      "155/3862\n",
      "156/3862\n",
      "157/3862\n",
      "158/3862\n",
      "159/3862\n",
      "160/3862\n",
      "161/3862\n",
      "162/3862\n",
      "163/3862\n",
      "164/3862\n",
      "165/3862\n",
      "166/3862\n",
      "167/3862\n",
      "168/3862\n",
      "169/3862\n",
      "170/3862\n",
      "171/3862\n",
      "172/3862\n",
      "173/3862\n",
      "174/3862\n",
      "175/3862\n",
      "176/3862\n",
      "177/3862\n",
      "178/3862\n",
      "179/3862\n",
      "180/3862\n",
      "181/3862\n",
      "182/3862\n",
      "183/3862\n",
      "184/3862\n",
      "185/3862\n",
      "186/3862\n",
      "187/3862\n",
      "188/3862\n",
      "189/3862\n",
      "190/3862\n",
      "191/3862\n",
      "192/3862\n",
      "193/3862\n",
      "194/3862\n",
      "195/3862\n",
      "196/3862\n",
      "197/3862\n",
      "198/3862\n",
      "199/3862\n",
      "200/3862\n",
      "201/3862\n",
      "202/3862\n",
      "203/3862\n",
      "204/3862\n",
      "205/3862\n",
      "206/3862\n",
      "207/3862\n",
      "208/3862\n",
      "209/3862\n",
      "210/3862\n",
      "211/3862\n",
      "212/3862\n",
      "213/3862\n",
      "214/3862\n",
      "215/3862\n",
      "216/3862\n",
      "217/3862\n",
      "218/3862\n",
      "219/3862\n",
      "220/3862\n",
      "221/3862\n",
      "222/3862\n",
      "223/3862\n",
      "224/3862\n",
      "225/3862\n",
      "226/3862\n",
      "227/3862\n",
      "228/3862\n",
      "229/3862\n",
      "230/3862\n",
      "231/3862\n",
      "232/3862\n",
      "233/3862\n",
      "234/3862\n",
      "235/3862\n",
      "236/3862\n",
      "237/3862\n",
      "238/3862\n",
      "239/3862\n",
      "240/3862\n",
      "241/3862\n",
      "242/3862\n",
      "243/3862\n",
      "244/3862\n",
      "245/3862\n",
      "246/3862\n",
      "247/3862\n",
      "248/3862\n",
      "249/3862\n",
      "250/3862\n",
      "251/3862\n",
      "252/3862\n",
      "253/3862\n",
      "254/3862\n",
      "255/3862\n",
      "256/3862\n",
      "257/3862\n",
      "258/3862\n",
      "259/3862\n",
      "260/3862\n",
      "261/3862\n",
      "262/3862\n",
      "263/3862\n",
      "264/3862\n",
      "265/3862\n",
      "266/3862\n",
      "267/3862\n",
      "268/3862\n",
      "269/3862\n",
      "270/3862\n",
      "271/3862\n",
      "272/3862\n",
      "273/3862\n",
      "274/3862\n",
      "275/3862\n",
      "276/3862\n",
      "277/3862\n",
      "278/3862\n",
      "279/3862\n",
      "280/3862\n",
      "281/3862\n",
      "282/3862\n",
      "283/3862\n",
      "284/3862\n",
      "285/3862\n",
      "286/3862\n",
      "287/3862\n",
      "288/3862\n",
      "289/3862\n",
      "290/3862\n",
      "291/3862\n",
      "292/3862\n",
      "293/3862\n",
      "294/3862\n",
      "295/3862\n",
      "296/3862\n",
      "297/3862\n",
      "298/3862\n",
      "299/3862\n",
      "300/3862\n",
      "301/3862\n",
      "302/3862\n",
      "303/3862\n",
      "304/3862\n",
      "305/3862\n",
      "306/3862\n",
      "307/3862\n",
      "308/3862\n",
      "309/3862\n",
      "310/3862\n",
      "311/3862\n",
      "312/3862\n",
      "313/3862\n",
      "314/3862\n",
      "315/3862\n",
      "316/3862\n",
      "317/3862\n",
      "318/3862\n",
      "319/3862\n",
      "320/3862\n",
      "321/3862\n",
      "322/3862\n",
      "323/3862\n",
      "324/3862\n",
      "325/3862\n",
      "326/3862\n",
      "327/3862\n",
      "328/3862\n",
      "329/3862\n",
      "330/3862\n",
      "331/3862\n",
      "332/3862\n",
      "333/3862\n",
      "334/3862\n",
      "335/3862\n",
      "336/3862\n",
      "337/3862\n",
      "338/3862\n",
      "339/3862\n",
      "340/3862\n",
      "341/3862\n",
      "342/3862\n",
      "343/3862\n",
      "344/3862\n",
      "345/3862\n",
      "346/3862\n",
      "347/3862\n",
      "348/3862\n",
      "349/3862\n",
      "350/3862\n",
      "351/3862\n",
      "352/3862\n",
      "353/3862\n",
      "354/3862\n",
      "355/3862\n",
      "356/3862\n",
      "357/3862\n",
      "358/3862\n",
      "359/3862\n",
      "360/3862\n",
      "361/3862\n",
      "362/3862\n",
      "363/3862\n",
      "364/3862\n",
      "365/3862\n",
      "366/3862\n",
      "367/3862\n",
      "368/3862\n",
      "369/3862\n",
      "370/3862\n",
      "371/3862\n",
      "372/3862\n",
      "373/3862\n",
      "374/3862\n",
      "375/3862\n",
      "376/3862\n",
      "377/3862\n",
      "378/3862\n",
      "379/3862\n",
      "380/3862\n",
      "381/3862\n",
      "382/3862\n",
      "383/3862\n",
      "384/3862\n",
      "385/3862\n",
      "386/3862\n",
      "387/3862\n",
      "388/3862\n",
      "389/3862\n",
      "390/3862\n",
      "391/3862\n",
      "392/3862\n",
      "393/3862\n",
      "394/3862\n",
      "395/3862\n",
      "396/3862\n",
      "397/3862\n",
      "398/3862\n",
      "399/3862\n",
      "400/3862\n",
      "401/3862\n",
      "402/3862\n",
      "403/3862\n",
      "404/3862\n",
      "405/3862\n",
      "406/3862\n",
      "407/3862\n",
      "408/3862\n",
      "409/3862\n",
      "410/3862\n",
      "411/3862\n",
      "412/3862\n",
      "413/3862\n",
      "414/3862\n",
      "415/3862\n",
      "416/3862\n",
      "417/3862\n",
      "418/3862\n",
      "419/3862\n",
      "420/3862\n",
      "421/3862\n",
      "422/3862\n",
      "423/3862\n",
      "424/3862\n",
      "425/3862\n",
      "426/3862\n",
      "427/3862\n",
      "428/3862\n",
      "429/3862\n",
      "430/3862\n",
      "431/3862\n",
      "432/3862\n",
      "433/3862\n",
      "434/3862\n",
      "435/3862\n",
      "436/3862\n",
      "437/3862\n",
      "438/3862\n",
      "439/3862\n",
      "440/3862\n",
      "441/3862\n",
      "442/3862\n",
      "443/3862\n",
      "444/3862\n",
      "445/3862\n",
      "446/3862\n",
      "447/3862\n",
      "448/3862\n",
      "449/3862\n",
      "450/3862\n",
      "451/3862\n",
      "452/3862\n",
      "453/3862\n",
      "454/3862\n",
      "455/3862\n",
      "456/3862\n",
      "457/3862\n",
      "458/3862\n",
      "459/3862\n",
      "460/3862\n",
      "461/3862\n",
      "462/3862\n",
      "463/3862\n",
      "464/3862\n",
      "465/3862\n",
      "466/3862\n",
      "467/3862\n",
      "468/3862\n",
      "469/3862\n",
      "470/3862\n",
      "471/3862\n",
      "472/3862\n",
      "473/3862\n",
      "474/3862\n",
      "475/3862\n",
      "476/3862\n",
      "477/3862\n",
      "478/3862\n",
      "479/3862\n",
      "480/3862\n",
      "481/3862\n",
      "482/3862\n",
      "483/3862\n",
      "484/3862\n",
      "485/3862\n",
      "486/3862\n",
      "487/3862\n",
      "488/3862\n",
      "489/3862\n",
      "490/3862\n",
      "491/3862\n",
      "492/3862\n",
      "493/3862\n",
      "494/3862\n",
      "495/3862\n",
      "496/3862\n",
      "497/3862\n",
      "498/3862\n",
      "499/3862\n",
      "500/3862\n",
      "501/3862\n",
      "502/3862\n",
      "503/3862\n",
      "504/3862\n",
      "505/3862\n",
      "506/3862\n",
      "507/3862\n",
      "508/3862\n",
      "509/3862\n",
      "510/3862\n",
      "511/3862\n",
      "512/3862\n",
      "513/3862\n",
      "514/3862\n",
      "515/3862\n",
      "516/3862\n",
      "517/3862\n",
      "518/3862\n",
      "519/3862\n",
      "520/3862\n",
      "521/3862\n",
      "522/3862\n",
      "523/3862\n",
      "524/3862\n",
      "525/3862\n",
      "526/3862\n",
      "527/3862\n",
      "528/3862\n",
      "529/3862\n",
      "530/3862\n",
      "531/3862\n",
      "532/3862\n",
      "533/3862\n",
      "534/3862\n",
      "535/3862\n",
      "536/3862\n",
      "537/3862\n",
      "538/3862\n",
      "539/3862\n",
      "540/3862\n",
      "541/3862\n",
      "542/3862\n",
      "543/3862\n",
      "544/3862\n",
      "545/3862\n",
      "546/3862\n",
      "547/3862\n",
      "548/3862\n",
      "549/3862\n",
      "550/3862\n",
      "551/3862\n",
      "552/3862\n",
      "553/3862\n",
      "554/3862\n",
      "555/3862\n",
      "556/3862\n",
      "557/3862\n",
      "558/3862\n",
      "559/3862\n",
      "560/3862\n",
      "561/3862\n",
      "562/3862\n",
      "563/3862\n",
      "564/3862\n",
      "565/3862\n",
      "566/3862\n",
      "567/3862\n",
      "568/3862\n",
      "569/3862\n",
      "570/3862\n",
      "571/3862\n",
      "572/3862\n",
      "573/3862\n",
      "574/3862\n",
      "575/3862\n",
      "576/3862\n",
      "577/3862\n",
      "578/3862\n",
      "579/3862\n",
      "580/3862\n",
      "581/3862\n",
      "582/3862\n",
      "583/3862\n",
      "584/3862\n",
      "585/3862\n",
      "586/3862\n",
      "587/3862\n",
      "588/3862\n",
      "589/3862\n",
      "590/3862\n",
      "591/3862\n",
      "592/3862\n",
      "593/3862\n",
      "594/3862\n",
      "595/3862\n",
      "596/3862\n",
      "597/3862\n",
      "598/3862\n",
      "599/3862\n",
      "600/3862\n",
      "601/3862\n",
      "602/3862\n",
      "603/3862\n",
      "604/3862\n",
      "605/3862\n",
      "606/3862\n",
      "607/3862\n",
      "608/3862\n",
      "609/3862\n",
      "610/3862\n",
      "611/3862\n",
      "612/3862\n",
      "613/3862\n",
      "614/3862\n",
      "615/3862\n",
      "616/3862\n",
      "617/3862\n",
      "618/3862\n",
      "619/3862\n",
      "620/3862\n",
      "621/3862\n",
      "622/3862\n",
      "623/3862\n",
      "624/3862\n",
      "625/3862\n",
      "626/3862\n",
      "627/3862\n",
      "628/3862\n",
      "629/3862\n",
      "630/3862\n",
      "631/3862\n",
      "632/3862\n",
      "633/3862\n",
      "634/3862\n",
      "635/3862\n",
      "636/3862\n",
      "637/3862\n",
      "638/3862\n",
      "639/3862\n",
      "640/3862\n",
      "641/3862\n",
      "642/3862\n",
      "643/3862\n",
      "644/3862\n",
      "645/3862\n",
      "646/3862\n",
      "647/3862\n",
      "648/3862\n",
      "649/3862\n",
      "650/3862\n",
      "651/3862\n",
      "652/3862\n",
      "653/3862\n",
      "654/3862\n",
      "655/3862\n",
      "656/3862\n",
      "657/3862\n",
      "658/3862\n",
      "659/3862\n",
      "660/3862\n",
      "661/3862\n",
      "662/3862\n",
      "663/3862\n",
      "664/3862\n",
      "665/3862\n",
      "666/3862\n",
      "667/3862\n",
      "668/3862\n",
      "669/3862\n",
      "670/3862\n",
      "671/3862\n",
      "672/3862\n",
      "673/3862\n",
      "674/3862\n",
      "675/3862\n",
      "676/3862\n",
      "677/3862\n",
      "678/3862\n",
      "679/3862\n",
      "680/3862\n",
      "681/3862\n",
      "682/3862\n",
      "683/3862\n",
      "684/3862\n",
      "685/3862\n",
      "686/3862\n",
      "687/3862\n",
      "688/3862\n",
      "689/3862\n",
      "690/3862\n",
      "691/3862\n",
      "692/3862\n",
      "693/3862\n",
      "694/3862\n",
      "695/3862\n",
      "696/3862\n",
      "697/3862\n",
      "698/3862\n",
      "699/3862\n",
      "700/3862\n",
      "701/3862\n",
      "702/3862\n",
      "703/3862\n",
      "704/3862\n",
      "705/3862\n",
      "706/3862\n",
      "707/3862\n",
      "708/3862\n",
      "709/3862\n",
      "710/3862\n",
      "711/3862\n",
      "712/3862\n",
      "713/3862\n",
      "714/3862\n",
      "715/3862\n",
      "716/3862\n",
      "717/3862\n",
      "718/3862\n",
      "719/3862\n",
      "720/3862\n",
      "721/3862\n",
      "722/3862\n",
      "723/3862\n",
      "724/3862\n",
      "725/3862\n",
      "726/3862\n",
      "727/3862\n",
      "728/3862\n",
      "729/3862\n",
      "730/3862\n",
      "731/3862\n",
      "732/3862\n",
      "733/3862\n",
      "734/3862\n",
      "735/3862\n",
      "736/3862\n",
      "737/3862\n",
      "738/3862\n",
      "739/3862\n",
      "740/3862\n",
      "741/3862\n",
      "742/3862\n",
      "743/3862\n",
      "744/3862\n",
      "745/3862\n",
      "746/3862\n",
      "747/3862\n",
      "748/3862\n",
      "749/3862\n",
      "750/3862\n",
      "751/3862\n",
      "752/3862\n",
      "753/3862\n",
      "754/3862\n",
      "755/3862\n",
      "756/3862\n",
      "757/3862\n",
      "758/3862\n",
      "759/3862\n",
      "760/3862\n",
      "761/3862\n",
      "762/3862\n",
      "763/3862\n",
      "764/3862\n",
      "765/3862\n",
      "766/3862\n",
      "767/3862\n",
      "768/3862\n",
      "769/3862\n",
      "770/3862\n",
      "771/3862\n",
      "772/3862\n",
      "773/3862\n",
      "774/3862\n",
      "775/3862\n",
      "776/3862\n",
      "777/3862\n",
      "778/3862\n",
      "779/3862\n",
      "780/3862\n",
      "781/3862\n",
      "782/3862\n",
      "783/3862\n",
      "784/3862\n",
      "785/3862\n",
      "786/3862\n",
      "787/3862\n",
      "788/3862\n",
      "789/3862\n",
      "790/3862\n",
      "791/3862\n",
      "792/3862\n",
      "793/3862\n",
      "794/3862\n",
      "795/3862\n",
      "796/3862\n",
      "797/3862\n",
      "798/3862\n",
      "799/3862\n",
      "800/3862\n",
      "801/3862\n",
      "802/3862\n",
      "803/3862\n",
      "804/3862\n",
      "805/3862\n",
      "806/3862\n",
      "807/3862\n",
      "808/3862\n",
      "809/3862\n",
      "810/3862\n",
      "811/3862\n",
      "812/3862\n",
      "813/3862\n",
      "814/3862\n",
      "815/3862\n",
      "816/3862\n",
      "817/3862\n",
      "818/3862\n",
      "819/3862\n",
      "820/3862\n",
      "821/3862\n",
      "822/3862\n",
      "823/3862\n",
      "824/3862\n",
      "825/3862\n",
      "826/3862\n",
      "827/3862\n",
      "828/3862\n",
      "829/3862\n",
      "830/3862\n",
      "831/3862\n",
      "832/3862\n",
      "833/3862\n",
      "834/3862\n",
      "835/3862\n",
      "836/3862\n",
      "837/3862\n",
      "838/3862\n",
      "839/3862\n",
      "840/3862\n",
      "841/3862\n",
      "842/3862\n",
      "843/3862\n",
      "844/3862\n",
      "845/3862\n",
      "846/3862\n",
      "847/3862\n",
      "848/3862\n",
      "849/3862\n",
      "850/3862\n",
      "851/3862\n",
      "852/3862\n",
      "853/3862\n",
      "854/3862\n",
      "855/3862\n",
      "856/3862\n",
      "857/3862\n",
      "858/3862\n",
      "859/3862\n",
      "860/3862\n",
      "861/3862\n",
      "862/3862\n",
      "863/3862\n",
      "864/3862\n",
      "865/3862\n",
      "866/3862\n",
      "867/3862\n",
      "868/3862\n",
      "869/3862\n",
      "870/3862\n",
      "871/3862\n",
      "872/3862\n",
      "873/3862\n",
      "874/3862\n",
      "875/3862\n",
      "876/3862\n",
      "877/3862\n",
      "878/3862\n",
      "879/3862\n",
      "880/3862\n",
      "881/3862\n",
      "882/3862\n",
      "883/3862\n",
      "884/3862\n",
      "885/3862\n",
      "886/3862\n",
      "887/3862\n",
      "888/3862\n",
      "889/3862\n",
      "890/3862\n",
      "891/3862\n",
      "892/3862\n",
      "893/3862\n",
      "894/3862\n",
      "895/3862\n",
      "896/3862\n",
      "897/3862\n",
      "898/3862\n",
      "899/3862\n",
      "900/3862\n",
      "901/3862\n",
      "902/3862\n",
      "903/3862\n",
      "904/3862\n",
      "905/3862\n",
      "906/3862\n",
      "907/3862\n",
      "908/3862\n",
      "909/3862\n",
      "910/3862\n",
      "911/3862\n",
      "912/3862\n",
      "913/3862\n",
      "914/3862\n",
      "915/3862\n",
      "916/3862\n",
      "917/3862\n",
      "918/3862\n",
      "919/3862\n",
      "920/3862\n",
      "921/3862\n",
      "922/3862\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "923/3862\n",
      "924/3862\n",
      "925/3862\n",
      "926/3862\n",
      "927/3862\n",
      "928/3862\n",
      "929/3862\n",
      "930/3862\n",
      "931/3862\n",
      "932/3862\n",
      "933/3862\n",
      "934/3862\n",
      "935/3862\n",
      "936/3862\n",
      "937/3862\n",
      "938/3862\n",
      "939/3862\n",
      "940/3862\n",
      "941/3862\n",
      "942/3862\n",
      "943/3862\n",
      "944/3862\n",
      "945/3862\n",
      "946/3862\n",
      "947/3862\n",
      "948/3862\n",
      "949/3862\n",
      "950/3862\n",
      "951/3862\n",
      "952/3862\n",
      "953/3862\n",
      "954/3862\n",
      "955/3862\n",
      "956/3862\n",
      "957/3862\n",
      "958/3862\n",
      "959/3862\n",
      "960/3862\n",
      "961/3862\n",
      "962/3862\n",
      "963/3862\n",
      "964/3862\n",
      "965/3862\n",
      "966/3862\n",
      "967/3862\n",
      "968/3862\n",
      "969/3862\n",
      "970/3862\n",
      "971/3862\n",
      "972/3862\n",
      "973/3862\n",
      "974/3862\n",
      "975/3862\n",
      "976/3862\n",
      "977/3862\n",
      "978/3862\n",
      "979/3862\n",
      "980/3862\n",
      "981/3862\n",
      "982/3862\n",
      "983/3862\n",
      "984/3862\n",
      "985/3862\n",
      "986/3862\n",
      "987/3862\n",
      "988/3862\n",
      "989/3862\n",
      "990/3862\n",
      "991/3862\n",
      "992/3862\n",
      "993/3862\n",
      "994/3862\n",
      "995/3862\n",
      "996/3862\n",
      "997/3862\n",
      "998/3862\n",
      "999/3862\n",
      "1000/3862\n",
      "1001/3862\n",
      "1002/3862\n",
      "1003/3862\n",
      "1004/3862\n",
      "1005/3862\n",
      "1006/3862\n",
      "1007/3862\n",
      "1008/3862\n",
      "1009/3862\n",
      "1010/3862\n",
      "1011/3862\n",
      "1012/3862\n",
      "1013/3862\n",
      "1014/3862\n",
      "1015/3862\n",
      "1016/3862\n",
      "1017/3862\n",
      "1018/3862\n",
      "1019/3862\n",
      "1020/3862\n",
      "1021/3862\n",
      "1022/3862\n",
      "1023/3862\n",
      "1024/3862\n",
      "1025/3862\n",
      "1026/3862\n",
      "1027/3862\n",
      "1028/3862\n",
      "1029/3862\n",
      "1030/3862\n",
      "1031/3862\n",
      "1032/3862\n",
      "1033/3862\n",
      "1034/3862\n",
      "1035/3862\n",
      "1036/3862\n",
      "1037/3862\n",
      "1038/3862\n",
      "1039/3862\n",
      "1040/3862\n",
      "1041/3862\n",
      "1042/3862\n",
      "1043/3862\n",
      "1044/3862\n",
      "1045/3862\n",
      "1046/3862\n",
      "1047/3862\n",
      "1048/3862\n",
      "1049/3862\n",
      "1050/3862\n",
      "1051/3862\n",
      "1052/3862\n",
      "1053/3862\n",
      "1054/3862\n",
      "1055/3862\n",
      "1056/3862\n",
      "1057/3862\n",
      "1058/3862\n",
      "1059/3862\n",
      "1060/3862\n",
      "1061/3862\n",
      "1062/3862\n",
      "1063/3862\n",
      "1064/3862\n",
      "1065/3862\n",
      "1066/3862\n",
      "1067/3862\n",
      "1068/3862\n",
      "1069/3862\n",
      "1070/3862\n",
      "1071/3862\n",
      "1072/3862\n",
      "1073/3862\n",
      "1074/3862\n",
      "1075/3862\n",
      "1076/3862\n",
      "1077/3862\n",
      "1078/3862\n",
      "1079/3862\n",
      "1080/3862\n",
      "1081/3862\n",
      "1082/3862\n",
      "1083/3862\n",
      "1084/3862\n",
      "1085/3862\n",
      "1086/3862\n",
      "1087/3862\n",
      "1088/3862\n",
      "1089/3862\n",
      "1090/3862\n",
      "1091/3862\n",
      "1092/3862\n",
      "1093/3862\n",
      "1094/3862\n",
      "1095/3862\n",
      "1096/3862\n",
      "1097/3862\n",
      "1098/3862\n",
      "1099/3862\n",
      "1100/3862\n",
      "1101/3862\n",
      "1102/3862\n",
      "1103/3862\n",
      "1104/3862\n",
      "1105/3862\n",
      "1106/3862\n",
      "1107/3862\n",
      "1108/3862\n",
      "1109/3862\n",
      "1110/3862\n",
      "1111/3862\n",
      "1112/3862\n",
      "1113/3862\n",
      "1114/3862\n",
      "1115/3862\n",
      "1116/3862\n",
      "1117/3862\n",
      "1118/3862\n",
      "1119/3862\n",
      "1120/3862\n",
      "1121/3862\n",
      "1122/3862\n",
      "1123/3862\n",
      "1124/3862\n",
      "1125/3862\n",
      "1126/3862\n",
      "1127/3862\n",
      "1128/3862\n",
      "1129/3862\n",
      "1130/3862\n",
      "1131/3862\n",
      "1132/3862\n",
      "1133/3862\n",
      "1134/3862\n",
      "1135/3862\n",
      "1136/3862\n",
      "1137/3862\n",
      "1138/3862\n",
      "1139/3862\n",
      "1140/3862\n",
      "1141/3862\n",
      "1142/3862\n",
      "1143/3862\n",
      "1144/3862\n",
      "1145/3862\n",
      "1146/3862\n",
      "1147/3862\n",
      "1148/3862\n",
      "1149/3862\n",
      "1150/3862\n",
      "1151/3862\n",
      "1152/3862\n",
      "1153/3862\n",
      "1154/3862\n",
      "1155/3862\n",
      "1156/3862\n",
      "1157/3862\n",
      "1158/3862\n",
      "1159/3862\n",
      "1160/3862\n",
      "1161/3862\n",
      "1162/3862\n",
      "1163/3862\n",
      "1164/3862\n",
      "1165/3862\n",
      "1166/3862\n",
      "1167/3862\n",
      "1168/3862\n",
      "1169/3862\n",
      "1170/3862\n",
      "1171/3862\n",
      "1172/3862\n",
      "1173/3862\n",
      "1174/3862\n",
      "1175/3862\n",
      "1176/3862\n",
      "1177/3862\n",
      "1178/3862\n",
      "1179/3862\n",
      "1180/3862\n",
      "1181/3862\n",
      "1182/3862\n",
      "1183/3862\n",
      "1184/3862\n",
      "1185/3862\n",
      "1186/3862\n",
      "1187/3862\n",
      "1188/3862\n",
      "1189/3862\n",
      "1190/3862\n",
      "1191/3862\n",
      "1192/3862\n",
      "1193/3862\n",
      "1194/3862\n",
      "1195/3862\n",
      "1196/3862\n",
      "1197/3862\n",
      "1198/3862\n",
      "1199/3862\n",
      "1200/3862\n",
      "1201/3862\n",
      "1202/3862\n",
      "1203/3862\n",
      "1204/3862\n",
      "1205/3862\n",
      "1206/3862\n",
      "1207/3862\n",
      "1208/3862\n",
      "1209/3862\n",
      "1210/3862\n",
      "1211/3862\n",
      "1212/3862\n",
      "1213/3862\n",
      "1214/3862\n",
      "1215/3862\n",
      "1216/3862\n",
      "1217/3862\n",
      "1218/3862\n",
      "1219/3862\n",
      "1220/3862\n",
      "1221/3862\n",
      "1222/3862\n",
      "1223/3862\n",
      "1224/3862\n",
      "1225/3862\n",
      "1226/3862\n",
      "1227/3862\n",
      "1228/3862\n",
      "1229/3862\n",
      "1230/3862\n",
      "1231/3862\n",
      "1232/3862\n",
      "1233/3862\n",
      "1234/3862\n",
      "1235/3862\n",
      "1236/3862\n",
      "1237/3862\n",
      "1238/3862\n",
      "1239/3862\n",
      "1240/3862\n",
      "1241/3862\n",
      "1242/3862\n",
      "1243/3862\n",
      "1244/3862\n",
      "1245/3862\n",
      "1246/3862\n",
      "1247/3862\n",
      "1248/3862\n",
      "1249/3862\n",
      "1250/3862\n",
      "1251/3862\n",
      "1252/3862\n",
      "1253/3862\n",
      "1254/3862\n",
      "1255/3862\n",
      "1256/3862\n",
      "1257/3862\n",
      "1258/3862\n",
      "1259/3862\n",
      "1260/3862\n",
      "1261/3862\n",
      "1262/3862\n",
      "1263/3862\n",
      "1264/3862\n",
      "1265/3862\n",
      "1266/3862\n",
      "1267/3862\n",
      "1268/3862\n",
      "1269/3862\n",
      "1270/3862\n",
      "1271/3862\n",
      "1272/3862\n",
      "1273/3862\n",
      "1274/3862\n",
      "1275/3862\n",
      "1276/3862\n",
      "1277/3862\n",
      "1278/3862\n",
      "1279/3862\n",
      "1280/3862\n",
      "1281/3862\n",
      "1282/3862\n",
      "1283/3862\n",
      "1284/3862\n",
      "1285/3862\n",
      "1286/3862\n",
      "1287/3862\n",
      "1288/3862\n",
      "1289/3862\n",
      "1290/3862\n",
      "1291/3862\n",
      "1292/3862\n",
      "1293/3862\n",
      "1294/3862\n",
      "1295/3862\n",
      "1296/3862\n",
      "1297/3862\n",
      "1298/3862\n",
      "1299/3862\n",
      "1300/3862\n",
      "1301/3862\n",
      "1302/3862\n",
      "1303/3862\n",
      "1304/3862\n",
      "1305/3862\n",
      "1306/3862\n",
      "1307/3862\n",
      "1308/3862\n",
      "1309/3862\n",
      "1310/3862\n",
      "1311/3862\n",
      "1312/3862\n",
      "1313/3862\n",
      "1314/3862\n",
      "1315/3862\n",
      "1316/3862\n",
      "1317/3862\n",
      "1318/3862\n",
      "1319/3862\n",
      "1320/3862\n",
      "1321/3862\n",
      "1322/3862\n",
      "1323/3862\n",
      "1324/3862\n",
      "1325/3862\n",
      "1326/3862\n",
      "1327/3862\n",
      "1328/3862\n",
      "1329/3862\n",
      "1330/3862\n",
      "1331/3862\n",
      "1332/3862\n",
      "1333/3862\n",
      "1334/3862\n",
      "1335/3862\n",
      "1336/3862\n",
      "1337/3862\n",
      "1338/3862\n",
      "1339/3862\n",
      "1340/3862\n",
      "1341/3862\n",
      "1342/3862\n",
      "1343/3862\n",
      "1344/3862\n",
      "1345/3862\n",
      "1346/3862\n",
      "1347/3862\n",
      "1348/3862\n",
      "1349/3862\n",
      "1350/3862\n",
      "1351/3862\n",
      "1352/3862\n",
      "1353/3862\n",
      "1354/3862\n",
      "1355/3862\n",
      "1356/3862\n",
      "1357/3862\n",
      "1358/3862\n",
      "1359/3862\n",
      "1360/3862\n",
      "1361/3862\n",
      "1362/3862\n",
      "1363/3862\n",
      "1364/3862\n",
      "1365/3862\n",
      "1366/3862\n",
      "1367/3862\n",
      "1368/3862\n",
      "1369/3862\n",
      "1370/3862\n",
      "1371/3862\n",
      "1372/3862\n",
      "1373/3862\n",
      "1374/3862\n",
      "1375/3862\n",
      "1376/3862\n",
      "1377/3862\n",
      "1378/3862\n",
      "1379/3862\n",
      "1380/3862\n",
      "1381/3862\n",
      "1382/3862\n",
      "1383/3862\n",
      "1384/3862\n",
      "1385/3862\n",
      "1386/3862\n",
      "1387/3862\n",
      "1388/3862\n",
      "1389/3862\n",
      "1390/3862\n",
      "1391/3862\n",
      "1392/3862\n",
      "1393/3862\n",
      "1394/3862\n",
      "1395/3862\n",
      "1396/3862\n",
      "1397/3862\n",
      "1398/3862\n",
      "1399/3862\n",
      "1400/3862\n",
      "1401/3862\n",
      "1402/3862\n",
      "1403/3862\n",
      "1404/3862\n",
      "1405/3862\n",
      "1406/3862\n",
      "1407/3862\n",
      "1408/3862\n",
      "1409/3862\n",
      "1410/3862\n",
      "1411/3862\n",
      "1412/3862\n",
      "1413/3862\n",
      "1414/3862\n",
      "1415/3862\n",
      "1416/3862\n",
      "1417/3862\n",
      "1418/3862\n",
      "1419/3862\n",
      "1420/3862\n",
      "1421/3862\n",
      "1422/3862\n",
      "1423/3862\n",
      "1424/3862\n",
      "1425/3862\n",
      "1426/3862\n",
      "1427/3862\n",
      "1428/3862\n",
      "1429/3862\n",
      "1430/3862\n",
      "1431/3862\n",
      "1432/3862\n",
      "1433/3862\n",
      "1434/3862\n",
      "1435/3862\n",
      "1436/3862\n",
      "1437/3862\n",
      "1438/3862\n",
      "1439/3862\n",
      "1440/3862\n",
      "1441/3862\n",
      "1442/3862\n",
      "1443/3862\n",
      "1444/3862\n",
      "1445/3862\n",
      "1446/3862\n",
      "1447/3862\n",
      "1448/3862\n",
      "1449/3862\n",
      "1450/3862\n",
      "1451/3862\n",
      "1452/3862\n",
      "1453/3862\n",
      "1454/3862\n",
      "1455/3862\n",
      "1456/3862\n",
      "1457/3862\n",
      "1458/3862\n",
      "1459/3862\n",
      "1460/3862\n",
      "1461/3862\n",
      "1462/3862\n",
      "1463/3862\n",
      "1464/3862\n",
      "1465/3862\n",
      "1466/3862\n",
      "1467/3862\n",
      "1468/3862\n",
      "1469/3862\n",
      "1470/3862\n",
      "1471/3862\n",
      "1472/3862\n",
      "1473/3862\n",
      "1474/3862\n",
      "1475/3862\n",
      "1476/3862\n",
      "1477/3862\n",
      "1478/3862\n",
      "1479/3862\n",
      "1480/3862\n",
      "1481/3862\n",
      "1482/3862\n",
      "1483/3862\n",
      "1484/3862\n",
      "1485/3862\n",
      "1486/3862\n",
      "1487/3862\n",
      "1488/3862\n",
      "1489/3862\n",
      "1490/3862\n",
      "1491/3862\n",
      "1492/3862\n",
      "1493/3862\n",
      "1494/3862\n",
      "1495/3862\n",
      "1496/3862\n",
      "1497/3862\n",
      "1498/3862\n",
      "1499/3862\n",
      "1500/3862\n",
      "1501/3862\n",
      "1502/3862\n",
      "1503/3862\n",
      "1504/3862\n",
      "1505/3862\n",
      "1506/3862\n",
      "1507/3862\n",
      "1508/3862\n",
      "1509/3862\n",
      "1510/3862\n",
      "1511/3862\n",
      "1512/3862\n",
      "1513/3862\n",
      "1514/3862\n",
      "1515/3862\n",
      "1516/3862\n",
      "1517/3862\n",
      "1518/3862\n",
      "1519/3862\n",
      "1520/3862\n",
      "1521/3862\n",
      "1522/3862\n",
      "1523/3862\n",
      "1524/3862\n",
      "1525/3862\n",
      "1526/3862\n",
      "1527/3862\n",
      "1528/3862\n",
      "1529/3862\n",
      "1530/3862\n",
      "1531/3862\n",
      "1532/3862\n",
      "1533/3862\n",
      "1534/3862\n",
      "1535/3862\n",
      "1536/3862\n",
      "1537/3862\n",
      "1538/3862\n",
      "1539/3862\n",
      "1540/3862\n",
      "1541/3862\n",
      "1542/3862\n",
      "1543/3862\n",
      "1544/3862\n",
      "1545/3862\n",
      "1546/3862\n",
      "1547/3862\n",
      "1548/3862\n",
      "1549/3862\n",
      "1550/3862\n",
      "1551/3862\n",
      "1552/3862\n",
      "1553/3862\n",
      "1554/3862\n",
      "1555/3862\n",
      "1556/3862\n",
      "1557/3862\n",
      "1558/3862\n",
      "1559/3862\n",
      "1560/3862\n",
      "1561/3862\n",
      "1562/3862\n",
      "1563/3862\n",
      "1564/3862\n",
      "1565/3862\n",
      "1566/3862\n",
      "1567/3862\n",
      "1568/3862\n",
      "1569/3862\n",
      "1570/3862\n",
      "1571/3862\n",
      "1572/3862\n",
      "1573/3862\n",
      "1574/3862\n",
      "1575/3862\n",
      "1576/3862\n",
      "1577/3862\n",
      "1578/3862\n",
      "1579/3862\n",
      "1580/3862\n",
      "1581/3862\n",
      "1582/3862\n",
      "1583/3862\n",
      "1584/3862\n",
      "1585/3862\n",
      "1586/3862\n",
      "1587/3862\n",
      "1588/3862\n",
      "1589/3862\n",
      "1590/3862\n",
      "1591/3862\n",
      "1592/3862\n",
      "1593/3862\n",
      "1594/3862\n",
      "1595/3862\n",
      "1596/3862\n",
      "1597/3862\n",
      "1598/3862\n",
      "1599/3862\n",
      "1600/3862\n",
      "1601/3862\n",
      "1602/3862\n",
      "1603/3862\n",
      "1604/3862\n",
      "1605/3862\n",
      "1606/3862\n",
      "1607/3862\n",
      "1608/3862\n",
      "1609/3862\n",
      "1610/3862\n",
      "1611/3862\n",
      "1612/3862\n",
      "1613/3862\n",
      "1614/3862\n",
      "1615/3862\n",
      "1616/3862\n",
      "1617/3862\n",
      "1618/3862\n",
      "1619/3862\n",
      "1620/3862\n",
      "1621/3862\n",
      "1622/3862\n",
      "1623/3862\n",
      "1624/3862\n",
      "1625/3862\n",
      "1626/3862\n",
      "1627/3862\n",
      "1628/3862\n",
      "1629/3862\n",
      "1630/3862\n",
      "1631/3862\n",
      "1632/3862\n",
      "1633/3862\n",
      "1634/3862\n",
      "1635/3862\n",
      "1636/3862\n",
      "1637/3862\n",
      "1638/3862\n",
      "1639/3862\n",
      "1640/3862\n",
      "1641/3862\n",
      "1642/3862\n",
      "1643/3862\n",
      "1644/3862\n",
      "1645/3862\n",
      "1646/3862\n",
      "1647/3862\n",
      "1648/3862\n",
      "1649/3862\n",
      "1650/3862\n",
      "1651/3862\n",
      "1652/3862\n",
      "1653/3862\n",
      "1654/3862\n",
      "1655/3862\n",
      "1656/3862\n",
      "1657/3862\n",
      "1658/3862\n",
      "1659/3862\n",
      "1660/3862\n",
      "1661/3862\n",
      "1662/3862\n",
      "1663/3862\n",
      "1664/3862\n",
      "1665/3862\n",
      "1666/3862\n",
      "1667/3862\n",
      "1668/3862\n",
      "1669/3862\n",
      "1670/3862\n",
      "1671/3862\n",
      "1672/3862\n",
      "1673/3862\n",
      "1674/3862\n",
      "1675/3862\n",
      "1676/3862\n",
      "1677/3862\n",
      "1678/3862\n",
      "1679/3862\n",
      "1680/3862\n",
      "1681/3862\n",
      "1682/3862\n",
      "1683/3862\n",
      "1684/3862\n",
      "1685/3862\n",
      "1686/3862\n",
      "1687/3862\n",
      "1688/3862\n",
      "1689/3862\n",
      "1690/3862\n",
      "1691/3862\n",
      "1692/3862\n",
      "1693/3862\n",
      "1694/3862\n",
      "1695/3862\n",
      "1696/3862\n",
      "1697/3862\n",
      "1698/3862\n",
      "1699/3862\n",
      "1700/3862\n",
      "1701/3862\n",
      "1702/3862\n",
      "1703/3862\n",
      "1704/3862\n",
      "1705/3862\n",
      "1706/3862\n",
      "1707/3862\n",
      "1708/3862\n",
      "1709/3862\n",
      "1710/3862\n",
      "1711/3862\n",
      "1712/3862\n",
      "1713/3862\n",
      "1714/3862\n",
      "1715/3862\n",
      "1716/3862\n",
      "1717/3862\n",
      "1718/3862\n",
      "1719/3862\n",
      "1720/3862\n",
      "1721/3862\n",
      "1722/3862\n",
      "1723/3862\n",
      "1724/3862\n",
      "1725/3862\n",
      "1726/3862\n",
      "1727/3862\n",
      "1728/3862\n",
      "1729/3862\n",
      "1730/3862\n",
      "1731/3862\n",
      "1732/3862\n",
      "1733/3862\n",
      "1734/3862\n",
      "1735/3862\n",
      "1736/3862\n",
      "1737/3862\n",
      "1738/3862\n",
      "1739/3862\n",
      "1740/3862\n",
      "1741/3862\n",
      "1742/3862\n",
      "1743/3862\n",
      "1744/3862\n",
      "1745/3862\n",
      "1746/3862\n",
      "1747/3862\n",
      "1748/3862\n",
      "1749/3862\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1750/3862\n",
      "1751/3862\n",
      "1752/3862\n",
      "1753/3862\n",
      "1754/3862\n",
      "1755/3862\n",
      "1756/3862\n",
      "1757/3862\n",
      "1758/3862\n",
      "1759/3862\n",
      "1760/3862\n",
      "1761/3862\n",
      "1762/3862\n",
      "1763/3862\n",
      "1764/3862\n",
      "1765/3862\n",
      "1766/3862\n",
      "1767/3862\n",
      "1768/3862\n",
      "1769/3862\n",
      "1770/3862\n",
      "1771/3862\n",
      "1772/3862\n",
      "1773/3862\n",
      "1774/3862\n",
      "1775/3862\n",
      "1776/3862\n",
      "1777/3862\n",
      "1778/3862\n",
      "1779/3862\n",
      "1780/3862\n",
      "1781/3862\n",
      "1782/3862\n",
      "1783/3862\n",
      "1784/3862\n",
      "1785/3862\n",
      "1786/3862\n",
      "1787/3862\n",
      "1788/3862\n",
      "1789/3862\n",
      "1790/3862\n",
      "1791/3862\n",
      "1792/3862\n",
      "1793/3862\n",
      "1794/3862\n",
      "1795/3862\n",
      "1796/3862\n",
      "1797/3862\n",
      "1798/3862\n",
      "1799/3862\n",
      "1800/3862\n",
      "1801/3862\n",
      "1802/3862\n",
      "1803/3862\n",
      "1804/3862\n",
      "1805/3862\n",
      "1806/3862\n",
      "1807/3862\n",
      "1808/3862\n",
      "1809/3862\n",
      "1810/3862\n",
      "1811/3862\n",
      "1812/3862\n",
      "1813/3862\n",
      "1814/3862\n",
      "1815/3862\n",
      "1816/3862\n",
      "1817/3862\n",
      "1818/3862\n",
      "1819/3862\n",
      "1820/3862\n",
      "1821/3862\n",
      "1822/3862\n",
      "1823/3862\n",
      "1824/3862\n",
      "1825/3862\n",
      "1826/3862\n",
      "1827/3862\n",
      "1828/3862\n",
      "1829/3862\n",
      "1830/3862\n",
      "1831/3862\n",
      "1832/3862\n",
      "1833/3862\n",
      "1834/3862\n",
      "1835/3862\n",
      "1836/3862\n",
      "1837/3862\n",
      "1838/3862\n",
      "1839/3862\n",
      "1840/3862\n",
      "1841/3862\n",
      "1842/3862\n",
      "1843/3862\n",
      "1844/3862\n",
      "1845/3862\n",
      "1846/3862\n",
      "1847/3862\n",
      "1848/3862\n",
      "1849/3862\n",
      "1850/3862\n",
      "1851/3862\n",
      "1852/3862\n",
      "1853/3862\n",
      "1854/3862\n",
      "1855/3862\n",
      "1856/3862\n",
      "1857/3862\n",
      "1858/3862\n",
      "1859/3862\n",
      "1860/3862\n",
      "1861/3862\n",
      "1862/3862\n",
      "1863/3862\n",
      "1864/3862\n",
      "1865/3862\n",
      "1866/3862\n",
      "1867/3862\n",
      "1868/3862\n",
      "1869/3862\n",
      "1870/3862\n",
      "1871/3862\n",
      "1872/3862\n",
      "1873/3862\n",
      "1874/3862\n",
      "1875/3862\n",
      "1876/3862\n",
      "1877/3862\n",
      "1878/3862\n",
      "1879/3862\n",
      "1880/3862\n",
      "1881/3862\n",
      "1882/3862\n",
      "1883/3862\n",
      "1884/3862\n",
      "1885/3862\n",
      "1886/3862\n",
      "1887/3862\n",
      "1888/3862\n",
      "1889/3862\n",
      "1890/3862\n",
      "1891/3862\n",
      "1892/3862\n",
      "1893/3862\n",
      "1894/3862\n",
      "1895/3862\n",
      "1896/3862\n",
      "1897/3862\n",
      "1898/3862\n",
      "1899/3862\n",
      "1900/3862\n",
      "1901/3862\n",
      "1902/3862\n",
      "1903/3862\n",
      "1904/3862\n",
      "1905/3862\n",
      "1906/3862\n",
      "1907/3862\n",
      "1908/3862\n",
      "1909/3862\n",
      "1910/3862\n",
      "1911/3862\n",
      "1912/3862\n",
      "1913/3862\n",
      "1914/3862\n",
      "1915/3862\n",
      "1916/3862\n",
      "1917/3862\n",
      "1918/3862\n",
      "1919/3862\n",
      "1920/3862\n",
      "1921/3862\n",
      "1922/3862\n",
      "1923/3862\n",
      "1924/3862\n",
      "1925/3862\n",
      "1926/3862\n",
      "1927/3862\n",
      "1928/3862\n",
      "1929/3862\n",
      "1930/3862\n",
      "1931/3862\n",
      "1932/3862\n",
      "1933/3862\n",
      "1934/3862\n",
      "1935/3862\n",
      "1936/3862\n",
      "1937/3862\n",
      "1938/3862\n",
      "1939/3862\n",
      "1940/3862\n",
      "1941/3862\n",
      "1942/3862\n",
      "1943/3862\n",
      "1944/3862\n",
      "1945/3862\n",
      "1946/3862\n",
      "1947/3862\n",
      "1948/3862\n",
      "1949/3862\n",
      "1950/3862\n",
      "1951/3862\n",
      "1952/3862\n",
      "1953/3862\n",
      "1954/3862\n",
      "1955/3862\n",
      "1956/3862\n",
      "1957/3862\n",
      "1958/3862\n",
      "1959/3862\n",
      "1960/3862\n",
      "1961/3862\n",
      "1962/3862\n",
      "1963/3862\n",
      "1964/3862\n",
      "1965/3862\n",
      "1966/3862\n",
      "1967/3862\n",
      "1968/3862\n",
      "1969/3862\n",
      "1970/3862\n",
      "1971/3862\n",
      "1972/3862\n",
      "1973/3862\n",
      "1974/3862\n",
      "1975/3862\n",
      "1976/3862\n",
      "1977/3862\n",
      "1978/3862\n",
      "1979/3862\n",
      "1980/3862\n",
      "1981/3862\n",
      "1982/3862\n",
      "1983/3862\n",
      "1984/3862\n",
      "1985/3862\n",
      "1986/3862\n",
      "1987/3862\n",
      "1988/3862\n",
      "1989/3862\n",
      "1990/3862\n",
      "1991/3862\n",
      "1992/3862\n",
      "1993/3862\n",
      "1994/3862\n",
      "1995/3862\n",
      "1996/3862\n",
      "1997/3862\n",
      "1998/3862\n",
      "1999/3862\n",
      "2000/3862\n",
      "2001/3862\n",
      "2002/3862\n",
      "2003/3862\n",
      "2004/3862\n",
      "2005/3862\n",
      "2006/3862\n",
      "2007/3862\n",
      "2008/3862\n",
      "2009/3862\n",
      "2010/3862\n",
      "2011/3862\n",
      "2012/3862\n",
      "2013/3862\n",
      "2014/3862\n",
      "2015/3862\n",
      "2016/3862\n",
      "2017/3862\n",
      "2018/3862\n",
      "2019/3862\n",
      "2020/3862\n",
      "2021/3862\n",
      "2022/3862\n",
      "2023/3862\n",
      "2024/3862\n",
      "2025/3862\n",
      "2026/3862\n",
      "2027/3862\n",
      "2028/3862\n",
      "2029/3862\n",
      "2030/3862\n",
      "2031/3862\n",
      "2032/3862\n",
      "2033/3862\n",
      "2034/3862\n",
      "2035/3862\n",
      "2036/3862\n",
      "2037/3862\n",
      "2038/3862\n",
      "2039/3862\n",
      "2040/3862\n",
      "2041/3862\n",
      "2042/3862\n",
      "2043/3862\n",
      "2044/3862\n",
      "2045/3862\n",
      "2046/3862\n",
      "2047/3862\n",
      "2048/3862\n",
      "2049/3862\n",
      "2050/3862\n",
      "2051/3862\n",
      "2052/3862\n",
      "2053/3862\n",
      "2054/3862\n",
      "2055/3862\n",
      "2056/3862\n",
      "2057/3862\n",
      "2058/3862\n",
      "2059/3862\n",
      "2060/3862\n",
      "2061/3862\n",
      "2062/3862\n",
      "2063/3862\n",
      "2064/3862\n",
      "2065/3862\n",
      "2066/3862\n",
      "2067/3862\n",
      "2068/3862\n",
      "2069/3862\n",
      "2070/3862\n",
      "2071/3862\n",
      "2072/3862\n",
      "2073/3862\n",
      "2074/3862\n",
      "2075/3862\n",
      "2076/3862\n",
      "2077/3862\n",
      "2078/3862\n",
      "2079/3862\n",
      "2080/3862\n",
      "2081/3862\n",
      "2082/3862\n",
      "2083/3862\n",
      "2084/3862\n",
      "2085/3862\n",
      "2086/3862\n",
      "2087/3862\n",
      "2088/3862\n",
      "2089/3862\n",
      "2090/3862\n",
      "2091/3862\n",
      "2092/3862\n",
      "2093/3862\n",
      "2094/3862\n",
      "2095/3862\n",
      "2096/3862\n",
      "2097/3862\n",
      "2098/3862\n",
      "2099/3862\n",
      "2100/3862\n",
      "2101/3862\n",
      "2102/3862\n",
      "2103/3862\n",
      "2104/3862\n",
      "2105/3862\n",
      "2106/3862\n",
      "2107/3862\n",
      "2108/3862\n",
      "2109/3862\n",
      "2110/3862\n",
      "2111/3862\n",
      "2112/3862\n",
      "2113/3862\n",
      "2114/3862\n",
      "2115/3862\n",
      "2116/3862\n",
      "2117/3862\n",
      "2118/3862\n",
      "2119/3862\n",
      "2120/3862\n",
      "2121/3862\n",
      "2122/3862\n",
      "2123/3862\n",
      "2124/3862\n",
      "2125/3862\n",
      "2126/3862\n",
      "2127/3862\n",
      "2128/3862\n",
      "2129/3862\n",
      "2130/3862\n",
      "2131/3862\n",
      "2132/3862\n",
      "2133/3862\n",
      "2134/3862\n",
      "2135/3862\n",
      "2136/3862\n",
      "2137/3862\n",
      "2138/3862\n",
      "2139/3862\n",
      "2140/3862\n",
      "2141/3862\n",
      "2142/3862\n",
      "2143/3862\n",
      "2144/3862\n",
      "2145/3862\n",
      "2146/3862\n",
      "2147/3862\n",
      "2148/3862\n",
      "2149/3862\n",
      "2150/3862\n",
      "2151/3862\n",
      "2152/3862\n",
      "2153/3862\n",
      "2154/3862\n",
      "2155/3862\n",
      "2156/3862\n",
      "2157/3862\n",
      "2158/3862\n",
      "2159/3862\n",
      "2160/3862\n",
      "2161/3862\n",
      "2162/3862\n",
      "2163/3862\n",
      "2164/3862\n",
      "2165/3862\n",
      "2166/3862\n",
      "2167/3862\n",
      "2168/3862\n",
      "2169/3862\n",
      "2170/3862\n",
      "2171/3862\n",
      "2172/3862\n",
      "2173/3862\n",
      "2174/3862\n",
      "2175/3862\n",
      "2176/3862\n",
      "2177/3862\n",
      "2178/3862\n",
      "2179/3862\n",
      "2180/3862\n",
      "2181/3862\n",
      "2182/3862\n",
      "2183/3862\n",
      "2184/3862\n",
      "2185/3862\n",
      "2186/3862\n",
      "2187/3862\n",
      "2188/3862\n",
      "2189/3862\n",
      "2190/3862\n",
      "2191/3862\n",
      "2192/3862\n",
      "2193/3862\n",
      "2194/3862\n",
      "2195/3862\n",
      "2196/3862\n",
      "2197/3862\n",
      "2198/3862\n",
      "2199/3862\n",
      "2200/3862\n",
      "2201/3862\n",
      "2202/3862\n",
      "2203/3862\n",
      "2204/3862\n",
      "2205/3862\n",
      "2206/3862\n",
      "2207/3862\n",
      "2208/3862\n",
      "2209/3862\n",
      "2210/3862\n",
      "2211/3862\n",
      "2212/3862\n",
      "2213/3862\n",
      "2214/3862\n",
      "2215/3862\n",
      "2216/3862\n",
      "2217/3862\n",
      "2218/3862\n",
      "2219/3862\n",
      "2220/3862\n",
      "2221/3862\n",
      "2222/3862\n",
      "2223/3862\n",
      "2224/3862\n",
      "2225/3862\n",
      "2226/3862\n",
      "2227/3862\n",
      "2228/3862\n",
      "2229/3862\n",
      "2230/3862\n",
      "2231/3862\n",
      "2232/3862\n",
      "2233/3862\n",
      "2234/3862\n",
      "2235/3862\n",
      "2236/3862\n",
      "2237/3862\n",
      "2238/3862\n",
      "2239/3862\n",
      "2240/3862\n",
      "2241/3862\n",
      "2242/3862\n",
      "2243/3862\n",
      "2244/3862\n",
      "2245/3862\n",
      "2246/3862\n",
      "2247/3862\n",
      "2248/3862\n",
      "2249/3862\n",
      "2250/3862\n",
      "2251/3862\n",
      "2252/3862\n",
      "2253/3862\n",
      "2254/3862\n",
      "2255/3862\n",
      "2256/3862\n",
      "2257/3862\n",
      "2258/3862\n",
      "2259/3862\n",
      "2260/3862\n",
      "2261/3862\n",
      "2262/3862\n",
      "2263/3862\n",
      "2264/3862\n",
      "2265/3862\n",
      "2266/3862\n",
      "2267/3862\n",
      "2268/3862\n",
      "2269/3862\n",
      "2270/3862\n",
      "2271/3862\n",
      "2272/3862\n",
      "2273/3862\n",
      "2274/3862\n",
      "2275/3862\n",
      "2276/3862\n",
      "2277/3862\n",
      "2278/3862\n",
      "2279/3862\n",
      "2280/3862\n",
      "2281/3862\n",
      "2282/3862\n",
      "2283/3862\n",
      "2284/3862\n",
      "2285/3862\n",
      "2286/3862\n",
      "2287/3862\n",
      "2288/3862\n",
      "2289/3862\n",
      "2290/3862\n",
      "2291/3862\n",
      "2292/3862\n",
      "2293/3862\n",
      "2294/3862\n",
      "2295/3862\n",
      "2296/3862\n",
      "2297/3862\n",
      "2298/3862\n",
      "2299/3862\n",
      "2300/3862\n",
      "2301/3862\n",
      "2302/3862\n",
      "2303/3862\n",
      "2304/3862\n",
      "2305/3862\n",
      "2306/3862\n",
      "2307/3862\n",
      "2308/3862\n",
      "2309/3862\n",
      "2310/3862\n",
      "2311/3862\n",
      "2312/3862\n",
      "2313/3862\n",
      "2314/3862\n",
      "2315/3862\n",
      "2316/3862\n",
      "2317/3862\n",
      "2318/3862\n",
      "2319/3862\n",
      "2320/3862\n",
      "2321/3862\n",
      "2322/3862\n",
      "2323/3862\n",
      "2324/3862\n",
      "2325/3862\n",
      "2326/3862\n",
      "2327/3862\n",
      "2328/3862\n",
      "2329/3862\n",
      "2330/3862\n",
      "2331/3862\n",
      "2332/3862\n",
      "2333/3862\n",
      "2334/3862\n",
      "2335/3862\n",
      "2336/3862\n",
      "2337/3862\n",
      "2338/3862\n",
      "2339/3862\n",
      "2340/3862\n",
      "2341/3862\n",
      "2342/3862\n",
      "2343/3862\n",
      "2344/3862\n",
      "2345/3862\n",
      "2346/3862\n",
      "2347/3862\n",
      "2348/3862\n",
      "2349/3862\n",
      "2350/3862\n",
      "2351/3862\n",
      "2352/3862\n",
      "2353/3862\n",
      "2354/3862\n",
      "2355/3862\n",
      "2356/3862\n",
      "2357/3862\n",
      "2358/3862\n",
      "2359/3862\n",
      "2360/3862\n",
      "2361/3862\n",
      "2362/3862\n",
      "2363/3862\n",
      "2364/3862\n",
      "2365/3862\n",
      "2366/3862\n",
      "2367/3862\n",
      "2368/3862\n",
      "2369/3862\n",
      "2370/3862\n",
      "2371/3862\n",
      "2372/3862\n",
      "2373/3862\n",
      "2374/3862\n",
      "2375/3862\n",
      "2376/3862\n",
      "2377/3862\n",
      "2378/3862\n",
      "2379/3862\n",
      "2380/3862\n",
      "2381/3862\n",
      "2382/3862\n",
      "2383/3862\n",
      "2384/3862\n",
      "2385/3862\n",
      "2386/3862\n",
      "2387/3862\n",
      "2388/3862\n",
      "2389/3862\n",
      "2390/3862\n",
      "2391/3862\n",
      "2392/3862\n",
      "2393/3862\n",
      "2394/3862\n",
      "2395/3862\n",
      "2396/3862\n",
      "2397/3862\n",
      "2398/3862\n",
      "2399/3862\n",
      "2400/3862\n",
      "2401/3862\n",
      "2402/3862\n",
      "2403/3862\n",
      "2404/3862\n",
      "2405/3862\n",
      "2406/3862\n",
      "2407/3862\n",
      "2408/3862\n",
      "2409/3862\n",
      "2410/3862\n",
      "2411/3862\n",
      "2412/3862\n",
      "2413/3862\n",
      "2414/3862\n",
      "2415/3862\n",
      "2416/3862\n",
      "2417/3862\n",
      "2418/3862\n",
      "2419/3862\n",
      "2420/3862\n",
      "2421/3862\n",
      "2422/3862\n",
      "2423/3862\n",
      "2424/3862\n",
      "2425/3862\n",
      "2426/3862\n",
      "2427/3862\n",
      "2428/3862\n",
      "2429/3862\n",
      "2430/3862\n",
      "2431/3862\n",
      "2432/3862\n",
      "2433/3862\n",
      "2434/3862\n",
      "2435/3862\n",
      "2436/3862\n",
      "2437/3862\n",
      "2438/3862\n",
      "2439/3862\n",
      "2440/3862\n",
      "2441/3862\n",
      "2442/3862\n",
      "2443/3862\n",
      "2444/3862\n",
      "2445/3862\n",
      "2446/3862\n",
      "2447/3862\n",
      "2448/3862\n",
      "2449/3862\n",
      "2450/3862\n",
      "2451/3862\n",
      "2452/3862\n",
      "2453/3862\n",
      "2454/3862\n",
      "2455/3862\n",
      "2456/3862\n",
      "2457/3862\n",
      "2458/3862\n",
      "2459/3862\n",
      "2460/3862\n",
      "2461/3862\n",
      "2462/3862\n",
      "2463/3862\n",
      "2464/3862\n",
      "2465/3862\n",
      "2466/3862\n",
      "2467/3862\n",
      "2468/3862\n",
      "2469/3862\n",
      "2470/3862\n",
      "2471/3862\n",
      "2472/3862\n",
      "2473/3862\n",
      "2474/3862\n",
      "2475/3862\n",
      "2476/3862\n",
      "2477/3862\n",
      "2478/3862\n",
      "2479/3862\n",
      "2480/3862\n",
      "2481/3862\n",
      "2482/3862\n",
      "2483/3862\n",
      "2484/3862\n",
      "2485/3862\n",
      "2486/3862\n",
      "2487/3862\n",
      "2488/3862\n",
      "2489/3862\n",
      "2490/3862\n",
      "2491/3862\n",
      "2492/3862\n",
      "2493/3862\n",
      "2494/3862\n",
      "2495/3862\n",
      "2496/3862\n",
      "2497/3862\n",
      "2498/3862\n",
      "2499/3862\n",
      "2500/3862\n",
      "2501/3862\n",
      "2502/3862\n",
      "2503/3862\n",
      "2504/3862\n",
      "2505/3862\n",
      "2506/3862\n",
      "2507/3862\n",
      "2508/3862\n",
      "2509/3862\n",
      "2510/3862\n",
      "2511/3862\n",
      "2512/3862\n",
      "2513/3862\n",
      "2514/3862\n",
      "2515/3862\n",
      "2516/3862\n",
      "2517/3862\n",
      "2518/3862\n",
      "2519/3862\n",
      "2520/3862\n",
      "2521/3862\n",
      "2522/3862\n",
      "2523/3862\n",
      "2524/3862\n",
      "2525/3862\n",
      "2526/3862\n",
      "2527/3862\n",
      "2528/3862\n",
      "2529/3862\n",
      "2530/3862\n",
      "2531/3862\n",
      "2532/3862\n",
      "2533/3862\n",
      "2534/3862\n",
      "2535/3862\n",
      "2536/3862\n",
      "2537/3862\n",
      "2538/3862\n",
      "2539/3862\n",
      "2540/3862\n",
      "2541/3862\n",
      "2542/3862\n",
      "2543/3862\n",
      "2544/3862\n",
      "2545/3862\n",
      "2546/3862\n",
      "2547/3862\n",
      "2548/3862\n",
      "2549/3862\n",
      "2550/3862\n",
      "2551/3862\n",
      "2552/3862\n",
      "2553/3862\n",
      "2554/3862\n",
      "2555/3862\n",
      "2556/3862\n",
      "2557/3862\n",
      "2558/3862\n",
      "2559/3862\n",
      "2560/3862\n",
      "2561/3862\n",
      "2562/3862\n",
      "2563/3862\n",
      "2564/3862\n",
      "2565/3862\n",
      "2566/3862\n",
      "2567/3862\n",
      "2568/3862\n",
      "2569/3862\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2570/3862\n",
      "2571/3862\n",
      "2572/3862\n",
      "2573/3862\n",
      "2574/3862\n",
      "2575/3862\n",
      "2576/3862\n",
      "2577/3862\n",
      "2578/3862\n",
      "2579/3862\n",
      "2580/3862\n",
      "2581/3862\n",
      "2582/3862\n",
      "2583/3862\n",
      "2584/3862\n",
      "2585/3862\n",
      "2586/3862\n",
      "2587/3862\n",
      "2588/3862\n",
      "2589/3862\n",
      "2590/3862\n",
      "2591/3862\n",
      "2592/3862\n",
      "2593/3862\n",
      "2594/3862\n",
      "2595/3862\n",
      "2596/3862\n",
      "2597/3862\n",
      "2598/3862\n",
      "2599/3862\n",
      "2600/3862\n",
      "2601/3862\n",
      "2602/3862\n",
      "2603/3862\n",
      "2604/3862\n",
      "2605/3862\n",
      "2606/3862\n",
      "2607/3862\n",
      "2608/3862\n",
      "2609/3862\n",
      "2610/3862\n",
      "2611/3862\n",
      "2612/3862\n",
      "2613/3862\n",
      "2614/3862\n",
      "2615/3862\n",
      "2616/3862\n",
      "2617/3862\n",
      "2618/3862\n",
      "2619/3862\n",
      "2620/3862\n",
      "2621/3862\n",
      "2622/3862\n",
      "2623/3862\n",
      "2624/3862\n",
      "2625/3862\n",
      "2626/3862\n",
      "2627/3862\n",
      "2628/3862\n",
      "2629/3862\n",
      "2630/3862\n",
      "2631/3862\n",
      "2632/3862\n",
      "2633/3862\n",
      "2634/3862\n",
      "2635/3862\n",
      "2636/3862\n",
      "2637/3862\n",
      "2638/3862\n",
      "2639/3862\n",
      "2640/3862\n",
      "2641/3862\n",
      "2642/3862\n",
      "2643/3862\n",
      "2644/3862\n",
      "2645/3862\n",
      "2646/3862\n",
      "2647/3862\n",
      "2648/3862\n",
      "2649/3862\n",
      "2650/3862\n",
      "2651/3862\n",
      "2652/3862\n",
      "2653/3862\n",
      "2654/3862\n",
      "2655/3862\n",
      "2656/3862\n",
      "2657/3862\n",
      "2658/3862\n",
      "2659/3862\n",
      "2660/3862\n",
      "2661/3862\n",
      "2662/3862\n",
      "2663/3862\n",
      "2664/3862\n",
      "2665/3862\n",
      "2666/3862\n",
      "2667/3862\n",
      "2668/3862\n",
      "2669/3862\n",
      "2670/3862\n",
      "2671/3862\n",
      "2672/3862\n",
      "2673/3862\n",
      "2674/3862\n",
      "2675/3862\n",
      "2676/3862\n",
      "2677/3862\n",
      "2678/3862\n",
      "2679/3862\n",
      "2680/3862\n",
      "2681/3862\n",
      "2682/3862\n",
      "2683/3862\n",
      "2684/3862\n",
      "2685/3862\n",
      "2686/3862\n",
      "2687/3862\n",
      "2688/3862\n",
      "2689/3862\n",
      "2690/3862\n",
      "2691/3862\n",
      "2692/3862\n",
      "2693/3862\n",
      "2694/3862\n",
      "2695/3862\n",
      "2696/3862\n",
      "2697/3862\n",
      "2698/3862\n",
      "2699/3862\n",
      "2700/3862\n",
      "2701/3862\n",
      "2702/3862\n",
      "2703/3862\n",
      "2704/3862\n",
      "2705/3862\n",
      "2706/3862\n",
      "2707/3862\n",
      "2708/3862\n",
      "2709/3862\n",
      "2710/3862\n",
      "2711/3862\n",
      "2712/3862\n",
      "2713/3862\n",
      "2714/3862\n",
      "2715/3862\n",
      "2716/3862\n",
      "2717/3862\n",
      "2718/3862\n",
      "2719/3862\n",
      "2720/3862\n",
      "2721/3862\n",
      "2722/3862\n",
      "2723/3862\n",
      "2724/3862\n",
      "2725/3862\n",
      "2726/3862\n",
      "2727/3862\n",
      "2728/3862\n",
      "2729/3862\n",
      "2730/3862\n",
      "2731/3862\n",
      "2732/3862\n",
      "2733/3862\n",
      "2734/3862\n",
      "2735/3862\n",
      "2736/3862\n",
      "2737/3862\n",
      "2738/3862\n",
      "2739/3862\n",
      "2740/3862\n",
      "2741/3862\n",
      "2742/3862\n",
      "2743/3862\n",
      "2744/3862\n",
      "2745/3862\n",
      "2746/3862\n",
      "2747/3862\n",
      "2748/3862\n",
      "2749/3862\n",
      "2750/3862\n",
      "2751/3862\n",
      "2752/3862\n",
      "2753/3862\n",
      "2754/3862\n",
      "2755/3862\n",
      "2756/3862\n",
      "2757/3862\n",
      "2758/3862\n",
      "2759/3862\n",
      "2760/3862\n",
      "2761/3862\n",
      "2762/3862\n",
      "2763/3862\n",
      "2764/3862\n",
      "2765/3862\n",
      "2766/3862\n",
      "2767/3862\n",
      "2768/3862\n",
      "2769/3862\n",
      "2770/3862\n",
      "2771/3862\n",
      "2772/3862\n",
      "2773/3862\n",
      "2774/3862\n",
      "2775/3862\n",
      "2776/3862\n",
      "2777/3862\n",
      "2778/3862\n",
      "2779/3862\n",
      "2780/3862\n",
      "2781/3862\n",
      "2782/3862\n",
      "2783/3862\n",
      "2784/3862\n",
      "2785/3862\n",
      "2786/3862\n",
      "2787/3862\n",
      "2788/3862\n",
      "2789/3862\n",
      "2790/3862\n",
      "2791/3862\n",
      "2792/3862\n",
      "2793/3862\n",
      "2794/3862\n",
      "2795/3862\n",
      "2796/3862\n",
      "2797/3862\n",
      "2798/3862\n",
      "2799/3862\n",
      "2800/3862\n",
      "2801/3862\n",
      "2802/3862\n",
      "2803/3862\n",
      "2804/3862\n",
      "2805/3862\n",
      "2806/3862\n",
      "2807/3862\n",
      "2808/3862\n",
      "2809/3862\n",
      "2810/3862\n",
      "2811/3862\n",
      "2812/3862\n",
      "2813/3862\n",
      "2814/3862\n",
      "2815/3862\n",
      "2816/3862\n",
      "2817/3862\n",
      "2818/3862\n",
      "2819/3862\n",
      "2820/3862\n",
      "2821/3862\n",
      "2822/3862\n",
      "2823/3862\n",
      "2824/3862\n",
      "2825/3862\n",
      "2826/3862\n",
      "2827/3862\n",
      "2828/3862\n",
      "2829/3862\n",
      "2830/3862\n",
      "2831/3862\n",
      "2832/3862\n",
      "2833/3862\n",
      "2834/3862\n",
      "2835/3862\n",
      "2836/3862\n",
      "2837/3862\n",
      "2838/3862\n",
      "2839/3862\n",
      "2840/3862\n",
      "2841/3862\n",
      "2842/3862\n",
      "2843/3862\n",
      "2844/3862\n",
      "2845/3862\n",
      "2846/3862\n",
      "2847/3862\n",
      "2848/3862\n",
      "2849/3862\n",
      "2850/3862\n",
      "2851/3862\n",
      "2852/3862\n",
      "2853/3862\n",
      "2854/3862\n",
      "2855/3862\n",
      "2856/3862\n",
      "2857/3862\n",
      "2858/3862\n",
      "2859/3862\n",
      "2860/3862\n",
      "2861/3862\n",
      "2862/3862\n",
      "2863/3862\n",
      "2864/3862\n",
      "2865/3862\n",
      "2866/3862\n",
      "2867/3862\n",
      "2868/3862\n",
      "2869/3862\n",
      "2870/3862\n",
      "2871/3862\n",
      "2872/3862\n",
      "2873/3862\n",
      "2874/3862\n",
      "2875/3862\n",
      "2876/3862\n",
      "2877/3862\n",
      "2878/3862\n",
      "2879/3862\n",
      "2880/3862\n",
      "2881/3862\n",
      "2882/3862\n",
      "2883/3862\n",
      "2884/3862\n",
      "2885/3862\n",
      "2886/3862\n",
      "2887/3862\n",
      "2888/3862\n",
      "2889/3862\n",
      "2890/3862\n",
      "2891/3862\n",
      "2892/3862\n",
      "2893/3862\n",
      "2894/3862\n",
      "2895/3862\n",
      "2896/3862\n",
      "2897/3862\n",
      "2898/3862\n",
      "2899/3862\n",
      "2900/3862\n",
      "2901/3862\n",
      "2902/3862\n",
      "2903/3862\n",
      "2904/3862\n",
      "2905/3862\n",
      "2906/3862\n",
      "2907/3862\n",
      "2908/3862\n",
      "2909/3862\n",
      "2910/3862\n",
      "2911/3862\n",
      "2912/3862\n",
      "2913/3862\n",
      "2914/3862\n",
      "2915/3862\n",
      "2916/3862\n",
      "2917/3862\n",
      "2918/3862\n",
      "2919/3862\n",
      "2920/3862\n",
      "2921/3862\n",
      "2922/3862\n",
      "2923/3862\n",
      "2924/3862\n",
      "2925/3862\n",
      "2926/3862\n",
      "2927/3862\n",
      "2928/3862\n",
      "2929/3862\n",
      "2930/3862\n",
      "2931/3862\n",
      "2932/3862\n",
      "2933/3862\n",
      "2934/3862\n",
      "2935/3862\n",
      "2936/3862\n",
      "2937/3862\n",
      "2938/3862\n",
      "2939/3862\n",
      "2940/3862\n",
      "2941/3862\n",
      "2942/3862\n",
      "2943/3862\n",
      "2944/3862\n",
      "2945/3862\n",
      "2946/3862\n",
      "2947/3862\n",
      "2948/3862\n",
      "2949/3862\n",
      "2950/3862\n",
      "2951/3862\n",
      "2952/3862\n",
      "2953/3862\n",
      "2954/3862\n",
      "2955/3862\n",
      "2956/3862\n",
      "2957/3862\n",
      "2958/3862\n",
      "2959/3862\n",
      "2960/3862\n",
      "2961/3862\n",
      "2962/3862\n",
      "2963/3862\n",
      "2964/3862\n",
      "2965/3862\n",
      "2966/3862\n",
      "2967/3862\n",
      "2968/3862\n",
      "2969/3862\n",
      "2970/3862\n",
      "2971/3862\n",
      "2972/3862\n",
      "2973/3862\n",
      "2974/3862\n",
      "2975/3862\n",
      "2976/3862\n",
      "2977/3862\n",
      "2978/3862\n",
      "2979/3862\n",
      "2980/3862\n",
      "2981/3862\n",
      "2982/3862\n",
      "2983/3862\n",
      "2984/3862\n",
      "2985/3862\n",
      "2986/3862\n",
      "2987/3862\n",
      "2988/3862\n",
      "2989/3862\n",
      "2990/3862\n",
      "2991/3862\n",
      "2992/3862\n",
      "2993/3862\n",
      "2994/3862\n",
      "2995/3862\n",
      "2996/3862\n",
      "2997/3862\n",
      "2998/3862\n",
      "2999/3862\n",
      "3000/3862\n",
      "3001/3862\n",
      "3002/3862\n",
      "3003/3862\n",
      "3004/3862\n",
      "3005/3862\n",
      "3006/3862\n",
      "3007/3862\n",
      "3008/3862\n",
      "3009/3862\n",
      "3010/3862\n",
      "3011/3862\n",
      "3012/3862\n",
      "3013/3862\n",
      "3014/3862\n",
      "3015/3862\n",
      "3016/3862\n",
      "3017/3862\n",
      "3018/3862\n",
      "3019/3862\n",
      "3020/3862\n",
      "3021/3862\n",
      "3022/3862\n",
      "3023/3862\n",
      "3024/3862\n",
      "3025/3862\n",
      "3026/3862\n",
      "3027/3862\n",
      "3028/3862\n",
      "3029/3862\n",
      "3030/3862\n",
      "3031/3862\n",
      "3032/3862\n",
      "3033/3862\n",
      "3034/3862\n",
      "3035/3862\n",
      "3036/3862\n",
      "3037/3862\n",
      "3038/3862\n",
      "3039/3862\n",
      "3040/3862\n",
      "3041/3862\n",
      "3042/3862\n",
      "3043/3862\n",
      "3044/3862\n",
      "3045/3862\n",
      "3046/3862\n",
      "3047/3862\n",
      "3048/3862\n",
      "3049/3862\n",
      "3050/3862\n",
      "3051/3862\n",
      "3052/3862\n",
      "3053/3862\n",
      "3054/3862\n",
      "3055/3862\n",
      "3056/3862\n",
      "3057/3862\n",
      "3058/3862\n",
      "3059/3862\n",
      "3060/3862\n",
      "3061/3862\n",
      "3062/3862\n",
      "3063/3862\n",
      "3064/3862\n",
      "3065/3862\n",
      "3066/3862\n",
      "3067/3862\n",
      "3068/3862\n",
      "3069/3862\n",
      "3070/3862\n",
      "3071/3862\n",
      "3072/3862\n",
      "3073/3862\n",
      "3074/3862\n",
      "3075/3862\n",
      "3076/3862\n",
      "3077/3862\n",
      "3078/3862\n",
      "3079/3862\n",
      "3080/3862\n",
      "3081/3862\n",
      "3082/3862\n",
      "3083/3862\n",
      "3084/3862\n",
      "3085/3862\n",
      "3086/3862\n",
      "3087/3862\n",
      "3088/3862\n",
      "3089/3862\n",
      "3090/3862\n",
      "3091/3862\n",
      "3092/3862\n",
      "3093/3862\n",
      "3094/3862\n",
      "3095/3862\n",
      "3096/3862\n",
      "3097/3862\n",
      "3098/3862\n",
      "3099/3862\n",
      "3100/3862\n",
      "3101/3862\n",
      "3102/3862\n",
      "3103/3862\n",
      "3104/3862\n",
      "3105/3862\n",
      "3106/3862\n",
      "3107/3862\n",
      "3108/3862\n",
      "3109/3862\n",
      "3110/3862\n",
      "3111/3862\n",
      "3112/3862\n",
      "3113/3862\n",
      "3114/3862\n",
      "3115/3862\n",
      "3116/3862\n",
      "3117/3862\n",
      "3118/3862\n",
      "3119/3862\n",
      "3120/3862\n",
      "3121/3862\n",
      "3122/3862\n",
      "3123/3862\n",
      "3124/3862\n",
      "3125/3862\n",
      "3126/3862\n",
      "3127/3862\n",
      "3128/3862\n",
      "3129/3862\n",
      "3130/3862\n",
      "3131/3862\n",
      "3132/3862\n",
      "3133/3862\n",
      "3134/3862\n",
      "3135/3862\n",
      "3136/3862\n",
      "3137/3862\n",
      "3138/3862\n",
      "3139/3862\n",
      "3140/3862\n",
      "3141/3862\n",
      "3142/3862\n",
      "3143/3862\n",
      "3144/3862\n",
      "3145/3862\n",
      "3146/3862\n",
      "3147/3862\n",
      "3148/3862\n",
      "3149/3862\n",
      "3150/3862\n",
      "3151/3862\n",
      "3152/3862\n",
      "3153/3862\n",
      "3154/3862\n",
      "3155/3862\n",
      "3156/3862\n",
      "3157/3862\n",
      "3158/3862\n",
      "3159/3862\n",
      "3160/3862\n",
      "3161/3862\n",
      "3162/3862\n",
      "3163/3862\n",
      "3164/3862\n",
      "3165/3862\n",
      "3166/3862\n",
      "3167/3862\n",
      "3168/3862\n",
      "3169/3862\n",
      "3170/3862\n",
      "3171/3862\n",
      "3172/3862\n",
      "3173/3862\n",
      "3174/3862\n",
      "3175/3862\n",
      "3176/3862\n",
      "3177/3862\n",
      "3178/3862\n",
      "3179/3862\n",
      "3180/3862\n",
      "3181/3862\n",
      "3182/3862\n",
      "3183/3862\n",
      "3184/3862\n",
      "3185/3862\n",
      "3186/3862\n",
      "3187/3862\n",
      "3188/3862\n",
      "3189/3862\n",
      "3190/3862\n",
      "3191/3862\n",
      "3192/3862\n",
      "3193/3862\n",
      "3194/3862\n",
      "3195/3862\n",
      "3196/3862\n",
      "3197/3862\n",
      "3198/3862\n",
      "3199/3862\n",
      "3200/3862\n",
      "3201/3862\n",
      "3202/3862\n",
      "3203/3862\n",
      "3204/3862\n",
      "3205/3862\n",
      "3206/3862\n",
      "3207/3862\n",
      "3208/3862\n",
      "3209/3862\n",
      "3210/3862\n",
      "3211/3862\n",
      "3212/3862\n",
      "3213/3862\n",
      "3214/3862\n",
      "3215/3862\n",
      "3216/3862\n",
      "3217/3862\n",
      "3218/3862\n",
      "3219/3862\n",
      "3220/3862\n",
      "3221/3862\n",
      "3222/3862\n",
      "3223/3862\n",
      "3224/3862\n",
      "3225/3862\n",
      "3226/3862\n",
      "3227/3862\n",
      "3228/3862\n",
      "3229/3862\n",
      "3230/3862\n",
      "3231/3862\n",
      "3232/3862\n",
      "3233/3862\n",
      "3234/3862\n",
      "3235/3862\n",
      "3236/3862\n",
      "3237/3862\n",
      "3238/3862\n",
      "3239/3862\n",
      "3240/3862\n",
      "3241/3862\n",
      "3242/3862\n",
      "3243/3862\n",
      "3244/3862\n",
      "3245/3862\n",
      "3246/3862\n",
      "3247/3862\n",
      "3248/3862\n",
      "3249/3862\n",
      "3250/3862\n",
      "3251/3862\n",
      "3252/3862\n",
      "3253/3862\n",
      "3254/3862\n",
      "3255/3862\n",
      "3256/3862\n",
      "3257/3862\n",
      "3258/3862\n",
      "3259/3862\n",
      "3260/3862\n",
      "3261/3862\n",
      "3262/3862\n",
      "3263/3862\n",
      "3264/3862\n",
      "3265/3862\n",
      "3266/3862\n",
      "3267/3862\n",
      "3268/3862\n",
      "3269/3862\n",
      "3270/3862\n",
      "3271/3862\n",
      "3272/3862\n",
      "3273/3862\n",
      "3274/3862\n",
      "3275/3862\n",
      "3276/3862\n",
      "3277/3862\n",
      "3278/3862\n",
      "3279/3862\n",
      "3280/3862\n",
      "3281/3862\n",
      "3282/3862\n",
      "3283/3862\n",
      "3284/3862\n",
      "3285/3862\n",
      "3286/3862\n",
      "3287/3862\n",
      "3288/3862\n",
      "3289/3862\n",
      "3290/3862\n",
      "3291/3862\n",
      "3292/3862\n",
      "3293/3862\n",
      "3294/3862\n",
      "3295/3862\n",
      "3296/3862\n",
      "3297/3862\n",
      "3298/3862\n",
      "3299/3862\n",
      "3300/3862\n",
      "3301/3862\n",
      "3302/3862\n",
      "3303/3862\n",
      "3304/3862\n",
      "3305/3862\n",
      "3306/3862\n",
      "3307/3862\n",
      "3308/3862\n",
      "3309/3862\n",
      "3310/3862\n",
      "3311/3862\n",
      "3312/3862\n",
      "3313/3862\n",
      "3314/3862\n",
      "3315/3862\n",
      "3316/3862\n",
      "3317/3862\n",
      "3318/3862\n",
      "3319/3862\n",
      "3320/3862\n",
      "3321/3862\n",
      "3322/3862\n",
      "3323/3862\n",
      "3324/3862\n",
      "3325/3862\n",
      "3326/3862\n",
      "3327/3862\n",
      "3328/3862\n",
      "3329/3862\n",
      "3330/3862\n",
      "3331/3862\n",
      "3332/3862\n",
      "3333/3862\n",
      "3334/3862\n",
      "3335/3862\n",
      "3336/3862\n",
      "3337/3862\n",
      "3338/3862\n",
      "3339/3862\n",
      "3340/3862\n",
      "3341/3862\n",
      "3342/3862\n",
      "3343/3862\n",
      "3344/3862\n",
      "3345/3862\n",
      "3346/3862\n",
      "3347/3862\n",
      "3348/3862\n",
      "3349/3862\n",
      "3350/3862\n",
      "3351/3862\n",
      "3352/3862\n",
      "3353/3862\n",
      "3354/3862\n",
      "3355/3862\n",
      "3356/3862\n",
      "3357/3862\n",
      "3358/3862\n",
      "3359/3862\n",
      "3360/3862\n",
      "3361/3862\n",
      "3362/3862\n",
      "3363/3862\n",
      "3364/3862\n",
      "3365/3862\n",
      "3366/3862\n",
      "3367/3862\n",
      "3368/3862\n",
      "3369/3862\n",
      "3370/3862\n",
      "3371/3862\n",
      "3372/3862\n",
      "3373/3862\n",
      "3374/3862\n",
      "3375/3862\n",
      "3376/3862\n",
      "3377/3862\n",
      "3378/3862\n",
      "3379/3862\n",
      "3380/3862\n",
      "3381/3862\n",
      "3382/3862\n",
      "3383/3862\n",
      "3384/3862\n",
      "3385/3862\n",
      "3386/3862\n",
      "3387/3862\n",
      "3388/3862\n",
      "3389/3862\n",
      "3390/3862\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3391/3862\n",
      "3392/3862\n",
      "3393/3862\n",
      "3394/3862\n",
      "3395/3862\n",
      "3396/3862\n",
      "3397/3862\n",
      "3398/3862\n",
      "3399/3862\n",
      "3400/3862\n",
      "3401/3862\n",
      "3402/3862\n",
      "3403/3862\n",
      "3404/3862\n",
      "3405/3862\n",
      "3406/3862\n",
      "3407/3862\n",
      "3408/3862\n",
      "3409/3862\n",
      "3410/3862\n",
      "3411/3862\n",
      "3412/3862\n",
      "3413/3862\n",
      "3414/3862\n",
      "3415/3862\n",
      "3416/3862\n",
      "3417/3862\n",
      "3418/3862\n",
      "3419/3862\n",
      "3420/3862\n",
      "3421/3862\n",
      "3422/3862\n",
      "3423/3862\n",
      "3424/3862\n",
      "3425/3862\n",
      "3426/3862\n",
      "3427/3862\n",
      "3428/3862\n",
      "3429/3862\n",
      "3430/3862\n",
      "3431/3862\n",
      "3432/3862\n",
      "3433/3862\n",
      "3434/3862\n",
      "3435/3862\n",
      "3436/3862\n",
      "3437/3862\n",
      "3438/3862\n",
      "3439/3862\n",
      "3440/3862\n",
      "3441/3862\n",
      "3442/3862\n",
      "3443/3862\n",
      "3444/3862\n",
      "3445/3862\n",
      "3446/3862\n",
      "3447/3862\n",
      "3448/3862\n",
      "3449/3862\n",
      "3450/3862\n",
      "3451/3862\n",
      "3452/3862\n",
      "3453/3862\n",
      "3454/3862\n",
      "3455/3862\n",
      "3456/3862\n",
      "3457/3862\n",
      "3458/3862\n",
      "3459/3862\n",
      "3460/3862\n",
      "3461/3862\n",
      "3462/3862\n",
      "3463/3862\n",
      "3464/3862\n",
      "3465/3862\n",
      "3466/3862\n",
      "3467/3862\n",
      "3468/3862\n",
      "3469/3862\n",
      "3470/3862\n",
      "3471/3862\n",
      "3472/3862\n",
      "3473/3862\n",
      "3474/3862\n",
      "3475/3862\n",
      "3476/3862\n",
      "3477/3862\n",
      "3478/3862\n",
      "3479/3862\n",
      "3480/3862\n",
      "3481/3862\n",
      "3482/3862\n",
      "3483/3862\n",
      "3484/3862\n",
      "3485/3862\n",
      "3486/3862\n",
      "3487/3862\n",
      "3488/3862\n",
      "3489/3862\n",
      "3490/3862\n",
      "3491/3862\n",
      "3492/3862\n",
      "3493/3862\n",
      "3494/3862\n",
      "3495/3862\n",
      "3496/3862\n",
      "3497/3862\n",
      "3498/3862\n",
      "3499/3862\n",
      "3500/3862\n",
      "3501/3862\n",
      "3502/3862\n",
      "3503/3862\n",
      "3504/3862\n",
      "3505/3862\n",
      "3506/3862\n",
      "3507/3862\n",
      "3508/3862\n",
      "3509/3862\n",
      "3510/3862\n",
      "3511/3862\n",
      "3512/3862\n",
      "3513/3862\n",
      "3514/3862\n",
      "3515/3862\n",
      "3516/3862\n",
      "3517/3862\n",
      "3518/3862\n",
      "3519/3862\n",
      "3520/3862\n",
      "3521/3862\n",
      "3522/3862\n",
      "3523/3862\n",
      "3524/3862\n",
      "3525/3862\n",
      "3526/3862\n",
      "3527/3862\n",
      "3528/3862\n",
      "3529/3862\n",
      "3530/3862\n",
      "3531/3862\n",
      "3532/3862\n",
      "3533/3862\n",
      "3534/3862\n",
      "3535/3862\n",
      "3536/3862\n",
      "3537/3862\n",
      "3538/3862\n",
      "3539/3862\n",
      "3540/3862\n",
      "3541/3862\n",
      "3542/3862\n",
      "3543/3862\n",
      "3544/3862\n",
      "3545/3862\n",
      "3546/3862\n",
      "3547/3862\n",
      "3548/3862\n",
      "3549/3862\n",
      "3550/3862\n",
      "3551/3862\n",
      "3552/3862\n",
      "3553/3862\n",
      "3554/3862\n",
      "3555/3862\n",
      "3556/3862\n",
      "3557/3862\n",
      "3558/3862\n",
      "3559/3862\n",
      "3560/3862\n",
      "3561/3862\n",
      "3562/3862\n",
      "3563/3862\n",
      "3564/3862\n",
      "3565/3862\n",
      "3566/3862\n",
      "3567/3862\n",
      "3568/3862\n",
      "3569/3862\n",
      "3570/3862\n",
      "3571/3862\n",
      "3572/3862\n",
      "3573/3862\n",
      "3574/3862\n",
      "3575/3862\n",
      "3576/3862\n",
      "3577/3862\n",
      "3578/3862\n",
      "3579/3862\n",
      "3580/3862\n",
      "3581/3862\n",
      "3582/3862\n",
      "3583/3862\n",
      "3584/3862\n",
      "3585/3862\n",
      "3586/3862\n",
      "3587/3862\n",
      "3588/3862\n",
      "3589/3862\n",
      "3590/3862\n",
      "3591/3862\n",
      "3592/3862\n",
      "3593/3862\n",
      "3594/3862\n",
      "3595/3862\n",
      "3596/3862\n",
      "3597/3862\n",
      "3598/3862\n",
      "3599/3862\n",
      "3600/3862\n",
      "3601/3862\n",
      "3602/3862\n",
      "3603/3862\n",
      "3604/3862\n",
      "3605/3862\n",
      "3606/3862\n",
      "3607/3862\n",
      "3608/3862\n",
      "3609/3862\n",
      "3610/3862\n",
      "3611/3862\n",
      "3612/3862\n",
      "3613/3862\n",
      "3614/3862\n",
      "3615/3862\n",
      "3616/3862\n",
      "3617/3862\n",
      "3618/3862\n",
      "3619/3862\n",
      "3620/3862\n",
      "3621/3862\n",
      "3622/3862\n",
      "3623/3862\n",
      "3624/3862\n",
      "3625/3862\n",
      "3626/3862\n",
      "3627/3862\n",
      "3628/3862\n",
      "3629/3862\n",
      "3630/3862\n",
      "3631/3862\n",
      "3632/3862\n",
      "3633/3862\n",
      "3634/3862\n",
      "3635/3862\n",
      "3636/3862\n",
      "3637/3862\n",
      "3638/3862\n",
      "3639/3862\n",
      "3640/3862\n",
      "3641/3862\n",
      "3642/3862\n",
      "3643/3862\n",
      "3644/3862\n",
      "3645/3862\n",
      "3646/3862\n",
      "3647/3862\n",
      "3648/3862\n",
      "3649/3862\n",
      "3650/3862\n",
      "3651/3862\n",
      "3652/3862\n",
      "3653/3862\n",
      "3654/3862\n",
      "3655/3862\n",
      "3656/3862\n",
      "3657/3862\n",
      "3658/3862\n",
      "3659/3862\n",
      "3660/3862\n",
      "3661/3862\n",
      "3662/3862\n",
      "3663/3862\n",
      "3664/3862\n",
      "3665/3862\n",
      "3666/3862\n",
      "3667/3862\n",
      "3668/3862\n",
      "3669/3862\n",
      "3670/3862\n",
      "3671/3862\n",
      "3672/3862\n",
      "3673/3862\n",
      "3674/3862\n",
      "3675/3862\n",
      "3676/3862\n",
      "3677/3862\n",
      "3678/3862\n",
      "3679/3862\n",
      "3680/3862\n",
      "3681/3862\n",
      "3682/3862\n",
      "3683/3862\n",
      "3684/3862\n",
      "3685/3862\n",
      "3686/3862\n",
      "3687/3862\n",
      "3688/3862\n",
      "3689/3862\n",
      "3690/3862\n",
      "3691/3862\n",
      "3692/3862\n",
      "3693/3862\n",
      "3694/3862\n",
      "3695/3862\n",
      "3696/3862\n",
      "3697/3862\n",
      "3698/3862\n",
      "3699/3862\n",
      "3700/3862\n",
      "3701/3862\n",
      "3702/3862\n",
      "3703/3862\n",
      "3704/3862\n",
      "3705/3862\n",
      "3706/3862\n",
      "3707/3862\n",
      "3708/3862\n",
      "3709/3862\n",
      "3710/3862\n",
      "3711/3862\n",
      "3712/3862\n",
      "3713/3862\n",
      "3714/3862\n",
      "3715/3862\n",
      "3716/3862\n",
      "3717/3862\n",
      "3718/3862\n",
      "3719/3862\n",
      "3720/3862\n",
      "3721/3862\n",
      "3722/3862\n",
      "3723/3862\n",
      "3724/3862\n",
      "3725/3862\n",
      "3726/3862\n",
      "3727/3862\n",
      "3728/3862\n",
      "3729/3862\n",
      "3730/3862\n",
      "3731/3862\n",
      "3732/3862\n",
      "3733/3862\n",
      "3734/3862\n",
      "3735/3862\n",
      "3736/3862\n",
      "3737/3862\n",
      "3738/3862\n",
      "3739/3862\n",
      "3740/3862\n",
      "3741/3862\n",
      "3742/3862\n",
      "3743/3862\n",
      "3744/3862\n",
      "3745/3862\n",
      "3746/3862\n",
      "3747/3862\n",
      "3748/3862\n",
      "3749/3862\n",
      "3750/3862\n",
      "3751/3862\n",
      "3752/3862\n",
      "3753/3862\n",
      "3754/3862\n",
      "3755/3862\n",
      "3756/3862\n",
      "3757/3862\n",
      "3758/3862\n",
      "3759/3862\n",
      "3760/3862\n",
      "3761/3862\n",
      "3762/3862\n",
      "3763/3862\n",
      "3764/3862\n",
      "3765/3862\n",
      "3766/3862\n",
      "3767/3862\n",
      "3768/3862\n",
      "3769/3862\n",
      "3770/3862\n",
      "3771/3862\n",
      "3772/3862\n",
      "3773/3862\n",
      "3774/3862\n",
      "3775/3862\n",
      "3776/3862\n",
      "3777/3862\n",
      "3778/3862\n",
      "3779/3862\n",
      "3780/3862\n",
      "3781/3862\n",
      "3782/3862\n",
      "3783/3862\n",
      "3784/3862\n",
      "3785/3862\n",
      "3786/3862\n",
      "3787/3862\n",
      "3788/3862\n",
      "3789/3862\n",
      "3790/3862\n",
      "3791/3862\n",
      "3792/3862\n",
      "3793/3862\n",
      "3794/3862\n",
      "3795/3862\n",
      "3796/3862\n",
      "3797/3862\n",
      "3798/3862\n",
      "3799/3862\n",
      "3800/3862\n",
      "3801/3862\n",
      "3802/3862\n",
      "3803/3862\n",
      "3804/3862\n",
      "3805/3862\n",
      "3806/3862\n",
      "3807/3862\n",
      "3808/3862\n",
      "3809/3862\n",
      "3810/3862\n",
      "3811/3862\n",
      "3812/3862\n",
      "3813/3862\n",
      "3814/3862\n",
      "3815/3862\n",
      "3816/3862\n",
      "3817/3862\n",
      "3818/3862\n",
      "3819/3862\n",
      "3820/3862\n",
      "3821/3862\n",
      "3822/3862\n",
      "3823/3862\n",
      "3824/3862\n",
      "3825/3862\n",
      "3826/3862\n",
      "3827/3862\n",
      "3828/3862\n",
      "3829/3862\n",
      "3830/3862\n",
      "3831/3862\n",
      "3832/3862\n",
      "3833/3862\n",
      "3834/3862\n",
      "3835/3862\n",
      "3836/3862\n",
      "3837/3862\n",
      "3838/3862\n",
      "3839/3862\n",
      "3840/3862\n",
      "3841/3862\n",
      "3842/3862\n",
      "3843/3862\n",
      "3844/3862\n",
      "3845/3862\n",
      "3846/3862\n",
      "3847/3862\n",
      "3848/3862\n",
      "3849/3862\n",
      "3850/3862\n",
      "3851/3862\n",
      "3852/3862\n",
      "3853/3862\n",
      "3854/3862\n",
      "3855/3862\n",
      "3856/3862\n",
      "3857/3862\n",
      "3858/3862\n",
      "3859/3862\n",
      "3860/3862\n",
      "3861/3862\n",
      "3862/3862\n"
     ]
    }
   ],
   "source": [
    "from osgeo import osr\n",
    "import wradlib as wrl\n",
    "import os\n",
    "\n",
    "#Data Paths\n",
    "wdir = \"/media/elisabeth/0A581C5C581C4939/master/radolan_data/out/bin/2009-2019/\"\n",
    "out2 = \"out/\"\n",
    "os.chdir(wdir)\n",
    "files = os.listdir()\n",
    "files = files[1:]\n",
    "\n",
    "proj_stereo = wrl.georef.create_osr(\"dwd-radolan\") #(osr.SpatialReference) – GDAL/OSR object defining projection\n",
    "proj_wgs = osr.SpatialReference() #Initialize a SpatialReference object\n",
    "proj_wgs.ImportFromEPSG(4326) # Spatial Referenz for EPSG 3426\n",
    "#print(proj_wgs)\n",
    "\n",
    "radolan_grid_xy = wrl.georef.get_radolan_grid(900,900) \n",
    "radolan_grid_ll = wrl.georef.projection.reproject(radolan_grid_xy, projection_source=proj_stereo, projection_target=proj_wgs)\n",
    "#print(\"{0}, ({1:.4f}, {2:.4f})\".format(radolan_grid_ll.shape, *radolan_grid_ll[0,0,:]))\n",
    "\n",
    "# create Gauss Krueger zone 3 projection osr object\n",
    "proj_gk3 = osr.SpatialReference()  #Initialize a SpatialReference object\n",
    "proj_gk3.ImportFromEPSG(31467) #Spatial Referenz for EPSG 31467\n",
    "#print(proj_gk3)\n",
    "\n",
    "# transform radolan polar stereographic projection to wgs84 and then to gk3\n",
    "radolan_grid_gk = wrl.georef.reproject(radolan_grid_ll,projection_source=proj_wgs,projection_target=proj_gk3)\n",
    "#print(\"{0}, ({1:.4f}, {2:.4f})\".format(radolan_grid_gk.shape, *radolan_grid_gk[0,0,:]))\n",
    "\n",
    "amount = len(files)\n",
    "print(amount)\n",
    "n=0\n",
    "\n",
    "for file in files:\n",
    "    n=n+1\n",
    "    data_raw, attributes = wrl.io.radolan.read_radolan_composite(wdir + file, missing=-9999, loaddata=True) #data : numpy.array of shape (number of rows, number of columns) | xarray.Dataset >// attrs : dictionary of metadata information from the file heade\n",
    "    data, xy = wrl.georef.set_raster_origin(data_raw, radolan_grid_xy, 'upper')\n",
    "\n",
    "    rad_stereo = wrl.georef.create_raster_dataset(data, xy, projection=proj_stereo, nodata=-9999) #dataset (gdal.Dataset) – In-Memory raster dataset\n",
    "\n",
    "    cols = rad_stereo.RasterXSize\n",
    "    rows = rad_stereo.RasterYSize\n",
    "\n",
    "    filename = file[15:25] + \".tif\"\n",
    "\n",
    "    wrl.io.write_raster_dataset(wdir + out2 + \"/\" + filename, rad_stereo, \"GTiff\")\n",
    "    print(str(n)+\"/\"+ str(amount))\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Import Modules\n",
    "import wradlib as wrl\n",
    "import numpy as np\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')\n",
    "\n",
    "#Working dir\n",
    "wdir = \"/home/elisabeth/Dokumente/master/radolan_data/raw/\"\n",
    "file = \"raa01-sf_10000-1701012350-dwd---bin\"\n",
    "#Read the original data: data = 2d Array(row,col) and metadata (productType, datetime,nodataflag)\n",
    "#Origin in lower left corner\n",
    "data_raw, meta  = wrl.io.read_radolan_composite(wdir + file)\n",
    "\n",
    "#Round data to 2 decimal places to save memory\n",
    "data_raw_s = np.around(data_raw, decimals = 2, out = None)\n",
    "\n",
    "# Create Spatail Reference (Radolan projection)\n",
    "proj_osr = wrl.georef.create_osr(\"dwd-radolan\")\n",
    "\n",
    "#Load EPSG 3857 \n",
    "epsg_3857 = wrl.georef.epsg_to_osr(3857)\n",
    "\n",
    "# Get projected RADOLAN coordinates for corner definition\n",
    "#origin in lower left corner\n",
    "xy_raw = wrl.georef.get_radolan_grid(900, 900)\n",
    "\n",
    "#flip the data and the coords up-down (from lower left to upper left corner)\n",
    "data, xy = wrl.georef.set_raster_origin(data_raw, xy_raw, 'upper')\n",
    "\n",
    "#Create raster dataset: Geotiff with radolan projection\n",
    "ds = wrl.georef.create_raster_dataset(data, xy, projection=proj_osr)\n",
    "\n",
    "#Reproject dataset to EPSG 3857\n",
    "ds2 = wrl.georef.reproject_raster_dataset(ds,projection_target=epsg_3857, spacing=1000)\n",
    "#Bug appears it will fixed in future https://groups.google.com/forum/#!topic/wradlib-users/-dvRhDCjgV0\n",
    "#This line is needed until fixed\n",
    "ds2.SetProjection(epsg_3857.ExportToWkt())\n",
    "\n",
    "#write raster dataset\n",
    "filename = file[15:25] + \".tif\"\n",
    "wrl.io.write_raster_dataset(wdir + filename, ds2, \"GTiff\")\n",
    "#"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from osgeo import gdal\n",
    "import numpy as np\n",
    "raster = None\n",
    "array = None\n",
    "raster = gdal.Open(\"/media/elisabeth/0A581C5C581C4939/master/radolan_data/out/bin/2009-2019/out/0901022350.tif\")\n",
    "#band = raster.GetRasterBand(1)\n",
    "array = raster.ReadAsArray()\n",
    "print(np.median(array))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# load radolan files\n",
    "print(os.getcwd())\n",
    "rw_filename = wrl.util.get_wradlib_data_file(\"raa01-sf_10000-0901012350-dwd---bin.gz\")\n",
    "ds, rwattrs = wrl.io.read_radolan_composite(rw_filename, loaddata='xarray')\n",
    "# print the xarray dataset\n",
    "print(ds)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}