source: LMDZ6/branches/Optimisation_LMDZ/libf/phymar/qsat.F90 @ 5428

Last change on this file since 5428 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 
1subroutine qsat (dq, q, e, p, t, r)
2
3#include "tsmbkind.h"
4
5implicit none
6
7REAL_B :: DQ, Q, E, P, T, R, TF, A, B
8
9tf=273.16_JPRB
10if (t.gt.tf) then
11  a = 17.269_JPRB
12  b = 35.86_JPRB
13else
14  a = 21.874_JPRB
15  b = 7.66_JPRB
16end if 
17e = 6.11_JPRB  * exp(a*(t-tf)/(t-b))
18q = 0.622_JPRB * r*e/(p-r*e)
19dq=a*p*q*(tf-b)
20dq=dq/((p-r*e)*(t-b)*(t-b))
21
22!print 9001, P,T,R,A,B,E,Q
239001 format (1x,'QSAT ',5F10.3,2E12.5)
24
25!----------------------------------------
26return
27end subroutine qsat
Note: See TracBrowser for help on using the repository browser.