Convertir la température en température potentielle.
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



 

 

 

 

POC (Pôle Océan & Couplages) - Observatoire Midi-Pyrénées - 14 Avenue Edouard Belin - 31400 Toulouse - France