Last change
on this file since 5377 was
2089,
checked in by Laurent Fairhead, 10 years ago
|
Inclusion de la physique de MAR
Integration of MAR physics
|
File size:
477 bytes
|
Line | |
---|
1 | subroutine qsat (dq, q, e, p, t, r) |
---|
2 | |
---|
3 | #include "tsmbkind.h" |
---|
4 | |
---|
5 | implicit none |
---|
6 | |
---|
7 | REAL_B :: DQ, Q, E, P, T, R, TF, A, B |
---|
8 | |
---|
9 | tf=273.16_JPRB |
---|
10 | if (t.gt.tf) then |
---|
11 | a = 17.269_JPRB |
---|
12 | b = 35.86_JPRB |
---|
13 | else |
---|
14 | a = 21.874_JPRB |
---|
15 | b = 7.66_JPRB |
---|
16 | end if |
---|
17 | e = 6.11_JPRB * exp(a*(t-tf)/(t-b)) |
---|
18 | q = 0.622_JPRB * r*e/(p-r*e) |
---|
19 | dq=a*p*q*(tf-b) |
---|
20 | dq=dq/((p-r*e)*(t-b)*(t-b)) |
---|
21 | |
---|
22 | !print 9001, P,T,R,A,B,E,Q |
---|
23 | 9001 format (1x,'QSAT ',5F10.3,2E12.5) |
---|
24 | |
---|
25 | !---------------------------------------- |
---|
26 | return |
---|
27 | end subroutine qsat |
---|
Note: See
TracBrowser
for help on using the repository browser.