Context Navigation

dynredem.F @ 1033

Last change on this file since 1033 was 960, checked in by lsce, 16 years ago

Ajoute du parametre config_inca dans conf_gcm.F config_inca='none'(sans INCA, par defaut) config_inca='chem'(avec INCA config chemie) config_inca='aero'(avec INCA config aerosol)
Menage parmis les cles CPP INCA
Enleve le calcul d'omega dans calfis.F et active le calcul correspondant dans advtrac.F(avant uniquement pour INCA).

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 20.7 KB

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

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

Download in other formats:

Original Format