1 | ! ADAPTATION GCM POUR CP(T) |
---|
2 | c====================================================================== |
---|
3 | c S. Lebonnois, 10/2007: |
---|
4 | c |
---|
5 | c VENUS: Cp(T) = cpp*(T/T0)^nu |
---|
6 | c avec T0=460. et nu=0.35 |
---|
7 | c cpp=RCPD=cp0 = 1000. |
---|
8 | c R/RCPD = RKAPPA |
---|
9 | c |
---|
10 | c La fonction d'Exner reste pk = RCPD*(play/pref)**RKAPPA |
---|
11 | c |
---|
12 | c T et teta (temperature potentielle) sont liees par: |
---|
13 | c |
---|
14 | c integrale[teta a T](cp/T dT) = integrale[pref a p](R/p dp) |
---|
15 | c |
---|
16 | c Dans le cas de l'expression pour Venus, ca donne: |
---|
17 | c |
---|
18 | c teta**nu = T**nu - nu * T0**nu * ln[ (p/pref)**RKAPPA ] |
---|
19 | c ou |
---|
20 | c teta**nu = T**nu - nu * T0**nu * ln[pk/RCPD] |
---|
21 | c |
---|
22 | c On passe de T a teta par t2tpot(t,teta,pk) |
---|
23 | c On passe de teta a T par tpot2t(teta,t,pk) |
---|
24 | c |
---|
25 | c Pour DT <-> Dteta, on utilise: dteta = dT *(T/teta)**(nu-1) |
---|
26 | c -> routine dt2dtpot(dt,dteta,t,teta) |
---|
27 | c (utilisee seulement pour le contregradient) |
---|
28 | c |
---|
29 | c====================================================================== |
---|
30 | |
---|
31 | FUNCTION cpdet(t) |
---|
32 | IMPLICIT none |
---|
33 | !----------------------------------------------------------------------- |
---|
34 | ! INCLUDE cpdet.h |
---|
35 | |
---|
36 | real cp0,t0,nu |
---|
37 | ! parameter (cp0 = 1000.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
38 | parameter (cp0 = 800.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
39 | ! parameter (t0 = 460.) |
---|
40 | parameter (t0 = 480.) |
---|
41 | ! parameter (t0 = 300.) |
---|
42 | parameter (nu = 0.35) |
---|
43 | ! parameter (nu = 1.0) |
---|
44 | !----------------------------------------------------------------------- |
---|
45 | |
---|
46 | real cpdet,t |
---|
47 | |
---|
48 | cpdet = cp0*(t/t0)**nu |
---|
49 | |
---|
50 | return |
---|
51 | end |
---|
52 | |
---|
53 | c====================================================================== |
---|
54 | c====================================================================== |
---|
55 | |
---|
56 | SUBROUTINE t2tpot(npoints,yt, yteta, ypk) |
---|
57 | IMPLICIT none |
---|
58 | c====================================================================== |
---|
59 | c Arguments: |
---|
60 | c |
---|
61 | c yt --------input-R- Temperature |
---|
62 | c yteta-------output-R- Temperature potentielle |
---|
63 | c ypk --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA |
---|
64 | c |
---|
65 | c====================================================================== |
---|
66 | !----------------------------------------------------------------------- |
---|
67 | ! INCLUDE cpdet.h |
---|
68 | |
---|
69 | real cp0,t0,nu |
---|
70 | ! parameter (cp0 = 1000.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
71 | parameter (cp0 = 800.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
72 | ! parameter (t0 = 460.) |
---|
73 | parameter (t0 = 480.) |
---|
74 | ! parameter (t0 = 300.) |
---|
75 | parameter (nu = 0.35) |
---|
76 | ! parameter (nu = 1.0) |
---|
77 | !----------------------------------------------------------------------- |
---|
78 | |
---|
79 | integer npoints |
---|
80 | REAL yt(npoints), yteta(npoints), ypk(npoints) |
---|
81 | |
---|
82 | yteta = yt**nu - nu * t0**nu * log(ypk/cp0) |
---|
83 | yteta = yteta**(1./nu) |
---|
84 | |
---|
85 | |
---|
86 | !yteta = yt/ypk |
---|
87 | |
---|
88 | |
---|
89 | return |
---|
90 | end |
---|
91 | |
---|
92 | c====================================================================== |
---|
93 | c====================================================================== |
---|
94 | |
---|
95 | SUBROUTINE tpot2t(npoints,yteta, yt, ypk) |
---|
96 | IMPLICIT none |
---|
97 | c====================================================================== |
---|
98 | c Arguments: |
---|
99 | c |
---|
100 | c yteta--------input-R- Temperature potentielle |
---|
101 | c yt -------output-R- Temperature |
---|
102 | c ypk --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA |
---|
103 | c |
---|
104 | c====================================================================== |
---|
105 | !----------------------------------------------------------------------- |
---|
106 | ! INCLUDE cpdet.h |
---|
107 | |
---|
108 | real cp0,t0,nu |
---|
109 | ! parameter (cp0 = 1000.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
110 | parameter (cp0 = 800.) !doit etre egal a cpp (dyn) et RCPD (phy) |
---|
111 | ! parameter (t0 = 460.) |
---|
112 | parameter (t0 = 480.) |
---|
113 | ! parameter (t0 = 300.) |
---|
114 | parameter (nu = 0.35) |
---|
115 | ! parameter (nu = 1.0) |
---|
116 | !----------------------------------------------------------------------- |
---|
117 | |
---|
118 | integer npoints |
---|
119 | REAL yt(npoints), yteta(npoints), ypk(npoints) |
---|
120 | |
---|
121 | yt = yteta**nu + nu * t0**nu * log(ypk/cp0) |
---|
122 | yt = yt**(1./nu) |
---|
123 | |
---|
124 | !yt = yteta*ypk |
---|
125 | |
---|
126 | return |
---|
127 | end |
---|
128 | |
---|
129 | c$$$c====================================================================== |
---|
130 | c$$$c====================================================================== |
---|
131 | c$$$ |
---|
132 | c$$$ SUBROUTINE dt2dtpot(npoints,ydt, ydteta, yt, yteta, ypk, ydpk) |
---|
133 | c$$$ IMPLICIT none |
---|
134 | c$$$c====================================================================== |
---|
135 | c$$$c Arguments: |
---|
136 | c$$$c |
---|
137 | c$$$c ydt --------input-R- D Temperature |
---|
138 | c$$$c ydteta------output-R- D Temperature potentielle |
---|
139 | c$$$c yt --------input-R- Temperature |
---|
140 | c$$$c yteta--------input-R- Temperature potentielle |
---|
141 | c$$$c ypk --------input-R- Fonction d'Exner: RCPD*(pplay/pref)**RKAPPA |
---|
142 | c$$$c ydpk -------input-R- D Fonction d'Exner |
---|
143 | c$$$c |
---|
144 | c$$$c====================================================================== |
---|
145 | c$$$#include "cpdet.h" |
---|
146 | c$$$ |
---|
147 | c$$$ integer npoints |
---|
148 | c$$$ REAL ydt(npoints), ydteta(npoints), ydpk(npoints) |
---|
149 | c$$$ REAL yt(npoints), yteta(npoints), ypk(npoints) |
---|
150 | c$$$ |
---|
151 | c$$$ ydteta = yteta**(1.-nu)* (ydt * yt**(nu-1.) |
---|
152 | c$$$ . - t0**nu * ydpk / ypk) |
---|
153 | c$$$ |
---|
154 | c$$$ return |
---|
155 | c$$$ end |
---|
156 | |
---|