BUGFIX:

wrong format in write_legend() resulted in the namefiles not being printed.
aamaricci · May 16, 2024 · 7d80916 · 7d80916
1 parent 47d9783
commit 7d80916
Show file tree

Hide file tree

Showing 3 changed files with 72 additions and 72 deletions.
diff --git a/src/ED_NONSU2/ED_OBSERVABLES_NONSU2.f90 b/src/ED_NONSU2/ED_OBSERVABLES_NONSU2.f90
@@ -932,47 +932,47 @@ subroutine write_legend()
     !
     unit = free_unit()
     open(unit,file="observables_info.ed")
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# docc_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# docc_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"docc_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_up_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_up_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_up_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_dw_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_dw_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_dw_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# magX_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# magX_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"magX_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# magY_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# magY_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"magY_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# magZ_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# magZ_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"magZ_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Sz2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Sz2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"Sz2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# n2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# n2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"n2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Z_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Z_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(ispin-1)*Norb))//"z_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# sig_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# sig_last.ed"
     write(unit,"(A1,90(A10,6X))") "#",((reg(txtfy(iorb+(ispin-1)*Norb))//"sig_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# phisc_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# phisc_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"phisc_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# imp_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# imp_last.ed"
     write(unit,"(A1,90(A10,6X))") "#", "1s2tot", "2egs"
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# exciton_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# exciton_last.ed"
     write(unit,"(A1,90(A10,6X))") "#","1S_0" , "2T_z", "3T_x", "4T_y"
     close(unit)
     !

diff --git a/src/ED_NORMAL/ED_OBSERVABLES_NORMAL.f90 b/src/ED_NORMAL/ED_OBSERVABLES_NORMAL.f90
@@ -721,38 +721,38 @@ subroutine write_legend()
 
     unit = free_unit()
     open(unit,file="observables_info.ed")
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# docc_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# docc_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"docc_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_up_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_up_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_up_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_dw_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_dw_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_dw_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# magZ_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# magZ_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"magZ_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Sz2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Sz2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"Sz2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# n2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# n2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"n2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Z_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Z_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(ispin-1)*Norb))//"z_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# sig_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# sig_last.ed"
     write(unit,"(A1,90(A10,6X))") "#",((reg(txtfy(iorb+(ispin-1)*Norb))//"sig_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# imp_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# imp_last.ed"
     write(unit,"(A1,90(A10,6X))") "#", "1s2tot", "2egs", "3nph", "4w_ph","5X_ph", "6X2_ph"
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# exciton_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# exciton_last.ed"
     write(unit,"(A1,90(A10,6X))") "#","1S_0" , "2T_z"
     close(unit)
     !

diff --git a/src/ED_SUPERC/ED_OBSERVABLES_SUPERC.f90 b/src/ED_SUPERC/ED_OBSERVABLES_SUPERC.f90
@@ -575,38 +575,38 @@ subroutine write_legend()
 
     unit = free_unit()
     open(unit,file="observables_info.ed")
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# docc_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# docc_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"docc_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_up_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_up_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_up_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# dens_dw_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# dens_dw_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"dens_dw_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# magz_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# magz_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",(reg(txtfy(iorb))//"magz_"//reg(txtfy(iorb)),iorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Sz2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Sz2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"Sz2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# n2_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# n2_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"n2_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# Z_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# Z_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(ispin-1)*Norb))//"z_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# sig_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# sig_last.ed"
     write(unit,"(A1,90(A10,6X))") "#",((reg(txtfy(iorb+(ispin-1)*Norb))//"sig_"//reg(txtfy(iorb))//"s"//reg(txtfy(ispin)),iorb=1,Norb),ispin=1,Nspin)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# phisc_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# phisc_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ",((reg(txtfy(iorb+(jorb-1)*Norb))//"phisc_"//reg(txtfy(iorb)),iorb=1,Norb),jorb=1,Norb)
-    write(unit,"(A1,90(A10,6X))") "# *****"
-    write(unit,"(A1,90(A10,6X))") "# imp_last.ed"
+    write(unit,"(A)") "# *****"
+    write(unit,"(A)") "# imp_last.ed"
     write(unit,"(A1,90(A10,6X))") "# ", "1s2tot", "2egs", "3nph", "4w_ph", "5X_ph", "6X2_ph"
     close(unit)