How Do You Read and Manipulate Complex Number Data in Fortran?

[Raziel13]

Hello, I need to read a fortran data with complex numbers and real numbers,
the first column is the real numbers, the second and third complex numbers (real, imaginary).

I need to read the first 64 lines and then the next 64 lines in separate ways and save in a variable. for example

read from 1 to 64 the first column that I will call k2 which is a real number and the second and third column that I will call E0 which is a complex number (part real and imaginary),

then read from 65 to 128 again the first column in which I can call time again k2 because the values are the same, but the second and third column will be my E1 because the values are no longer repeated.
93
then do the same for 129 to 193 and the same for 194 to 256.

Can someone help me ?
thanks

   0.3475999158309906D-03   0.1000000000000000D+01  -0.0000000000000000D+00
   0.1833296439467811D-02   0.9999949615756608D+00   0.1829832969435655D-14
   0.4513595263596993D-02   0.9999685090269346D+00   0.8610978679697199D-15
   0.8401876425510740D-02   0.9998904236081211D+00   0.8610978679697199D-15
   0.1351685480299403D-01   0.9997161185401198D+00   0.4305489339848600D-15
   0.1988323578908195D-01   0.9993856054848251D+00  -0.1076372334962150D-15
   0.2753201395850637D-01   0.9988221275112609D+00   0.5381861674810750D-15
   0.3650080488705306D-01   0.9979304502530658D+00   0.4305489339848600D-15
   0.4683424556300787D-01   0.9965948061656795D+00   0.1076372334962150D-14
   0.5858447526034363D-01   0.9946764945490979D+00   0.1184009568458365D-14
   0.7181170689468865D-01   0.9920111543230785D+00   0.3659665938871310D-14
   0.8658490010557232D-01   0.9884057503391426D+00   0.4520763806841030D-14
   0.1029825493929753D+00   0.9836353513267678D+00   0.6565871243269115D-14
   0.1210936032577395D+00   0.9774398339009408D+00   0.8288066979208554D-14
   0.1410185334589130D+00   0.9695207297376155D+00   0.1130190951710258D-13
   0.1628705773018296D+00   0.9595385515367878D+00   0.1350847280377498D-13
   0.1867771804509159D+00   0.9471110994656604D+00   0.1684522704215765D-13
   0.2128816733065060D+00   0.9318134768763994D+00   0.2066634883127328D-13
   0.2413452207319886D+00   0.9131808456567779D+00   0.2437983338689270D-13
   0.2723490931675873D+00   0.8907153367325286D+00   0.2830859240950455D-13
   0.3060973173045934D+00   0.8638989960823151D+00   0.3137625356414667D-13
   0.3428197770124591D+00   0.8322151579545033D+00   0.3299081206658989D-13
   0.3827758507475470D+00   0.7951811047074102D+00   0.3369045408431530D-13
   0.4262586910527104D+00   0.7523951072774169D+00   0.3035369984593263D-13
   0.4736002760595692D+00   0.7036005955085440D+00   0.2357255413567108D-13
   0.5251773935465452D+00   0.6487687219107892D+00   0.1087136058311772D-13
   0.5814187569572651D+00   0.5881971611468716D+00  -0.7426969111238834D-14
   0.6428135023366771D+00   0.5226167217783453D+00  -0.3153770941439099D-13
   0.7099213787340648D+00   0.4532876377370986D+00  -0.6172995341007930D-13
   0.7833850267593251D+00   0.3820549050592221D+00  -0.9439785377618055D-13
   0.8639448467923458D+00   0.3113201577456072D+00  -0.1265006586664267D-12
   0.9524570982449548D+00   0.2438841128848114D+00  -0.1546747045340610D-12
   0.1049916055896533D+01   0.1826305284653794D+00  -0.1718697525850813D-12
   0.1157481294914600D+01   0.1300708966838944D+00  -0.1748028671978531D-12
   0.1276511505634412D+01   0.8784606175054865D-01  -0.1617047613467825D-12
   0.1408606685127873D+01   0.5635473357345396D-01  -0.1363360108771433D-12
   0.1555661162008767D+01   0.3468869480994791D-01  -0.1043240249026596D-12
   0.1719930752205679D+01   0.2095036471103247D-01  -0.7363732236559808D-13
   0.1904118517453612D+01   0.1283919806400845D-01  -0.4880003073634647D-13
   0.2111485255710072D+01   0.8246343265954895D-02  -0.3031165405159817D-13
   0.2345993221933724D+01   0.5616414055055371D-02  -0.1962277221589201D-13
   0.2612495009930844D+01   0.4004079448995915D-02  -0.1401218141990180D-13
   0.2916984570477848D+01   0.2924767359531206D-02  -0.9237671244252310D-14
   0.3266934871581749D+01   0.2155876043675108D-02  -0.6236021986715284D-14
   0.3671758140882298D+01   0.1592209468999256D-02  -0.5399310679459550D-14
   0.4143442314569344D+01   0.1174791948100784D-02  -0.3846244151230472D-14
   0.4697445226110210D+01   0.8645900884451562D-03  -0.2418106189417654D-14
   0.5353973100920618D+01   0.6333910886683653D-03  -0.1268942073013972D-14
   0.6139844387895875D+01   0.4603362803121361D-03  -0.1171290715672191D-14
   0.7091266484425319D+01   0.3306447016578899D-03  -0.5865388309657028D-15
   0.8258074523321996D+01   0.2341796781250851D-03  -0.1013829215889642D-14
   0.9710383029886540D+01   0.1632957168524006D-03   0.4687974675184881D-15
   0.1154935789936476D+02   0.1117266529545774D-03   0.1878067379123717D-15
   0.1392530609899710D+02   0.7451295627757056D-04  -0.2119830696546856D-15
   0.1706936855807103D+02   0.4811281408853923D-04  -0.3162566396041439D-15
   0.2135189727044204D+02   0.2982137989322365D-04   0.2274611776940583D-15
   0.2739665613112830D+02   0.1755172699612436D-04  -0.4665079426238695D-15
   0.3632135308034874D+02   0.9663125194392776D-05  -0.7393432663015580D-15
   0.5029362477052696D+02   0.4864767650024737D-05   0.9990635048538857D-15
   0.7398170762169443D+02   0.2161879609180998D-05   0.1299916766248431D-14
   0.1190210316547486D+03   0.7982433375114770D-06   0.4217067937690845D-15
   0.2215528733967777D+03   0.2169925441877510D-06   0.3852726656010782D-15
   0.5454655223627067D+03   0.3157102436351680D-07   0.3253355408323738D-15
   0.2876870662092502D+04   0.3494143826275012D-09   0.4140718487780305D-15
   0.3475999158309906D-03  -0.8733997719885090D-15  -0.1455056665833896D-08
   0.1833296439467811D-02   0.2416115304760093D-16  -0.7674057388064735D-08
   0.4513595263596993D-02   0.6313988500916080D-15  -0.1889182202030143D-07
   0.8401876425510740D-02   0.2652459378271201D-15  -0.3515653925677366D-07
   0.1351685480299403D-01   0.5974413568903742D-15  -0.5652419730346947D-07
   0.1988323578908195D-01  -0.3659643586645440D-15  -0.8304847156857286D-07
   0.2753201395850637D-01   0.4654433055924298D-15  -0.1147639369810866D-06
   0.3650080488705306D-01  -0.4810980821862747D-15  -0.1516632990964897D-06
   0.4683424556300787D-01   0.7793848409895145D-15  -0.1936677117106545D-06
   0.5858447526034363D-01   0.1815349753580583D-15  -0.2405890593921787D-06
   0.7181170689468865D-01  -0.1375700461128853D-15  -0.2920832320133602D-06
   0.8658490010557232D-01   0.1595081373327718D-14  -0.3475938060222055D-06
   0.1029825493929753D+00   0.5296778222158051D-15  -0.4062864285807373D-06
   0.1210936032577395D+00   0.7863744025492185D-15  -0.4669757665713327D-06
   0.1410185334589130D+00   0.1390888323847228D-14  -0.5280494923715648D-06
   0.1628705773018296D+00   0.8988982791518660D-15  -0.5873977311507574D-06
   0.1867771804509159D+00   0.3379035513801605D-15  -0.6423620783025709D-06
   0.2128816733065060D+00   0.1628514914058735D-14  -0.6897260164978229D-06
   0.2413452207319886D+00   0.2129162532714469D-14  -0.7257778952315798D-06
   0.2723490931675873D+00   0.1912735245181821D-14  -0.7464876736207860D-06
   0.3060973173045934D+00   0.1970019714412867D-14  -0.7478461711433848D-06
   0.3428197770124591D+00   0.1440059935933419D-14  -0.7264150608421358D-06
   0.3827758507475470D+00   0.2694062150719036D-14  -0.6801172368370289D-06
   0.4262586910527104D+00   0.1983053808005684D-14  -0.6092463604097870D-06
   0.4736002760595692D+00   0.2287396274954972D-14  -0.5175747587457052D-06
   0.5251773935465452D+00   0.1723009398021452D-14  -0.4132801411110632D-06
   0.5814187569572651D+00   0.1830203949662034D-14  -0.3092083034950250D-06
   0.6428135023366771D+00   0.2242053437485190D-14  -0.2218128026173135D-06
   0.7099213787340648D+00   0.1665706964889018D-14  -0.1681275892419838D-06
   0.7833850267593251D+00   0.1588318131262786D-14  -0.1605937019942488D-06
   0.8639448467923458D+00   0.1609886964224556D-14  -0.2007131787033556D-06
   0.9524570982449548D+00   0.2385090463370206D-14  -0.2742290087894907D-06
   0.1049916055896533D+01   0.1251399861379684D-14  -0.3518728354072258D-06
   0.1157481294914600D+01   0.3092226943766533D-15  -0.3987447853792095D-06
   0.1276511505634412D+01  -0.3657043568546687D-15  -0.3904456501467158D-06
   0.1408606685127873D+01   0.5260705681662269D-15  -0.3269097207443796D-06
   0.1555661162008767D+01   0.4356172177970416D-15  -0.2325158794195402D-06
   0.1719930752205679D+01   0.1127584428112931D-14  -0.1401791169097913D-06
   0.1904118517453612D+01   0.3333114970094542D-15  -0.7200028112835087D-07
   0.2111485255710072D+01   0.3875241775972240D-16  -0.3205642042665038D-07
   0.2345993221933724D+01   0.3131162184320873D-15  -0.1289824011960974D-07
   0.2612495009930844D+01  -0.7765031918218519D-15  -0.5130539217639537D-08
   0.2916984570477848D+01   0.2022355383834758D-15  -0.2307679181447039D-08
   0.3266934871581749D+01  -0.5420592942883447D-15  -0.1267576387470081D-08
   0.3671758140882298D+01  -0.5826579323652031D-15  -0.8061997052306844D-09
   0.4143442314569344D+01   0.3342887842195372D-15  -0.5427584922350266D-09
   0.4697445226110210D+01   0.2348953258393024D-15  -0.3666879024144684D-09
   0.5353973100920618D+01  -0.5469195170770313D-15  -0.2431829525939626D-09
   0.6139844387895875D+01   0.5377557829284162D-16  -0.1569965510493959D-09
   0.7091266484425319D+01  -0.2910341105947241D-15  -0.9802805615852434D-10
   0.8258074523321996D+01  -0.9860788121527835D-16  -0.5841877824494821D-10
   0.9710383029886540D+01  -0.3003046174034079D-15  -0.3219681333822845D-10
   0.1154935789936476D+02  -0.2483339064857475D-16  -0.1897158684902246D-11
   0.1392530609899710D+02   0.4901567447468588D-15  -0.1172031033380893D-10
   0.1706936855807103D+02  -0.2398416562072209D-15  -0.4818932621194005D-11
   0.2135189727044204D+02  -0.1775856923938541D-15  -0.2138845788398870D-11
   0.2739665613112830D+02  -0.2686764687542991D-15  -0.9889042059520092D-12
   0.3632135308034874D+02  -0.1206920071946245D-15  -0.4194568948885728D-12
   0.5029362477052696D+02  -0.2447518635624472D-15  -0.1541611363581582D-12
   0.7398170762169443D+02   0.3010623066066283D-15  -0.4623266143420919D-13
   0.1190210316547486D+03  -0.3369955857880252D-15  -0.1076229507803411D-13
   0.2215528733967777D+03   0.7619316117631777D-15  -0.2021459215536639D-14
   0.5454655223627067D+03  -0.5864664272748860D-15   0.6148311046782166D-16
   0.2876870662092502D+04  -0.4056929326897661D-15   0.3342501892659557D-16
   0.3475999158309906D-03   0.1848656039410625D-03   0.2028840969744941D-15
   0.1833296439467811D-02   0.1848562032489020D-03   0.1956443367625295D-15
   0.4513595263596993D-02   0.1848068478938808D-03   0.8616628583310893D-15
   0.8401876425510740D-02   0.1846611537137981D-03   0.2429589946995253D-15
   0.1351685480299403D-01   0.1843359200091747D-03  -0.1387905392359032D-15
   0.1988323578908195D-01   0.1837191821660309D-03   0.8022994524551250D-15
   0.2753201395850637D-01   0.1826676188326764D-03   0.3041356254561632D-15
   0.3650080488705306D-01   0.1810032828268234D-03   0.6824557992911874D-15
   0.4683424556300787D-01   0.1785096270908737D-03   0.1183457717407725D-14
   0.5858447526034363D-01   0.1749268103198742D-03   0.2254140730823994D-14
   0.7181170689468865D-01   0.1699462996607832D-03   0.3438754707575916D-14
   0.8658490010557232D-01   0.1632048544058181D-03   0.5300266555180440D-14
   0.1029825493929753D+00   0.1542780972306471D-03   0.7308413110526391D-14
   0.1210936032577395D+00   0.1426740947165766D-03   0.9834642687937251D-14
   0.1410185334589130D+00   0.1278277305276408D-03   0.1385933181818531D-13
   0.1628705773018296D+00   0.1090972425316932D-03   0.1790127943307660D-13
   0.1867771804509159D+00   0.8576521725961464D-04   0.2319975247234246D-13
   0.2128816733065060D+00   0.5704772630100787D-04   0.2932827847315187D-13
   0.2413452207319886D+00   0.2211728831791783D-04   0.3754280801503499D-13
   0.2723490931675873D+00  -0.1985197093365777D-04   0.4675922313935462D-13
   0.3060973173045934D+00  -0.6960532203289708D-04   0.5794068787982738D-13
   0.3428197770124591D+00  -0.1276796825900629D-03   0.7037730170187581D-13
   0.3827758507475470D+00  -0.1942103614243251D-03   0.8322096480066935D-13
   0.4262586910527104D+00  -0.2686580848945300D-03   0.9945717165120308D-13
   0.4736002760595692D+00  -0.3494540529878669D-03   0.1165993733944693D-12
   0.5251773935465452D+00  -0.4335893355066667D-03   0.1328448434590674D-12
   0.5814187569572651D+00  -0.5162276627590251D-03   0.1482343925441231D-12
   0.6428135023366771D+00  -0.5904998653809707D-03   0.1633163240863794D-12
   0.7099213787340648D+00  -0.6477289755677672D-03   0.1753175602777499D-12
   0.7833850267593251D+00  -0.6783879846336999D-03   0.1803122327439831D-12
   0.8639448467923458D+00  -0.6740118634383764D-03   0.1805934665508378D-12
   0.9524570982449548D+00  -0.6299527648066868D-03   0.1712384349831494D-12
   0.1049916055896533D+01  -0.5482406528706048D-03   0.1540266737288811D-12
   0.1157481294914600D+01  -0.4391112160383313D-03   0.1293711250456402D-12
   0.1276511505634412D+01  -0.3195992939993111D-03   0.1002451361156523D-12
   0.1408606685127873D+01  -0.2086956646787709D-03   0.7079785162039551D-13
   0.1555661162008767D+01  -0.1208668528681256D-03   0.4577660964093063D-13
   0.1719930752205679D+01  -0.6157351244253159D-04   0.2619904322360354D-13
   0.1904118517453612D+01  -0.2752100904109864D-04   0.1336022688805535D-13
   0.2111485255710072D+01  -0.1084855352911824D-04   0.8074977743831828D-14
   0.2345993221933724D+01  -0.3758654643861246D-05   0.2632538098109554D-14
   0.2612495009930844D+01  -0.8701264203233716D-06   0.1384795276449009D-14
   0.2916984570477848D+01   0.7830013643097617D-06   0.7572418817393842D-15
   0.3266934871581749D+01   0.2637989678420009D-05   0.5378260682463941D-15
   0.3671758140882298D+01   0.5135337628825519D-05  -0.8818160973694087D-16
   0.4143442314569344D+01   0.7883682020575292D-05   0.3273010514792248D-15
   0.4697445226110210D+01   0.1013876606363818D-04   0.5272977849209071D-15
   0.5353973100920618D+01   0.1139565774038222D-04  -0.6779810427808658D-15
   0.6139844387895875D+01   0.1136919396889468D-04  -0.9653369472285132D-16
   0.7091266484425319D+01   0.9837991686548136D-05  -0.1410537855090627D-15
   0.8258074523321996D+01   0.7130470326935327D-05  -0.9197353269278784D-16
   0.9710383029886540D+01   0.4285450977624413D-05  -0.8261756330679615D-15
   0.1154935789936476D+02   0.2246177877841349D-05  -0.1258099267437897D-15
   0.1392530609899710D+02   0.1031545452112179D-05   0.1769071664079108D-15
   0.1706936855807103D+02   0.4378818457430551D-06  -0.1343729002331806D-14
   0.2135189727044204D+02   0.1700480665169558D-06  -0.2906121695268202D-15
   0.2739665613112830D+02   0.6219438509173145D-07   0.3363119497171824D-16
   0.3632135308034874D+02   0.2605711365528449D-07   0.5381881050161532D-15
   0.5029362477052696D+02   0.9682208178426647D-08  -0.4781738140635490D-15
   0.7398170762169443D+02   0.2788834406768809D-08  -0.1045431857546606D-14
   0.1190210316547486D+03   0.1027600894319595D-08   0.4626019379675053D-15
   0.2215528733967777D+03   0.3075858809608018D-09   0.4030756312836121D-16
   0.5454655223627067D+03   0.5439778738986877D-10   0.1889392503974009D-15
   0.2876870662092502D+04   0.9608248579180385D-12  -0.7279525934514845D-15
   0.3475999158309906D-03   0.3711665442749685D-04  -0.6841116809258749D-14
   0.1833296439467811D-02   0.2657461358558079D-01  -0.5312134611318257D-11
   0.4513595263596993D-02   0.2919882139978147D+00  -0.5742989975052226D-10
   0.8401876425510740D-02   0.4604174422851450D+00  -0.9182274060201710D-10
   0.1351685480299403D-01   0.5130959599539389D+00  -0.1021401999815633D-09
   0.1988323578908195D-01   0.5245352482643203D+00  -0.1039035131406983D-09
   0.2753201395850637D-01   0.5266736940021662D+00  -0.1043129651769179D-09
   0.3650080488705306D-01   0.5266923025744201D+00  -0.1049570125635425D-09
   0.4683424556300787D-01   0.5259373636630101D+00  -0.1046986293845348D-09
   0.5858447526034363D-01   0.5246365566402391D+00  -0.1051149702036982D-09
   0.7181170689468865D-01   0.5227607562534307D+00  -0.1057278027926089D-09
   0.8658490010557232D-01   0.5202001078776072D+00  -0.1065177524492376D-09
   0.1029825493929753D+00   0.5168050174542463D+00  -0.1069862973266466D-09
   0.1210936032577395D+00   0.5123961761748981D+00  -0.1091742931905410D-09
   0.1410185334589130D+00   0.5067667073985429D+00  -0.1096434300727342D-09
   0.1628705773018296D+00   0.4996825609151804D+00  -0.1111616532511983D-09
   0.1867771804509159D+00   0.4908833146653196D+00  -0.1135675068756890D-09
   0.2128816733065060D+00   0.4800846467916449D+00  -0.1151234838138019D-09
   0.2413452207319886D+00   0.4669837008363743D+00  -0.1196067629540778D-09
   0.2723490931675873D+00   0.4512688225276932D+00  -0.1224104437935705D-09
   0.3060973173045934D+00   0.4326355095482476D+00  -0.1272717179885768D-09
   0.3428197770124591D+00   0.4108107437213624D+00  -0.1318178841825229D-09
   0.3827758507475470D+00   0.3855879452074180D+00  -0.1375479253983856D-09
   0.4262586910527104D+00   0.3568742131955406D+00  -0.1443716416344949D-09
   0.4736002760595692D+00   0.3247497005494051D+00  -0.1531417888732246D-09
   0.5251773935465452D+00   0.2895351810455067D+00  -0.1614128491695407D-09
   0.5814187569572651D+00   0.2518575961321720D+00  -0.1688927488903513D-09
   0.6428135023366771D+00   0.2126950854520832D+00  -0.1765449892766436D-09
   0.7099213787340648D+00   0.1733752076935676D+00  -0.1844437862009131D-09
   0.7833850267593251D+00   0.1354978447851120D+00  -0.1939325330651177D-09
   0.8639448467923458D+00   0.1007642207693468D+00  -0.1984606363858178D-09
   0.9524570982449548D+00   0.7071997015524834D-01  -0.2016128196780628D-09
   0.1049916055896533D+01   0.4646041215981891D-01  -0.2026235904828726D-09
   0.1157481294914600D+01   0.2838615694883695D-01  -0.1992640524694559D-09
   0.1276511505634412D+01   0.1611230685311180D-01  -0.1970765544277665D-09
   0.1408606685127873D+01   0.8598484894333474D-02  -0.1911255819037417D-09
   0.1555661162008767D+01   0.4473672244359425D-02  -0.1866136409863838D-09
   0.1719930752205679D+01   0.2421494775253672D-02  -0.1815605730084602D-09
   0.1904118517453612D+01   0.1452906340458963D-02  -0.1758845226919879D-09
   0.2111485255710072D+01   0.9748601873694136D-03  -0.1709177520251039D-09
   0.2345993221933724D+01   0.7007806571062660D-03  -0.1651563153197262D-09
   0.2612495009930844D+01   0.5167979811052178D-03  -0.1630088822424430D-09
   0.2916984570477848D+01   0.3828065733824403D-03  -0.1557593854333271D-09
   0.3266934871581749D+01   0.2830854629256306D-03  -0.1549588621528963D-09
   0.3671758140882298D+01   0.2087729455436860D-03  -0.1512318383784860D-09
   0.4143442314569344D+01   0.1534769531120022D-03  -0.1486769193920228D-09
   0.4697445226110210D+01   0.1123562040326752D-03  -0.1472571032257998D-09
   0.5353973100920618D+01   0.8175221985657395D-04  -0.1451370215552787D-09
   0.6139844387895875D+01   0.5894429819511890D-04  -0.1436889794330563D-09
   0.7091266484425319D+01   0.4196875125118042D-04  -0.1426188147621635D-09
   0.8258074523321996D+01   0.2942150885807123D-04  -0.1417458532709309D-09
   0.9710383029886540D+01   0.2025054334971717D-04  -0.1417073402463571D-09
   0.1154935789936476D+02   0.1363300080959531D-04  -0.1407792926769395D-09
   0.1392530609899710D+02   0.8927910808379660D-05  -0.1410592706999196D-09
   0.1706936855807103D+02   0.5650897506135906D-05  -0.1414467350637302D-09
   0.2135189727044204D+02   0.3428296379756167D-05  -0.1397833656893875D-09
   0.2739665613112830D+02   0.1971677079401429D-05  -0.1402732323249772D-09
   0.3632135308034874D+02   0.1058455302418004D-05  -0.1402163795036603D-09
   0.5029362477052696D+02   0.5183089093494448D-06  -0.1389509663166628D-09
   0.7398170762169443D+02   0.2232992863388601D-06  -0.1386070954279062D-09
   0.1190210316547486D+03   0.7952919191617421D-07  -0.1383200189277319D-09
   0.2215528733967777D+03   0.2067350176573896D-07  -0.1388069459183895D-09
   0.5454655223627067D+03   0.2787861460038876D-08  -0.1390822755609575D-09
   0.2876870662092502D+04  -0.4627967448074155D-11  -0.6726700201895319D-10

 

[jedishrfu]

This looks like a homework problem.

Whether it is or it isn't, we can't help you until you make an attempt at a solution.

Also it would help if you tell us what fortran you're using and what this program is for.

 

[Raziel13]

is not a homework problem ,
i using fortran 90 , the program is to fit the first four moments of Chebyschev from BSA of pion.

 

[jedishrfu]

Great so what have you written so far?

We have [ code ] tags to bracket your sample code:

[ code=fortran ]
  DO I=1,20
    PRINT(...)
  END
[ /code ]

When posting don't forget to remove the spaces around "code=fortran" and "/code"

 

[Raziel13]

i try, but the variables are not being allocated , so when I do the last DO to write, it just repeats the first read value

PROGRAM chesb_tay

      IMPLICIT None
 
      call chebtay
 
      STOP
      END

 
 
      subroutine chebtay 
      complex*8 EH0,EH1,EH2,EH3
      double precision k2, EH0R ,EH1R, EH2R,EH3R
      double precision  EH0IM,EH1IM,EH2IM,EH3IM,
      integer  i0,i1,i2,i3,j
      open (1,file='data.out',status='old')
      open (12,file='chesb.out',status='unknown')
 
 
      do i0 =1,64
      read(1,*) k2,EH0R,EH0IM
      EH0=CMPLX(EH0R,EH0IM)
      end do
      do i1= 65, 129
      read(1,*) k2,EH1R,EH1IM
      EH1=CMPLX(EH1R,EH1IM)
      end do

      do j =1,64
      write (12,*) k2,EH0,EH1
      end do
      close(unit=12)
      close(unit=1)
      return
      end

 

[jedishrfu]

So you want to save the 64 values into a variable then you need to look at arrays in fortran.

Here's something to help here:

https://www.tutorialspoint.com/fortran/fortran_arrays.htm

https://web.stanford.edu/class/me200c/tutorial_90/07_arrays.html

 

[FactChecker]

I could understand it printing the LAST value in the input file, but not the first. It looks like EH0 and EH1 will have only the final value in their respective loops when they get to the last write loop. You should make them arrays and increment their index. If you have a debugger, step through to see what is happening, otherwise, just put some writes to the terminal (or a text file) inside the loops to see what is happening at each step of each loop. It should soon become obvious.

 

[sysprog]

As @jedishrfu and @FactChecker suggested, you'll need to use arrays -- you're not referencing your counter variables (i0, i1, j) to subscript the variables that are being assigned the input values inside the loops, so you'll just re-use the same non-arrayed variables each time through the loops.

Also (just a sideline observation), your input file has 256 lines, but your code reads only the first half of them.

I shudder to consider how you got to the point of implementing moment methods for pion amplitude distribution analysis without even a rudimentary knowledge of programming -- whether in fortran or another computer language -- but I did encounter something like that a long time ago --

I had a buddy who'd done well as an undergrad in his homeland in Civil Engineering, and was still trying to do all his Fourier transforms and such by hand in graduate school in the US in the early '80s. He was failing. He said "too many variables".

When I learned that he was trying to do it all by hand, I asked him how fast can you count to ten million, and he shook his head and smiled.

I showed him how to use TSO (time sharing option) on the university's IBM mainframe, gave him a copy of a textbook entitled "Ten Statement Fortran", and he did the rest -- he aced all his classes, and wound up writing a faster 'fast Fourier transform' subroutine that got accepted into the university's fortlib as an alternate.

Best wishes to you, @Raziel13.

 

[Dr Transport]

try something like this, fill in the rest of your data. I think I have my FORTRAN syntax correct, it has been a while

ic = complex(0,1)
do ii = 1,256
read(1,*) K(ii),RealPart(ii),ImagPart(ii)
enddo

k2 = K(1:64)
E0 = RealPart(1:64) + ic*ImagPart(1:64)

...