Last change
on this file since 325 was
325,
checked in by pinsard, 13 years ago
|
rehab MSG, EPSAT tools (to be cont.)
|
-
Property svn:keywords set to
Id
|
File size:
979 bytes
|
Line | |
---|
1 | function rmix=hspe2rmix(hspe,p) |
---|
2 | |
---|
3 | %HSPE2RMIX Calcul du rapport de melange a partir de l'humidite specifique |
---|
4 | |
---|
5 | %+ |
---|
6 | % |
---|
7 | % .. _hspe2rmix.m: |
---|
8 | % |
---|
9 | % =========== |
---|
10 | % hspe2rmix.m |
---|
11 | % =========== |
---|
12 | % |
---|
13 | % DESCRIPTION |
---|
14 | % =========== |
---|
15 | % |
---|
16 | % Calcul du rapport de melange a partir de |
---|
17 | % l'humidite specifique et la pression en utilisant : |
---|
18 | % |
---|
19 | % - rmix = 0.622*e/(p-e) |
---|
20 | % - hspe = 0.622*e/(p-0.378*e)) avec e la pression partielle de vapeur d'eau |
---|
21 | % |
---|
22 | % Entrees : |
---|
23 | % |
---|
24 | % hspe |
---|
25 | % l'humidite specifique en g/kg ou kg/kg |
---|
26 | % p |
---|
27 | % la pression de l'air en hPa ou Pa |
---|
28 | % |
---|
29 | % Sortie : |
---|
30 | % |
---|
31 | % rmix |
---|
32 | % rapport de melange en kg/kg |
---|
33 | % |
---|
34 | % EXAMPLES |
---|
35 | % ======== |
---|
36 | % |
---|
37 | % voir :ref:`cal_hcl.m` pour exemple d'utilisation |
---|
38 | % |
---|
39 | % SEE ALSO |
---|
40 | % ======== |
---|
41 | % |
---|
42 | % TODO |
---|
43 | % ==== |
---|
44 | % |
---|
45 | % EVOLUTIONS |
---|
46 | % ========== |
---|
47 | % |
---|
48 | % $Id$ |
---|
49 | % |
---|
50 | % $URL$ |
---|
51 | % |
---|
52 | %- |
---|
53 | |
---|
54 | %test si hspe est en g/kg ou kg/kg |
---|
55 | ind_nan=find(~isnan(hspe)); |
---|
56 | if hspe(ind_nan(1))>1 |
---|
57 | hspe=hspe/1000; |
---|
58 | end |
---|
59 | |
---|
60 | % Pression partielle de vapeur d'eau |
---|
61 | e=hspe.*p./(0.622+0.378*hspe); |
---|
62 | |
---|
63 | % Rapport de melange |
---|
64 | rmix=hspe.*((p-0.378*e)./(p-e)); |
---|
65 | |
---|
Note: See
TracBrowser
for help on using the repository browser.