Changeset 36 for trunk/examples/fermi_fortran/fermi_example.f
- Timestamp:
- Aug 11, 2005, 11:15:45 AM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/fermi_fortran/fermi_example.f
r30 r36 8 8 9 9 CHARACTER*100 inFile, outFile 10 character*100 T_path, Ef_path, outPath 10 character*100 T_path, Ef_path, outPath, outPath2,outPath3 11 11 character*40 xy 12 12 integer handle … … 50 50 15 continue 51 51 52 outPath = "output.curve(f13).about.group" 53 call rp_lib_put_str(handle,outPath,"group1",1) 52 54 outPath = "output.curve(f13).xaxis.label" 53 55 call rp_lib_put_str(handle,outPath,"f13 x label",1) … … 57 59 outPath = "output.curve(f13).component.xy" 58 60 61 outPath2 = "output.curve(f14).about.group" 62 call rp_lib_put_str(handle,outPath2,"group1",1) 63 outPath2 = "output.curve(f14).component.xy" 64 65 outPath3 = "output.curve(f15).about.group" 66 call rp_lib_put_str(handle,outPath3,"group1",1) 67 outPath3 = "output.curve(f15).component.xy" 68 59 69 do 16,E = Emin,Emax,dE 60 70 f = -1.0/(1.0+exp((E-Ef)/kT)) 61 71 write(xy,'(E20.12,F13.9,A)')f,E,char(10) 62 72 call rp_lib_put_str(handle,outPath,xy,1) 73 f = -4.0/(1.0+exp((E-Ef)/kT)) 74 write(xy,'(E20.12,F13.9,A)')f,E,char(10) 75 call rp_lib_put_str(handle,outPath2,xy,1) 76 f = 1.0/(1.0+exp((E-Ef)/kT)) 77 write(xy,'(E20.12,F13.9,A)')f,E,char(10) 78 call rp_lib_put_str(handle,outPath3,xy,1) 63 79 16 continue 64 80
Note: See TracChangeset
for help on using the changeset viewer.