Commit af19620a authored by Kévin Duquette's avatar Kévin Duquette

Add all files

parent d68836a9
@article{Wenegrat2017,
author = {Wenegrat, Jacobo O. and Thomas, Leif N.},
doi = {10.1175/JPO-D-16-0239.1},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/wenegrat(2017){\_}Ekman Transport in Balanced Currents with Curvature.pdf:pdf},
number = {1965},
pages = {1189--1203},
title = {{Ekman Transport in Balanced Currents with Curvature}},
year = {2017}
}
@article{Thomas2005,
author = {Thomas, Leif N. and Craig, M. Lee},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/thomas2005{\_}Intensification of Ocean Fronts by Down-Front Winds.pdf:pdf},
number = {1972},
pages = {1086--1102},
title = {{Intensification of Ocean Fronts by Down-Front Winds}},
year = {2005}
}
@article{Thomas2002,
author = {Thomas, Leif N. and Rhines, Peter B.},
doi = {10.1017/S0022112002002367},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/thomas2002{\_}Nonlinear stratified spin-up.pdf:pdf},
pages = {211--244},
title = {{Nonlinear stratified spin-up}},
volume = {473},
year = {2002}
}
@article{Tan2001,
author = {Tan, Zhe-min},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/Tan2001{\_}Article{\_}AnApproximateAnalyticalSolutio.pdf:pdf},
keywords = {baroclinity,boundary-layer winds,ekman pumping,semi-geostrophic dynamics,variable eddy viscosity},
number = {1994},
pages = {361--385},
title = {{An approximate analytical solution for the baroclinic and variable eddy diffusivity semi-geostrophic ekman boundary layer}},
year = {2001}
}
@article{,
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/R{\'{e}}sum{\'{e}}{\_}lecture.pdf:pdf},
number = {1994},
title = {{R{\'{e}}sum{\'{e}} litt{\'{e}}rature sur l ' asym{\'{e}}trie d ' Ekman A Comprehensive Analysis of Nonlinear Corrections to the Classical Ekman Pumping}},
year = {2006}
}
@article{Polvani1994,
author = {Polvani, L M and Mcwilliams, J C and Spall, M A and Ford, R},
doi = {10.1063/1.166002},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/Polvani et al(1994). The coherent structures of shallowwater turbulence$\backslash$: Deformationradius effects, cyclone$\backslash$:anticyclone asymmetry and gravitywave generation.pdf:pdf},
title = {{The coherent structures of shallowwater turbulence : Deformationradius effects , cyclone / anticyclone asymmetry and gravitywave generation The coherent structures of shallow-water turbulence : Deformation-radius effects , cyclone / anticyclone asymmetry and gravity-wave generation}},
volume = {177},
year = {1994}
}
@article{Mahadevan2006,
author = {Mahadevan, Amala and Tandon, Amit},
doi = {10.1016/j.ocemod.2006.05.006},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/mahadevan2006{\_}An analysis of mechanisms for submesoscale vertical motion at ocean fronts.pdf:pdf},
pages = {241--256},
title = {{An analysis of mechanisms for submesoscale vertical motion at ocean fronts}},
volume = {14},
year = {2006}
}
@article{Kuo2000,
author = {Kuo, Allen C and Polvani, Lorenzo M and Kuo, Allen C and Polvani, Lorenzo M},
doi = {10.1063/1.870363},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/kuo(2000){\_}Nonlinear geostrophic adjustment, cycone,anticyclone asymetry, and potential vorticity rearangement.pdf:pdf},
title = {{vorticity rearrangement Nonlinear geostrophic adjustment , cyclone {\~{O}} anticyclone asymmetry , and potential vorticity rearrangement}},
volume = {1087},
year = {2000}
}
@article{Ishida2006,
author = {Ishida, Shin-ichi and Iwayama, Takahiro},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/ishida2006{\_}A Comprehensive Analysis of Nonlinear Corrections to the Classical Ekman Pumping.pdf:pdf},
number = {5},
pages = {839--851},
title = {{A Comprehensive Analysis of Nonlinear Corrections to the Classical Ekman Pumping}},
volume = {84},
year = {2006}
}
@article{Hart2001,
author = {Hart, J E},
doi = {10.1063/1.870300},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/hart 2000{\_}A note on nonlinear corrections to the Ekman layer pumping velocitypdf.pdf:pdf},
number = {2000},
title = {{A note on nonlinear corrections to the Ekman layer pumping velocity}},
volume = {131},
year = {2001}
}
@article{Bannon1998,
author = {Bannon, Peter R.},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/bannon(1998){\_}A comparison of ekman pumping in approximate models of the aceelerating planetary boundary layer.pdf:pdf},
pages = {1446--1451},
title = {{A Comparison of Ekman Pumping in Approximate Models of the Accelerating Planetary Boundary Layer}},
year = {1998}
}
@article{Arai2013,
author = {Arai, Masazumi and Yamagata, Toshio and Arai, Masazumi},
doi = {10.1063/1.166001},
file = {:Users/Administrator/Documents/uqar/Project/Litterature/Assymetry/Arai et Yamagata (1994). Assymetry evolution of eddiesin rotating shallow water.pdf:pdf},
number = {1994},
title = {{Asymmetric evolution of eddies in rotating shallow water}},
volume = {163},
year = {2013}
}
% ============================= Introduction =================================
\subsection{}
\relax
\select@language{french}
\@writefile{toc}{\select@language{french}}
\@writefile{lof}{\select@language{french}}
\@writefile{lot}{\select@language{french}}
\bibstyle{te-fr2.bst}
\bibdata{assymetry.bib}
\begin{thebibliography}{0}
\providecommand{\natexlab}[1]{#1}
\providecommand{\url}[1]{\texttt{#1}}
\providecommand{\urlprefix}{URL }
\providecommand{\bibAnnoteFile}[1]{%
\IfFileExists{#1}{\begin{quotation}\noindent\textsc{Key:} #1\\
\textsc{Annotation:}\ \input{#1}\end{quotation}}{}}
\providecommand{\bibAnnote}[2]{%
\begin{quotation}\noindent\textsc{Key:} #1\\
\textsc{Annotation:}\ #2\end{quotation}}
\end{thebibliography}
This is BibTeX, Version 0.99d (TeX Live 2017)
Capacity: max_strings=100000, hash_size=100000, hash_prime=85009
The top-level auxiliary file: main.aux
The style file: te-fr2.bst.bst
I found no \citation commands---while reading file main.aux
Database file #1: assymetry.bib.bib
You can't pop an empty literal stack
while executing---line 3003 of file te-fr2.bst.bst
"" can't start a style-file command---line 3025 of file te-fr2.bst.bst
:
: 
(Error may have been on previous line)
You've used 0 entries,
2890 wiz_defined-function locations,
600 strings with 4896 characters,
and the built_in function-call counts, 41 in all, are:
= -- 0
> -- 0
< -- 0
+ -- 0
- -- 0
* -- 2
:= -- 8
add.period$ -- 0
call.type$ -- 0
change.case$ -- 0
chr.to.int$ -- 0
cite$ -- 0
duplicate$ -- 0
empty$ -- 1
format.name$ -- 0
if$ -- 1
int.to.chr$ -- 1
int.to.str$ -- 1
missing$ -- 0
newline$ -- 13
num.names$ -- 0
pop$ -- 0
preamble$ -- 1
purify$ -- 0
quote$ -- 0
skip$ -- 1
stack$ -- 0
substring$ -- 0
swap$ -- 0
text.length$ -- 0
text.prefix$ -- 0
top$ -- 0
type$ -- 0
warning$ -- 0
while$ -- 0
width$ -- 0
write$ -- 12
(There were 3 error messages)
This diff is collapsed.
This diff is collapsed.
% ---------------------------------- Package ------------------------------
% Document
\documentclass[titlepage,12pt]{article} % Document type 'article', avec une page titre
\usepackage{natbib} % Bibliographie
\usepackage[margin=2.5cm]{geometry} % Élargir les marges
\usepackage{setspace} % Gêré les espaces dans le texte
\usepackage{rotating} % Rotationer les pages en horizontales
\usepackage{comment} % Pour les commentaires
\usepackage{setspace} % Double iterligne
\usepackage{titlesec} % Space the title
% Écriture francais et langue
\usepackage[french]{babel} % Conventions canadien-francais
\usepackage[utf8]{inputenc} % UFT8, supporte les caractères accentués
\usepackage[T1]{fontenc} % Le PDF supporte les caractères accentués
\usepackage{csquotes} % Pour les guillemets
% Math et sciences
\usepackage{amsmath, amsfonts, amsthm, bm} % Tous les packages de math
\usepackage{chemist} % Pour les formules chimiques
\usepackage{listings} % Pour inclure le code source en annexe
\usepackage{color} % Pour les couleurs dans le code source
% Graphiques,tableau et figures
\usepackage{tikz} % Image en tikz
\usepackage{float} % Gêré les floats (\FLOATBARRIER)
\usepackage{graphicx}
\usepackage[section]{placeins} % Géré la position des graphiques précisement
\usepackage{multirow} % Tableau multiple
\usepackage{tabulary} % Tableau avec case en Y auto
% Package de pdf et output
\usepackage{pgfplots}
\usepackage{pgfgantt}
\pgfplotsset{compat=newest}
\pgfplotsset{plot coordinates/math parser=false}
\usepackage{pdflscape} % Mettre la page vertical
\usepackage[toc,page]{appendix} % Appendix
\usepackage{booktabs}
%\usepackage{arydshln}
% ------------------------------- Setting special ------------------------------
% On set le titre annexe
%Indique le nom de la tables des annexes dans la toc
\renewcommand{\appendixtocname}{Tables des Annexes}
\renewcommand{\appendixpagename}{Annexes} %Nom du titre de la page des annexes
% On set les spacing
\onehalfspacing
\setlength{\textfloatsep}{5pt}
% Pour tabulary
\newcolumntype{P}[1]{>{\centering\arraybackslash}p{#1}}
\setlength{\parskip}{0.76em}
\titlespacing{\section}{2pt}{2pt plus 2pt minus 2pt}{0pt plus 1pt minus 1pt}
\titlespacing{\subsection}{2pt}{0pt plus 1pt minus 1pt}{0pt plus 1pt minus 1pt}
%\titlespacing{\subsubsection}{0pt}{\parskip}{\parskip}
\ No newline at end of file
This diff is collapsed.
File added
This diff is collapsed.
MUDPACK Version 5.0.1
----------------------
Multigrid Software for Solving Elliptic Partial Differential Equations
12/06/2011
Author: John Adams
CAVEAT
------
MUDPACK 5.0.1 is written in Fortran 77 and 90 but it may not fully comply
with either Fortran standard. We have not run the code through a rigorous
standards checker, other than the various compilers we have used in building
the code. Users whose applications require strict adherence to the standard
must provide this assurance themselves.
Documentation
-------------
Documentation for this software package is provided in file MUDPACK5.0.1.html
and companion image file Stackedgrid.gif. We intend the document for browsing
on your local computer.
Recent bugfixes
---------------
None. This library is regarded as legacy.
Compiling the Library and Test Programs
---------------------------------------
Our Makefile requires the following: a Unix or Linux or Mac operating system, gmake,
and a Fortran 90 compiler. If successful, this Makefile builds a static library
libmudpack.a and several binary executables.
Examine file make.inc to see if your OS and compiler are represented. If they
are not, you should modify file make.inc and the Makefile in each directory
(main, src and test) so that they are represented.
If you desire double precision floating-point arithmetic, you should modify Makefile
and make.inc with compiler options for promoting single to double precision. Your
compiler's user guide will have information on the required compiler options.
This diff is collapsed.
program main
implicit none
integer mainVar
integer para
parameter (para=3)
real mat1(10,10), mat2(10,10)
mat1(:,:) = 0
mat2 = mat1
!mainVar =1
!
! call sub1(mainVar)
write(*,*) 'mainVar in Main' ,mainVar
write(*,*) 'para in Main', para
end program main
subroutine sub1(mainVar)
implicit none
integer mainVar, para
real var1
parameter (var1=3.1416)
write(*,*) 'para in sub1 ', para
write(*,*) 'mainVar In sub1', mainVar
call sub2(mainVar)
end
subroutine sub2(mainVar)
implicit none
integer mainVar, para
real var1
mainVar = var1
write(*,*) 'mainVar in sub2', mainVar
write(*,*) 'para in sub2 ', para
end
include ../make.inc
SRC= main.f90
#tcud2.f tcud24.f tcud24cr.f tcud24sp.f tcud2cr.f tcud2sp.f tcud3.f \
# tcud34.f tcud34sp.f tcud3cr.f tcud3sp.f tcuh2.f tcuh24.f tcuh24cr.f \
# tcuh2cr.f tcuh3.f tcuh34.f tmud2.f tmud24.f tmud24cr.f tmud24sp.f \
# tmud2cr.f tmud2sa.f tmud2sp.f tmud3.f tmud34.f tmud34sp.f tmud3cr.f \
# tmud3sa.f tmud3sp.f tmuh2.f tmuh24.f tmuh24cr.f tmuh2cr.f tmuh3.f \
# tmuh34.f
EXES=$(subst .f90,.exe, $(SRC))
all : $(EXES)
$(EXES) : $(SRC)
%.exe : %.f90
rm -f $@
$(F90) $< -o $@ -L../lib -l mudpack
./$@
clean :
rm -f $(EXES)
This diff is collapsed.
This diff is collapsed.
%% Open
nx=3585;
ny=65;
home = pwd;
if strcmp(fgetMachin, 'mingan')
path2Sim = '/share/work/duqk0001/run/2DAh/Ro0.003';
elseif strcmp(fgetMachin, 'mac')
path2Sim = '/Users/Administrator/Documents/uqar/Leads/result/data/varRoEqAh_v2/Ro0.003';
else
error('Wrong location in genVreal.m ')
end
cd(path2Sim)
dynDiag = rdmds('dynDiag', Inf);
U = squeeze(dynDiag(:,1,:,1));
V = squeeze(dynDiag(:,1,:,2));
W = squeeze(dynDiag(:,1,:,3));
Vmean = mean(V(floor(nx/2),:));
V2Write = resizem(V, [nx ny]);
U2Write = resizem(U,[nx ny]);
cd(home)
%% Write
%fid=fopen('Vini.bin','w','b'); fwrite(fid,V2Write,'real*4'); fclose(fid);
%V
fileID = fopen('Vini.txt','w');
fprintf(fileID,[repmat('%6.5f ',1,nx) '\n'],fliplr(V2Write));
fclose(fileID);
% U
fileID = fopen('Uini.txt','w');
fprintf(fileID,[repmat('%6.5f ',1,nx) '\n'],fliplr(U2Write));
fclose(fileID);
! Initilisation
zeta(:,:) = 0.0
psi(:,:) = 0.0
v(:,:) = 0.0
rhs(:,:) = 0.0
jacoV(:,:)= 0.0
jacoZ(:,:)= 0.0
call readField(realV,'Vini.txt')
call readField(realU,'Uini.txt')
call writeField("realV.bin", realV)
call writeField("realU.bin",realU)
! set initialisation and iguess to 0
intlD = 0
iguessD = 0
intlN = 0
iguessN = 0
write(*,*) 'Initialize Dirichlet'
! Dirichlet
call mud2(iprmD,fprm,workD,cof,bndc,rhs,zeta(1:nx,:),mgopt,ierror)
call checkIerror(ierror,'Dirichlet')
write(*,*) 'Initialize Neumman'
! Neumman
call mud2(iprmN,fprm,workN,cof,bndc,rhs,v(1:nx,:),mgopt,ierror)
call checkIerror(ierror,'Neumman')
! set initialisation and iguess to 0
intlD = 1
intlN = 1
iguessD = 1
iguessN = 1
!
! -------------------------------------------------------------------
! main.f90
! This code do:
! 1) generate à initial function psi
! 2) Find the solution of psi with dirichlet condition
! 3) Find the solution of V with neuman condition
!
! In the code all equation refer of psi, RHS as
! Laplacien(psi) = rhs
!
! zeta, psi, v are the field to be compute
! stage 1,2,3 refer as computing psi,zeta,v
!
! psi = cos(2*pi*x) * sin(2*pi*y)
! with periodic boundary at xa and xb and
! diriclet bndy at yc,yd -> psi = 0
!
! V = cos(2*pi*x) * cos(pi*y)
! with periodic boundary at xa and xb and
! neuman bndy at yc,yd -> dV/dx = 0
! V(yc) = -tau/po and V(yd) = 0
!
!
! * * * * * * * * * * * * y(nx)=yd = 0
! * kbdy=4 *
! * *
! * *
! * *
! * kbdy=1 kbdy=2 *
! * *
! * *
! * *
! * kbdy=3 *
! * * * * * * * * * * * * y(1)=yc = Ly
!
! x(0)=xa= 0 x(nx)=xb = Lx
!
! see tmud2.m from the mudpack5.0.1 to have more information
!
! -------------------------------------------------------------------
module gridInfo
! Declare constante
real, parameter :: Pi = 3.1415927, nu = 0.016, f = 1.E-4, tau= 1E-3, Vmean = 0.0031
! LOW RES TEST
!integer, parameter :: nx=97, ny=49
!real, parameter :: Lx=1075, Ly=195, dx = (Lx-0) / float(nx-1), dy = (Ly-0)/float(ny-1)
! HIGH RES 10KM
integer, parameter :: nx=3585, ny=65
real, parameter :: Lx=10755, Ly=195, dx = (Lx-0) / float(nx-1), dy = (Ly-0)/float(ny-1)
end module gridInfo
program main
use gridInfo
implicit none
! Time parameter
integer niter, it
parameter (niter =3)
! Field
real readV(nx,ny)
real zeta(nx+1,ny), psi(nx+1,ny), v(nx+1,ny), rhs(nx,ny), laplace(nx+1,ny), jacoZ(nx+1,ny), jacoV(nx+1,ny), realV(nx+1,ny)
real realU(nx+1,ny)
real A, epsi, sumErr, psi0(nx+1,ny)
! MudPack Multigrid parameter
integer iixp,jjyq,iiex,jjey,nnx,nny,isx,jsy,llwork
parameter (iixp = 3 , jjyq = 3 , iiex = 6, jjey = 5) ! LOW RES
!parameter (iixp = 112 , jjyq = 2 , iiex = 6, jjey = 6) ! HIGH RES D parameter (nnx=iixp*2**(iiex-1)+1, nny=jjyq*2**(jjey-1)+1)
include 'mudPackVar.f90' ! Variable needed for mudPack package
include 'initialize.f90' ! Initialising field and discretization of mud2
!include 'rhsRealV.f90' ! Create a V field as calcul in MIT
write(*,*) '------ Starting Main.f90 ------'
! Set V ini as realV
v = realV
epsi=0.0 ! Time dependante contante Y = A + epsi(it) * B
!--------------- ITERATIVE SOLVER --------------------
do it=1,niter
!if (it.gt.1) then ! Set Guessing matrice to TRUE
iguessD = 1
iguessN = 1
!endif
!if (it>5) then ! Set epsilon as function of time
!epsi = float((it-5))/float(niter-5)
!end if
psi0=psi
print*, 'Itertation:' , it , epsi
! Executing solver for a diri bndy
print*, 'zeta'
include 'rhsZeta.f90'
if (it.eq.1) then
call writeField('zetaini.bin',zeta)
endif
! set initialisation and iguess to 0
call mud2(iprmD,fprm,workD,cof,bndc,rhs,zeta(1:nx,:),mgopt,ierror)
call checkIerror(ierror, 'zeta')
! Executing sovler for a dirichlet bndy
print*, 'phi'
include 'rhsPsi.f90'
call mud2(iprmD,fprm,workD,cof,bndc,rhs,psi(1:nx,:),mgopt,ierror)
call checkIerror(ierror, 'psi')
! Executing solver for a neumman bndy
!print*, 'v'
!include 'rhsV.f90'
!call mud2(iprmN,fprm,workN,cof,bndc,rhs,v(1:nx,:),mgopt,ierror)
!call checkIerror(ierror, 'v')
! If we wante to force V
!print*, 'V sumdiff',sum(realV-v)
v = realV
print*, 'Sum Err',sum(psi-psi0)
! print*, 'Mean V', sum(v)/(nx*ny), 'Mean Jacov', sum(jacoV)/(nx*ny)