%% Template file for running NGspice
%  Designed to work with CPPsim and the Sue2 editor.
%  This Matlab script can be placed in any directory.
%  The script is divided into three sections:
%     Section 1: Input the location of your CppSim directory, library name,
%                and schematic name.
%     Section 2: Input the simulation and any device parameters.
%     Section 3: Load simulation results and analyze data.

%% Copyright (c) 2021 by David E. Kotecki. All rights reserved.

% Redistribution and use in source and binary forms, with or without
% modification, are permitted provided that the following conditions are
% met:

% 1. Redistributions of source code must retain the above copyright notice,
% this list of conditions and the following disclaimer.

% 2. Redistributions in binary form must reproduce the above copyright
% notice, this list of conditions and the following disclaimer in the
% documentation and/or other materials provided with the distribution.

% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
% IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
% THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
% PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
% CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
% EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
% PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
% PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
% LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
% NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
% SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

%% Section 1: Define: CPPSim location, library, and schematic
clear variables;  
CppSim_Location = sprintf('/Users/Kotecki/CppSim'); % location of CppSim directory
Design_Library = sprintf('ECE444_2021'); % name of design library
Schematic_Name = sprintf('PSRR_test_tran'); % name of schematic

%% Section 2: Generate HSPC file and run NGspice
addpath(sprintf('%s/CppSimShared/HspiceToolbox', CppSim_Location)); % add ngspice matlab toolbox to the path
Working_Dir = sprintf('%s/SimRuns/%s/%s', CppSim_Location, Design_Library, Schematic_Name);
if ~exist(Working_Dir, 'dir')  
    mkdir(Working_Dir)  % create working directory if it does not exist
end
cd(Working_Dir) % set current folder to the working directory
hspc_filename = sprintf('%s.hspc', Schematic_Name);   % define hspc filename

hspcfile = fopen(hspc_filename, 'w');    % open file 'hspc_filename' for write
fprintf(hspcfile, '**** NGspice HSPC file **** \n');
fprintf(hspcfile, '**** File: %s/%s **** \n', pwd, hspc_filename);
fprintf(hspcfile, '**** Date: %s **** \n\n', datestr(datetime('now')));

fprintf(hspcfile, '**** Simulation Statement ****\n');
fprintf(hspcfile, '.tran 1e-2 30 0 0 \n\n');
%fprintf(hspcfile, '.ac dec 1000 1 1e6 \n\n');

fprintf(hspcfile, '**** Paramenter Statements ****\n');
fprintf(hspcfile, '.param r_null = 500 \n');
fprintf(hspcfile, '.param c_miller = 1e-12 \n\n');
OpAmpParameters

fprintf(hspcfile, '**** Include Statements ****\n');
fprintf(hspcfile, '.include ../../../SpiceModels/ECE214_models.mod \n');
fprintf(hspcfile, '.include ../../../SpiceModels/ECE444_models.mod \n\n');

fprintf(hspcfile, '**** Define Vdd Supply Voltage ****\n');
fprintf(hspcfile, 'Vsupply vdd gnd 3.3 \n\n');

fprintf(hspcfile, '**** Simulation Options ****\n');
fprintf(hspcfile, '.options post=1 delmax=5p relv=1e-6 reli=1e-6 relmos=1e-6 method=gear \n');
fprintf(hspcfile, '.temp 27\n');
fprintf(hspcfile, '.global gnd vdd\n');
fprintf(hspcfile, '.op \n\n');
fprintf(hspcfile, '**** End of NGspice hspc file \n');
fclose(hspcfile);
 
ngsim(hspc_filename); % run ngspice  

%% Section 3: Load simulation results and analyze data

data = loadsig('simrun.raw');  % load data from simulation

%frequency = evalsig(data, 'FREQUENCY');  % create vector of time values
time = evalsig(data, 'TIME');
Vpsrr = evalsig(data, 'vpsrr');
%Vdm = evalsig(data, 'vdm');

fs = 16;   % define font size (fs)
lw = 1.5;  % define linewidth (lw)
FigHandle = figure('Position', [200, 75, 850, 600]);  % set figure size and location

%subplot(2,1,1) % plot A+ vs. frequency
plot(time, Vpsrr, 'linewidth', 1.7)
hold on
set(gca, 'fontsize', 16);
xlabel('Time (s)', 'fontsize', 16);
ylabel('Vout (V)', 'fontsize', 16);
title('Vout as power supply voltage is varied');
%axis([1e5 1e8 -10, 10]);
grid on

return
subplot(2,1,2) % plot PSRR + vs. frequency
semilogx(frequency, 20.*log10(abs(Vdm./Vpsrr)), 'linewidth', 1.7)
grid on
set(gca, 'fontsize', 16);
xlabel('Frequency (Hz)', 'fontsize', 16);
ylabel('PSRR+ (dB)', 'fontsize', 16);
%axis([1e5 1e8 -200 200]);

return
FigHandle2 = figure('Position', [300, 75, 850, 600]);
semilogx(frequency, 20*log10(abs(Vdm) ./ abs(Vcm)), 'linewidth', 1.7)
grid on
set(gca, 'fontsize', 16);
xlabel('Frequency (Hz)', 'fontsize', 16);
ylabel('CMMR', 'fontsize', 16);
title('CMMR vs. Frequency');

return
FigHandle3 = figure('Position', [400, 75, 850, 600]);
semilogx(frequency, 20*log10(abs(Vcm)), 'linewidth', 1.7)
grid on
set(gca, 'fontsize', 16);
xlabel('Frequency (Hz)', 'fontsize', 16);
ylabel('Common Mode Gain', 'fontsize', 16);
title('Common Mode Gain vs. Frequency');


%% end of .m file