source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3d_common/dynredem.F @ 298

Last change on this file since 298 was 222, checked in by ymipsl, 10 years ago

Creating temporary dynamico/lmdz/saturn branche

YM
File size: 25.4 KB
Line 
1!
2! $Id: dynredem.F 1635 2012-07-12 11:37:16Z lguez $
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 control_mod, only : planet_type
12
13      IMPLICIT NONE
14c=======================================================================
15c Ecriture du fichier de redemarrage sous format NetCDF (initialisation)
16c=======================================================================
17c   Declarations:
18c   -------------
19#include "dimensions.h"
20#include "paramet.h"
21#include "comconst.h"
22#include "comvert.h"
23#include "comgeom2.h"
24#include "temps.h"
25#include "ener.h"
26#include "logic.h"
27#include "netcdf.inc"
28#include "description.h"
29#include "serre.h"
30#include "iniprint.h"
31
32c   Arguments:
33c   ----------
34      INTEGER iday_end
35      REAL phis(iip1, jjp1)
36      CHARACTER*(*) fichnom
37
38c   Local:
39c   ------
40      INTEGER iq,l
41      INTEGER length
42      PARAMETER (length = 100)
43      REAL tab_cntrl(length) ! tableau des parametres du run
44      INTEGER ierr
45      character*20 modname
46      character*80 abort_message
47
48c   Variables locales pour NetCDF:
49c
50      INTEGER dims2(2), dims3(3), dims4(4)
51      INTEGER idim_index
52      INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
53      INTEGER idim_s, idim_sig
54      INTEGER idim_tim
55      INTEGER nid,nvarid
56
57      REAL zan0,zjulian,hours
58      INTEGER yyears0,jjour0, mmois0
59      character*30 unites
60
61      character(len=12) :: start_file_type="earth" ! default start file type
62      INTEGER idecal
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      !!! AS: idecal is a hack to be able to read planeto starts...
78      !!!     .... while keeping everything OK for LMDZ EARTH
79      if ((planet_type.eq."generic").or.(planet_type.eq."mars")) then
80          write(lunout,*) trim(modname),' : Planeto-like start file'
81          start_file_type="planeto"
82          idecal = 4
83      else
84          write(lunout,*) trim(modname),' : Earth-like start file'
85          idecal = 5
86      endif
87
88      DO l=1,length
89       tab_cntrl(l) = 0.
90      ENDDO
91       tab_cntrl(1)  = REAL(iim)
92       tab_cntrl(2)  = REAL(jjm)
93       tab_cntrl(3)  = REAL(llm)
94       if (start_file_type.eq."earth") then
95         tab_cntrl(4)=REAL(day_ref)
96       else
97         !tab_cntrl(4)=REAL(day_end)
98         tab_cntrl(4)=REAL(iday_end)
99       endif
100       tab_cntrl(5)  = REAL(annee_ref)
101       tab_cntrl(idecal+1)  = rad
102       tab_cntrl(idecal+2)  = omeg
103       tab_cntrl(idecal+3)  = g
104       tab_cntrl(idecal+4)  = cpp
105       tab_cntrl(idecal+5) = kappa
106       tab_cntrl(idecal+6) = daysec
107       tab_cntrl(idecal+7) = dtvr
108       tab_cntrl(idecal+8) = etot0
109       tab_cntrl(idecal+9) = ptot0
110       tab_cntrl(idecal+10) = ztot0
111       tab_cntrl(idecal+11) = stot0
112       tab_cntrl(idecal+12) = ang0
113       tab_cntrl(idecal+13) = pa
114       tab_cntrl(idecal+14) = preff
115c
116c    .....    parametres  pour le zoom      ......   
117
118       tab_cntrl(idecal+15)  = clon
119       tab_cntrl(idecal+16)  = clat
120       tab_cntrl(idecal+17)  = grossismx
121       tab_cntrl(idecal+18)  = grossismy
122c
123      IF ( fxyhypb )   THEN
124       tab_cntrl(idecal+19) = 1.
125       tab_cntrl(idecal+20) = dzoomx
126       tab_cntrl(idecal+21) = dzoomy
127       tab_cntrl(idecal+22) = 0.
128       tab_cntrl(idecal+23) = taux
129       tab_cntrl(idecal+24) = tauy
130      ELSE
131       tab_cntrl(idecal+19) = 0.
132       tab_cntrl(idecal+20) = dzoomx
133       tab_cntrl(idecal+21) = dzoomy
134       tab_cntrl(idecal+22) = 0.
135       tab_cntrl(idecal+23) = 0.
136       tab_cntrl(idecal+24) = 0.
137       IF( ysinus )  tab_cntrl(idecal+22) = 1.
138      ENDIF
139
140      if (start_file_type.eq."earth") then
141       tab_cntrl(idecal+25) = REAL(iday_end)
142       tab_cntrl(idecal+26) = REAL(itau_dyn + itaufin)
143c start_time: start_time of simulation (not necessarily 0.)
144       tab_cntrl(idecal+27) = start_time
145      endif
146
147      if (planet_type=="mars") then ! For Mars only
148        tab_cntrl(29)=hour_ini
149      endif
150c
151c    .........................................................
152c
153c Creation du fichier:
154c
155      ierr = NF_CREATE(fichnom, NF_CLOBBER, nid)
156      IF (ierr.NE.NF_NOERR) THEN
157         write(lunout,*)"dynredem0: Pb d ouverture du fichier "
158     &                  //trim(fichnom)
159         write(lunout,*)' ierr = ', ierr
160         CALL ABORT_GCM("DYNREDEM0", "", 1)
161      ENDIF
162c
163c Preciser quelques attributs globaux:
164c
165      ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 27,
166     .                       "Fichier demmarage dynamique")
167c
168c Definir les dimensions du fichiers:
169c
170      if (start_file_type.eq."earth") then
171        ierr = NF_DEF_DIM (nid, "index", length, idim_index)
172        ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
173        ierr = NF_DEF_DIM (nid, "rlatu", jjp1, idim_rlatu)
174        ierr = NF_DEF_DIM (nid, "rlonv", iip1, idim_rlonv)
175        ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
176        ierr = NF_DEF_DIM (nid, "sigs", llm, idim_s)
177        ierr = NF_DEF_DIM (nid, "sig", llmp1, idim_sig)
178        ierr = NF_DEF_DIM (nid, "temps", NF_UNLIMITED, idim_tim)
179      else
180        ierr = NF_DEF_DIM (nid, "index", length, idim_index)
181        ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
182        ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_rlatu)
183        ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_rlonv)
184        ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
185        ierr = NF_DEF_DIM (nid, "altitude", llm, idim_s)
186        ierr = NF_DEF_DIM (nid, "interlayer", llmp1, idim_sig)
187        ierr = NF_DEF_DIM (nid, "Time", NF_UNLIMITED, idim_tim)
188      endif
189c
190      ierr = NF_ENDDEF(nid) ! sortir du mode de definition
191c
192c Definir et enregistrer certains champs invariants:
193c
194      ierr = NF_REDEF (nid)
195cIM 220306 BEG
196#ifdef NC_DOUBLE
197      ierr = NF_DEF_VAR (nid,"controle",NF_DOUBLE,1,idim_index,nvarid)
198#else
199      ierr = NF_DEF_VAR (nid,"controle",NF_FLOAT,1,idim_index,nvarid)
200#endif
201cIM 220306 END
202      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
203     .                       "Parametres de controle")
204      ierr = NF_ENDDEF(nid)
205      call NF95_PUT_VAR(nid,nvarid,tab_cntrl)
206c
207      ierr = NF_REDEF (nid)
208cIM 220306 BEG
209#ifdef NC_DOUBLE
210      ierr = NF_DEF_VAR (nid,"rlonu",NF_DOUBLE,1,idim_rlonu,nvarid)
211#else
212      ierr = NF_DEF_VAR (nid,"rlonu",NF_FLOAT,1,idim_rlonu,nvarid)
213#endif
214cIM 220306 END
215      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
216     .                       "Longitudes des points U")
217      ierr = NF_ENDDEF(nid)
218      call NF95_PUT_VAR(nid,nvarid,rlonu)
219c
220      ierr = NF_REDEF (nid)
221cIM 220306 BEG
222#ifdef NC_DOUBLE
223      ierr = NF_DEF_VAR (nid,"rlatu",NF_DOUBLE,1,idim_rlatu,nvarid)
224#else
225      ierr = NF_DEF_VAR (nid,"rlatu",NF_FLOAT,1,idim_rlatu,nvarid)
226#endif
227cIM 220306 END
228      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
229     .                       "Latitudes des points U")
230      ierr = NF_ENDDEF(nid)
231      call NF95_PUT_VAR (nid,nvarid,rlatu)
232c
233      ierr = NF_REDEF (nid)
234cIM 220306 BEG
235#ifdef NC_DOUBLE
236      ierr = NF_DEF_VAR (nid,"rlonv",NF_DOUBLE,1,idim_rlonv,nvarid)
237#else
238      ierr = NF_DEF_VAR (nid,"rlonv",NF_FLOAT,1,idim_rlonv,nvarid)
239#endif
240cIM 220306 END
241      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23,
242     .                       "Longitudes des points V")
243      ierr = NF_ENDDEF(nid)
244      call NF95_PUT_VAR(nid,nvarid,rlonv)
245c
246      ierr = NF_REDEF (nid)
247cIM 220306 BEG
248#ifdef NC_DOUBLE
249      ierr = NF_DEF_VAR (nid,"rlatv",NF_DOUBLE,1,idim_rlatv,nvarid)
250#else
251      ierr = NF_DEF_VAR (nid,"rlatv",NF_FLOAT,1,idim_rlatv,nvarid)
252#endif
253cIM 220306 END
254      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
255     .                       "Latitudes des points V")
256      ierr = NF_ENDDEF(nid)
257      call NF95_PUT_VAR(nid,nvarid,rlatv)
258c
259      if (start_file_type.eq."earth") then
260        ierr = NF_REDEF (nid)
261cIM 220306 BEG
262#ifdef NC_DOUBLE
263        ierr = NF_DEF_VAR (nid,"nivsigs",NF_DOUBLE,1,idim_s,nvarid)
264#else
265        ierr = NF_DEF_VAR (nid,"nivsigs",NF_FLOAT,1,idim_s,nvarid)
266#endif
267cIM 220306 END
268        ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 28,
269     .                       "Numero naturel des couches s")
270        ierr = NF_ENDDEF(nid)
271        call NF95_PUT_VAR(nid,nvarid,nivsigs)
272c
273        ierr = NF_REDEF (nid)
274cIM 220306 BEG
275#ifdef NC_DOUBLE
276        ierr = NF_DEF_VAR (nid,"nivsig",NF_DOUBLE,1,idim_sig,nvarid)
277#else
278        ierr = NF_DEF_VAR (nid,"nivsig",NF_FLOAT,1,idim_sig,nvarid)
279#endif
280cIM 220306 END
281        ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 32,
282     .                       "Numero naturel des couches sigma")
283        ierr = NF_ENDDEF(nid)
284        call NF95_PUT_VAR(nid,nvarid,nivsig)
285      endif ! of if (start_file_type.eq."earth")
286c
287      ierr = NF_REDEF (nid)
288cIM 220306 BEG
289#ifdef NC_DOUBLE
290      ierr = NF_DEF_VAR (nid,"ap",NF_DOUBLE,1,idim_sig,nvarid)
291#else
292      ierr = NF_DEF_VAR (nid,"ap",NF_FLOAT,1,idim_sig,nvarid)
293#endif
294cIM 220306 END
295      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
296     .                       "Coefficient A pour hybride")
297      ierr = NF_ENDDEF(nid)
298      call NF95_PUT_VAR(nid,nvarid,ap)
299c
300      ierr = NF_REDEF (nid)
301cIM 220306 BEG
302#ifdef NC_DOUBLE
303      ierr = NF_DEF_VAR (nid,"bp",NF_DOUBLE,1,idim_sig,nvarid)
304#else
305      ierr = NF_DEF_VAR (nid,"bp",NF_FLOAT,1,idim_sig,nvarid)
306#endif
307cIM 220306 END
308      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26,
309     .                       "Coefficient B pour hybride")
310      ierr = NF_ENDDEF(nid)
311      call NF95_PUT_VAR(nid,nvarid,bp)
312c
313      if (start_file_type.ne."earth") then
314        ierr = NF_REDEF (nid)
315cIM 220306 BEG
316#ifdef NC_DOUBLE
317        ierr = NF_DEF_VAR (nid,"aps",NF_DOUBLE,1,idim_s,nvarid)
318#else
319        ierr = NF_DEF_VAR (nid,"aps",NF_FLOAT,1,idim_s,nvarid)
320#endif
321cIM 220306 END
322        ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 37,
323     .                       "Coef AS: hybrid pressure at midlayers")
324        ierr = NF_ENDDEF(nid)
325        call NF95_PUT_VAR(nid,nvarid,aps)
326c
327        ierr = NF_REDEF (nid)
328cIM 220306 BEG
329#ifdef NC_DOUBLE
330        ierr = NF_DEF_VAR (nid,"bps",NF_DOUBLE,1,idim_s,nvarid)
331#else
332        ierr = NF_DEF_VAR (nid,"bps",NF_FLOAT,1,idim_s,nvarid)
333#endif
334cIM 220306 END
335        ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 34,
336     .                       "Coef BS: hybrid sigma at midlayers")
337        ierr = NF_ENDDEF(nid)
338        call NF95_PUT_VAR(nid,nvarid,bps)
339      endif ! of if (start_file_type.ne."earth")
340c
341      ierr = NF_REDEF (nid)
342cIM 220306 BEG
343#ifdef NC_DOUBLE
344      ierr = NF_DEF_VAR (nid,"presnivs",NF_DOUBLE,1,idim_s,nvarid)
345#else
346      ierr = NF_DEF_VAR (nid,"presnivs",NF_FLOAT,1,idim_s,nvarid)
347#endif
348cIM 220306 END
349      ierr = NF_ENDDEF(nid)
350      call NF95_PUT_VAR(nid,nvarid,presnivs)
351c
352      if (start_file_type.ne."earth") then
353        ierr = NF_REDEF (nid)
354#ifdef NC_DOUBLE
355        ierr = NF_DEF_VAR(nid,"latitude",NF_DOUBLE,1,idim_rlatu,nvarid)
356#else
357        ierr = NF_DEF_VAR(nid,"latitude",NF_FLOAT,1,idim_rlatu,nvarid)
358#endif
359        ierr =NF_PUT_ATT_TEXT(nid,nvarid,'units',13,"degrees_north")
360        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
361     .        "North latitude")
362        ierr = NF_ENDDEF(nid)
363        call NF95_PUT_VAR(nid,nvarid,rlatu*180/pi)
364c
365        ierr = NF_REDEF (nid)
366#ifdef NC_DOUBLE
367        ierr=NF_DEF_VAR(nid,"longitude",NF_DOUBLE,1,idim_rlonv,nvarid)
368#else
369        ierr=NF_DEF_VAR(nid,"longitude",NF_FLOAT,1,idim_rlonv,nvarid)
370#endif
371        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 14,
372     .        "East longitude")
373        ierr = NF_PUT_ATT_TEXT(nid,nvarid,'units',12,"degrees_east")
374        ierr = NF_ENDDEF(nid)
375        call NF95_PUT_VAR(nid,nvarid,rlonv*180/pi)
376c
377        ierr = NF_REDEF (nid)
378#ifdef NC_DOUBLE
379        ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
380     .       idim_s,nvarid)
381#else
382        ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
383     .       idim_s,nvarid)
384#endif
385        ierr = NF_PUT_ATT_TEXT(nid,nvarid,"long_name",10,"pseudo-alt")
386        ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
387        ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
388        ierr = NF_ENDDEF(nid)
389        call NF95_PUT_VAR(nid,nvarid,pseudoalt)
390      endif ! of if (start_file_type.ne."earth")
391c
392c Coefficients de passage cov. <-> contra. <--> naturel
393c
394      ierr = NF_REDEF (nid)
395      dims2(1) = idim_rlonu
396      dims2(2) = idim_rlatu
397cIM 220306 BEG
398#ifdef NC_DOUBLE
399      ierr = NF_DEF_VAR (nid,"cu",NF_DOUBLE,2,dims2,nvarid)
400#else
401      ierr = NF_DEF_VAR (nid,"cu",NF_FLOAT,2,dims2,nvarid)
402#endif
403cIM 220306 END
404      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
405     .                       "Coefficient de passage pour U")
406      ierr = NF_ENDDEF(nid)
407      call NF95_PUT_VAR(nid,nvarid,cu)
408c
409      ierr = NF_REDEF (nid)
410      dims2(1) = idim_rlonv
411      dims2(2) = idim_rlatv
412cIM 220306 BEG
413#ifdef NC_DOUBLE
414      ierr = NF_DEF_VAR (nid,"cv",NF_DOUBLE,2,dims2,nvarid)
415#else
416      ierr = NF_DEF_VAR (nid,"cv",NF_FLOAT,2,dims2,nvarid)
417#endif
418cIM 220306 END
419      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29,
420     .                       "Coefficient de passage pour V")
421      ierr = NF_ENDDEF(nid)
422      call NF95_PUT_VAR(nid,nvarid,cv)
423c
424c Aire de chaque maille:
425c
426      ierr = NF_REDEF (nid)
427      dims2(1) = idim_rlonv
428      dims2(2) = idim_rlatu
429cIM 220306 BEG
430#ifdef NC_DOUBLE
431      ierr = NF_DEF_VAR (nid,"aire",NF_DOUBLE,2,dims2,nvarid)
432#else
433      ierr = NF_DEF_VAR (nid,"aire",NF_FLOAT,2,dims2,nvarid)
434#endif
435cIM 220306 END
436      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22,
437     .                       "Aires de chaque maille")
438      ierr = NF_ENDDEF(nid)
439      call NF95_PUT_VAR(nid,nvarid,aire)
440c
441c Geopentiel au sol:
442c
443      ierr = NF_REDEF (nid)
444      dims2(1) = idim_rlonv
445      dims2(2) = idim_rlatu
446cIM 220306 BEG
447#ifdef NC_DOUBLE
448      ierr = NF_DEF_VAR (nid,"phisinit",NF_DOUBLE,2,dims2,nvarid)
449#else
450      ierr = NF_DEF_VAR (nid,"phisinit",NF_FLOAT,2,dims2,nvarid)
451#endif
452cIM 220306 END
453      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
454     .                       "Geopotentiel au sol")
455      ierr = NF_ENDDEF(nid)
456      call NF95_PUT_VAR(nid,nvarid,phis)
457c
458c Definir les variables pour pouvoir les enregistrer plus tard:
459c
460      ierr = NF_REDEF (nid) ! entrer dans le mode de definition
461c
462      if (start_file_type.eq."earth") then
463cIM 220306 BEG
464#ifdef NC_DOUBLE
465        ierr = NF_DEF_VAR (nid,"temps",NF_DOUBLE,1,idim_tim,nvarid)
466#else
467        ierr = NF_DEF_VAR (nid,"temps",NF_FLOAT,1,idim_tim,nvarid)
468#endif
469cIM 220306 END
470      else ! start_file_type=="planeto"
471#ifdef NC_DOUBLE
472        ierr = NF_DEF_VAR (nid,"Time",NF_DOUBLE,1,idim_tim,nvarid)
473#else
474        ierr = NF_DEF_VAR (nid,"Time",NF_FLOAT,1,idim_tim,nvarid)
475#endif
476      endif ! of if (start_file_type.eq."earth")
477      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19,
478     .                       "Temps de simulation")
479      write(unites,200)yyears0,mmois0,jjour0
480200   format('days since ',i4,'-',i2.2,'-',i2.2,' 00:00:00')
481      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "units", 30,
482     .                         unites)
483
484c
485      dims4(1) = idim_rlonu
486      dims4(2) = idim_rlatu
487      dims4(3) = idim_s
488      dims4(4) = idim_tim
489cIM 220306 BEG
490#ifdef NC_DOUBLE
491      ierr = NF_DEF_VAR (nid,"ucov",NF_DOUBLE,4,dims4,nvarid)
492#else
493      ierr = NF_DEF_VAR (nid,"ucov",NF_FLOAT,4,dims4,nvarid)
494#endif
495cIM 220306 END
496      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
497     .                       "Vitesse U")
498c
499      dims4(1) = idim_rlonv
500      dims4(2) = idim_rlatv
501      dims4(3) = idim_s
502      dims4(4) = idim_tim
503cIM 220306 BEG
504#ifdef NC_DOUBLE
505      ierr = NF_DEF_VAR (nid,"vcov",NF_DOUBLE,4,dims4,nvarid)
506#else
507      ierr = NF_DEF_VAR (nid,"vcov",NF_FLOAT,4,dims4,nvarid)
508#endif
509cIM 220306 END
510      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9,
511     .                       "Vitesse V")
512c
513      dims4(1) = idim_rlonv
514      dims4(2) = idim_rlatu
515      dims4(3) = idim_s
516      dims4(4) = idim_tim
517cIM 220306 BEG
518#ifdef NC_DOUBLE
519      ierr = NF_DEF_VAR (nid,"teta",NF_DOUBLE,4,dims4,nvarid)
520#else
521      ierr = NF_DEF_VAR (nid,"teta",NF_FLOAT,4,dims4,nvarid)
522#endif
523cIM 220306 END
524      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 11,
525     .                       "Temperature")
526c
527      dims4(1) = idim_rlonv
528      dims4(2) = idim_rlatu
529      dims4(3) = idim_s
530      dims4(4) = idim_tim
531      IF(nqtot.GE.1) THEN
532      DO iq=1,nqtot
533cIM 220306 BEG
534#ifdef NC_DOUBLE
535      ierr = NF_DEF_VAR (nid,tname(iq),NF_DOUBLE,4,dims4,nvarid)
536#else
537      ierr = NF_DEF_VAR (nid,tname(iq),NF_FLOAT,4,dims4,nvarid)
538#endif
539cIM 220306 END
540      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))
541      ENDDO
542      ENDIF
543c
544      dims4(1) = idim_rlonv
545      dims4(2) = idim_rlatu
546      dims4(3) = idim_s
547      dims4(4) = idim_tim
548cIM 220306 BEG
549#ifdef NC_DOUBLE
550      ierr = NF_DEF_VAR (nid,"masse",NF_DOUBLE,4,dims4,nvarid)
551#else
552      ierr = NF_DEF_VAR (nid,"masse",NF_FLOAT,4,dims4,nvarid)
553#endif
554cIM 220306 END
555      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,
556     .                       "C est quoi ?")
557c
558      dims3(1) = idim_rlonv
559      dims3(2) = idim_rlatu
560      dims3(3) = idim_tim
561cIM 220306 BEG
562#ifdef NC_DOUBLE
563      ierr = NF_DEF_VAR (nid,"ps",NF_DOUBLE,3,dims3,nvarid)
564#else
565      ierr = NF_DEF_VAR (nid,"ps",NF_FLOAT,3,dims3,nvarid)
566#endif
567cIM 220306 END
568      ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 15,
569     .                       "Pression au sol")
570c
571      ierr = NF_ENDDEF(nid) ! sortir du mode de definition
572      ierr = NF_CLOSE(nid) ! fermer le fichier
573
574      write(lunout,*)'dynredem0: iim,jjm,llm,iday_end',
575     &               iim,jjm,llm,iday_end
576      write(lunout,*)'dynredem0: rad,omeg,g,cpp,kappa',
577     &        rad,omeg,g,cpp,kappa
578
579      RETURN
580      END
581
582!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
583
584      SUBROUTINE dynredem1(fichnom,time,
585     .                     vcov,ucov,teta,q,masse,ps)
586      USE infotrac
587      USE control_mod, only : planet_type
588      use netcdf, only: NF90_get_VAR
589      use netcdf95, only: NF95_PUT_VAR
590 
591      IMPLICIT NONE
592c=================================================================
593c  Ecriture du fichier de redemarrage sous format NetCDF
594c=================================================================
595#include "dimensions.h"
596#include "paramet.h"
597#include "description.h"
598#include "netcdf.inc"
599#include "comvert.h"
600#include "comgeom.h"
601#include "temps.h"
602#include "iniprint.h"
603
604
605      INTEGER l
606      REAL vcov(iip1,jjm,llm),ucov(iip1, jjp1,llm) 
607      REAL teta(iip1, jjp1,llm)                   
608      REAL ps(iip1, jjp1),masse(iip1, jjp1,llm)                   
609      REAL q(iip1, jjp1, llm, nqtot)
610      CHARACTER*(*) fichnom
611     
612      REAL time
613      INTEGER nid, nvarid, nid_trac, nvarid_trac
614      REAL trac_tmp(ip1jmp1,llm)     
615      INTEGER ierr, ierr_file 
616      INTEGER iq
617      INTEGER length
618      PARAMETER (length = 100)
619      REAL tab_cntrl(length) ! tableau des parametres du run
620      character(len=*),parameter :: modname='dynredem1'
621      character*80 abort_message
622c
623      INTEGER nb
624      SAVE nb
625      DATA nb / 0 /
626      character(len=12) :: start_file_type="earth" ! default start file type
627
628      if ((planet_type.eq."generic").or.(planet_type.eq."mars")) then
629          write(lunout,*) trim(modname),' : Planeto-like start file'
630          start_file_type="planeto"
631      else
632          write(lunout,*) trim(modname),' : Earth-like start file'
633      endif
634
635      ierr = NF_OPEN(fichnom, NF_WRITE, nid)
636      IF (ierr .NE. NF_NOERR) THEN
637         PRINT*, "dynredem1: Pb. d ouverture "//trim(fichnom)
638         CALL abort_gcm("dynredem1", "", 1)
639      ENDIF
640
641c  Ecriture/extension de la coordonnee temps
642
643      nb = nb + 1
644      if (start_file_type.eq."earth") then
645        ierr = NF_INQ_VARID(nid, "temps", nvarid)
646        IF (ierr .NE. NF_NOERR) THEN
647          write(lunout,*) NF_STRERROR(ierr)
648          abort_message='Variable temps n est pas definie'
649          CALL abort_gcm(modname,abort_message,ierr)
650        ENDIF
651      else
652        ierr = NF_INQ_VARID(nid,"Time", nvarid)
653        IF (ierr .NE. NF_NOERR) THEN
654          write(lunout,*) NF_STRERROR(ierr)
655          abort_message='Variable Time not found'
656          CALL abort_gcm(modname,abort_message,ierr)
657        ENDIF
658      endif ! of if (start_file_type.eq."earth")
659      call NF95_PUT_VAR(nid,nvarid,time,start=(/nb/))
660      write(lunout,*) "dynredem1: Enregistrement pour ", nb, time
661
662c
663c  Re-ecriture du tableau de controle, itaufin n'est plus defini quand
664c  on passe dans dynredem0
665      ierr = NF_INQ_VARID (nid, "controle", nvarid)
666      IF (ierr .NE. NF_NOERR) THEN
667         abort_message="dynredem1: Le champ <controle> est absent"
668         ierr = 1
669         CALL abort_gcm(modname,abort_message,ierr)
670      ENDIF
671      ierr = NF90_GET_VAR(nid, nvarid, tab_cntrl)
672      if (start_file_type=="earth") then
673        tab_cntrl(31) = REAL(itau_dyn + itaufin)
674      else
675        tab_cntrl(31) = 0
676      endif
677      call NF95_PUT_VAR(nid,nvarid,tab_cntrl)
678
679c  Ecriture des champs
680c
681      ierr = NF_INQ_VARID(nid, "ucov", nvarid)
682      IF (ierr .NE. NF_NOERR) THEN
683         abort_message="Variable ucov n est pas definie"
684         ierr=1
685         CALL abort_gcm(modname,abort_message,ierr)
686      ENDIF
687      call NF95_PUT_VAR(nid,nvarid,ucov,start=(/1,1,1,nb/))
688
689      ierr = NF_INQ_VARID(nid, "vcov", nvarid)
690      IF (ierr .NE. NF_NOERR) THEN
691         abort_message="Variable vcov n est pas definie"
692         ierr=1
693         CALL abort_gcm(modname,abort_message,ierr)
694      ENDIF
695      call NF95_PUT_VAR(nid,nvarid,vcov,start=(/1,1,1,nb/))
696
697      ierr = NF_INQ_VARID(nid, "teta", nvarid)
698      IF (ierr .NE. NF_NOERR) THEN
699         abort_message="Variable teta n est pas definie"
700         ierr=1
701         CALL abort_gcm(modname,abort_message,ierr)
702      ENDIF
703      call NF95_PUT_VAR(nid,nvarid,teta,start=(/1,1,1,nb/))
704
705      IF (type_trac == 'inca') THEN
706! Ajout Anne pour lecture valeurs traceurs dans un fichier start_trac.nc
707         ierr_file = NF_OPEN ("start_trac.nc", NF_NOWRITE,nid_trac)
708         IF (ierr_file .NE.NF_NOERR) THEN
709            write(lunout,*)'dynredem1: Pb d''ouverture du fichier',
710     &                     ' start_trac.nc'
711            write(lunout,*)' ierr = ', ierr_file
712         ENDIF
713      END IF
714
715      IF(nqtot.GE.1) THEN
716      do iq=1,nqtot
717
718         IF (type_trac /= 'inca') THEN
719            ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
720            IF (ierr .NE. NF_NOERR) THEN
721               abort_message="Variable  tname(iq) n est pas definie"
722               ierr=1
723               CALL abort_gcm(modname,abort_message,ierr)
724            ENDIF
725            call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq),start=(/1,1,1,nb/))
726        ELSE ! type_trac = inca
727! lecture de la valeur du traceur dans start_trac.nc
728           IF (ierr_file .ne. 2) THEN
729             ierr = NF_INQ_VARID (nid_trac, tname(iq), nvarid_trac)
730             IF (ierr .NE. NF_NOERR) THEN
731                write(lunout,*) "dynredem1: ",trim(tname(iq)),
732     &                          " est absent de start_trac.nc"
733                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
734                IF (ierr .NE. NF_NOERR) THEN
735                   abort_message="dynredem1: Variable "//
736     &                     trim(tname(iq))//" n est pas definie"
737                   ierr=1
738                   CALL abort_gcm(modname,abort_message,ierr)
739                ENDIF
740                call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq))
741               
742             ELSE
743                write(lunout,*) "dynredem1: ",trim(tname(iq)),
744     &              " est present dans start_trac.nc"
745               ierr = NF90_GET_VAR(nid_trac, nvarid_trac, trac_tmp)
746                IF (ierr .NE. NF_NOERR) THEN
747                   abort_message="dynredem1: Lecture echouee pour"//
748     &                    trim(tname(iq))
749                   ierr=1
750                   CALL abort_gcm(modname,abort_message,ierr)
751                ENDIF
752                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
753                IF (ierr .NE. NF_NOERR) THEN
754                   abort_message="dynredem1: Variable "//
755     &                trim(tname(iq))//" n est pas definie"
756                   ierr=1
757                   CALL abort_gcm(modname,abort_message,ierr)
758                ENDIF
759                call NF95_PUT_VAR(nid, nvarid, trac_tmp)
760               
761             ENDIF ! IF (ierr .NE. NF_NOERR)
762! fin lecture du traceur
763          ELSE                  ! si il n'y a pas de fichier start_trac.nc
764!             print *, 'il n y a pas de fichier start_trac'
765             ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
766             IF (ierr .NE. NF_NOERR) THEN
767                abort_message="dynredem1: Variable "//
768     &                trim(tname(iq))//" n est pas definie"
769                   ierr=1
770                   CALL abort_gcm(modname,abort_message,ierr)
771             ENDIF
772             call NF95_PUT_VAR(nid,nvarid,q(:,:,:,iq),
773     &                  start=(/1,1,1,nb/))
774          ENDIF ! (ierr_file .ne. 2)
775       END IF   !type_trac
776     
777      ENDDO
778      ENDIF
779c
780      ierr = NF_INQ_VARID(nid, "masse", nvarid)
781      IF (ierr .NE. NF_NOERR) THEN
782         abort_message="dynredem1: Variable masse n est pas definie"
783         ierr=1
784         CALL abort_gcm(modname,abort_message,ierr)
785      ENDIF
786      call NF95_PUT_VAR(nid,nvarid,masse,start=(/1,1,1,nb/))
787c
788      ierr = NF_INQ_VARID(nid, "ps", nvarid)
789      IF (ierr .NE. NF_NOERR) THEN
790         abort_message="dynredem1: Variable ps n est pas definie"
791         ierr=1
792         CALL abort_gcm(modname,abort_message,ierr)
793      ENDIF
794      call NF95_PUT_VAR(nid,nvarid,ps,start=(/1,1,nb/))
795
796      ierr = NF_CLOSE(nid)
797c
798      RETURN
799      END
800
Note: See TracBrowser for help on using the repository browser.