Changeset 142 for trunk/SRC/ToBeReviewed/CALCULS/hdyn.pro
- Timestamp:
- 07/21/06 14:47:49 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/CALCULS/hdyn.pro
r134 r142 3 3 ;------------------------------------------------------------ 4 4 ;+ 5 ; NAME:hdyn6 5 ; 7 ; PURPOSE:calcule la hauteur dynamique par rapport a un etat de 8 ; reference pour une profondeur de reference. Cf les mots cles pour 9 ; les differentes possibilites. 10 ; Par defaut l''etat de reference est 11 ; rho=1020 et la profondeur de reference est gdepw[ka] avec ka le 12 ; premier niveau W directement au dessus de 1000 m. 6 ; @file_comments 7 ; Calculate the height by rapport to a reference state for a depth reference. 8 ; See keywords for different possibilities. By default, the state reference 9 ; is rho=1020 and the depth reference is gdepw[ka] with ka the first W level 10 ; directly above 1000 m. 13 11 ; 14 ; CATEGORY: calculs de post-traitement 12 ; @categories 13 ; Post-traiment calculation 15 14 ; 16 ; CALLING SEQUENCE:res=hdyn(sn,tn) 15 ; @param TABSN {in}{required} 16 ; array representing the salinity 17 17 ; 18 ; INPUTS:sn et tn sont des tableaux de meme taille representant la19 ; salinite et la temperature.18 ; @param TABTN {in}{required} 19 ; array representing the temperature.Has the same size than SN 20 20 ; 21 ; KEYWORD PARAMETERS: 21 ; @keyword GILL We activate this key if we want to calculate the dinamic height 22 ; like in the GILL page 215, which means by rapport to a reference state which 23 ; vary in depth and which is determinated by a reference temperature tref at 0°C 24 ; and a reference salinity sref at 35psu 22 25 ; 23 ; GILL: activer cette cle si on veut faire le calcul de la 24 ; hauteur dynamique comme ds le GILL page 215 cad par rapport a 25 ; un etat de reference qui varie en profondeur et qui est 26 ; determine par une temperature de reference tref a 0 degre et 27 ; une salinite de reference sref a 35psu 26 ; @keyword LEVEL 27 ; It is the same reference level to take. This level is defined like 28 ; that gdepw[level] is the reference depth 29 ; 30 ; @keyword SREF 31 ; Give a value to this keyword to change the reference salinity used in the 32 ; calculation when GILL is activated. 33 ; 34 ; @keyword TREF 35 ; Give a value to this keyword to change the reference temperature used in the 36 ; calculation when GILL is activated. 37 ; 38 ; @keyword PROFREF 39 ; Give a depth to this keyword which will be considerated as the refeence depth 40 ; (in this case, LEVEL has not any effect). the calculation will be effectuated 41 ; until this depth effecting an interpolation between the the last W level above 42 ; PROFREF and PROFREF. 28 43 ; 29 ; LEVEL: C''est le niveau de reference a prendre. Ce niveau est 30 ; definit tel que gdepw[level] est la profondeur de reference 44 ; @keyword SURFACE_LEVEL 45 ; It is the level where we wan to calculate the dynamic height. By default, it is 46 ; the level 0. 31 47 ; 32 ; SREF: donner une valeur a ce mot cle pour changer la salinite 33 ; de reference utiliser ds le calcul lorsque GILL est active. 34 ; 35 ; TREF: donner une valeur a ce mot cle pour changer la temperature 36 ; de reference utiliser ds le calcul lorsque GILL est active. 48 ; @returns 49 ; An array of the same size of sn and tn representing the dynamic height calculated 50 ; from a reference depth nd by rapport to a reference state. 37 51 ; 38 ; PROFREF: donner a ce mot cle une profondeur qui sera prise comme 39 ; la profondeur de reference (ds ce cas LEVEL n'' a aucun 40 ; effet). le calcul sera alors effectue jusqu''a cette 41 ; profondeur en effectuant une interpolation entre le dernier 42 ; niveau W au dessus de PROFREF et PROFREF. 52 ; @uses 53 ; common.pro 43 54 ; 44 ; SURFACE_LEVEL: C''est le niveau auquel on veut calculer la 45 ; hauteur dynamique. Par defaut c''est le niveau 0. 55 ; @restrictions 56 ; Points for which we can not calculate the dynamic height (whose the batymetry 57 ; is less deep than the reference depth) are put at the value !values.f_nan 46 58 ; 47 ; OUTPUTS:un tableau de la meme taille que sn et tn representant la 48 ; hauteur dynamique calculee a partir d''une profondeur de reference 49 ; et par rapport a un etat de reference. 59 ; @restrictions 60 ; approximation: The pressure in decibars is equal to the depth in meters (the pressure increase of 1bar all 10m) 50 61 ; 51 ; COMMON BLOCKS:52 ; common.pro62 ; @history 63 ; Sebastien Masson (smasson@lodyc.jussieu.fr) 53 64 ; 54 ; SIDE EFFECTS: 55 ; les points pour lesquels on nje peut calcule la hauteur dynamique 56 ; (dont la batymetrie est moins profonde que la profondeur de 57 ; reference) sont mis a la valeur !values.f_nan 58 ; 59 ; RESTRICTIONS:approximation: la pression en decibars est egale a la 60 ; profondeur en metres (la pression augmente de 1bar tous les 10m) 61 ; 62 ; EXAMPLE: 63 ; 64 ; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr) 65 ; @version 66 ; $Id$ 65 67 ; 66 68 ;- … … 72 74 compile_opt idl2, strictarrsubs 73 75 ; 74 tempsun = systime(1) ; pour key_performance76 tempsun = systime(1) ; for key_performance 75 77 @common 76 78 77 79 if NOT keyword_set(surface_level) then surface_level = 0 78 ; u tile si GILL est active80 ; useful if GILL is activated 79 81 if NOT keyword_set(tref) then tref = 0. 80 82 if NOT keyword_set(sref) then sref = 35. 81 ; on determine si besoin est la profondeur de reference et le niveau 82 ; W situe directement au dessus. 83 ; If needed, we determinate the reference depth and the W level situated directly above. 83 84 if keyword_set(profref) then begin 84 85 rien = where(gdepw LE profref, level)
Note: See TracChangeset
for help on using the changeset viewer.