Context Navigation

dynredem_p.F @ 2598

Last change on this file since 2598 was 2598, checked in by Ehouarn Millour, 8 years ago
Cleanup in the dynamics: turn serre.h into module serre_mod.F90 EM
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 19.5 KB

Line
1	!
2	! $Id: dynredem_p.F 2598 2016-07-22 09:28:39Z emillour $
3	!
4	c
5	SUBROUTINE dynredem0_p(fichnom,iday_end,phis)
6	#ifdef CPP_IOIPSL
7	USE IOIPSL
8	#endif
9	USE parallel_lmdz
10	USE infotrac
11	use netcdf95, only: NF95_PUT_VAR
12	USE comconst_mod, ONLY: cpp, daysec, dtvr, g, kappa, omeg, rad
13	USE serre_mod, ONLY: clon,clat,grossismx,grossismy,dzoomx,dzoomy,
14	& taux,tauy
15
16	IMPLICIT NONE
17	c=======================================================================
18	c Ecriture du fichier de redemarrage sous format NetCDF (initialisation)
19	c=======================================================================
20	c Declarations:
21	c -------------
22	#include "dimensions.h"
23	#include "paramet.h"
24	#include "comvert.h"
25	#include "comgeom2.h"
26	#include "temps.h"
27	#include "ener.h"
28	#include "logic.h"
29	#include "netcdf.inc"
30	#include "description.h"
31
32	c Arguments:
33	c ----------
34	INTEGER iday_end
35	REAL phis(iip1, jjp1)
36	CHARACTER() fichnom
37
38	c Local:
39	c ------
40	INTEGER iq,l
41	INTEGER length
42	PARAMETER (length = 100)
43	REAL tab_cntrl(length) ! tableau des parametres du run
44	INTEGER ierr
45	character*20 modname
46	character*80 abort_message
47
48	c Variables locales pour NetCDF:
49	c
50	INTEGER dims2(2), dims3(3), dims4(4)
51	INTEGER idim_index
52	INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
53	INTEGER idim_s, idim_sig
54	INTEGER idim_tim
55	INTEGER nid,nvarid
56
57	REAL zan0,zjulian,hours
58	INTEGER yyears0,jjour0, mmois0
59	character*30 unites
60
61
62	c-----------------------------------------------------------------------
63	if (mpi_rank==0) then
64
65	modname='dynredem0_p'
66
67	#ifdef CPP_IOIPSL
68	call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
69	call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
70	#else
71	! set yyears0, mmois0, jjour0 to 0,1,1 (hours is not used)
72	yyears0=0
73	mmois0=1
74	jjour0=1
75	#endif
76
77	DO l=1,length
78	tab_cntrl(l) = 0.
79	ENDDO
80	tab_cntrl(1) = REAL(iim)
81	tab_cntrl(2) = REAL(jjm)
82	tab_cntrl(3) = REAL(llm)
83	tab_cntrl(4) = REAL(day_ref)
84	tab_cntrl(5) = REAL(annee_ref)
85	tab_cntrl(6) = rad
86	tab_cntrl(7) = omeg
87	tab_cntrl(8) = g
88	tab_cntrl(9) = cpp
89	tab_cntrl(10) = kappa
90	tab_cntrl(11) = daysec
91	tab_cntrl(12) = dtvr
92	tab_cntrl(13) = etot0
93	tab_cntrl(14) = ptot0
94	tab_cntrl(15) = ztot0
95	tab_cntrl(16) = stot0
96	tab_cntrl(17) = ang0
97	tab_cntrl(18) = pa
98	tab_cntrl(19) = preff
99	c
100	c ..... parametres pour le zoom ......
101
102	tab_cntrl(20) = clon
103	tab_cntrl(21) = clat
104	tab_cntrl(22) = grossismx
105	tab_cntrl(23) = grossismy
106	c
107	IF ( fxyhypb ) THEN
108	tab_cntrl(24) = 1.
109	tab_cntrl(25) = dzoomx
110	tab_cntrl(26) = dzoomy
111	tab_cntrl(27) = 0.
112	tab_cntrl(28) = taux
113	tab_cntrl(29) = tauy
114	ELSE
115	tab_cntrl(24) = 0.
116	tab_cntrl(25) = dzoomx
117	tab_cntrl(26) = dzoomy
118	tab_cntrl(27) = 0.
119	tab_cntrl(28) = 0.
120	tab_cntrl(29) = 0.
121	IF( ysinus ) tab_cntrl(27) = 1.
122	ENDIF
123
124	tab_cntrl(30) = REAL(iday_end)
125	tab_cntrl(31) = REAL(itau_dyn + itaufin)
126	c start_time: start_time of simulation (not necessarily 0.)
127	tab_cntrl(32) = start_time
128	c
129	c .........................................................
130	c
131	c Creation du fichier:
132	c
133	ierr = NF_CREATE(fichnom, NF_CLOBBER, nid)
134	IF (ierr.NE.NF_NOERR) THEN
135	WRITE(6,*)" Pb d ouverture du fichier "//fichnom
136	WRITE(6,*)' ierr = ', ierr
137	CALL ABORT
138	ENDIF
139	c
140	c Preciser quelques attributs globaux:
141	c
142	ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 27,
143	. "Fichier demmarage dynamique")
144	c
145	c Definir les dimensions du fichiers:
146	c
147	ierr = NF_DEF_DIM (nid, "index", length, idim_index)
148	ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
149	ierr = NF_DEF_DIM (nid, "rlatu", jjp1, idim_rlatu)
150	ierr = NF_DEF_DIM (nid, "rlonv", iip1, idim_rlonv)
151	ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
152	ierr = NF_DEF_DIM (nid, "sigs", llm, idim_s)
153	ierr = NF_DEF_DIM (nid, "sig", llmp1, idim_sig)
154	ierr = NF_DEF_DIM (nid, "temps", NF_UNLIMITED, idim_tim)
155	c
156	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
157	c
158	c Definir et enregistrer certains champs invariants:
159	c
160	ierr = NF_REDEF (nid)
161	cIM 220306 BEG
162	#ifdef NC_DOUBLE
163	ierr = NF_DEF_VAR (nid,"controle",NF_DOUBLE,1,idim_index,nvarid)
164	#else
165	ierr = NF_DEF_VAR (nid,"controle",NF_FLOAT,1,idim_index,nvarid)
166	#endif
167	cIM 220306 END
168	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
169	. "Parametres de controle")
170	ierr = NF_ENDDEF(nid)
171	call NF95_PUT_VAR(nid,nvarid,tab_cntrl)
172	c
173	ierr = NF_REDEF (nid)
174	cIM 220306 BEG
175	#ifdef NC_DOUBLE
176	ierr = NF_DEF_VAR (nid,"rlonu",NF_DOUBLE,1,idim_rlonu,nvarid)
177	#else
178	ierr = NF_DEF_VAR (nid,"rlonu",NF_FLOAT,1,idim_rlonu,nvarid)
179	#endif
180	cIM 220306 END
181	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
182	. "Longitudes des points U")
183	ierr = NF_ENDDEF(nid)
184	call NF95_PUT_VAR(nid,nvarid,rlonu)
185	c
186	ierr = NF_REDEF (nid)
187	cIM 220306 BEG
188	#ifdef NC_DOUBLE
189	ierr = NF_DEF_VAR (nid,"rlatu",NF_DOUBLE,1,idim_rlatu,nvarid)
190	#else
191	ierr = NF_DEF_VAR (nid,"rlatu",NF_FLOAT,1,idim_rlatu,nvarid)
192	#endif
193	cIM 220306 END
194	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
195	. "Latitudes des points U")
196	ierr = NF_ENDDEF(nid)
197	call NF95_PUT_VAR (nid,nvarid,rlatu)
198	c
199	ierr = NF_REDEF (nid)
200	cIM 220306 BEG
201	#ifdef NC_DOUBLE
202	ierr = NF_DEF_VAR (nid,"rlonv",NF_DOUBLE,1,idim_rlonv,nvarid)
203	#else
204	ierr = NF_DEF_VAR (nid,"rlonv",NF_FLOAT,1,idim_rlonv,nvarid)
205	#endif
206	cIM 220306 END
207	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
208	. "Longitudes des points V")
209	ierr = NF_ENDDEF(nid)
210	call NF95_PUT_VAR(nid,nvarid,rlonv)
211	c
212	ierr = NF_REDEF (nid)
213	cIM 220306 BEG
214	#ifdef NC_DOUBLE
215	ierr = NF_DEF_VAR (nid,"rlatv",NF_DOUBLE,1,idim_rlatv,nvarid)
216	#else
217	ierr = NF_DEF_VAR (nid,"rlatv",NF_FLOAT,1,idim_rlatv,nvarid)
218	#endif
219	cIM 220306 END
220	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
221	. "Latitudes des points V")
222	ierr = NF_ENDDEF(nid)
223	call NF95_PUT_VAR(nid,nvarid,rlatv)
224	c
225	ierr = NF_REDEF (nid)
226	cIM 220306 BEG
227	#ifdef NC_DOUBLE
228	ierr = NF_DEF_VAR (nid,"nivsigs",NF_DOUBLE,1,idim_s,nvarid)
229	#else
230	ierr = NF_DEF_VAR (nid,"nivsigs",NF_FLOAT,1,idim_s,nvarid)
231	#endif
232	cIM 220306 END
233	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 28,
234	. "Numero naturel des couches s")
235	ierr = NF_ENDDEF(nid)
236	call NF95_PUT_VAR(nid,nvarid,nivsigs)
237	c
238	ierr = NF_REDEF (nid)
239	cIM 220306 BEG
240	#ifdef NC_DOUBLE
241	ierr = NF_DEF_VAR (nid,"nivsig",NF_DOUBLE,1,idim_sig,nvarid)
242	#else
243	ierr = NF_DEF_VAR (nid,"nivsig",NF_FLOAT,1,idim_sig,nvarid)
244	#endif
245	cIM 220306 END
246	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 32,
247	. "Numero naturel des couches sigma")
248	ierr = NF_ENDDEF(nid)
249	call NF95_PUT_VAR(nid,nvarid,nivsig)
250	c
251	ierr = NF_REDEF (nid)
252	cIM 220306 BEG
253	#ifdef NC_DOUBLE
254	ierr = NF_DEF_VAR (nid,"ap",NF_DOUBLE,1,idim_sig,nvarid)
255	#else
256	ierr = NF_DEF_VAR (nid,"ap",NF_FLOAT,1,idim_sig,nvarid)
257	#endif
258	cIM 220306 END
259	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
260	. "Coefficient A pour hybride")
261	ierr = NF_ENDDEF(nid)
262	call NF95_PUT_VAR(nid,nvarid,ap)
263	c
264	ierr = NF_REDEF (nid)
265	cIM 220306 BEG
266	#ifdef NC_DOUBLE
267	ierr = NF_DEF_VAR (nid,"bp",NF_DOUBLE,1,idim_sig,nvarid)
268	#else
269	ierr = NF_DEF_VAR (nid,"bp",NF_FLOAT,1,idim_sig,nvarid)
270	#endif
271	cIM 220306 END
272	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
273	. "Coefficient B pour hybride")
274	ierr = NF_ENDDEF(nid)
275	call NF95_PUT_VAR(nid,nvarid,bp)
276	c
277	ierr = NF_REDEF (nid)
278	cIM 220306 BEG
279	#ifdef NC_DOUBLE
280	ierr = NF_DEF_VAR (nid,"presnivs",NF_DOUBLE,1,idim_s,nvarid)
281	#else
282	ierr = NF_DEF_VAR (nid,"presnivs",NF_FLOAT,1,idim_s,nvarid)
283	#endif
284	cIM 220306 END
285	ierr = NF_ENDDEF(nid)
286	call NF95_PUT_VAR(nid,nvarid,presnivs)
287	c
288	c Coefficients de passage cov. <-> contra. <--> naturel
289	c
290	ierr = NF_REDEF (nid)
291	dims2(1) = idim_rlonu
292	dims2(2) = idim_rlatu
293	cIM 220306 BEG
294	#ifdef NC_DOUBLE
295	ierr = NF_DEF_VAR (nid,"cu",NF_DOUBLE,2,dims2,nvarid)
296	#else
297	ierr = NF_DEF_VAR (nid,"cu",NF_FLOAT,2,dims2,nvarid)
298	#endif
299	cIM 220306 END
300	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
301	. "Coefficient de passage pour U")
302	ierr = NF_ENDDEF(nid)
303	call NF95_PUT_VAR(nid,nvarid,cu)
304	c
305	ierr = NF_REDEF (nid)
306	dims2(1) = idim_rlonv
307	dims2(2) = idim_rlatv
308	cIM 220306 BEG
309	#ifdef NC_DOUBLE
310	ierr = NF_DEF_VAR (nid,"cv",NF_DOUBLE,2,dims2,nvarid)
311	#else
312	ierr = NF_DEF_VAR (nid,"cv",NF_FLOAT,2,dims2,nvarid)
313	#endif
314	cIM 220306 END
315	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
316	. "Coefficient de passage pour V")
317	ierr = NF_ENDDEF(nid)
318	call NF95_PUT_VAR(nid,nvarid,cv)
319	c
320	c Aire de chaque maille:
321	c
322	ierr = NF_REDEF (nid)
323	dims2(1) = idim_rlonv
324	dims2(2) = idim_rlatu
325	cIM 220306 BEG
326	#ifdef NC_DOUBLE
327	ierr = NF_DEF_VAR (nid,"aire",NF_DOUBLE,2,dims2,nvarid)
328	#else
329	ierr = NF_DEF_VAR (nid,"aire",NF_FLOAT,2,dims2,nvarid)
330	#endif
331	cIM 220306 END
332	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
333	. "Aires de chaque maille")
334	ierr = NF_ENDDEF(nid)
335	call NF95_PUT_VAR(nid,nvarid,aire)
336	c
337	c Geopentiel au sol:
338	c
339	ierr = NF_REDEF (nid)
340	dims2(1) = idim_rlonv
341	dims2(2) = idim_rlatu
342	cIM 220306 BEG
343	#ifdef NC_DOUBLE
344	ierr = NF_DEF_VAR (nid,"phisinit",NF_DOUBLE,2,dims2,nvarid)
345	#else
346	ierr = NF_DEF_VAR (nid,"phisinit",NF_FLOAT,2,dims2,nvarid)
347	#endif
348	cIM 220306 END
349	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
350	. "Geopotentiel au sol")
351	ierr = NF_ENDDEF(nid)
352	call NF95_PUT_VAR(nid,nvarid,phis)
353	c
354	c Definir les variables pour pouvoir les enregistrer plus tard:
355	c
356	ierr = NF_REDEF (nid) ! entrer dans le mode de definition
357	c
358	cIM 220306 BEG
359	#ifdef NC_DOUBLE
360	ierr = NF_DEF_VAR (nid,"temps",NF_DOUBLE,1,idim_tim,nvarid)
361	#else
362	ierr = NF_DEF_VAR (nid,"temps",NF_FLOAT,1,idim_tim,nvarid)
363	#endif
364	cIM 220306 END
365	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
366	. "Temps de simulation")
367	write(unites,200)yyears0,mmois0,jjour0
368	200 format('days since ',i4,'-',i2.2,'-',i2.2,' 00:00:00')
369	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "units", 30,
370	. unites)
371
372	c
373	dims4(1) = idim_rlonu
374	dims4(2) = idim_rlatu
375	dims4(3) = idim_s
376	dims4(4) = idim_tim
377	cIM 220306 BEG
378	#ifdef NC_DOUBLE
379	ierr = NF_DEF_VAR (nid,"ucov",NF_DOUBLE,4,dims4,nvarid)
380	#else
381	ierr = NF_DEF_VAR (nid,"ucov",NF_FLOAT,4,dims4,nvarid)
382	#endif
383	cIM 220306 END
384	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
385	. "Vitesse U")
386	c
387	dims4(1) = idim_rlonv
388	dims4(2) = idim_rlatv
389	dims4(3) = idim_s
390	dims4(4) = idim_tim
391	cIM 220306 BEG
392	#ifdef NC_DOUBLE
393	ierr = NF_DEF_VAR (nid,"vcov",NF_DOUBLE,4,dims4,nvarid)
394	#else
395	ierr = NF_DEF_VAR (nid,"vcov",NF_FLOAT,4,dims4,nvarid)
396	#endif
397	cIM 220306 END
398	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
399	. "Vitesse V")
400	c
401	dims4(1) = idim_rlonv
402	dims4(2) = idim_rlatu
403	dims4(3) = idim_s
404	dims4(4) = idim_tim
405	cIM 220306 BEG
406	#ifdef NC_DOUBLE
407	ierr = NF_DEF_VAR (nid,"teta",NF_DOUBLE,4,dims4,nvarid)
408	#else
409	ierr = NF_DEF_VAR (nid,"teta",NF_FLOAT,4,dims4,nvarid)
410	#endif
411	cIM 220306 END
412	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 11,
413	. "Temperature")
414	c
415	dims4(1) = idim_rlonv
416	dims4(2) = idim_rlatu
417	dims4(3) = idim_s
418	dims4(4) = idim_tim
419
420	DO iq=1,nqtot
421	cIM 220306 BEG
422	#ifdef NC_DOUBLE
423	ierr = NF_DEF_VAR (nid,tname(iq),NF_DOUBLE,4,dims4,nvarid)
424	#else
425	ierr = NF_DEF_VAR (nid,tname(iq),NF_FLOAT,4,dims4,nvarid)
426	#endif
427	cIM 220306 END
428	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))
429	ENDDO
430	c
431	dims4(1) = idim_rlonv
432	dims4(2) = idim_rlatu
433	dims4(3) = idim_s
434	dims4(4) = idim_tim
435	cIM 220306 BEG
436	#ifdef NC_DOUBLE
437	ierr = NF_DEF_VAR (nid,"masse",NF_DOUBLE,4,dims4,nvarid)
438	#else
439	ierr = NF_DEF_VAR (nid,"masse",NF_FLOAT,4,dims4,nvarid)
440	#endif
441	cIM 220306 END
442	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,
443	. "C est quoi ?")
444	c
445	dims3(1) = idim_rlonv
446	dims3(2) = idim_rlatu
447	dims3(3) = idim_tim
448	cIM 220306 BEG
449	#ifdef NC_DOUBLE
450	ierr = NF_DEF_VAR (nid,"ps",NF_DOUBLE,3,dims3,nvarid)
451	#else
452	ierr = NF_DEF_VAR (nid,"ps",NF_FLOAT,3,dims3,nvarid)
453	#endif
454	cIM 220306 END
455	ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 15,
456	. "Pression au sol")
457	c
458	ierr = NF_ENDDEF(nid) ! sortir du mode de definition
459	ierr = NF_CLOSE(nid) ! fermer le fichier
460
461	PRINT*,'iim,jjm,llm,iday_end',iim,jjm,llm,iday_end
462	PRINT*,'rad,omeg,g,cpp,kappa',
463	, rad,omeg,g,cpp,kappa
464
465	endif ! mpi_rank==0
466	RETURN
467	END
468	SUBROUTINE dynredem1_p(fichnom,time,
469	. vcov,ucov,teta,q,masse,ps)
470	USE parallel_lmdz
471	USE infotrac
472	USE control_mod
473	use netcdf, only: NF90_get_VAR
474	use netcdf95, only: NF95_PUT_VAR
475
476	IMPLICIT NONE
477	c=================================================================
478	c Ecriture du fichier de redemarrage sous format NetCDF
479	c=================================================================
480	#include "dimensions.h"
481	#include "paramet.h"
482	#include "description.h"
483	#include "netcdf.inc"
484	#include "comvert.h"
485	#include "comgeom.h"
486	#include "temps.h"
487
488
489	INTEGER l
490	REAL vcov(iip1,jjm,llm),ucov(iip1, jjp1,llm)
491	REAL teta(iip1, jjp1,llm)
492	REAL ps(iip1, jjp1),masse(iip1, jjp1,llm)
493	REAL q(iip1, jjp1, llm, nqtot)
494	CHARACTER() fichnom
495
496	REAL time
497	INTEGER nid, nvarid, nid_trac, nvarid_trac
498	REAL trac_tmp(ip1jmp1,llm)
499	INTEGER ierr, ierr_file
500	INTEGER iq
501	INTEGER length
502	PARAMETER (length = 100)
503	REAL tab_cntrl(length) ! tableau des parametres du run
504	character*20 modname
505	character*80 abort_message
506	c
507	INTEGER nb
508	SAVE nb
509	DATA nb / 0 /
510
511	logical exist_file
512
513	call Gather_Field(ucov,ip1jmp1,llm,0)
514	call Gather_Field(vcov,ip1jm,llm,0)
515	call Gather_Field(teta,ip1jmp1,llm,0)
516	call Gather_Field(masse,ip1jmp1,llm,0)
517	call Gather_Field(ps,ip1jmp1,1,0)
518
519	do iq=1,nqtot
520	call Gather_Field(q(:,:,:,iq),ip1jmp1,llm,0)
521	enddo
522
523
524	if (mpi_rank==0) then
525
526	modname = 'dynredem1'
527	ierr = NF_OPEN(fichnom, NF_WRITE, nid)
528	IF (ierr .NE. NF_NOERR) THEN
529	PRINT*, "Pb. d ouverture "//fichnom
530	CALL abort
531	ENDIF
532
533	c Ecriture/extension de la coordonnee temps
534
535	nb = nb + 1
536	ierr = NF_INQ_VARID(nid, "temps", nvarid)
537	IF (ierr .NE. NF_NOERR) THEN
538	print *, NF_STRERROR(ierr)
539	abort_message='Variable temps n est pas definie'
540	CALL abort_gcm(modname,abort_message,ierr)
541	ENDIF
542	call NF95_PUT_VAR(nid,nvarid,time,start=(/nb/))
543	PRINT*, "Enregistrement pour ", nb, time
544
545	c
546	c Re-ecriture du tableau de controle, itaufin n'est plus defini quand
547	c on passe dans dynredem0
548	ierr = NF_INQ_VARID (nid, "controle", nvarid)
549	IF (ierr .NE. NF_NOERR) THEN
550	abort_message="dynredem1: Le champ <controle> est absent"
551	ierr = 1
552	CALL abort_gcm(modname,abort_message,ierr)
553	ENDIF
554	ierr = NF90_GET_VAR(nid, nvarid, tab_cntrl)
555	tab_cntrl(31) = REAL(itau_dyn + itaufin)
556	call NF95_PUT_VAR(nid,nvarid,tab_cntrl)
557
558	c Ecriture des champs
559	c
560	ierr = NF_INQ_VARID(nid, "ucov", nvarid)
561	IF (ierr .NE. NF_NOERR) THEN
562	PRINT*, "Variable ucov n est pas definie"
563	CALL abort
564	ENDIF
565	call NF95_PUT_VAR(nid,nvarid,ucov)
566
567	ierr = NF_INQ_VARID(nid, "vcov", nvarid)
568	IF (ierr .NE. NF_NOERR) THEN
569	PRINT*, "Variable vcov n est pas definie"
570	CALL abort
571	ENDIF
572	call NF95_PUT_VAR(nid,nvarid,vcov)
573
574	ierr = NF_INQ_VARID(nid, "teta", nvarid)
575	IF (ierr .NE. NF_NOERR) THEN
576	PRINT*, "Variable teta n est pas definie"
577	CALL abort
578	ENDIF
579	call NF95_PUT_VAR(nid,nvarid,teta)
580
581	IF (type_trac == 'inca') THEN
582	! Ajout Anne pour lecture valeurs traceurs dans un fichier start_trac.nc
583	inquire(FILE="start_trac.nc", EXIST=exist_file)
584	print *, "EXIST", exist_file
585	if (exist_file) then
586	ierr_file = NF_OPEN ("start_trac.nc", NF_NOWRITE,nid_trac)
587	IF (ierr_file .NE.NF_NOERR) THEN
588	write(6,*)' Pb d''ouverture du fichier start_trac.nc'
589	write(6,*)' ierr = ', ierr_file
590	ENDIF
591	else
592	ierr_file = 2
593	endif
594	END IF
595
596	do iq=1,nqtot
597
598	IF (type_trac /= 'inca') THEN
599	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
600	IF (ierr .NE. NF_NOERR) THEN
601	PRINT*, "Variable tname(iq) n est pas definie"
602	CALL abort
603	ENDIF
604	call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq))
605	ELSE ! type_trac = inca
606	! lecture de la valeur du traceur dans start_trac.nc
607	IF (ierr_file .ne. 2) THEN
608	ierr = NF_INQ_VARID (nid_trac, tname(iq), nvarid_trac)
609	IF (ierr .NE. NF_NOERR) THEN
610	PRINT*, tname(iq),"est absent de start_trac.nc"
611	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
612	IF (ierr .NE. NF_NOERR) THEN
613	PRINT*, "Variable ", tname(iq)," n est pas definie"
614	CALL abort
615	ENDIF
616	call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq))
617
618	ELSE
619	PRINT*, tname(iq), "est present dans start_trac.nc"
620	ierr = NF90_GET_VAR(nid_trac, nvarid_trac, trac_tmp)
621	IF (ierr .NE. NF_NOERR) THEN
622	PRINT*, "Lecture echouee pour", tname(iq)
623	CALL abort
624	ENDIF
625	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
626	IF (ierr .NE. NF_NOERR) THEN
627	PRINT*, "Variable ", tname(iq)," n est pas definie"
628	CALL abort
629	ENDIF
630	call NF95_PUT_VAR(nid, nvarid, trac_tmp)
631
632	ENDIF ! IF (ierr .NE. NF_NOERR)
633	! fin lecture du traceur
634	ELSE ! si il n'y a pas de fichier start_trac.nc
635	! print *, 'il n y a pas de fichier start_trac'
636	ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
637	IF (ierr .NE. NF_NOERR) THEN
638	PRINT*, "Variable tname(iq) n est pas definie"
639	CALL abort
640	ENDIF
641	call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq))
642	ENDIF ! (ierr_file .ne. 2)
643	END IF !type_trac
644
645	ENDDO
646
647
648
649	c
650	ierr = NF_INQ_VARID(nid, "masse", nvarid)
651	IF (ierr .NE. NF_NOERR) THEN
652	PRINT*, "Variable masse n est pas definie"
653	CALL abort
654	ENDIF
655	call NF95_PUT_VAR(nid,nvarid,masse)
656	c
657	ierr = NF_INQ_VARID(nid, "ps", nvarid)
658	IF (ierr .NE. NF_NOERR) THEN
659	PRINT*, "Variable ps n est pas definie"
660	CALL abort
661	ENDIF
662	call NF95_PUT_VAR(nid,nvarid,ps)
663
664	ierr = NF_CLOSE(nid)
665	c
666	endif ! mpi_rank==0
667
668	RETURN
669	END
670

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

Download in other formats:

Original Format