Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

profile.F @ 3436

Last change on this file since 3436 was 3175, checked in by emillour, 11 months ago
Pluto PCM: Add the old Pluto LMDZ for reference (required prior step to making an LMDZ.PLUTO using the same framework as the other physics packages). TB+EM
File size: 7.2 KB

Line
1	SUBROUTINE profile(nlev,zkm,temp)
2	! to use 'getin'
3	USE ioipsl_getincom
4	IMPLICIT NONE
5	c=======================================================================
6	c Subroutine utilisee dans le modele 1-D "testphys1d"
7	c pour l'initialisation du profil atmospherique
8	c=======================================================================
9	c
10	c differents profils d'atmospheres. T=f(z)
11	c entree:
12	c nlev nombre de niveaux
13	c zkm alititudes en km
14	c ichoice choix de l'atmosphere:
15	c 1 Temperature constante
16	c 2 profil Savidjari
17	c 3 Lindner (profil polaire)
18	c 4 Inversion pour Francois
19	c 5 Seiff (moyen)
20	c 6 T constante + perturbation gauss (level) (christophe 10/98)
21	c 7 T constante + perturbation gauss (km) (christophe 10/98)
22	c 8 Lecture du profile dans un fichier ASCII (profile)
23	c tref temperature de reference
24	c isin ajout d'une perturbation (isin=1)
25	c pic pic perturbation gauss pour ichoice = 6 ou 7
26	c largeur largeur de la perturbation gauss pour ichoice = 6 ou 7
27	c hauteur hauteur de la perturbation gauss pour ichoice = 6 ou 7
28	c
29	c sortie:
30	c temp temperatures en K
31	c
32	c=======================================================================
33	c-----------------------------------------------------------------------
34	c declarations:
35	c -------------
36
37	c arguments:
38	c ----------
39
40	INTEGER nlev, unit
41	REAL zkm(nlev),temp(nlev)
42
43	c local:
44	c ------
45
46	INTEGER il,ichoice,nseiff,iseiff,isin,iter
47	REAL pi
48	PARAMETER(nseiff=37)
49	REAL tref,t1,t2,t3,ww
50	REAL tseiff(nseiff)
51	DATA tseiff/214.,213.8,213.4,212.4,209.3,205.,201.4,197.8,
52	$ 194.6,191.4,188.2,185.2,182.5,180.,177.5,175,
53	$ 172.5,170.,167.5,164.8,162.4,160.,158.,156.,
54	$ 154.1,152.2,150.3,148.7,147.2,145.7,144.2,143.,
55	$ 142.,141.,140,139.5,139./
56	REAL pic,largeur
57	REAL hauteur,tmp
58	REAL phitop,ttop
59	REAL tsurf,r,rcp,g,s(nlev)
60
61	c-----------------------------------------------------------------------
62	c read input profile type:
63	c--------------------------
64
65	ichoice=1 ! default value for ichoice
66	call getin("ichoice",ichoice)
67	tref=200 ! default value for tref
68	call getin("tref",tref)
69	isin=0 ! default value for isin (=0 means no perturbation)
70	call getin("isin",isin)
71	pic=26.522 ! default value for pic
72	call getin("pic",pic)
73	largeur=10 ! default value for largeur
74	call getin("largeur",largeur)
75	hauteur=30 ! default value for hauteur
76	call getin("hauteur",hauteur)
77
78	c-----------------------------------------------------------------------
79	c ichoice=1 temperature constante:
80	c --------------------------------
81
82	IF(ichoice.EQ.1) THEN
83	DO il=1,nlev
84	temp(il)=tref
85	ENDDO
86
87	c-----------------------------------------------------------------------
88	c ichoice=2 savidjari:
89	c --------------------
90
91	ELSE IF(ichoice.EQ.2) THEN
92	ccc DO il=1,7
93	cc temp(il)=AMAX1(219.-2.5*zkm(il),140.)
94	c pour pluton on prend en compte un profil lineaire qui va de 40 K au sol à 100 K à 50km
95	ccc temp(il)=38
96	cc temp(il)=6*zkm(il)+40
97	ccc ENDDO
98	ccc DO il=8,nlev
99	ccc temp(il)=100
100	ccc ENDDO
101
102
103	temp(1)=44.
104	temp(2)=44.
105	temp(3)=44.
106	temp(4)=44.
107	temp(5)=44.
108	temp(6)=44.
109	temp(7)=44.
110	temp(8)=44.
111	temp(9)=44.
112	temp(10)=44.
113	temp(11)=44.
114	temp(12)=44.
115	temp(13)=44.
116	temp(14)=44.
117	temp(15)=44.
118	temp(16)=44.
119	temp(17)=44.
120	temp(18)=45.85745
121	temp(19)=56.81386
122	temp(20)=71.71157
123	temp(21)=79.29989
124	temp(22)=82.10749
125	temp(23)=83.93003
126	temp(24)=85.58517
127	temp(25)=85.60072
128	c-----------------------------------------------------------------------
129	c ichoice=3 Lindner:
130	c ------------------
131
132	ELSE IF(ichoice.EQ.3) THEN
133	DO il=1,nlev
134	IF(zkm(il).LT.2.5) THEN
135	temp(il)=150.+30.*zkm(il)/2.5
136	ELSE IF(zkm(il).LT.5.) THEN
137	temp(il)=180.
138	ELSE
139	temp(il)=AMAX1(180.-2.*(zkm(il)-5.),130.)
140	ENDIF
141	ENDDO
142
143	c-----------------------------------------------------------------------
144	c ichoice=4 Inversion pour Francois:
145	c ------------------
146
147	ELSE IF(ichoice.EQ.4) THEN
148	DO il=1,nlev
149	IF(zkm(il).LT.20.) THEN
150	temp(il)=135.
151	ELSE
152	temp(il)=AMIN1(135.+5.*(zkm(il)-20.),200.)
153	ENDIF
154	ENDDO
155
156
157	c-----------------------------------------------------------------------
158	c ichoice=5 Seiff:
159	c ----------------
160
161	ELSE IF(ichoice.EQ.5) THEN
162	DO il=1,nlev
163	iseiff=INT(zkm(il)/2.)+1
164	IF(iseiff.LT.nseiff-1) THEN
165	temp(il)=tseiff(iseiff)+(zkm(il)-2.(iseiff-1))
166	$ (tseiff(iseiff+1)-tseiff(iseiff))/2.
167	ELSE
168	temp(il)=tseiff(nseiff)
169	ENDIF
170	ENDDO
171	c IF(ichoice.EQ.6) THEN
172	c DO iter=1,6
173	c t2=temp(1)
174	c t3=temp(2)
175	c DO il=2,nlev-1
176	c t1=t2
177	c t2=t3
178	c t3=temp(il+1)
179	c ww=tanh(zkm(il)/20.)
180	c ww=wwwwww
181	c temp(il)=t2+ww.5(t1-2.*t2+t3)
182	c ENDDO
183	c ENDDO
184	c ENDIF
185
186	c-----------------------------------------------------------------------
187	c ichoice=6
188	c ---------
189
190	ELSE IF(ichoice.EQ.6) THEN
191	DO il=1,nlev
192	tmp=il-pic
193	temp(il)=tref + hauteurexp(-tmptmp/largeur/largeur)
194	ENDDO
195
196
197	c-----------------------------------------------------------------------
198	c ichoice=7
199	c ---------
200
201	ELSE IF(ichoice.EQ.7) THEN
202	DO il=1,nlev
203	tmp=zkm(il)-pic
204	temp(il)=tref + hauteurexp(-tmptmp*4/largeur/largeur)
205	ENDDO
206
207	c-----------------------------------------------------------------------
208	c ichoice=8
209	c ---------
210
211	ELSE IF(ichoice.GE.8) THEN
212	OPEN(11,file='profile',status='old',form='formatted',err=101)
213	DO il=1,nlev
214	READ (11,*) temp(il)
215	ENDDO
216
217	GOTO 201
218	101 STOP'fichier profile inexistant'
219	201 CONTINUE
220	CLOSE(10)
221
222	c-----------------------------------------------------------------------
223
224	ENDIF
225
226	c-----------------------------------------------------------------------
227	c rajout eventuel d'une perturbation:
228	c -----------------------------------
229
230	IF(isin.EQ.1) THEN
231	pi=2.*ASIN(1.)
232	DO il=1,nlev
233	c if (nlev.EQ.501) then
234	c if (zkm(il).LE.70.5) then
235	c temp(il)=temp(il)+(1.-1000./(1000+zkm(il)zkm(il)))(
236	c s 6.SIN(zkm(il)pi/6.)+9.SIN(zkm(il)pi/10.3) )
237	c endif
238	c else
239	temp(il)=temp(il)+(1.-1000./(1000+zkm(il)zkm(il)))(
240	s 6.SIN(zkm(il)pi/6.)+9.SIN(zkm(il)pi/10.3) )
241	c endif
242	ENDDO
243	ENDIF
244
245
246	c-----------------------------------------------------------------------
247	c Ecriture du profil de temperature dans un fichier profile.out
248	c -------------------------------------------------------------
249
250
251	OPEN(12,file='profile.out',form='formatted')
252	DO il=1,nlev
253	write(12,*) temp(il)
254	ENDDO
255	CLOSE(12)
256
257	RETURN
258	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format