From f7dda9b2a0797210b8273e349c37ea48b42e43d4 Mon Sep 17 00:00:00 2001 From: Yeongdo Park Date: Fri, 15 Mar 2019 14:03:41 -0400 Subject: [PATCH] makefile change compiler (ifort->gfortran) / fix array boundary bug / added sample file --- code/makefile | 13 +- code/ysolve.f90 | 2 +- sample/4pi-IC1/itape | 19 ++ sample/4pi-IC1/sfield.diff | 512 +++++++++++++++++++++++++++++++++++++ 4 files changed, 539 insertions(+), 7 deletions(-) create mode 100644 sample/4pi-IC1/itape create mode 100644 sample/4pi-IC1/sfield.diff diff --git a/code/makefile b/code/makefile index 6d1afd7..2692fd1 100644 --- a/code/makefile +++ b/code/makefile @@ -1,17 +1,18 @@ -flags = -O3 -r8 -i8 -msse3 -axT -m64 -fPIC -i_dynamic -compiler = ifort +flags = -Wall -O3 -fdefault-integer-8 -fdefault-double-8 -fdefault-real-8 -march=native -ex : test.o ysolve.mod Compact.mod +compiler = gfortran + +ex : test.o ysolve.mod compact.mod ${compiler} -o ex test.o ysolve.o Compact.o test.o : test.f90 ysolve.mod ${compiler} -c ${flags} test.f90 -ysolve.mod: ysolve.f90 Compact.mod +ysolve.mod: ysolve.f90 compact.mod ${compiler} -c ${flags} ysolve.f90 -Compact.mod: Compact.f90 +compact.mod: Compact.f90 ${compiler} -c ${flags} Compact.f90 clean: - rm *.o *.mod ex + rm -f *.o *.mod ex diff --git a/code/ysolve.f90 b/code/ysolve.f90 index a96504e..284bddb 100644 --- a/code/ysolve.f90 +++ b/code/ysolve.f90 @@ -32,7 +32,7 @@ CALL READ_INTRO - CALL ludcmp(nx,1,1,1,0,0) + CALL ludcmp(nx,5,5,1,0,0) CALL SET_IC diff --git a/sample/4pi-IC1/itape b/sample/4pi-IC1/itape new file mode 100644 index 0000000..3e343cc --- /dev/null +++ b/sample/4pi-IC1/itape @@ -0,0 +1,19 @@ + nx 512 + l_0 4. !3[pi] + int_pr 400 + tar_lo 0.60 + dt 0.001 + sc 0.75 + vis 0.020 !2pi(162)ref !2pi(256)mid !2pi(256)Max !2p(162)Low + pre 2.10E+4 !185.75 !185.75 !316.88 !440.5 !102.69 + ac 26.7 !22.250 !22.25 !22.56 !22.3 !22.39 + bc 3. !7.000 !7.0 !7.0 !7.0 !7.0 + u0 0.21 + tf 80. + dt_uf 0.4 + ctmp 0. !1.0e-14 + c_cut 0.001 !0.012 ! c < c_cut -> wrate = 0. + c_ref 0.01 !0.003 ! + min_wr 0. ! 5.0e-14 + prof_wr 1.0 + diff --git a/sample/4pi-IC1/sfield.diff b/sample/4pi-IC1/sfield.diff new file mode 100644 index 0000000..cb4598c --- /dev/null +++ b/sample/4pi-IC1/sfield.diff @@ -0,0 +1,512 @@ + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.10000000000000000000E+01 + 0.99999999999999955591E+00 + 0.99999999999999900080E+00 + 0.99999999999999833467E+00 + 0.99999999999999777955E+00 + 0.99999999999999722444E+00 + 0.99999999999999666933E+00 + 0.99999999999999600320E+00 + 0.99999999999999544809E+00 + 0.99999999999999489297E+00 + 0.99999999999999433786E+00 + 0.99999999999999367173E+00 + 0.99999999999999311662E+00 + 0.99999999999999256151E+00 + 0.99999999999999200639E+00 + 0.99999999999999145128E+00 + 0.99999999999999078515E+00 + 0.99999999999999023004E+00 + 0.99999999999998967493E+00 + 0.99999999999998911981E+00 + 0.99999999999998845368E+00 + 0.99999999999998789857E+00 + 0.99999999999998734346E+00 + 0.99999999999998678835E+00 + 0.99999999999998612221E+00 + 0.99999999999998556710E+00 + 0.99999999999998501199E+00 + 0.99999999999998445688E+00 + 0.99999999999998356870E+00 + 0.99999999999998190336E+00 + 0.99999999999997968292E+00 + 0.99999999999997735145E+00 + 0.99999999999997501998E+00 + 0.99999999999997268851E+00 + 0.99999999999997035705E+00 + 0.99999999999996802558E+00 + 0.99999999999996569411E+00 + 0.99999999999996336264E+00 + 0.99999999999996103117E+00 + 0.99999999999995869970E+00 + 0.99999999999995636824E+00 + 0.99999999999995381472E+00 + 0.99999999999994970690E+00 + 0.99999999999994371169E+00 + 0.99999999999993716138E+00 + 0.99999999999992905675E+00 + 0.99999999999991795452E+00 + 0.99999999999990329957E+00 + 0.99999999999988420374E+00 + 0.99999999999985988985E+00 + 0.99999999999983113508E+00 + 0.99999999999979627407E+00 + 0.99999999999975286435E+00 + 0.99999999999969890752E+00 + 0.99999999999963207209E+00 + 0.99999999999954969354E+00 + 0.99999999999944855222E+00 + 0.99999999999932664974E+00 + 0.99999999999917976723E+00 + 0.99999999999900135439E+00 + 0.99999999999878363965E+00 + 0.99999999999851818533E+00 + 0.99999999999819488838E+00 + 0.99999999999780142534E+00 + 0.99999999999732425149E+00 + 0.99999999999674615836E+00 + 0.99999999999604549661E+00 + 0.99999999999519473270E+00 + 0.99999999999416089302E+00 + 0.99999999999290478669E+00 + 0.99999999999137878515E+00 + 0.99999999998952582292E+00 + 0.99999999998727751027E+00 + 0.99999999998454935923E+00 + 0.99999999998123856315E+00 + 0.99999999997721944478E+00 + 0.99999999997233945948E+00 + 0.99999999996641442124E+00 + 0.99999999995922084217E+00 + 0.99999999995048860502E+00 + 0.99999999993988941682E+00 + 0.99999999992702381935E+00 + 0.99999999991140720024E+00 + 0.99999999989244925391E+00 + 0.99999999986943466368E+00 + 0.99999999984149545718E+00 + 0.99999999980757814377E+00 + 0.99999999976640552291E+00 + 0.99999999971642616892E+00 + 0.99999999965575603333E+00 + 0.99999999958210794571E+00 + 0.99999999949270379496E+00 + 0.99999999938417261092E+00 + 0.99999999925242266663E+00 + 0.99999999909248704633E+00 + 0.99999999889833801614E+00 + 0.99999999866265665283E+00 + 0.99999999837655806356E+00 + 0.99999999802925620962E+00 + 0.99999999760765656553E+00 + 0.99999999709586462338E+00 + 0.99999999647458592822E+00 + 0.99999999572039888207E+00 + 0.99999999480487267167E+00 + 0.99999999369349279466E+00 + 0.99999999234436087736E+00 + 0.99999999070661538259E+00 + 0.99999998871851258908E+00 + 0.99999998630510178455E+00 + 0.99999998337539874882E+00 + 0.99999997981895571986E+00 + 0.99999997550169827942E+00 + 0.99999997026086850926E+00 + 0.99999996389889023174E+00 + 0.99999995617591974639E+00 + 0.99999994680080428466E+00 + 0.99999993542010701120E+00 + 0.99999992160478190506E+00 + 0.99999990483400025276E+00 + 0.99999988447551402260E+00 + 0.99999985976182004244E+00 + 0.99999982976122636646E+00 + 0.99999979334272626197E+00 + 0.99999974913336220350E+00 + 0.99999969546646616525E+00 + 0.99999963031882654896E+00 + 0.99999955123441575111E+00 + 0.99999945523179711859E+00 + 0.99999933869172663581E+00 + 0.99999919722070995665E+00 + 0.99999902548536856450E+00 + 0.99999881701137705026E+00 + 0.99999856393938424404E+00 + 0.99999825672871978099E+00 + 0.99999788379771803459E+00 + 0.99999743108709260309E+00 + 0.99999688152991039747E+00 + 0.99999621440817809681E+00 + 0.99999540457178559549E+00 + 0.99999442149036110727E+00 + 0.99999322810229362979E+00 + 0.99999177941753114585E+00 + 0.99999002082148102932E+00 + 0.99998788601607191406E+00 + 0.99998529452035334053E+00 + 0.99998214863640944028E+00 + 0.99997832976620304457E+00 + 0.99997369394049728797E+00 + 0.99996806639129753869E+00 + 0.99996123496319622159E+00 + 0.99995294211523566563E+00 + 0.99994287521176139499E+00 + 0.99993065473623277128E+00 + 0.99991581998366618755E+00 + 0.99989781169231983782E+00 + 0.99987595095984271865E+00 + 0.99984941364904711403E+00 + 0.99981719931839607707E+00 + 0.99977809350593149773E+00 + 0.99973062194474593412E+00 + 0.99967299498395911783E+00 + 0.99960304011990008810E+00 + 0.99951812009394858105E+00 + 0.99941503346938664887E+00 + 0.99928989393905676142E+00 + 0.99913798381370888535E+00 + 0.99895357616905511300E+00 + 0.99872971896017059823E+00 + 0.99845797296804317256E+00 + 0.99812809362601828678E+00 + 0.99772764480696751743E+00 + 0.99724152999601212155E+00 + 0.99665142316932142652E+00 + 0.99593507792424551539E+00 + 0.99506548881571632581E+00 + 0.99400987328225554052E+00 + 0.99272843578152536459E+00 + 0.99117286754502942348E+00 + 0.98928452538414335393E+00 + 0.98699222082667914524E+00 + 0.98420953639897412746E+00 + 0.98083156792273673563E+00 + 0.97673096989153151970E+00 + 0.97175315495585878889E+00 + 0.96571046669015281338E+00 + 0.95837510625746358262E+00 + 0.94947054699732003069E+00 + 0.93866111494419091077E+00 + 0.92553934677775184792E+00 + 0.90961066049177308201E+00 + 0.89027479554303234899E+00 + 0.86680343040412299871E+00 + 0.83831348371527447672E+00 + 0.80373627840195971395E+00 + 0.76178549678038254278E+00 + 0.71093677231427665486E+00 + 0.64946438797693928358E+00 + 0.57567439607637960197E+00 + 0.48867726080226464580E+00 + 0.39021372632551271087E+00 + 0.28725565231896382778E+00 + 0.19226393360764873486E+00 + 0.11726829144399875116E+00 + 0.66388592405444438183E-01 + 0.35735052559352491197E-01 + 0.18663494688699135138E-01 + 0.95861739814404323123E-02 + 0.48804801297117677769E-02 + 0.24734156922963800203E-02 + 0.12506014489007202545E-02 + 0.63157754912839003149E-03 + 0.31876768266224391911E-03 + 0.16083866804920426041E-03 + 0.81140988522889208432E-04 + 0.40931404387194565914E-04 + 0.20646959308359734089E-04 + 0.10414706605319981654E-04 + 0.52533178365350360556E-05 + 0.26498308272822090778E-05 + 0.13366002077599419203E-05 + 0.67419317163842360068E-06 + 0.34006888383009713071E-06 + 0.17153363956392798050E-06 + 0.86523013663131540347E-07 + 0.43642933836094156241E-07 + 0.22013860976036643648E-07 + 0.11103975431753144437E-07 + 0.56009379360483640550E-08 + 0.28251598532943308997E-08 + 0.14250342113031681796E-08 + 0.71879914987794337672E-09 + 0.36256829020668265959E-09 + 0.18288247147590120560E-09 + 0.92247444880166095687E-10 + 0.46530381056859186912E-10 + 0.23470312527053937351E-10 + 0.11838621511391101341E-10 + 0.59714994900649730370E-11 + 0.30120741781970911299E-11 + 0.15193153528658277907E-11 + 0.76635534349336187954E-12 + 0.38655603090620881529E-12 + 0.19498208800774406035E-12 + 0.98350592421664555797E-13 + 0.49608859606213857086E-13 + 0.25023122798035452088E-13 + 0.12621871970778032775E-13 + 0.63665775583831450198E-14 + 0.32113548529697610713E-14 + 0.16198341883247153929E-14 + 0.81705788298032221448E-15 + 0.41213081496362010617E-15 + 0.20788222252141758030E-15 + 0.10485752792898019950E-15 + 0.52891012179959855713E-16 + 0.26678667947572384215E-16 + 0.13456942760617543878E-16 + 0.67877942338950767493E-17 + 0.34238200593777083377E-17 + 0.17270034115736916252E-17 + 0.87111493357196944495E-18 + 0.43939764241728761359E-18 + 0.22163583784540955660E-18 + 0.11179496627973696902E-18 + 0.56390313980742287540E-19 + 0.28443744979446820243E-19 + 0.14347262346013967342E-19 + 0.72368788629658713366E-20 + 0.36503420941342606531E-20 + 0.18412630163533656484E-20 + 0.92874843178081658893E-21 + 0.46846845989643572762E-21 + 0.23629940079353131369E-21 + 0.11919138980610546423E-21 + 0.60121131734583241139E-22 + 0.30325600590169130271E-22 + 0.15296486021161559889E-22 + 0.77156752064934491220E-23 + 0.38918509656231018204E-23 + 0.19630821066540337005E-23 + 0.99019499757443654281E-24 + 0.49946262048744589289E-24 + 0.25193311405861048767E-24 + 0.12707716524877405253E-24 + 0.64098782678911125352E-25 + 0.32331960922128872140E-25 + 0.16308510916760300752E-25 + 0.82261490097266145613E-26 + 0.41493382121528909961E-26 + 0.20929608226734914759E-26 + 0.10557069058425297508E-26 + 0.53250737375005118717E-27 + 0.26860116338054305619E-27 + 0.13548466843061103341E-27 + 0.68339597448974827218E-28 + 0.34471063284033737477E-28 + 0.17387492000067136361E-28 + 0.87703960728252109386E-29 + 0.44238609728148886954E-29 + 0.22314323941917978960E-29 + 0.11255531221362557662E-29 + 0.56773838815292403712E-30 + 0.28637198106715474570E-30 + 0.14444841717878903835E-30 + 0.72860987124869990643E-31 + 0.36751689969988859164E-31 + 0.18537859133521737079E-31 + 0.93506508553734342104E-32 + 0.47165463276710835845E-32 + 0.23790653297982815327E-32 + 0.12000204069325784104E-32 + 0.60530030807381434816E-33 + 0.30531852694971670665E-33 + 0.15400521304108946875E-33 + 0.77681514714426336863E-34 + 0.39183204316061801542E-34 + 0.19764335261978557938E-34 + 0.99692956501713025506E-35 + 0.50285959250913387658E-35 + 0.25364657509591150488E-35 + 0.12794144929574403520E-35 + 0.64534734765118917729E-36 + 0.32551858792667241402E-36 + 0.16419429237826124546E-36 + 0.82820971365453889449E-37 + 0.41775589142375583560E-37 + 0.21071955803219700232E-37 + 0.10628870363971477251E-37 + 0.53612909152386015623E-38 + 0.27042798805083926584E-38 + 0.13640613404014510106E-38 + 0.68804392392550865578E-39 + 0.34705509732537829799E-39 + 0.17505748745274577409E-39 + 0.88300457620241207641E-40 + 0.44539487738556115796E-40 + 0.22466089321356465703E-40 + 0.11332082945315190675E-40 + 0.57159972099563394814E-41 + 0.28831966957792057037E-41 + 0.14543084751812942997E-41 + 0.73356533186943837181E-42 + 0.37001647539297827657E-42 + 0.18663939817511813594E-42 + 0.94142470045892432517E-43 + 0.47486247560796649892E-43 + 0.23952459567993361850E-43 + 0.12081820502279570138E-43 + 0.60941710906531297171E-44 + 0.30739507572675039900E-44 + 0.15505264157414346065E-44 + 0.78209846407854063681E-45 + 0.39449699231439324568E-45 + 0.19898757521338620792E-45 + 0.10037099359618692226E-45 + 0.50627966819950823847E-46 + 0.25537168981651044448E-46 + 0.12881161155782589064E-46 + 0.64973651872085706073E-47 + 0.32773252244420842013E-47 + 0.16531101942538169569E-47 + 0.83384257807709688750E-48 + 0.42059715524777974068E-48 + 0.21215271521691098301E-48 + 0.10701160008415271968E-48 + 0.53977544151920490809E-49 + 0.27226723741924465133E-49 + 0.13733386677109485153E-49 + 0.69272348524470945224E-50 + 0.34941550710824094672E-50 + 0.17624809784611188611E-50 + 0.88901011438939563904E-51 + 0.44842412096652436907E-51 + 0.22618886895647043808E-51 + 0.11409155316965344840E-51 + 0.57548731574270907696E-52 + 0.29028060481241547292E-52 + 0.14641995961546805913E-52 + 0.73855449583514009377E-53 + 0.37253305133453533327E-53 + 0.18790878008221824082E-53 + 0.94782756873507371397E-54 + 0.47809213580180706200E-54 + 0.24115366323488939138E-54 + 0.12163992029279474536E-54 + 0.61356190945977684487E-55 + 0.30948574762674016239E-55 + 0.15610719390530740726E-55 + 0.78741771349099561397E-56 + 0.39718006477978031807E-56 + 0.20034093614777023210E-56 + 0.10105363244275728726E-56 + 0.50972276964749757689E-57 + 0.25710797284390360111E-57 + 0.12968633062211458189E-57 + 0.65412284373552106435E-58 + 0.32988266120339507051E-58 + 0.16624599030477054159E-58 + 0.83494650946974640202E-59 + 0.41249148223150098961E-59 + 0.18716116359996794534E-59 + 0.43902014526217839815E-60