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