% construction='L'; %'L','M','H'
% glz_pct=32; %15,30,45,60,75
% glazing_type='double';%'triple'

%loop to go through all combinations
for z=1:1:2 %2 %glazing type
    for y=1:1:5 %5%glazing area
        for x=1:1:3 %3%construction
            switch x
                case 1
                    construction='L';
                case 2
                    construction='M';
                case 3
                    construction='H';
            end
            switch z
                case 1
                    glz='D';
                case 2
                    glz='T';
            end
        glz_pct=15*y;
        
                 
%construction specifier
switch construction
    case 'L'
        int_wall='AInternalWal';
        int_wall_layers=3;
        floor='AFloor_no_sl';
        floor_bottom='AFloor_no_bo';
        floor_layers=3;
        int_wall_air_gap=1;
    case 'M'
        int_wall='AInternalWal';
        int_wall_layers=3;
        floor='AFloor_sl2in';
        floor_bottom='AFloor_sl2_b';
        floor_layers=4;
        int_wall_air_gap=1;
    case 'H'
        int_wall='8in_concrete';
        int_wall_layers=1;
        floor='AFloor_sl8in';
        floor_bottom='AFloor_sl8_b';
        floor_layers=4;
        int_wall_air_gap=0;
end

switch construction
    case 'L'
        gap_location=2;
    case 'M'
        gap_location=3;
    case 'H'
        gap_location=3;
end

%glazing specifier
switch glz
    case 'D'
        glz_name='DblglzLowArg';
        glazing_layers=3; %including air gaps
        num_air_gaps=1;
    case 'T'
        glz_name='TrpglzLowArg';
        glazing_layers=5;  %including air gaps
        num_air_gaps=2;
end

%window vertex specifier
frame_area=50*glz_pct/100*0.416/2;
frame_width=(frame_area/2.5)^.5;
frame_height=2.5*frame_width;
glazed_area=glz_pct/100*50;
glazed_width=(glazed_area*8/5)^.5; 
glazed_height=5/8*glazed_width;
x1frame=1-frame_width/2;
x2frame=1+frame_width/2;
y1frame=2.5-frame_height/2;
y2frame=2.5+frame_height/2;
x1glazing=6-glazed_width/2;
x2glazing=6+glazed_width/2;
y1glazing=2.5-glazed_height/2;
y2glazing=2.5+glazed_height/2;        
        
a=['C:\Esru\Models\SDD\' construction glz num2str(glz_pct)];
b=[construction glz num2str(glz_pct)];
copyfile('C:\Esru\Models\SDD\Master',[a]);
movefile(['C:\Esru\Models\SDD\' b '\cfg\Master.cfg'],[a '\cfg\' b '.cfg']);

fid=fopen([a '\cfg\' b '.cfg'], 'wt');
fprintf(fid,'* CONFIGURATION3.0\n');
fprintf(fid,'# ESRU system configuration defined by file\n'); 
fprintf(fid,'# %s.cfg\n',b);
fprintf(fid,'*date Mon Mar 24 15:20:55 2008  # latest file modification\n'); 
fprintf(fid,'*root %s \n',b);
fprintf(fid,'*zonpth ..\\zones                  # path to zones\n');
fprintf(fid,'*netpth ..\\nets                   # path to networks\n');
fprintf(fid,'*ctlpth ..\\ctl                    # path to controls\n');
fprintf(fid,'*radpth ..\\rad                    # path to radiance files\n');
fprintf(fid,'*imgpth ..\\images                 # path to project images\n');
fprintf(fid,'*tmppth ..\\temp                   # path to project scratch folder\n');
fprintf(fid,'*docpth ..\\doc                    # path to project documents\n');
fprintf(fid,'*dbspth ..\\dbs                    # path to local databases\n');
fprintf(fid,'*indx    1 # Building only\n');
fprintf(fid,' 43.66  -9.400   # Latitude & Longitude (diff from meridian)\n');
fprintf(fid,'      1   0.200   # Site exposure & ground reflectivity\n');
fprintf(fid,'* DATABASES\n');
fprintf(fid,'*prm  C:/Esru/esp-r/databases/constr.db2\n');
fprintf(fid,'*mlc  C:/Esru/esp-r/databases/.Alouette\n');
fprintf(fid,'*opt  C:/Esru/esp-r/databases/optics.db2\n');
fprintf(fid,'*prs  C:/Esru/esp-r/databases/pressc.db1\n');
fprintf(fid,'*evn  C:/Esru/esp-r/databases/profiles.db1\n');
fprintf(fid,'*clm  C:/Esru/esp-r/climate/Toronto_SDD\n');
fprintf(fid,'*pdb  C:/Esru/esp-r/databases/plantc.db1\n');
fprintf(fid,'*slr_half_hr   0  # solar timing hour centred\n');
fprintf(fid,'*ctl  ..\\ctl\\ctlfile.ctl\n');
fprintf(fid,'*calename standard weekday Sat & Sun\n');
fprintf(fid,'*calentag weekdays,weekdays (all year),261\n');
fprintf(fid,'*calentag saturday,Saturdays (all year), 52\n');
fprintf(fid,'*calentag sunday,Sundays (all year), 52\n');
fprintf(fid,'*list\n');
fprintf(fid,' 1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,\n'); 
fprintf(fid,' 2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,\n'); 
fprintf(fid,' 1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,\n'); 
fprintf(fid,' 1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,\n'); 
fprintf(fid,' 3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,\n'); 
fprintf(fid,' 1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,\n'); 
fprintf(fid,' 1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,\n'); 
fprintf(fid,' 1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1,1,1,1,1,2,3,1\n');
fprintf(fid,'*end_list\n');
fprintf(fid,'*year  2001 # assessment year\n');
fprintf(fid,'# sim setup: no. sets startup zone_ts plant_ts save_lv @ each ts\n');
fprintf(fid,'*sps    1   4   1  10   3   0\n');
fprintf(fid,'   1   1  31  12  std  # period & name\n');
fprintf(fid,'*sblr %s.res\n', b);
fprintf(fid,'*end_set\n');
fprintf(fid,'*end_sps\n');
fprintf(fid,'* PROJ LOG\n');
fprintf(fid,'logfile.log\n');
fprintf(fid,'* Building\n');
fprintf(fid,'model for paper\n');
fprintf(fid,'      3  # no of zones\n');
fprintf(fid,'*zon   1   # reference for South\n');       
fprintf(fid,'*opr ..\\zones\\South.opr  # schedules\n');
fprintf(fid,'*geo ..\\zones\\South.geo  # geometry\n');
fprintf(fid,'*con ..\\zones\\South.con  # construction\n');
fprintf(fid,'*tmc ..\\zones/South.tmc  # transparent constr\n');
fprintf(fid,'*zend\n'); 
fprintf(fid,'*zon   2   # reference for North\n');       
fprintf(fid,'*opr ..\\zones\\North.opr  # schedules\n');
fprintf(fid,'*geo ..\\zones\\North.geo  # geometry\n');
fprintf(fid,'*con ..\\zones/North.con  # construction\n');
fprintf(fid,'*zend\n'); 
fprintf(fid,'*zon   3   # reference for Basement\n');    
fprintf(fid,'*opr ..\\zones\\Basement.opr  # schedules\n');
fprintf(fid,'*geo ..\\zones\\Basement.geo  # geometry\n');
fprintf(fid,'*con ..\\zones\\Basement.con  # construction\n');
fprintf(fid,'*zend\n'); 
fprintf(fid,'*cnn  connfile.cnn  # connections\n'); 
fprintf(fid,'      0   # no fluid flow network\n');
fclose(fid);

%South Zone
fid=fopen([a '\zones\South.geo'], 'wt');
fprintf(fid,'# geometry of South defined in: ..\\zones\\South.geo\n');
fprintf(fid,'GEN  South  South Zone  # type, name, descr\n');
fprintf(fid,'      18       9   0.000    # vertices, surfaces, rotation angle\n');
fprintf(fid,'#  X co-ord, Y co-ord, Z co-ord\n');
fprintf(fid,'      0.00000     0.00000     0.00000  # vert   1\n');
fprintf(fid,'     10.00000     0.00000     0.00000  # vert   2\n');
fprintf(fid,'     10.00000     5.00000     0.00000  # vert   3\n');
fprintf(fid,'      0.00000     5.00000     0.00000  # vert   4\n');
fprintf(fid,'      0.00000     0.00000     5.00000  # vert   5\n');
fprintf(fid,'     10.00000     0.00000     5.00000  # vert   6\n');
fprintf(fid,'     10.00000     5.00000     5.00000  # vert   7\n');
fprintf(fid,'      0.00000     5.00000     5.00000  # vert   8\n');
fprintf(fid,'      2.00000     0.00000     0.00000  # vert   9\n');
fprintf(fid,'      2.00000     0.00000     5.00000  # vert  10\n');
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  11\n',x1glazing, y1glazing);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  12\n',x2glazing, y1glazing);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  13\n',x2glazing, y2glazing);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  14\n',x1glazing, y2glazing);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  15\n',x1frame, y1frame);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  16\n',x2frame, y1frame);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  17\n',x2frame, y2frame);
fprintf(fid,'      %-11.5f 0.00000     %0.5f  # vert  18\n',x1frame, y2frame);
fprintf(fid,'# no of vertices followed by list of associated vert\n');
fprintf(fid,'   4,  2,  3,  7,  6,\n');
fprintf(fid,'   4,  3,  4,  8,  7,\n');
fprintf(fid,'   4,  4,  1,  5,  8,\n');
fprintf(fid,'   5,  5, 10,  6,  7,  8,\n');
fprintf(fid,'   5,  1,  4,  3,  2,  9,\n');
fprintf(fid,'  10,  1,  9, 10,  5,  1, 15, 18, 17, 16, 15,\n');
fprintf(fid,'  10,  9,  2,  6, 10,  9, 11, 14, 13, 12, 11,\n');
fprintf(fid,'   4, 11, 12, 13, 14,\n');
fprintf(fid,'   4, 15, 16, 17, 18,\n');
fprintf(fid,'# unused index\n');
fprintf(fid,' 0,0,0,0,0,0,0,0,0\n');
fprintf(fid,'# surfaces indentation (m)\n');
fprintf(fid,' 0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00\n');
fprintf(fid,'    3   0   0   0    # default insolation distribution\n');
fprintf(fid,'# surface attributes follow: \n');
fprintf(fid,'# id  surface      geom  loc/  construction environment\n');
fprintf(fid,'# no  name         type  posn  name         other side\n');
fprintf(fid,'  1, Surf-2        OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  2, Surf-3        OPAQ  VERT  %s North   \n',int_wall);       
fprintf(fid,'  3, Surf-4        OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  4, Surf-5        OPAQ  CEIL  Top_ceiling_ EXTERIOR\n');       
fprintf(fid,'  5, Surf-6        OPAQ  FLOR  %s Basement\n',floor);       
fprintf(fid,'  6, Frame_Wall    OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  7, Glazing_Wall  OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  8, Glazing       TRAN  VERT  %s EXTERIOR\n',glz_name);       
fprintf(fid,'  9, Frame         OPAQ  VERT  Frame_window EXTERIOR\n');       
fprintf(fid,'# base\n');
fprintf(fid,'  5  0  0  0  0  0    50.00 0\n');
fclose(fid);

%South construction
fid=fopen([a '\zones\South.con'], 'wt');
fprintf(fid,'# thermophysical properties of South defined in ..\\zones\\South.con\n');
fprintf(fid,'# no of |air |surface(from geo)| multilayer construction\n');
fprintf(fid,'# layers|gaps|  no.  name      | database name \n');
fprintf(fid,'     4,     0  #  1 Surf-2       AExteriorWal \n');
fprintf(fid,'     %i,     %i  #  2 Surf-3       %s \n',int_wall_layers,int_wall_air_gap, int_wall);
fprintf(fid,'     4,     0  #  3 Surf-4       AExteriorWal \n');
fprintf(fid,'     5,     0  #  4 Surf-5       Top_ceiling_ \n');
fprintf(fid,'     %i,     1  #  5 Surf-6       %s  \n',floor_layers,floor);
fprintf(fid,'     4,     0  #  6 Frame_Wall   AExteriorWal \n');
fprintf(fid,'     4,     0  #  7 Glazing_Wall AExteriorWal \n');
fprintf(fid,'     %i,     %i  #  8 Glazing      %s    \n',glazing_layers,num_air_gaps,glz_name);
fprintf(fid,'     1,     0  #  9 Frame        Frame_window \n');
switch construction 
    case 'L'  
        fprintf(fid,'# air gap position & resistance for surface  2\n');
        fprintf(fid,'  2,   0.170, \n');
    case 'M'
        fprintf(fid,'# air gap position & resistance for surface  2\n');
        fprintf(fid,'  2,   0.170, \n');
    case 'H'
        l=1;
end
fprintf(fid,'# air gap position & resistance for surface  5\n');
fprintf(fid,'  2,   0.170,\n');  %The same regardless of floor type
switch glz
    case 'D'
fprintf(fid,'# air gap position & resistance for surface  8\n');
fprintf(fid,'  2,   0.440, \n');
    case 'T'
fprintf(fid,'# air gap position & resistance for surface  8\n');
fprintf(fid,'  2,   0.220, 4,   0.210,\n');
end
fprintf(fid,'# conduc-  |  density | specific | thick- |dpnd|  ref. |  temp. |moisture| surf|lyr\n');
fprintf(fid,'# tivity   |          | heat     |ness(mm)|type|  temp | factor | factor |     |  \n');
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  1  1\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
switch construction
    case 'L'
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #  2  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.1000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     3\n');
    case 'M'
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #  2  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.1000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     3\n');
    case 'H'
fprintf(fid,'     1.0600,    1950.0,    1000.0,  0.2032,   0,   0.00, 0.00000, 0.00000  #  2  1\n');
end
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  3  1\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0020,   0,   0.00, 0.00000, 0.00000  #  4  1\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0328,       1.0,       1.0,  0.2790,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
switch construction
    case 'L'
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
    case 'M'
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.0580,   0,   0.00, 0.00000, 0.00000  #     4\n');
    case 'H'
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2030,   0,   0.00, 0.00000, 0.00000  #     4\n');
end
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  6  1\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  7  1\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
switch glz
    case 'D'
fprintf(fid,'     0.7600,    2710.0,     837.0,  0.0030,   0,   0.00, 0.00000, 0.00000  #  8  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0120,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.7600,    2710.0,     837.0,  0.0030,   0,   0.00, 0.00000, 0.00000  #     3\n');
    case 'T'
fprintf(fid,'     0.7600,    2710.0,     837.0,  0.0030,   0,   0.00, 0.00000, 0.00000  #  8  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0120,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.7600,    2710.0,     837.0,  0.0030,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0120,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.7600,    2710.0,     837.0,  0.0030,   0,   0.00, 0.00000, 0.00000  #     5\n');
end
fprintf(fid,'     0.1000,     800.0,     900.0,  0.0200,   0,   0.00, 0.00000, 0.00000  #  9  1\n');
fprintf(fid,'# for each surface: inside face emissivity\n');
fprintf(fid,' 0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.83000,0.90000\n');
fprintf(fid,'# for each surface: outside face emissivity\n');
fprintf(fid,' 0.90000,0.90000,0.90000,0.12000,0.90000,0.90000,0.90000,0.83000,0.90000\n');
fprintf(fid,'# for each surface: inside face solar absorptivity\n');
fprintf(fid,' 0.50000,0.40000,0.50000,0.50000,0.65000,0.50000,0.50000,0.0500000,0.20000\n');
fprintf(fid,'# for each surface: outside face solar absorptivity\n');
fprintf(fid,' 0.65000,0.40000,0.65000,0.20000,0.65000,0.65000,0.65000,0.0500000,0.20000\n');
fclose(fid);
 
%North Zone
fid=fopen([a '\zones\North.geo'], 'wt');
fprintf(fid,'# geometry of North defined in: ..\\zones\\North.geo\n');
fprintf(fid,'GEN  North  North describes a...  # type, name, descr\n');
fprintf(fid,'       8       6   0.000    # vertices, surfaces, rotation angle\n');
fprintf(fid,'#  X co-ord, Y co-ord, Z co-ord\n');
fprintf(fid,'      0.00000     5.00000     0.00000  # vert   1\n');
fprintf(fid,'     10.00000     5.00000     0.00000  # vert   2\n');
fprintf(fid,'     10.00000    10.00000     0.00000  # vert   3\n');
fprintf(fid,'      0.00000    10.00000     0.00000  # vert   4\n');
fprintf(fid,'      0.00000     5.00000     5.00000  # vert   5\n');
fprintf(fid,'     10.00000     5.00000     5.00000  # vert   6\n');
fprintf(fid,'     10.00000    10.00000     5.00000  # vert   7\n');
fprintf(fid,'      0.00000    10.00000     5.00000  # vert   8\n');
fprintf(fid,'# no of vertices followed by list of associated vert\n');
fprintf(fid,'   4,  1,  2,  6,  5,\n');
fprintf(fid,'   4,  2,  3,  7,  6,\n');
fprintf(fid,'   4,  3,  4,  8,  7,\n');
fprintf(fid,'   4,  4,  1,  5,  8,\n');
fprintf(fid,'   4,  5,  6,  7,  8,\n');
fprintf(fid,'   4,  1,  4,  3,  2,\n');
fprintf(fid,'# unused index\n');
fprintf(fid,' 0,0,0,0,0,0\n');
fprintf(fid,'# surfaces indentation (m)\n');
fprintf(fid,' 0.00,0.00,0.00,0.00,0.00,0.00\n');
fprintf(fid,'    3   0   0   0    # default insolation distribution\n');
fprintf(fid,'# surface attributes follow:\n'); 
fprintf(fid,'# id  surface      geom  loc/  construction environment\n');
fprintf(fid,'# no  name         type  posn  name         other side\n');
fprintf(fid,'  1, Surf-1        OPAQ  VERT  %s South\n', int_wall);          
fprintf(fid,'  2, Surf-2        OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  3, Surf-3        OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  4, Surf-4        OPAQ  VERT  AExteriorWal EXTERIOR\n');       
fprintf(fid,'  5, Surf-5        OPAQ  CEIL  Top_ceiling_ EXTERIOR\n');       
fprintf(fid,'  6, Surf-6        OPAQ  FLOR  AFloor_no_sl Basement\n');       
fprintf(fid,'# base\n');
fprintf(fid,'  6  0  0  0  0  0    50.00 0\n');
fclose(fid);

%North Construction
fid=fopen([a '\zones\North.con'], 'wt');
fprintf(fid,'# thermophysical properties of North defined in ..\\zones/North.con\n'); 
fprintf(fid,'# no of |air |surface(from geo)| multilayer construction\n'); 
fprintf(fid,'# layers|gaps|  no.  name      | database name \n'); 
fprintf(fid,'     %i,     %i  #  1 Surf-1       %s\n',int_wall_layers,int_wall_air_gap,int_wall);  
fprintf(fid,'     4,     0  #  2 Surf-2       AExteriorWal\n'); 
fprintf(fid,'     4,     0  #  3 Surf-3       AExteriorWal\n'); 
fprintf(fid,'     4,     0  #  4 Surf-4       AExteriorWal\n'); 
fprintf(fid,'     5,     0  #  5 Surf-5       Top_ceiling_\n'); 
fprintf(fid,'     3,     1  #  6 Surf-6       AFloor_no_sl\n');
switch construction 
    case 'L'  
        fprintf(fid,'# air gap position & resistance for surface  1\n');
        fprintf(fid,'  2,   0.170, \n');
    case 'M'
        fprintf(fid,'# air gap position & resistance for surface  1\n');
        fprintf(fid,'  2,   0.170, \n');
    case 'H'
        l=1;
end
fprintf(fid,'# air gap position & resistance for surface  6\n'); 
fprintf(fid,'  2,   0.170,\n'); 
fprintf(fid,'# conduc-  |  density | specific | thick- |dpnd|  ref. |  temp. |moisture| surf|lyr\n'); 
fprintf(fid,'# tivity   |          | heat     |ness(mm)|type|  temp | factor | factor |     |  \n'); 
switch construction
    case 'L'
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #  1  1\n'); 
fprintf(fid,'     0.0000,       0.0,       0.0,  0.1000,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
    case 'M'
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #  1  1\n'); 
fprintf(fid,'     0.0000,       0.0,       0.0,  0.1000,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
    case 'H'
fprintf(fid,'     1.0600,    1950.0,    1000.0,  0.2032,   0,   0.00, 0.00000, 0.00000  #  1  1\n');
end
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  2  1\n'); 
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n'); 
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  3  1\n'); 
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n'); 
fprintf(fid,'     0.0690,     630.0,    2200.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #  4  1\n'); 
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.0390,      79.0,     345.0,  0.1524,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n'); 
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0020,   0,   0.00, 0.00000, 0.00000  #  5  1\n'); 
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n'); 
fprintf(fid,'     0.0328,       1.0,       1.0,  0.2790,   0,   0.00, 0.00000, 0.00000  #     3\n'); 
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n'); 
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n'); 
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  6  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'# for each surface: inside face emissivity\n'); 
fprintf(fid,' 0.90000,0.90000,0.90000,0.90000,0.90000,0.90000\n'); 
fprintf(fid,'# for each surface: outside face emissivity\n'); 
fprintf(fid,' 0.90000,0.90000,0.90000,0.90000,0.12000,0.90000\n'); 
fprintf(fid,'# for each surface: inside face solar absorptivity\n'); 
fprintf(fid,' 0.40000,0.50000,0.50000,0.50000,0.50000,0.65000\n'); 
fprintf(fid,'# for each surface: outside face solar absorptivity\n'); 
fprintf(fid,' 0.40000,0.65000,0.65000,0.65000,0.20000,0.65000\n'); 
fclose(fid);

%Basement
fid=fopen([a '\zones\Basement.geo'], 'wt');
fprintf(fid,'# geometry of Basement defined in: ..\\zones\\Basement.geo\n'); 
fprintf(fid,'GEN  Basement  Basement describes a...  # type, name, descr\n'); 
fprintf(fid,'      14      11   0.000    # vertices, surfaces, rotation angle\n'); 
fprintf(fid,'#  X co-ord, Y co-ord, Z co-ord\n'); 
fprintf(fid,'      0.00000     0.00000    -2.50000  # vert   1\n'); 
fprintf(fid,'     10.00000     0.00000    -2.50000  # vert   2\n'); 
fprintf(fid,'     10.00000    10.00000    -2.50000  # vert   3\n'); 
fprintf(fid,'      0.00000    10.00000    -2.50000  # vert   4\n'); 
fprintf(fid,'      0.00000     0.00000     0.00000  # vert   5\n'); 
fprintf(fid,'     10.00000     0.00000     0.00000  # vert   6\n'); 
fprintf(fid,'     10.00000    10.00000     0.00000  # vert   7\n'); 
fprintf(fid,'      0.00000    10.00000     0.00000  # vert   8\n'); 
fprintf(fid,'      0.00000     5.00000     0.00000  # vert   9\n'); 
fprintf(fid,'     10.00000     5.00000     0.00000  # vert  10\n'); 
fprintf(fid,'      0.00000     0.00000    -1.00000  # vert  11\n'); 
fprintf(fid,'     10.00000     0.00000    -1.00000  # vert  12\n'); 
fprintf(fid,'     10.00000    10.00000    -1.00000  # vert  13\n'); 
fprintf(fid,'      0.00000    10.00000    -1.00000  # vert  14\n'); 
fprintf(fid,'# no of vertices followed by list of associated vert\n'); 
fprintf(fid,'   4,  2,  3, 13, 12,\n'); 
fprintf(fid,'   4,  3,  4, 14, 13,\n'); 
fprintf(fid,'   4,  4,  1, 11, 14,\n'); 
fprintf(fid,'   4,  1,  4,  3,  2,\n'); 
fprintf(fid,'   4,  9,  5,  6, 10,\n'); 
fprintf(fid,'   4,  8,  9, 10,  7,\n'); 
fprintf(fid,'   4,  1,  2, 12, 11,\n'); 
fprintf(fid,'   4, 11, 12,  6,  5,\n'); 
fprintf(fid,'   5, 14, 11,  5,  9,  8,\n'); 
fprintf(fid,'   4, 13, 14,  8,  7,\n'); 
fprintf(fid,'   5, 12, 13,  7, 10,  6,\n'); 
fprintf(fid,'# unused index\n'); 
fprintf(fid,' 0,0,0,0,0,0,0,0,0,0,0\n'); 
fprintf(fid,'# surfaces indentation (m)\n'); 
fprintf(fid,' 0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00\n'); 
fprintf(fid,'    3   0   0   0    # default insolation distribution\n'); 
fprintf(fid,'# surface attributes follow:\n');  
fprintf(fid,'# id  surface      geom  loc/  construction environment\n'); 
fprintf(fid,'# no  name         type  posn  name         other side\n'); 
fprintf(fid,'  1, East BG       OPAQ  VERT  ABasementBG  GROUND\n');          
fprintf(fid,'  2, North BG      OPAQ  VERT  ABasementBG  GROUND\n');          
fprintf(fid,'  3, West BG       OPAQ  VERT  ABasementBG  GROUND\n');          
fprintf(fid,'  4, Basement Flo  OPAQ  FLOR  ABasement_Fl GROUND\n');          
fprintf(fid,'  5, South Floor   OPAQ  CEIL  %s South\n', floor_bottom);           
fprintf(fid,'  6, North Floor   OPAQ  CEIL  AFloor_no_bo North\n');           
fprintf(fid,'  7, South_BG      OPAQ  VERT  ABasementBG  GROUND\n');          
fprintf(fid,'  8, South_AG      OPAQ  VERT  ABasementAG  EXTERIOR\n');        
fprintf(fid,'  9, West_AG       OPAQ  VERT  ABasementAG  EXTERIOR\n');        
fprintf(fid,' 10, North_AG      OPAQ  VERT  ABasementAG  EXTERIOR\n');        
fprintf(fid,' 11, East_AG       OPAQ  VERT  ABasementAG  EXTERIOR\n');        
fprintf(fid,'# base\n'); 
fprintf(fid,'  5  4  0  0  0  0   150.00 0\n'); 
fclose(fid);

fid=fopen([a '\zones\Basement.con'], 'wt');
fprintf(fid,'# thermophysical properties of Basement defined in ..\\zones\\Basement.con\n');
fprintf(fid,'# no of |air |surface(from geo)| multilayer construction\n');
fprintf(fid,'# layers|gaps|  no.  name      | database name \n');
fprintf(fid,'     5,     1  #  1 East BG      ABasementBG \n');
fprintf(fid,'     5,     1  #  2 North BG     ABasementBG \n');
fprintf(fid,'     5,     1  #  3 West BG      ABasementBG \n');
fprintf(fid,'     2,     0  #  4 Basement Flo ABasement_Fl\n');
fprintf(fid,'     %i,     1  #  5 South Floor  %s\n',floor_layers,floor_bottom);
fprintf(fid,'     3,     1  #  6 North Floor  AFloor_no_bo\n');
fprintf(fid,'     5,     1  #  7 South_BG     ABasementBG \n');
fprintf(fid,'     5,     0  #  8 South_AG     ABasementAG \n');
fprintf(fid,'     5,     0  #  9 West_AG      ABasementAG \n');
fprintf(fid,'     5,     0  # 10 North_AG     ABasementAG \n');
fprintf(fid,'     5,     0  # 11 East_AG      ABasementAG \n');
fprintf(fid,'# air gap position & resistance for surface  1\n');
fprintf(fid,'  4,   0.170,\n');
fprintf(fid,'# air gap position & resistance for surface  2\n');
fprintf(fid,'  4,   0.170,\n');
fprintf(fid,'# air gap position & resistance for surface  3\n');
fprintf(fid,'  4,   0.170,\n');
fprintf(fid,'# air gap position & resistance for surface  5\n');
switch construction
    case 'L'
        inverse_gap_location=2;
    case 'M'
        inverse_gap_location=3;
    case 'H'
        inverse_gap_location=3;
end
fprintf(fid,'  %i,   0.170,\n',inverse_gap_location);  %gap location is opposite that of the south construction.  If it's 3 for south, it's 2 for basement.  So, layers+1-location
fprintf(fid,'# air gap position & resistance for surface  6\n');
fprintf(fid,'  2,   0.170,\n');
fprintf(fid,'# air gap position & resistance for surface  7\n');
fprintf(fid,'  4,   0.170,\n');
fprintf(fid,'# conduc-  |  density | specific | thick- |dpnd|  ref. |  temp. |moisture| surf|lyr\n');
fprintf(fid,'# tivity   |          | heat     |ness(mm)|type|  temp | factor | factor |     |  \n');
fprintf(fid,'     0.0240,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #  1  1\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'     0.0240,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #  2  1\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'     0.0240,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #  3  1\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'     0.0350,       1.0,       1.0,  0.0470,   0,   0.00, 0.00000, 0.00000  #  4  1\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #     2\n');
switch construction
    case 'L'
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
    case 'M'
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.0580,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     4\n');
    case 'H'
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2030,   0,   0.00, 0.00000, 0.00000  #  5  1\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     4\n');
end
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #  6  1\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.3000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.1500,     700.0,    1420.0,  0.0160,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0240,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #  7  1\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0000,       0.0,       0.0,  0.0508,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0010,   0,   0.00, 0.00000, 0.00000  #  8  1\n');
fprintf(fid,'     0.0360,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0010,   0,   0.00, 0.00000, 0.00000  #  9  1\n');
fprintf(fid,'     0.0360,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0010,   0,   0.00, 0.00000, 0.00000  # 10  1\n');
fprintf(fid,'     0.0360,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'    50.0000,    7800.0,     502.0,  0.0010,   0,   0.00, 0.00000, 0.00000  # 11  1\n');
fprintf(fid,'     0.0360,       1.0,       1.0,  0.0254,   0,   0.00, 0.00000, 0.00000  #     2\n');
fprintf(fid,'     1.6930,    2200.0,     800.0,  0.2000,   0,   0.00, 0.00000, 0.00000  #     3\n');
fprintf(fid,'     0.0236,       1.0,       1.0,  0.0762,   0,   0.00, 0.00000, 0.00000  #     4\n');
fprintf(fid,'     0.1610,     800.0,    1100.0,  0.0127,   0,   0.00, 0.00000, 0.00000  #     5\n');
fprintf(fid,'# for each surface: inside face emissivity\n');
fprintf(fid,' 0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000\n');
fprintf(fid,'# for each surface: outside face emissivity\n');
fprintf(fid,' 0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.90000,0.12000,0.12000,0.12000,0.12000\n');
fprintf(fid,'# for each surface: inside face solar absorptivity\n');
fprintf(fid,' 0.50000,0.50000,0.50000,0.65000,0.65000,0.65000,0.50000,0.50000,0.50000,0.50000,0.50000\n');
fprintf(fid,'# for each surface: outside face solar absorptivity\n');
fprintf(fid,' 0.50000,0.50000,0.50000,0.65000,0.65000,0.65000,0.50000,0.20000,0.20000,0.20000,0.20000\n');
fclose(fid);

%Optical Properties
fid=fopen([a '\zones\South.tmc'], 'wt');
fprintf(fid,'# transparent properties of South defined in ..\\zones/South.tmc\n');
fprintf(fid,'   9   # surfaces\n');
fprintf(fid,'# tmc index for each surface\n');
fprintf(fid,'0,0,0,0,0,0,0,1,0\n');
switch glz
    case 'D'
fprintf(fid,'   3  DblGlazLoAr  # layers in tmc type 1\n');
fprintf(fid,'# Transmission @ 5 angles & visible tr.\n');
fprintf(fid,'   0.586   0.566   0.513   0.363   0.177   0.740\n');
fprintf(fid,'# For each layer absorption @ 5 angles\n');
fprintf(fid,'   0.100   0.110   0.131   0.133   0.000\n');
fprintf(fid,'   0.001   0.001   0.001   0.001   0.001\n');
fprintf(fid,'   0.151   0.151   0.143   0.118   0.068   \n');   
    case 'T'
fprintf(fid,'   5  Trip_glaz_3  # layers in tmc type 1\n');
fprintf(fid,'# Transmission @ 5 angles & visible tr.\n');
fprintf(fid,'   0.598   0.568   0.490   0.315   0.124   0.740\n');
fprintf(fid,'# For each layer absorption @ 5 angles\n');
fprintf(fid,'   0.100   0.110   0.120   0.135   0.136\n');
fprintf(fid,'   0.001   0.001   0.001   0.001   0.001\n');
fprintf(fid,'   0.077   0.082   0.086   0.083   0.066\n');
fprintf(fid,'   0.001   0.001   0.001   0.001   0.001\n');
fprintf(fid,'   0.058   0.060   0.058   0.046   0.026\n');
end
fprintf(fid,'   0  # optical control flag\n');
fclose(fid);

%connections file
fid=fopen([a '\cfg\connfile.cnn'], 'wt');
fprintf(fid,'*connections  for %s\n',b);
fprintf(fid,'*Date Wed Mar 26 19:32:07 2008\n');
fprintf(fid,'     26   # number of connections\n');
fprintf(fid,'   1   1   0   0   0  #   1  Surf-2 in South is External\n');
fprintf(fid,'   1   2   3   2   1  #   2  Surf-3 in South >|< Surf-1 in North\n');
fprintf(fid,'   1   3   0   0   0  #   3  Surf-4 in South is External\n');
fprintf(fid,'   1   4   0   0   0  #   4  Surf-5 in South is External\n');
fprintf(fid,'   1   5   3   3   5  #   5  Surf-6 in South >|< South Floor in Basement\n');
fprintf(fid,'   1   6   0   0   0  #   6  Frame_Wall in South is External\n');
fprintf(fid,'   1   7   0   0   0  #   7  Glazing_Wall in South is External\n');
fprintf(fid,'   1   8   0   0   0  #   8  Glazing in South is External\n');
fprintf(fid,'   1   9   0   0   0  #   9  Frame in South is External\n');
fprintf(fid,'   2   1   3   1   2  #  10  Surf-1 in North >|< Surf-3 in South\n');
fprintf(fid,'   2   2   0   0   0  #  11  Surf-2 in North is External\n');
fprintf(fid,'   2   3   0   0   0  #  12  Surf-3 in North is External\n');
fprintf(fid,'   2   4   0   0   0  #  13  Surf-4 in North is External\n');
fprintf(fid,'   2   5   0   0   0  #  14  Surf-5 in North is External\n');
fprintf(fid,'   2   6   3   3   6  #  15  Surf-6 in North >|< North Floor in Basement\n');
fprintf(fid,'   3   1   4   1   0  #  16  East BG in Basement >|< ground profile 1\n');
fprintf(fid,'   3   2   4   1   0  #  17  North BG in Basement >|< ground profile 1\n');
fprintf(fid,'   3   3   4   1   0  #  18  West BG in Basement >|< ground profile 1\n');
fprintf(fid,'   3   4   4   1   0  #  19  Basement Flo in Basement >|< ground profile 1\n');
fprintf(fid,'   3   5   3   1   5  #  20  South Floor in Basement >|< Surf-6 in South\n');
fprintf(fid,'   3   6   3   2   6  #  21  North Floor in Basement >|< Surf-6 in North\n');
fprintf(fid,'   3   7   4   1   0  #  22  South_BG in Basement >|< ground profile 1\n');
fprintf(fid,'   3   8   0   0   0  #  23  South_AG in Basement is External\n');
fprintf(fid,'   3   9   0   0   0  #  24  West_AG in Basement is External\n');
fprintf(fid,'   3  10   0   0   0  #  25  North_AG in Basement is External\n');
fprintf(fid,'   3  11   0   0   0  #  26  East_AG in Basement is External\n');
fclose(fid);    

dos(['bps -file c:\esru\models\sdd\' b '\cfg\' b '.cfg -mode text -p std silent']);  %run simulation with above input files
%['bps -file c:\esru\models\sdd\' b '\cfg\' b '.cfg -mode text -p std silent']
%movefile(['c:\Matlab7\work\' b '.res'],['c:\esru\models\sdd\' b]);
copyfile(['c:\esru\models\sdd\' b],['g:\sdd\' b]);
movefile(['c:\Matlab7\work\' b '.res'],['g:\sdd\' b]);
%copyfile(['c:\esru\models\sdd\' b],['g:\sdd\' b]);
        end
    end
end


%%%%%THIS WORKS --- DO NOT TOUCH%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%dos('bps -file c:\esru\models\sdd\hd15\cfg\hd15.cfg -mode text -p std silent');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%