Last change
on this file since 3080 was
2089,
checked in by Laurent Fairhead, 10 years ago
|
Inclusion de la physique de MAR
Integration of MAR physics
|
File size:
477 bytes
|
Rev | Line | |
---|
[2089] | 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.