Context Navigation

dynredem_p.F @ 5054

Last change on this file since 5054 was 2641, checked in by Laurent Fairhead, 8 years ago
Merged trunk changes r2593:2640 into testing branch
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.8 KB

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

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

Download in other formats:

Original Format