Skip to content
Snippets Groups Projects
Commit 2022078d authored by Victor Stepanenko's avatar Victor Stepanenko
Browse files

Some scripts have been improved

parent f41ed22d
No related branches found
No related tags found
No related merge requests found
......@@ -7,6 +7,14 @@ mkdir -p -v results/$1/hourly
mkdir -p -v results/$1/monthly
mkdir -p -v results/$1/daily
# Modifying driver file
sed -i '2d' driver_file.dat
sed -i "\$a setup/$1_driver.dat" driver_file.dat
# Modifying setup file
sed -i '2d' setup_file.dat
sed -i "\$a setup/$1_setup.dat" setup_file.dat
file=./setup/$1_setup.dat
if [ ! -f ${file} ];
then
......
......@@ -48,7 +48,7 @@ Short Users's guide on Linux}
\maketitle
\section{Compiling the model and creating a project}
The main prerequisite to install the LAKE model in the Fortran compiler.
The main prerequisite to install the LAKE model is the Fortran compiler.
It is recommended that you use Intel compiler, that is free under Linux
for non-commercial purposes. Intel compiler is set by default in the model install configuration (\emph{./Makefile}, \emph{./source/Makefile}).
So, if you want to use another compiler, you have to change these makefiles
......@@ -68,25 +68,9 @@ user@computer:~/LAKE/lake$ ./crproj <project_name>
\end{lstlisting}
This creates the folder structure in \emph{./results}, to which the output
files will be written. Now, you have to direct the model to files including
all parameters of the numerical experiment. To do this, modify the file
\emph{./driver\_file.dat} contents as follows
\begin{lstlisting}
# The file with parameters for driver of the Lake model
setup/<project_name>_driver.dat
\end{lstlisting}
and the file \emph{./setup\_file.dat} :
\begin{lstlisting}
# The file with parameters for driver of the Lake model
setup/<project_name>_setup.dat
\end{lstlisting}
After that, create \emph{./setup/project\_name\_driver.dat} and \emph{./setup/project\_name\_setup.dat} files. These files must have
a certain structure, so that we recommend that you simply copy them from
sample files already existing in \emph{./setup}. Once copied, you can modify them. This is what the next sections devoted to.
files will be written. After that, create \emph{./setup/project\_name\_driver.dat} and \emph{./setup/project\_name\_setup.dat} files. These files must have
a certain structure, so that we recommend that you create them by copying from
sample files already existing in \emph{./setup}. Once copied, you can modify them appropriately. This is what the next sections devoted to.
\section{Configuring the driver file}
Driver file is used to set the outer environment of the LAKE model
......
function series
rootdir='/home/victor/Files/main/models/LAKE/2014/January/lake/';
rootdir='/home/victor/Files/main/models/LAKE/2014/06/lake/';
%expname='essea1948-2011_nospinup/';
expname='Seida1_2007-2008/';
dir=strcat(rootdir,'results/',expname,'time_series/');
monthdays0 = 30.42;
work = cell(2,1);
[work{1},m,n]=getlmipascii(dir,'methane_series 1 1.dat');
......@@ -63,7 +65,7 @@ tsurfmeas(6,4) = 11.90;
timemeas = zeros(nmeas,1);
for i=1:nmeas
timemeas(i) = (tsurfmeas(i,1) - 2007)*12 + ...
tsurfmeas(i,2)-1 + tsurfmeas(i,3)/30.42;
tsurfmeas(i,2)-1 + tsurfmeas(i,3)/monthdays0; %
end
% Methane time series
......@@ -71,12 +73,14 @@ h1 = figure;
line(work{1}.data(:,ntime(1))/hr2mnth,work{1}.data(:,nsurfconc)*molm3tonM,'Color','b')
line(work{1}.data(:,ntime(1))/hr2mnth,work{1}.data(:,nbotconc)*molm3tonM,'Color','r')
legend('Surface concentration','Bottom concentration')
xlabel('Time, months'); ylabel('Concentration, nM')
h2 = figure;
line(work{1}.data(:,ntime(1))/hr2mnth,work{1}.data(:,nbotebul),'Color','g')
line(work{1}.data(:,ntime(1))/hr2mnth,work{1}.data(:,nbotdiff),'Color','r')
line(work{1}.data(:,ntime(1))/hr2mnth,work{1}.data(:,nsurfdiff),'Color','b')
legend('Bottom ebullition','Bottom diffusion','Surface diffusion')
xlabel('Time, months'); ylabel('Flux, mg/(m^2 day)')
% Temperature time series
h3 = figure;
......@@ -85,6 +89,7 @@ line(work{2}.data(:,ntime(2))/hr2mnth,work{2}.data(:,nbottemp),'Color','b')
line(timemeas,tsurfmeas(:,4),'Color','g','Marker','o')
% scatter(timemeas,tsurfmeas(:,4),'Marker','o')
legend('Surface temperature','Bottom temperature','Surface temperature measured')
xlabel('Time, months'); ylabel('Temperature, Celsius')
%line(4+work{2}.data(:,ntime(2))/hr2mnth,work{2}.data(:,nsignwaveheight),'Color','b') % 4 - starting month
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment