|
Last change
on this file since 5225 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.
