Context Navigation

startvar.F

Last change on this file was 633, checked in by (none), 20 years ago
This commit was manufactured by cvs2svn to create branch 'LMDZ4_par_0'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 39.7 KB

Line
1	!
2	! $Header$
3	!
4	C
5	C
6	MODULE startvar
7	!
8	!
9	! There are three ways to access data from the database of atmospheric data which
10	! can be used to initialize the model. This depends on the type of field which needs
11	! to be extracted. In any case the call should come after a restget and should be of the type :
12	! CALL startget(...)
13	!
14	! We will details the possible arguments to startget here :
15	!
16	! - A 2D variable on the dynamical grid :
17	! CALL startget(varname, iml, jml, lon_in, lat_in, champ, val_ex, jml2, lon_in2, lat_in2, interbar )
18	!
19	! - A 1D variable on the physical grid :
20	! CALL startget(varname, iml, jml, lon_in, lat_in, nbindex, champ, val_exp, jml2, lon_in2, lat_in2, interbar )
21	!
22	!
23	! - A 3D variable on the dynamical grid :
24	! CALL startget(varname, iml, jml, lon_in, lat_in, lml, pls, workvar, champ, val_exp, jml2, lon_in2, lat_in2, interbar )
25	!
26	!
27	! There is special constraint on the atmospheric data base except that the
28	! the data needs to be in netCDF and the variables should have the the following
29	! names in the file :
30	!
31	! 'RELIEF' : High resolution orography
32	! 'ST' : Surface temperature
33	! 'CDSW' : Soil moisture
34	! 'Z' : Surface geopotential
35	! 'SP' : Surface pressure
36	! 'U' : East ward wind
37	! 'V' : Northward wind
38	! 'TEMP' : Temperature
39	! 'R' : Relative humidity
40	!
41	USE ioipsl
42	!
43	!
44	IMPLICIT NONE
45	!
46	!
47	PRIVATE
48	PUBLIC startget
49	!
50	!
51	INTERFACE startget
52	MODULE PROCEDURE startget_phys2d, startget_phys1d, startget_dyn
53	END INTERFACE
54	!
55	INTEGER, SAVE :: fid_phys, fid_dyn
56	INTEGER, SAVE :: iml_phys, iml_rel, iml_dyn
57	INTEGER, SAVE :: jml_phys, jml_rel, jml_dyn
58	INTEGER, SAVE :: llm_dyn, ttm_dyn
59	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:) :: lon_phys, lon_rug,
60	. lon_alb, lon_rel, lon_dyn
61	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:) :: lat_phys, lat_rug,
62	. lat_alb, lat_rel, lat_dyn
63	REAL, ALLOCATABLE, SAVE, DIMENSION (:) :: levdyn_ini
64	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:) :: relief, zstd, zsig,
65	. zgam, zthe, zpic, zval
66	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:) :: rugo, masque, phis
67	!
68	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:) :: tsol, qsol, psol_dyn
69	REAL, ALLOCATABLE, SAVE, DIMENSION (:,:,:) :: var_ana3d
70	!
71	CONTAINS
72	!
73	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
74	!
75	SUBROUTINE startget_phys2d(varname, iml, jml, lon_in, lat_in,
76	. champ, val_exp, jml2, lon_in2, lat_in2 , interbar, masque_lu )
77	!
78	! There is a big mess with the size in logitude, should it be iml or iml+1.
79	! I have chosen to use the iml+1 as an argument to this routine and we declare
80	! internaly smaler fields when needed. This needs to be cleared once and for all in LMDZ.
81	! A convention is required.
82	!
83	!
84	CHARACTER(), INTENT(in) :: varname
85	INTEGER, INTENT(in) :: iml, jml ,jml2
86	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
87	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
88	REAL, INTENT(inout) :: champ(iml,jml)
89	REAL, INTENT(in) :: val_exp
90	REAL, INTENT(in), optional :: masque_lu(iml,jml)
91	LOGICAL interbar
92	!
93	! This routine only works if the variable does not exist or is constant
94	!
95	IF ( MINVAL(champ(:,:)).EQ.MAXVAL(champ(:,:)) .AND.
96	.MINVAL(champ(:,:)).EQ.val_exp ) THEN
97	!
98	SELECTCASE(varname)
99	!
100	CASE ('relief')
101	!
102	! If we do not have the orography we need to get it
103	!
104	IF ( .NOT.ALLOCATED(relief)) THEN
105	!
106	if (present(masque_lu)) then
107	CALL start_init_orog( iml, jml, lon_in, lat_in,
108	. jml2,lon_in2,lat_in2, interbar, masque_lu )
109	else
110	CALL start_init_orog( iml, jml, lon_in, lat_in,
111	. jml2,lon_in2,lat_in2, interbar)
112	endif
113	!
114	ENDIF
115	!
116	IF (SIZE(relief) .NE. SIZE(champ)) THEN
117	!
118	WRITE(,) 'STARTVAR module has been',
119	.' initialized to the wrong size'
120	STOP
121	!
122	ENDIF
123	!
124	champ(:,:) = relief(:,:)
125	!
126	CASE ('rugosite')
127	!
128	! If we do not have the orography we need to get it
129	!
130	IF ( .NOT.ALLOCATED(rugo)) THEN
131	!
132	CALL start_init_orog( iml, jml, lon_in, lat_in,
133	. jml2,lon_in2,lat_in2 , interbar )
134	!
135	ENDIF
136	!
137	IF (SIZE(rugo) .NE. SIZE(champ)) THEN
138	!
139	WRITE(,)
140	. 'STARTVAR module has been initialized to the wrong size'
141	STOP
142	!
143	ENDIF
144	!
145	champ(:,:) = rugo(:,:)
146	!
147	CASE ('masque')
148	!
149	! If we do not have the orography we need to get it
150	!
151	IF ( .NOT.ALLOCATED(masque)) THEN
152	!
153	CALL start_init_orog( iml, jml, lon_in, lat_in,
154	. jml2,lon_in2,lat_in2 , interbar )
155	!
156	ENDIF
157	!
158	IF (SIZE(masque) .NE. SIZE(champ)) THEN
159	!
160	WRITE(,)
161	. 'STARTVAR module has been initialized to the wrong size'
162	STOP
163	!
164	ENDIF
165	!
166	champ(:,:) = masque(:,:)
167	!
168	CASE ('surfgeo')
169	!
170	! If we do not have the orography we need to get it
171	!
172	IF ( .NOT.ALLOCATED(phis)) THEN
173	!
174	CALL start_init_orog( iml, jml, lon_in, lat_in,
175	. jml2,lon_in2, lat_in2 , interbar )
176	!
177	ENDIF
178	!
179	IF (SIZE(phis) .NE. SIZE(champ)) THEN
180	!
181	WRITE(,)
182	. 'STARTVAR module has been initialized to the wrong size'
183	STOP
184	!
185	ENDIF
186	!
187	champ(:,:) = phis(:,:)
188	!
189	CASE ('psol')
190	!
191	! If we do not have the orography we need to get it
192	!
193	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
194	!
195	CALL start_init_dyn( iml, jml, lon_in, lat_in,
196	. jml2,lon_in2, lat_in2 , interbar )
197	!
198	ENDIF
199	!
200	IF (SIZE(psol_dyn) .NE. SIZE(champ)) THEN
201	!
202	WRITE(,)
203	. 'STARTVAR module has been initialized to the wrong size'
204	STOP
205	!
206	ENDIF
207	!
208	champ(:,:) = psol_dyn(:,:)
209	!
210	CASE DEFAULT
211	!
212	WRITE(,) 'startget_phys2d'
213	WRITE(,) 'No rule is present to extract variable',
214	. varname(:LEN_TRIM(varname)),' from any data set'
215	STOP
216	!
217	END SELECT
218	!
219	ELSE
220	!
221	! There are a few fields we might need if we need to interpolate 3D filed. Thus if they come through here we
222	! will catch them
223	!
224	SELECTCASE(varname)
225	!
226	CASE ('surfgeo')
227	!
228	IF ( .NOT.ALLOCATED(phis)) THEN
229	ALLOCATE(phis(iml,jml))
230	ENDIF
231	!
232	IF (SIZE(phis) .NE. SIZE(champ)) THEN
233	!
234	WRITE(,)
235	. 'STARTVAR module has been initialized to the wrong size'
236	STOP
237	!
238	ENDIF
239	!
240	phis(:,:) = champ(:,:)
241	!
242	END SELECT
243	!
244	ENDIF
245	!
246	END SUBROUTINE startget_phys2d
247	!
248	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
249	!
250	SUBROUTINE start_init_orog ( iml,jml,lon_in, lat_in,jml2,lon_in2 ,
251	, lat_in2 , interbar, masque_lu )
252	!
253	INTEGER, INTENT(in) :: iml, jml, jml2
254	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
255	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
256	REAL, intent(in), optional :: masque_lu(iml,jml)
257	LOGICAL interbar
258	!
259	! LOCAL
260	!
261	LOGICAL interbar2
262	REAL :: lev(1), date, dt,chmin,chmax
263	INTEGER :: itau(1), fid
264	INTEGER :: llm_tmp, ttm_tmp
265	INTEGER :: i, j
266	INTEGER :: iret
267	CHARACTER*25 title
268	REAL, ALLOCATABLE :: relief_hi(:,:)
269	REAL, ALLOCATABLE :: lon_rad(:), lat_rad(:)
270	REAL, ALLOCATABLE :: lon_ini(:), lat_ini(:)
271	REAL, ALLOCATABLE :: tmp_var(:,:)
272	INTEGER, ALLOCATABLE :: tmp_int(:,:)
273	!
274	CHARACTER*120 :: orogfname
275	LOGICAL :: check=.TRUE.
276	!
277	!
278	orogfname = 'Relief.nc'
279	!
280	IF ( check ) WRITE(,) 'Reading the high resolution orography'
281	!
282	CALL flininfo(orogfname,iml_rel, jml_rel, llm_tmp, ttm_tmp, fid)
283	!
284	ALLOCATE (lat_rel(iml_rel,jml_rel), stat=iret)
285	ALLOCATE (lon_rel(iml_rel,jml_rel), stat=iret)
286	ALLOCATE (relief_hi(iml_rel,jml_rel), stat=iret)
287	!
288	CALL flinopen(orogfname, .FALSE., iml_rel, jml_rel,
289	.llm_tmp, lon_rel, lat_rel, lev, ttm_tmp,
290	. itau, date, dt, fid)
291	!
292	CALL flinget(fid, 'RELIEF', iml_rel, jml_rel, llm_tmp,
293	. ttm_tmp, 1, 1, relief_hi)
294	!
295	CALL flinclo(fid)
296	!
297	! In case we have a file which is in degrees we do the transformation
298	!
299	ALLOCATE(lon_rad(iml_rel))
300	ALLOCATE(lon_ini(iml_rel))
301
302	IF ( MAXVAL(lon_rel(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
303	lon_ini(:) = lon_rel(:,1) * 2.0 * ASIN(1.0) / 180.0
304	ELSE
305	lon_ini(:) = lon_rel(:,1)
306	ENDIF
307
308	ALLOCATE(lat_rad(jml_rel))
309	ALLOCATE(lat_ini(jml_rel))
310
311	IF ( MAXVAL(lat_rel(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
312	lat_ini(:) = lat_rel(1,:) * 2.0 * ASIN(1.0) / 180.0
313	ELSE
314	lat_ini(:) = lat_rel(1,:)
315	ENDIF
316	!
317	!
318
319	title='RELIEF'
320
321	interbar2 = .FALSE.
322	CALL conf_dat2d(title,iml_rel, jml_rel, lon_ini, lat_ini,
323	. lon_rad, lat_rad, relief_hi , interbar2 )
324
325	IF ( check ) WRITE(,) 'Computes all the parameters needed',
326	.' for the gravity wave drag code'
327	!
328	! Allocate the data we need to put in the interpolated fields
329	!
330	! RELIEF: orographie moyenne
331	ALLOCATE(relief(iml,jml))
332	! zphi : orographie moyenne
333	ALLOCATE(phis(iml,jml))
334	! zstd: deviation standard de l'orographie sous-maille
335	ALLOCATE(zstd(iml,jml))
336	! zsig: pente de l'orographie sous-maille
337	ALLOCATE(zsig(iml,jml))
338	! zgam: anisotropy de l'orographie sous maille
339	ALLOCATE(zgam(iml,jml))
340	! zthe: orientation de l'axe oriente dans la direction
341	! de plus grande pente de l'orographie sous maille
342	ALLOCATE(zthe(iml,jml))
343	! zpic: hauteur pics de la SSO
344	ALLOCATE(zpic(iml,jml))
345	! zval: hauteur vallees de la SSO
346	ALLOCATE(zval(iml,jml))
347	! masque : Masque terre ocean
348	ALLOCATE(tmp_int(iml,jml))
349	ALLOCATE(masque(iml,jml))
350
351	masque = -99999.
352	if (present(masque_lu)) then
353	masque = masque_lu
354	endif
355	!
356	CALL grid_noro(iml_rel, jml_rel, lon_rad, lat_rad, relief_hi,
357	. iml-1, jml, lon_in, lat_in,
358	. phis, relief, zstd, zsig, zgam, zthe, zpic, zval, masque)
359	phis = phis * 9.81
360	!
361	! masque(:,:) = FLOAT(tmp_int(:,:))
362	!
363	! Compute surface roughness
364	!
365	IF ( check ) WRITE(,)
366	.'Compute surface roughness induced by the orography'
367	!
368	ALLOCATE(rugo(iml,jml))
369	ALLOCATE(tmp_var(iml-1,jml))
370	!
371	CALL rugsoro(iml_rel, jml_rel, lon_rad, lat_rad, relief_hi,
372	. iml-1, jml, lon_in, lat_in, tmp_var)
373	!
374	DO j = 1, jml
375	DO i = 1, iml-1
376	rugo(i,j) = tmp_var(i,j)
377	ENDDO
378	rugo(iml,j) = tmp_var(1,j)
379	ENDDO
380	c
381	cc * rugo n'est pas utilise pour l'instant ****
382	!
383	! Build land-sea mask
384	!
385	!
386	RETURN
387	!
388	END SUBROUTINE start_init_orog
389	!
390	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
391	!
392	SUBROUTINE startget_phys1d(varname, iml, jml, lon_in,
393	.lat_in, nbindex, champ, val_exp ,jml2, lon_in2, lat_in2,interbar)
394	!
395	CHARACTER(), INTENT(in) :: varname
396	INTEGER, INTENT(in) :: iml, jml, nbindex, jml2
397	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
398	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
399	REAL, INTENT(inout) :: champ(nbindex)
400	REAL, INTENT(in) :: val_exp
401	LOGICAL interbar
402	!
403	!
404	! This routine only works if the variable does not exist or is constant
405	!
406	IF ( MINVAL(champ(:)).EQ.MAXVAL(champ(:)) .AND.
407	.MINVAL(champ(:)).EQ.val_exp ) THEN
408	SELECTCASE(varname)
409	CASE ('tsol')
410	IF ( .NOT.ALLOCATED(tsol)) THEN
411	CALL start_init_phys( iml, jml, lon_in, lat_in,
412	. jml2, lon_in2, lat_in2, interbar )
413	ENDIF
414	IF ( SIZE(tsol) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
415	WRITE(,)
416	. 'STARTVAR module has been initialized to the wrong size'
417	STOP
418	ENDIF
419	CALL gr_dyn_fi(1, iml, jml, nbindex, tsol, champ)
420	CASE ('qsol')
421	IF ( .NOT.ALLOCATED(qsol)) THEN
422	CALL start_init_phys( iml, jml, lon_in, lat_in,
423	. jml2, lon_in2,lat_in2 , interbar )
424	ENDIF
425	IF ( SIZE(qsol) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
426	WRITE(,)
427	. 'STARTVAR module has been initialized to the wrong size'
428	STOP
429	ENDIF
430	CALL gr_dyn_fi(1, iml, jml, nbindex, qsol, champ)
431	CASE ('psol')
432	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
433	CALL start_init_dyn( iml, jml, lon_in, lat_in,
434	. jml2, lon_in2,lat_in2 , interbar )
435	ENDIF
436	IF (SIZE(psol_dyn) .NE. SIZE(lon_in)*SIZE(lat_in)) THEN
437	WRITE(,)
438	. 'STARTVAR module has been initialized to the wrong size'
439	STOP
440	ENDIF
441	CALL gr_dyn_fi(1, iml, jml, nbindex, psol_dyn, champ)
442	CASE ('zmea')
443	IF ( .NOT.ALLOCATED(relief)) THEN
444	CALL start_init_orog( iml, jml, lon_in, lat_in,
445	. jml2, lon_in2,lat_in2 , interbar )
446	ENDIF
447	IF ( SIZE(relief) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
448	WRITE(,)
449	. 'STARTVAR module has been initialized to the wrong size'
450	STOP
451	ENDIF
452	CALL gr_dyn_fi(1, iml, jml, nbindex, relief, champ)
453	CASE ('zstd')
454	IF ( .NOT.ALLOCATED(zstd)) THEN
455	CALL start_init_orog( iml, jml, lon_in, lat_in,
456	. jml2, lon_in2,lat_in2 , interbar )
457	ENDIF
458	IF ( SIZE(zstd) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
459	WRITE(,)
460	. 'STARTVAR module has been initialized to the wrong size'
461	STOP
462	ENDIF
463	CALL gr_dyn_fi(1, iml, jml, nbindex,zstd, champ)
464	CASE ('zsig')
465	IF ( .NOT.ALLOCATED(zsig)) THEN
466	CALL start_init_orog( iml, jml, lon_in, lat_in,
467	. jml2, lon_in2,lat_in2 , interbar )
468	ENDIF
469	IF ( SIZE(zsig) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
470	WRITE(,)
471	. 'STARTVAR module has been initialized to the wrong size'
472	STOP
473	ENDIF
474	CALL gr_dyn_fi(1, iml, jml, nbindex,zsig, champ)
475	CASE ('zgam')
476	IF ( .NOT.ALLOCATED(zgam)) THEN
477	CALL start_init_orog( iml, jml, lon_in, lat_in,
478	. jml2, lon_in2,lat_in2 , interbar )
479	ENDIF
480	IF ( SIZE(zgam) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
481	WRITE(,)
482	. 'STARTVAR module has been initialized to the wrong size'
483	STOP
484	ENDIF
485	CALL gr_dyn_fi(1, iml, jml, nbindex,zgam, champ)
486	CASE ('zthe')
487	IF ( .NOT.ALLOCATED(zthe)) THEN
488	CALL start_init_orog( iml, jml, lon_in, lat_in,
489	. jml2, lon_in2,lat_in2 , interbar )
490	ENDIF
491	IF ( SIZE(zthe) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
492	WRITE(,)
493	. 'STARTVAR module has been initialized to the wrong size'
494	STOP
495	ENDIF
496	CALL gr_dyn_fi(1, iml, jml, nbindex,zthe, champ)
497	CASE ('zpic')
498	IF ( .NOT.ALLOCATED(zpic)) THEN
499	CALL start_init_orog( iml, jml, lon_in, lat_in,
500	. jml2, lon_in2,lat_in2 , interbar )
501	ENDIF
502	IF ( SIZE(zpic) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
503	WRITE(,)
504	. 'STARTVAR module has been initialized to the wrong size'
505	STOP
506	ENDIF
507	CALL gr_dyn_fi(1, iml, jml, nbindex,zpic, champ)
508	CASE ('zval')
509	IF ( .NOT.ALLOCATED(zval)) THEN
510	CALL start_init_orog( iml, jml, lon_in, lat_in,
511	. jml2, lon_in2,lat_in2 , interbar )
512	ENDIF
513	IF ( SIZE(zval) .NE. SIZE(lon_in)*SIZE(lat_in) ) THEN
514	WRITE(,)
515	. 'STARTVAR module has been initialized to the wrong size'
516	STOP
517	ENDIF
518	CALL gr_dyn_fi(1, iml, jml, nbindex,zval, champ)
519	CASE ('rads')
520	champ(:) = 0.0
521	CASE ('snow')
522	champ(:) = 0.0
523	CASE ('deltat')
524	champ(:) = 0.0
525	CASE ('rugmer')
526	champ(:) = 0.001
527	CASE ('agsno')
528	champ(:) = 50.0
529	CASE DEFAULT
530	WRITE(,) 'startget_phys1d'
531	WRITE(,) 'No rule is present to extract variable ',
532	. varname(:LEN_TRIM(varname)),' from any data set'
533	STOP
534	END SELECT
535	ELSE
536	!
537	! If we see tsol we catch it as we may need it for a 3D interpolation
538	!
539	SELECTCASE(varname)
540	CASE ('tsol')
541	IF ( .NOT.ALLOCATED(tsol)) THEN
542	ALLOCATE(tsol(SIZE(lon_in),SIZE(lat_in) ))
543	ENDIF
544	CALL gr_fi_dyn(1, iml, jml, nbindex, champ, tsol)
545	END SELECT
546	ENDIF
547	END SUBROUTINE startget_phys1d
548	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
549	!
550	SUBROUTINE start_init_phys( iml, jml, lon_in, lat_in, jml2,
551	. lon_in2, lat_in2 , interbar )
552	!
553	INTEGER, INTENT(in) :: iml, jml ,jml2
554	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
555	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
556	LOGICAL interbar
557	!
558	! LOCAL
559	!
560	!ac REAL :: lev(1), date, dt
561	REAL :: date, dt
562	REAL, DIMENSION(:), ALLOCATABLE :: levphys_ini
563	!ac
564	INTEGER :: itau(1)
565	INTEGER :: llm_tmp, ttm_tmp
566	INTEGER :: i, j
567	!
568	CHARACTER*25 title
569	CHARACTER*120 :: physfname
570	LOGICAL :: check=.TRUE.
571	!
572	REAL, ALLOCATABLE :: lon_rad(:), lat_rad(:)
573	REAL, ALLOCATABLE :: lon_ini(:), lat_ini(:)
574	REAL, ALLOCATABLE :: var_ana(:,:), tmp_var(:,:)
575	!
576	physfname = 'ECPHY.nc'
577	!
578	IF ( check ) WRITE(,) 'Opening the surface analysis'
579	!
580	CALL flininfo(physfname, iml_phys, jml_phys, llm_tmp,
581	. ttm_tmp, fid_phys)
582	!
583	ALLOCATE (lat_phys(iml_phys,jml_phys))
584	ALLOCATE (lon_phys(iml_phys,jml_phys))
585	!ac
586	ALLOCATE (levphys_ini(llm_tmp))
587	!
588	! CALL flinopen(physfname, .FALSE., iml_phys, jml_phys,
589	! . llm_tmp, lon_phys, lat_phys, lev, ttm_tmp,
590	! . itau, date, dt, fid_phys)
591	!
592	CALL flinopen(physfname, .FALSE., iml_phys, jml_phys,
593	. llm_tmp, lon_phys, lat_phys, levphys_ini, ttm_tmp,
594	. itau, date, dt, fid_phys)
595	!
596	DEALLOCATE (levphys_ini)
597	!ac
598	!
599	! Allocate the space we will need to get the data out of this file
600	!
601	ALLOCATE(var_ana(iml_phys, jml_phys))
602	!
603	! In case we have a file which is in degrees we do the transformation
604	!
605	ALLOCATE(lon_rad(iml_phys))
606	ALLOCATE(lon_ini(iml_phys))
607
608	IF ( MAXVAL(lon_phys(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
609	lon_ini(:) = lon_phys(:,1) * 2.0 * ASIN(1.0) / 180.0
610	ELSE
611	lon_ini(:) = lon_phys(:,1)
612	ENDIF
613
614	ALLOCATE(lat_rad(jml_phys))
615	ALLOCATE(lat_ini(jml_phys))
616
617	IF ( MAXVAL(lat_phys(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
618	lat_ini(:) = lat_phys(1,:) * 2.0 * ASIN(1.0) / 180.0
619	ELSE
620	lat_ini(:) = lat_phys(1,:)
621	ENDIF
622
623
624	!
625	! We get the two standard varibales
626	! Surface temperature
627	!
628	ALLOCATE(tsol(iml,jml))
629	ALLOCATE(tmp_var(iml-1,jml))
630	!
631	!
632
633	CALL flinget(fid_phys, 'ST', iml_phys, jml_phys,
634	.llm_tmp, ttm_tmp, 1, 1, var_ana)
635
636	title='ST'
637	CALL conf_dat2d(title,iml_phys, jml_phys, lon_ini, lat_ini,
638	. lon_rad, lat_rad, var_ana , interbar )
639
640	IF ( interbar ) THEN
641	WRITE(6,*) '-------------------------------------------------',
642	,'--------------'
643	WRITE(6,*) '$$$ Utilisation de l interpolation barycentrique ',
644	, ' pour ST $$$ '
645	WRITE(6,*) '-------------------------------------------------',
646	,'--------------'
647	CALL inter_barxy ( iml_phys,jml_phys -1,lon_rad,lat_rad ,
648	, var_ana, iml-1, jml-1, lon_in2, lat_in2, jml, tmp_var )
649	ELSE
650	CALL grille_m(iml_phys, jml_phys, lon_rad, lat_rad,
651	. var_ana, iml-1, jml, lon_in, lat_in, tmp_var )
652	ENDIF
653
654	CALL gr_int_dyn(tmp_var, tsol, iml-1, jml)
655	!
656	! Soil moisture
657	!
658	ALLOCATE(qsol(iml,jml))
659	CALL flinget(fid_phys, 'CDSW', iml_phys, jml_phys,
660	. llm_tmp, ttm_tmp, 1, 1, var_ana)
661
662	title='CDSW'
663	CALL conf_dat2d(title,iml_phys, jml_phys, lon_ini, lat_ini,
664	. lon_rad, lat_rad, var_ana, interbar )
665
666	IF ( interbar ) THEN
667	WRITE(6,*) '-------------------------------------------------',
668	,'--------------'
669	WRITE(6,*) '$$$ Utilisation de l interpolation barycentrique ',
670	, ' pour CDSW $$$ '
671	WRITE(6,*) '-------------------------------------------------',
672	,'--------------'
673	CALL inter_barxy ( iml_phys,jml_phys -1,lon_rad,lat_rad ,
674	, var_ana, iml-1, jml-1, lon_in2, lat_in2, jml, tmp_var )
675	ELSE
676	CALL grille_m(iml_phys, jml_phys, lon_rad, lat_rad,
677	. var_ana, iml-1, jml, lon_in, lat_in, tmp_var )
678	ENDIF
679	c
680	CALL gr_int_dyn(tmp_var, qsol, iml-1, jml)
681	!
682	CALL flinclo(fid_phys)
683	!
684	END SUBROUTINE start_init_phys
685	!
686	!
687	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
688	!
689	!
690	SUBROUTINE startget_dyn(varname, iml, jml, lon_in, lat_in,
691	. lml, pls, workvar, champ, val_exp,jml2, lon_in2, lat_in2 ,
692	, interbar )
693	!
694	! ARGUMENTS
695	!
696	CHARACTER(), INTENT(in) :: varname
697	INTEGER, INTENT(in) :: iml, jml, lml, jml2
698	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
699	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
700	REAL, INTENT(in) :: pls(iml, jml, lml)
701	REAL, INTENT(in) :: workvar(iml, jml, lml)
702	REAL, INTENT(inout) :: champ(iml, jml, lml)
703	REAL, INTENT(in) :: val_exp
704	LOGICAL interbar
705	!
706	! LOCAL
707	!
708	INTEGER :: il, ij, ii
709	REAL :: xppn, xpps
710	!
711	#include "dimensions.h"
712	#include "paramet.h"
713	#include "comgeom2.h"
714	#include "comconst.h"
715	!
716	! This routine only works if the variable does not exist or is constant
717	!
718	IF ( MINVAL(champ(:,:,:)).EQ.MAXVAL(champ(:,:,:)) .AND.
719	. MINVAL(champ(:,:,:)).EQ.val_exp ) THEN
720	!
721	SELECTCASE(varname)
722	CASE ('u')
723	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
724	CALL start_init_dyn( iml, jml, lon_in, lat_in, jml2 ,
725	. lon_in2,lat_in2 , interbar )
726	ENDIF
727	CALL start_inter_3d('U', iml, jml, lml, lon_in,
728	. lat_in, jml2, lon_in2, lat_in2, pls, champ,interbar )
729	DO il=1,lml
730	DO ij=1,jml
731	DO ii=1,iml-1
732	champ(ii,ij,il) = champ(ii,ij,il) * cu(ii,ij)
733	ENDDO
734	champ(iml,ij, il) = champ(1,ij, il)
735	ENDDO
736	ENDDO
737	CASE ('v')
738	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
739	CALL start_init_dyn( iml, jml, lon_in, lat_in , jml2,
740	. lon_in2, lat_in2 , interbar )
741	ENDIF
742	CALL start_inter_3d('V', iml, jml, lml, lon_in,
743	. lat_in, jml2, lon_in2, lat_in2, pls, champ, interbar )
744	DO il=1,lml
745	DO ij=1,jml
746	DO ii=1,iml-1
747	champ(ii,ij,il) = champ(ii,ij,il) * cv(ii,ij)
748	ENDDO
749	champ(iml,ij, il) = champ(1,ij, il)
750	ENDDO
751	ENDDO
752	CASE ('t')
753	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
754	CALL start_init_dyn( iml, jml, lon_in, lat_in, jml2 ,
755	. lon_in2, lat_in2 ,interbar )
756	ENDIF
757	CALL start_inter_3d('TEMP', iml, jml, lml, lon_in,
758	. lat_in, jml2, lon_in2, lat_in2, pls, champ, interbar )
759
760	CASE ('tpot')
761	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
762	CALL start_init_dyn( iml, jml, lon_in, lat_in , jml2 ,
763	. lon_in2, lat_in2 , interbar )
764	ENDIF
765	CALL start_inter_3d('TEMP', iml, jml, lml, lon_in,
766	. lat_in, jml2, lon_in2, lat_in2, pls, champ, interbar )
767	IF ( MINVAL(workvar(:,:,:)) .NE. MAXVAL(workvar(:,:,:)) )
768	. THEN
769	DO il=1,lml
770	DO ij=1,jml
771	DO ii=1,iml-1
772	champ(ii,ij,il) = champ(ii,ij,il) * cpp
773	. / workvar(ii,ij,il)
774	ENDDO
775	champ(iml,ij,il) = champ(1,ij,il)
776	ENDDO
777	ENDDO
778	DO il=1,lml
779	xppn = SUM(aire(:,1)*champ(:,1,il))/apoln
780	xpps = SUM(aire(:,jml)*champ(:,jml,il))/apols
781	champ(:,1,il) = xppn
782	champ(:,jml,il) = xpps
783	ENDDO
784	ELSE
785	WRITE(,)'Could not compute potential temperature as the'
786	WRITE(,)'Exner function is missing or constant.'
787	STOP
788	ENDIF
789	CASE ('q')
790	IF ( .NOT.ALLOCATED(psol_dyn)) THEN
791	CALL start_init_dyn( iml, jml, lon_in, lat_in, jml2 ,
792	. lon_in2, lat_in2 , interbar )
793	ENDIF
794	CALL start_inter_3d('R', iml, jml, lml, lon_in, lat_in,
795	. jml2, lon_in2, lat_in2, pls, champ, interbar )
796	IF ( MINVAL(workvar(:,:,:)) .NE. MAXVAL(workvar(:,:,:)) )
797	. THEN
798	DO il=1,lml
799	DO ij=1,jml
800	DO ii=1,iml-1
801	champ(ii,ij,il) = 0.01 * champ(ii,ij,il) *
802	. workvar(ii,ij,il)
803	ENDDO
804	champ(iml,ij,il) = champ(1,ij,il)
805	ENDDO
806	ENDDO
807	WHERE ( champ .LT. 0.) champ = 1.0E-10
808	DO il=1,lml
809	xppn = SUM(aire(:,1)*champ(:,1,il))/apoln
810	xpps = SUM(aire(:,jml)*champ(:,jml,il))/apols
811	champ(:,1,il) = xppn
812	champ(:,jml,il) = xpps
813	ENDDO
814	ELSE
815	WRITE(,)'Could not compute specific humidity as the'
816	WRITE(,)'saturated humidity is missing or constant.'
817	STOP
818	ENDIF
819	CASE DEFAULT
820	WRITE(,) 'startget_dyn'
821	WRITE(,) 'No rule is present to extract variable ',
822	. varname(:LEN_TRIM(varname)),' from any data set'
823	STOP
824	END SELECT
825	ENDIF
826	END SUBROUTINE startget_dyn
827	!
828	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
829	!
830	SUBROUTINE start_init_dyn( iml, jml, lon_in, lat_in,jml2,lon_in2 ,
831	, lat_in2 , interbar )
832	!
833	INTEGER, INTENT(in) :: iml, jml, jml2
834	REAL, INTENT(in) :: lon_in(iml), lat_in(jml)
835	REAL, INTENT(in) :: lon_in2(iml), lat_in2(jml2)
836	LOGICAL interbar
837	!
838	! LOCAL
839	!
840	REAL :: lev(1), date, dt
841	INTEGER :: itau(1)
842	INTEGER :: i, j
843	integer :: iret
844	!
845	CHARACTER*120 :: physfname
846	LOGICAL :: check=.TRUE.
847	!
848	REAL, ALLOCATABLE :: lon_rad(:), lat_rad(:)
849	REAL, ALLOCATABLE :: lon_ini(:), lat_ini(:)
850	REAL, ALLOCATABLE :: var_ana(:,:), tmp_var(:,:), z(:,:)
851	REAL, ALLOCATABLE :: xppn(:), xpps(:)
852	LOGICAL :: allo
853	!
854	!
855	#include "dimensions.h"
856	#include "paramet.h"
857	#include "comgeom2.h"
858
859	CHARACTER*25 title
860
861	!
862	physfname = 'ECDYN.nc'
863	!
864	IF ( check ) WRITE(,) 'Opening the surface analysis'
865	!
866	CALL flininfo(physfname, iml_dyn, jml_dyn, llm_dyn,
867	. ttm_dyn, fid_dyn)
868	IF ( check ) WRITE(,) 'Values read: ', iml_dyn, jml_dyn,
869	. llm_dyn, ttm_dyn
870	!
871	ALLOCATE (lat_dyn(iml_dyn,jml_dyn), stat=iret)
872	ALLOCATE (lon_dyn(iml_dyn,jml_dyn), stat=iret)
873	ALLOCATE (levdyn_ini(llm_dyn), stat=iret)
874	!
875	CALL flinopen(physfname, .FALSE., iml_dyn, jml_dyn, llm_dyn,
876	. lon_dyn, lat_dyn, levdyn_ini, ttm_dyn,
877	. itau, date, dt, fid_dyn)
878	!
879
880	allo = allocated (var_ana)
881	if (allo) then
882	DEALLOCATE(var_ana, stat=iret)
883	endif
884	ALLOCATE(var_ana(iml_dyn, jml_dyn), stat=iret)
885
886	allo = allocated (lon_rad)
887	if (allo) then
888	DEALLOCATE(lon_rad, stat=iret)
889	endif
890
891	ALLOCATE(lon_rad(iml_dyn), stat=iret)
892	ALLOCATE(lon_ini(iml_dyn))
893
894	IF ( MAXVAL(lon_dyn(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
895	lon_ini(:) = lon_dyn(:,1) * 2.0 * ASIN(1.0) / 180.0
896	ELSE
897	lon_ini(:) = lon_dyn(:,1)
898	ENDIF
899
900	ALLOCATE(lat_rad(jml_dyn))
901	ALLOCATE(lat_ini(jml_dyn))
902
903	IF ( MAXVAL(lat_dyn(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
904	lat_ini(:) = lat_dyn(1,:) * 2.0 * ASIN(1.0) / 180.0
905	ELSE
906	lat_ini(:) = lat_dyn(1,:)
907	ENDIF
908	!
909
910
911	ALLOCATE(z(iml, jml))
912	ALLOCATE(tmp_var(iml-1,jml))
913	!
914	CALL flinget(fid_dyn, 'Z', iml_dyn, jml_dyn, 0, ttm_dyn,
915	. 1, 1, var_ana)
916	c
917	title='Z'
918	CALL conf_dat2d( title,iml_dyn, jml_dyn,lon_ini, lat_ini,
919	. lon_rad, lat_rad, var_ana, interbar )
920	c
921	IF ( interbar ) THEN
922	WRITE(6,*) '-------------------------------------------------',
923	,'--------------'
924	WRITE(6,*) '$$$ Utilisation de l interpolation barycentrique ',
925	, ' pour Z $$$ '
926	WRITE(6,*) '-------------------------------------------------',
927	,'--------------'
928	CALL inter_barxy ( iml_dyn,jml_dyn -1,lon_rad,lat_rad ,
929	, var_ana, iml-1, jml-1, lon_in2, lat_in2, jml, tmp_var)
930	ELSE
931	CALL grille_m(iml_dyn, jml_dyn , lon_rad, lat_rad, var_ana,
932	. iml-1, jml, lon_in, lat_in, tmp_var)
933	ENDIF
934
935	CALL gr_int_dyn(tmp_var, z, iml-1, jml)
936	!
937	ALLOCATE(psol_dyn(iml, jml))
938	!
939	CALL flinget(fid_dyn, 'SP', iml_dyn, jml_dyn, 0, ttm_dyn,
940	. 1, 1, var_ana)
941
942	title='SP'
943	CALL conf_dat2d( title,iml_dyn, jml_dyn,lon_ini, lat_ini,
944	. lon_rad, lat_rad, var_ana, interbar )
945
946	IF ( interbar ) THEN
947	WRITE(6,*) '-------------------------------------------------',
948	,'--------------'
949	WRITE(6,*) '$$$ Utilisation de l interpolation barycentrique ',
950	, ' pour SP $$$ '
951	WRITE(6,*) '-------------------------------------------------',
952	,'--------------'
953	CALL inter_barxy ( iml_dyn,jml_dyn -1,lon_rad,lat_rad ,
954	, var_ana, iml-1, jml-1, lon_in2, lat_in2, jml, tmp_var)
955	ELSE
956	CALL grille_m(iml_dyn, jml_dyn , lon_rad, lat_rad, var_ana,
957	. iml-1, jml, lon_in, lat_in, tmp_var )
958	ENDIF
959
960	CALL gr_int_dyn(tmp_var, psol_dyn, iml-1, jml)
961	!
962	IF ( .NOT.ALLOCATED(tsol)) THEN
963	! These variables may have been allocated by the need to
964	! create a start field for them or by the varibale
965	! coming out of the restart file. In case we dor have it we will initialize it.
966	!
967	CALL start_init_phys( iml, jml, lon_in, lat_in,jml2,lon_in2,
968	. lat_in2 , interbar )
969	ELSE
970	IF ( SIZE(tsol) .NE. SIZE(psol_dyn) ) THEN
971	WRITE(,) 'start_init_dyn :'
972	WRITE(,) 'The temperature field we have does not ',
973	. 'have the right size'
974	STOP
975	ENDIF
976	ENDIF
977	IF ( .NOT.ALLOCATED(phis)) THEN
978	!
979	! These variables may have been allocated by the need to create a start field for them or by the varibale
980	! coming out of the restart file. In case we dor have it we will initialize it.
981	!
982	CALL start_init_orog( iml, jml, lon_in, lat_in, jml2, lon_in2 ,
983	. lat_in2 , interbar )
984	!
985	ELSE
986	!
987	IF (SIZE(phis) .NE. SIZE(psol_dyn)) THEN
988	!
989	WRITE(,) 'start_init_dyn :'
990	WRITE(,) 'The orography field we have does not ',
991	. ' have the right size'
992	STOP
993	ENDIF
994	!
995	ENDIF
996	!
997	! PSOL is computed in Pascals
998	!
999	!
1000	DO j = 1, jml
1001	DO i = 1, iml-1
1002	psol_dyn(i,j) = psol_dyn(i,j)*(1.0+(z(i,j)-phis(i,j))
1003	. /287.0/tsol(i,j))
1004	ENDDO
1005	psol_dyn(iml,j) = psol_dyn(1,j)
1006	ENDDO
1007	!
1008	!
1009	ALLOCATE(xppn(iml-1))
1010	ALLOCATE(xpps(iml-1))
1011	!
1012	DO i = 1, iml-1
1013	xppn(i) = aire( i,1) * psol_dyn( i,1)
1014	xpps(i) = aire( i,jml) * psol_dyn( i,jml)
1015	ENDDO
1016	!
1017	DO i = 1, iml
1018	psol_dyn(i,1 ) = SUM(xppn)/apoln
1019	psol_dyn(i,jml) = SUM(xpps)/apols
1020	ENDDO
1021	!
1022	RETURN
1023	!
1024	END SUBROUTINE start_init_dyn
1025	!
1026	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1027	!
1028	SUBROUTINE start_inter_3d(varname, iml, jml, lml, lon_in,
1029	. lat_in, jml2, lon_in2, lat_in2, pls_in, var3d, interbar )
1030	!
1031	! This subroutine gets a variables from a 3D file and does the interpolations needed
1032	!
1033	!
1034	! ARGUMENTS
1035	!
1036	CHARACTER() :: varname
1037	INTEGER :: iml, jml, lml, jml2
1038	REAL :: lon_in(iml), lat_in(jml), pls_in(iml, jml, lml)
1039	REAL :: lon_in2(iml) , lat_in2(jml2)
1040	REAL :: var3d(iml, jml, lml)
1041	LOGICAL interbar
1042	real chmin,chmax
1043	!
1044	! LOCAL
1045	!
1046	CHARACTER*25 title
1047	INTEGER :: ii, ij, il, jsort,i,j,l
1048	REAL :: bx, by
1049	REAL, ALLOCATABLE :: lon_rad(:), lat_rad(:)
1050	REAL, ALLOCATABLE :: lon_ini(:), lat_ini(:) , lev_dyn(:)
1051	REAL, ALLOCATABLE :: var_tmp2d(:,:), var_tmp3d(:,:,:)
1052	REAL, ALLOCATABLE :: ax(:), ay(:), yder(:)
1053	REAL, ALLOCATABLE :: varrr(:,:,:)
1054	INTEGER, ALLOCATABLE :: lind(:)
1055	!
1056	LOGICAL :: check = .TRUE.
1057	!
1058	IF ( .NOT. ALLOCATED(var_ana3d)) THEN
1059	ALLOCATE(var_ana3d(iml_dyn, jml_dyn, llm_dyn))
1060	ENDIF
1061	ALLOCATE(varrr(iml_dyn, jml_dyn, llm_dyn))
1062	!
1063	!
1064	IF ( check) WRITE(,) 'Going into flinget to extract the 3D ',
1065	. ' field.', fid_dyn
1066	IF ( check) WRITE(,) fid_dyn, varname, iml_dyn, jml_dyn,
1067	. llm_dyn,ttm_dyn
1068	!
1069	CALL flinget(fid_dyn, varname, iml_dyn, jml_dyn, llm_dyn,
1070	. ttm_dyn, 1, 1, var_ana3d)
1071	!
1072	IF ( check) WRITE(,) 'Allocating space for the interpolation',
1073	. iml, jml, llm_dyn
1074	!
1075	ALLOCATE(lon_rad(iml_dyn))
1076	ALLOCATE(lon_ini(iml_dyn))
1077
1078	IF ( MAXVAL(lon_dyn(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
1079	lon_ini(:) = lon_dyn(:,1) * 2.0 * ASIN(1.0) / 180.0
1080	ELSE
1081	lon_ini(:) = lon_dyn(:,1)
1082	ENDIF
1083
1084	ALLOCATE(lat_rad(jml_dyn))
1085	ALLOCATE(lat_ini(jml_dyn))
1086
1087	ALLOCATE(lev_dyn(llm_dyn))
1088
1089	IF ( MAXVAL(lat_dyn(:,:)) .GT. 2.0 * ASIN(1.0) ) THEN
1090	lat_ini(:) = lat_dyn(1,:) * 2.0 * ASIN(1.0) / 180.0
1091	ELSE
1092	lat_ini(:) = lat_dyn(1,:)
1093	ENDIF
1094	!
1095
1096	CALL conf_dat3d ( varname,iml_dyn, jml_dyn, llm_dyn, lon_ini,
1097	. lat_ini, levdyn_ini, lon_rad, lat_rad, lev_dyn, var_ana3d ,
1098	, interbar )
1099
1100	ALLOCATE(var_tmp2d(iml-1, jml))
1101	ALLOCATE(var_tmp3d(iml, jml, llm_dyn))
1102	ALLOCATE(ax(llm_dyn))
1103	ALLOCATE(ay(llm_dyn))
1104	ALLOCATE(yder(llm_dyn))
1105	ALLOCATE(lind(llm_dyn))
1106	!
1107
1108	DO il=1,llm_dyn
1109	!
1110	IF( interbar ) THEN
1111	IF( il.EQ.1 ) THEN
1112	WRITE(6,*) '-------------------------------------------------',
1113	,'--------------'
1114	WRITE(6,*) '$$$ Utilisation de l interpolation barycentrique ',
1115	, ' pour ', varname
1116	WRITE(6,*) '-------------------------------------------------',
1117	,'--------------'
1118	ENDIF
1119	CALL inter_barxy ( iml_dyn, jml_dyn -1,lon_rad, lat_rad,
1120	, var_ana3d(:,:,il),iml-1, jml2, lon_in2, lat_in2,jml,var_tmp2d )
1121	ELSE
1122	CALL grille_m(iml_dyn, jml_dyn, lon_rad, lat_rad,
1123	. var_ana3d(:,:,il), iml-1, jml, lon_in, lat_in, var_tmp2d )
1124	ENDIF
1125	!
1126	CALL gr_int_dyn(var_tmp2d, var_tmp3d(:,:,il), iml-1, jml)
1127	!
1128	ENDDO
1129	!
1130	DO il=1,llm_dyn
1131	lind(il) = llm_dyn-il+1
1132	ENDDO
1133	!
1134	c
1135	c ... Pour l'interpolation verticale ,on interpole du haut de l'atmosphere
1136	c vers le sol ...
1137	c
1138	DO ij=1,jml
1139	DO ii=1,iml-1
1140	!
1141	ax(:) = lev_dyn(lind(:))
1142	ay(:) = var_tmp3d(ii, ij, lind(:))
1143	!
1144
1145	CALL SPLINE(ax, ay, llm_dyn, 1.e30, 1.e30, yder)
1146	!
1147	DO il=1,lml
1148	bx = pls_in(ii, ij, il)
1149	CALL SPLINT(ax, ay, yder, llm_dyn, bx, by)
1150	var3d(ii, ij, il) = by
1151	ENDDO
1152	!
1153	ENDDO
1154	var3d(iml, ij, :) = var3d(1, ij, :)
1155	ENDDO
1156
1157	do il=1,lml
1158	call minmax(iml*jml,var3d(1,1,il),chmin,chmax)
1159	SELECTCASE(varname)
1160	CASE('U')
1161	WRITE(,) ' U min max l ',il,chmin,chmax
1162	CASE('V')
1163	WRITE(,) ' V min max l ',il,chmin,chmax
1164	CASE('TEMP')
1165	WRITE(,) ' TEMP min max l ',il,chmin,chmax
1166	CASE('R')
1167	WRITE(,) ' R min max l ',il,chmin,chmax
1168	END SELECT
1169	enddo
1170
1171	DEALLOCATE(lon_rad)
1172	DEALLOCATE(lat_rad)
1173	DEALLOCATE(var_tmp2d)
1174	DEALLOCATE(var_tmp3d)
1175	DEALLOCATE(ax)
1176	DEALLOCATE(ay)
1177	DEALLOCATE(yder)
1178	DEALLOCATE(lind)
1179
1180	!
1181	RETURN
1182	!
1183	END SUBROUTINE start_inter_3d
1184	!
1185	END MODULE startvar

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

Download in other formats:

Original Format