source: LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/dynredem_p.F @ 1410

Last change on this file since 1410 was 764, checked in by Laurent Fairhead, 17 years ago

Merge entre la version V3_conv et le HEAD
YM, JG, LF

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