Fortran complex number read

  • Fortran
  • Thread starter Raziel13
  • Start date
  • #1
3
0

Main Question or Discussion Point

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


Code:
   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
 

Answers and Replies

  • #2
11,473
4,996
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.
 
  • #3
3
0
is not a homework problem ,
i using fortran 90 , the program is to fit the first four moments of Chebyschev from BSA of pion.
 
  • #4
11,473
4,996
Great so what have you written so far?

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

Fortran:
[ code=fortran ]
  DO I=1,20
    PRINT(...)
  END
[ /code ]
When posting don't forget to remove the spaces around "code=fortran" and "/code"
 
  • #5
3
0
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
Fortran:
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
 
  • #7
FactChecker
Science Advisor
Gold Member
5,375
1,947
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.
 
  • #8
1,101
555
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.
 
Last edited:
  • #9
Dr Transport
Science Advisor
Gold Member
2,327
442
try something like this, fill in the rest of your data. I think I have my FORTRAN syntax correct, it has been a while

Code:
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)

......
 

Related Threads for: Fortran complex number read

Replies
2
Views
587
Replies
3
Views
4K
Replies
6
Views
9K
Replies
1
Views
8K
  • Last Post
Replies
9
Views
30K
  • Last Post
Replies
5
Views
9K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
3
Views
3K
Top