source: LMDZ4/trunk/libf/dyn3d/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).

JG

  • 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!
4c
5      SUBROUTINE dynredem0(fichnom,iday_end,phis,nq)
6      USE IOIPSL
7      IMPLICIT NONE
8c=======================================================================
9c Ecriture du fichier de redemarrage sous format NetCDF (initialisation)
10c=======================================================================
11c   Declarations:
12c   -------------
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
26c   Arguments:
27c   ----------
28      INTEGER iday_end
29      REAL phis(ip1jmp1)
30      CHARACTER*(*) fichnom
31      INTEGER nq
32
33c   Local:
34c   ------
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
43c   Variables locales pour NetCDF:
44c
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
57c-----------------------------------------------------------------------
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
86c
87c    .....    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
93c
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)
113c
114c    .........................................................
115c
116c Creation du fichier:
117c
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
124c
125c Preciser quelques attributs globaux:
126c
127      ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 27,
128     .                       "Fichier demmarage dynamique")
129c
130c Definir les dimensions du fichiers:
131c
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)
140c
141      ierr = NF_ENDDEF(nid) ! sortir du mode de definition
142c
143c Definir et enregistrer certains champs invariants:
144c
145      ierr = NF_REDEF (nid)
146cIM 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
152cIM 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
161c
162      ierr = NF_REDEF (nid)
163cIM 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
169cIM 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
178c
179      ierr = NF_REDEF (nid)
180cIM 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
186cIM 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
195c
196      ierr = NF_REDEF (nid)
197cIM 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
203cIM 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
212c
213      ierr = NF_REDEF (nid)
214cIM 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
220cIM 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
229c
230      ierr = NF_REDEF (nid)
231cIM 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
237cIM 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
246c
247      ierr = NF_REDEF (nid)
248cIM 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
254cIM 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
263c
264      ierr = NF_REDEF (nid)
265cIM 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
271cIM 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
280c
281      ierr = NF_REDEF (nid)
282cIM 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
288cIM 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
297c
298      ierr = NF_REDEF (nid)
299cIM 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
305cIM 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
312c
313c Coefficients de passage cov. <-> contra. <--> naturel
314c
315      ierr = NF_REDEF (nid)
316      dims2(1) = idim_rlonu
317      dims2(2) = idim_rlatu
318cIM 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
324cIM 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
333c
334      ierr = NF_REDEF (nid)
335      dims2(1) = idim_rlonv
336      dims2(2) = idim_rlatv
337cIM 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
343cIM 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
352c
353c Aire de chaque maille:
354c
355      ierr = NF_REDEF (nid)
356      dims2(1) = idim_rlonv
357      dims2(2) = idim_rlatu
358cIM 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
364cIM 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
373c
374c Geopentiel au sol:
375c
376      ierr = NF_REDEF (nid)
377      dims2(1) = idim_rlonv
378      dims2(2) = idim_rlatu
379cIM 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
385cIM 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
394c
395c Definir les variables pour pouvoir les enregistrer plus tard:
396c
397      ierr = NF_REDEF (nid) ! entrer dans le mode de definition
398c
399cIM 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
405cIM 220306 END
406      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
407     .                       "Temps de simulation")
408      write(unites,200)yyears0,mmois0,jjour0
409200   format('days since ',i4,'-',i2.2,'-',i2.2,' 00:00:00')
410      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "units", 30,
411     .                         unites)
412
413c
414      dims4(1) = idim_rlonu
415      dims4(2) = idim_rlatu
416      dims4(3) = idim_s
417      dims4(4) = idim_tim
418cIM 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
424cIM 220306 END
425      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
426     .                       "Vitesse U")
427c
428      dims4(1) = idim_rlonv
429      dims4(2) = idim_rlatv
430      dims4(3) = idim_s
431      dims4(4) = idim_tim
432cIM 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
438cIM 220306 END
439      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
440     .                       "Vitesse V")
441c
442      dims4(1) = idim_rlonv
443      dims4(2) = idim_rlatu
444      dims4(3) = idim_s
445      dims4(4) = idim_tim
446cIM 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
452cIM 220306 END
453      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 11,
454     .                       "Temperature")
455c
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
462cIM 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
468cIM 220306 END
469      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))
470      ENDDO
471      ENDIF
472c
473      dims4(1) = idim_rlonv
474      dims4(2) = idim_rlatu
475      dims4(3) = idim_s
476      dims4(4) = idim_tim
477cIM 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
483cIM 220306 END
484      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,
485     .                       "C est quoi ?")
486c
487      dims3(1) = idim_rlonv
488      dims3(2) = idim_rlatu
489      dims3(3) = idim_tim
490cIM 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
496cIM 220306 END
497      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 15,
498     .                       "Pression au sol")
499c
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
512c=================================================================
513c  Ecriture du fichier de redemarrage sous format NetCDF
514c=================================================================
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
542c
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
554c  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
570c
571c  Re-ecriture du tableau de controle, itaufin n'est plus defini quand
572c  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
591c  Ecriture des champs
592c
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
707c
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
718c
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)
731c
732      RETURN
733      END
734
Note: See TracBrowser for help on using the repository browser.