source: LMDZ5/trunk/libf/dyn3dpar/dynredem.F @ 4735

Last change on this file since 4735 was 2622, checked in by Ehouarn Millour, 8 years ago

Some code tidying: turn ener.h into ener_mod.F90
EM

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