source: LMDZ4/trunk/libf/dyn3dpar/dynredem.F @ 1172

Last change on this file since 1172 was 1146, checked in by Laurent Fairhead, 16 years ago

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