| Ce programme montre comment à partir d'un profil vertical
de salinité et de température, on peut obtenir
un profil de salinité et de température potentielle.
Ce type de programme peut être utile lorsque l'on souhaite
initialiser un modèle aux équations primitives.
En effet, la variable d'état du modèle est la
température potentielle. IL faut donc convertir la
température en température potentielle. Malheureusement,
les banques de données se contentent parfois de fournir
la température sur des niveaux z (cas de la banque
de données MEDATLAS). Le programme que je vous propose
est constitué d'un programme principal réduit
à la lecture d'un banal profil (1D z,T,S). Les sousprogrammes
associés permettent de faire une estimation à
priori de la pression (supposant une densité constante),
puis une estimation de la température potentielle,
puis une estimation de la densité potentielle et enfin
une nouvelle estimation de la densité. Cette dernière
étant plus précise que la densité initiale
à priori, on l'utilise pour refaire le calcul, jusqu'à
convergence du résulat final.
|
| |
| En pratique, télécharger
pack_eqstate.tar.gz, decompresser (gunzip) et redéployer
(tar xvf) les fichiers. |
| |
On obtient:
le programme: init_with_ts1dv.F
Deux fichiers d'entrée 1DV pour pouvoir faire une application
à titre d'exemple: S1DV_winter.ascii T1DV_winter.ascii
dimension.h: contenant les
dimensions spatiales du problème. Ici:
INTEGER IMAX,JMAX,KMAX
PARAMETER(IMAX=1,JMAX=1,KMAX=22) |
| |
L'exemple se contente de traiter un profil unique, donc les
dimensions horizontales IMAX et JMAX sont réduites à
1, KMAX=22 étant la dimension verticale, soit le nombre
de niveaux verticaux des fichiers d'entrée S1DV_winter.ascii
et T1DV_winter.ascii.
Le programme est trés facilement transposable à
un véritable cas 3D. |
| |
| La compilation est simplement: g77
init_with_ts1dv.F |
| Exécutable: a.out |
| Pour cet exemple, le résultat est un fichier ascii
(file.out) fournissant,
en plus des informations de départ, la pression, la température
potentielle, la densité potentielle et la densité. |
|
Toulouse 14 mars 2005
|
|