// inas2.cpp  Aug 6 2002   Q versus. 1/R
// inas.cpp  Aug 5 2002

#include <iostream.h>
#include <conio.h>
#include <stdlib.h>
#include <stdio.h>
#include <graphics.h>

void p(float diamA,float freqwavenum,float taud);

float bx0=30,by0=380,sfx=3000,sfy=20;
void makescr(void);

void main(void)
{ //
int gdriver=VGA;
int gmode=VGAHI;
initgraph(&gdriver,&gmode,"c:\\tc\\bgi");
cleardevice();


setcolor(14);
line(0+20,165,360+20,450);
setcolor(11);
line(bx0,by0,bx0+300,by0);
line(bx0,by0,bx0,by0-300);
int iy,iDA;
for (iy=1;iy<3;iy++)
{ //
line(bx0,by0-sfy*iy*10,bx0+11,by0-sfy*iy*10);
} //
setcolor(12);
for (iy=1;iy<20;iy++)
{ //
if ((iy%10)!=0)
line(bx0+1,by0-sfy*iy,bx0+7,by0-sfy*iy);
} //

setcolor(13);
for (iDA=11;iDA<=19;iDA=iDA+1)
{ //
line(bx0+sfx/iDA,by0,bx0+sfx/iDA,by0-4);
} //
setcolor(12);
for (iDA=10;iDA<=90;iDA=iDA+10)
{ //
line(bx0+sfx/iDA,by0,bx0+sfx/iDA,by0-7);
} //
setcolor(11);
for (iDA=100;iDA<=400;iDA=iDA+100)
{ //
line(bx0+sfx/iDA,by0,bx0+sfx/iDA,by0-11);
} //


// From: Cerullo de Silvestri, Banin

//Diam  f    damp time
// A   cm-1   (ps)
//p(24 ,28.52, 0.75 ); // in toluene solution
p(24 ,28.49, 1.12 );
p(29 ,24.7, 1.36 );
p(36 ,22.9, 1.47 );
p(42 ,22.07, 1.78 );
p(47 ,20.54, 2.14 );
p(55 ,18.31, 2.34 );

gotoxy(2,8);cout<<"Q";
gotoxy(2,12);cout<<"10";
gotoxy(3,24);cout<<"0";
gotoxy(11,9 );cout<<"InAs Nanoparticles";
gotoxy(8,10);cout<<"in polyvinyl butyral film";
gotoxy(22,12);cout<<"Data from:";
gotoxy(22,13);cout<<"G. Cerullo, ";
gotoxy(22,14);cout<<"S. De Silvestri,";
gotoxy(22,15);cout<<"U. Banin ";
gotoxy(22,16);cout<<"Phys. Rev. B";
gotoxy(22,17);cout<<"Vol. 60 (1999)";
gotoxy(10,25);cout<<"Inverse diameter";
makescr();
//getch();
} //

void p(float diamA,float freqwavenum, float tauD)
{ // p
float omegaR;
// lambda = 0.01/freqwavenum;
// freqHz = 3.0e8*freqwavenum/0.01;
// omega = 2.*3.14159*3.0e8*freqwavenum/0.01;
//omegaR = 2.*3.14159*3.0e8*(freqwavenum/0.01)*diamA*1e-10*0.5;
omegaR = 2.*3.14159*3.0e8*(freqwavenum/0.01)*(tauD*1e-12);

putpixel(bx0+sfx/diamA,by0-sfy*omegaR,13);
setcolor(13);
circle(bx0+sfx/diamA,by0-sfy*omegaR,4);
return;
} // p


void makescr(void)
{ // makescr
int scrwid=300;
int scrhgt=300;
int scrx0=0;
int scry0=110;
setcolor(15);
line(scrx0+scrwid,scry0-1,scrx0+scrwid,scry0+scrhgt);
line(scrx0,scry0-1,scrwid+scrwid-1,scry0-1);
line(scrx0,scry0+scrhgt,scrwid+scrwid-1,scry0+scrhgt);
gotoxy(70,5);cout<<"Make .scr ";
gotoxy(70,6);cout<<"file? ";
gotoxy(70,7);cout<<"(y/n/q)  ";
char ch678;
ch678=getch();
if (ch678=='q') {closegraph();exit(0);}
if (ch678!='y') return;
char scrfilename[40];
FILE *scrfout;
gotoxy(65,8);cout<<"Filename";
gotoxy(65,9);cout<<"? ";cin>>scrfilename;

scrfout=fopen(scrfilename,"wb");
unsigned char ch1b,ch2b,ch3b,ch4b,ch5,ch6;
ch3b=scrwid/256;
ch1b=scrwid%256;
ch4b=scrhgt/256;
ch2b=scrhgt%256;
fprintf(scrfout,"%c%c%c%c%c%c%c%c%c%c%c%c%c",'s','c','r',ch3b,ch1b,ch4b,ch2b,'0','0','0','0',13,10);
int ipixcol,ix,iy;
gotoxy(65,10);cout<<"saving ";
for (iy=scry0;iy<scry0+scrhgt;iy++)
{ // iy
gotoxy(65,11);cout<<iy;
for (ix=scrx0;ix<scrx0+scrwid;ix++)
{ //
ipixcol=getpixel(ix,iy);
//putpixel(ix-100,iy-100,ipixcol);
char ch1;
ch1=ipixcol+'A';
fprintf(scrfout,"%c",ch1);
} //
fprintf(scrfout,"%c%c",13,10);
} // iy
fclose(scrfout);
return;
} //