source: LMDZ4/branches/V3_test/libf/phylmd/phyredem.F @ 1195

Last change on this file since 1195 was 735, checked in by Laurent Fairhead, 18 years ago

Modifications version parallele
YM/LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 31.8 KB
Line 
1!
2! $Header$
3!
4c
5      SUBROUTINE phyredem (fichnom,dtime,radpas,
6     .           rlat_p,rlon_p, pctsrf_p,tsol_p,tsoil_p,
7cIM "slab" ocean
8     .           tslab,seaice,
9     .           qsurf_p,qsol_p,snow_p,albedo_p, alblw_p, evap_p,
10     .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,fder_p,
11     .           radsol_p,frugs_p,agesno_p,zmea_p,zstd_p,zsig_p,
12     .           zgam_p,zthe_p,zpic_p,zval_p,rugsrel_p,
13     .           t_ancien_p, q_ancien_p, rnebcon_p, ratqs_p, clwcon_p,
14     .           run_off_lic_0_p)
15      USE dimphy
16      IMPLICIT none
17c======================================================================
18c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
19c Objet: Ecriture de l'etat de redemarrage pour la physique
20c======================================================================
21cym#include "dimensions.h"
22cym#include "dimphy.h"
23#include "netcdf.inc"
24#include "indicesol.h"
25#include "dimsoil.h"
26#include "clesphys.h"
27#include "control.h"
28#include "temps.h"
29c======================================================================
30      CHARACTER*(*) fichnom
31      REAL dtime
32      INTEGER radpas
33      REAL rlat_p(klon), rlon_p(klon)
34      REAL tsol_p(klon,nbsrf)
35      REAL tsoil_p(klon,nsoilmx,nbsrf)
36cIM "slab" ocean
37      REAL tslab_p(klon), seaice_p(klon)
38      REAL qsurf_p(klon,nbsrf)
39      REAL qsol_p(klon)
40      REAL snow_p(klon,nbsrf)
41      REAL albedo_p(klon,nbsrf)
42cIM BEG
43      REAL alblw_p(klon,nbsrf)
44cIM END
45      REAL evap_p(klon,nbsrf)
46      REAL rain_fall_p(klon)
47      REAL snow_fall_p(klon)
48      real solsw_p(klon)
49      real sollw_p(klon)
50      real fder_p(klon)
51      REAL radsol_p(klon)
52      REAL frugs_p(klon,nbsrf)
53      REAL agesno_p(klon,nbsrf)
54      REAL zmea_p(klon)
55      REAL zstd_p(klon)
56      REAL zsig_p(klon)
57      REAL zgam_p(klon)
58      REAL zthe_p(klon)
59      REAL zpic_p(klon)
60      REAL zval_p(klon)
61      REAL rugsrel_p(klon)
62      REAL pctsrf_p(klon, nbsrf)
63      REAL t_ancien_p(klon,klev), q_ancien_p(klon,klev)
64      real clwcon_p(klon,klev),rnebcon_p(klon,klev),ratqs_p(klon,klev)
65      REAL run_off_lic_0_p(klon)
66     
67      REAL rlat(klon2), rlon(klon2)
68      REAL tsol(klon2,nbsrf)
69      REAL tsoil(klon2,nsoilmx,nbsrf)
70      REAL tslab(klon2), seaice(klon2)
71      REAL qsurf(klon2,nbsrf)
72      REAL qsol(klon2)
73      REAL snow(klon2,nbsrf)
74      REAL albedo(klon2,nbsrf)
75      REAL alblw(klon2,nbsrf)
76      REAL evap(klon2,nbsrf)
77      REAL rain_fall(klon2)
78      REAL snow_fall(klon2)
79      real solsw(klon2)
80      real sollw(klon2)
81      real fder(klon2)
82      REAL radsol(klon2)
83      REAL frugs(klon2,nbsrf)
84      REAL agesno(klon2,nbsrf)
85      REAL zmea(klon2)
86      REAL zstd(klon2)
87      REAL zsig(klon2)
88      REAL zgam(klon2)
89      REAL zthe(klon2)
90      REAL zpic(klon2)
91      REAL zval(klon2)
92      REAL rugsrel(klon2)
93      REAL pctsrf(klon2, nbsrf)
94      REAL t_ancien(klon2,klev), q_ancien(klon2,klev)
95      real clwcon(klon2,klev),rnebcon(klon2,klev),ratqs(klon2,klev)
96      REAL run_off_lic_0(klon2)
97      REAL masq(klon2)
98c
99      INTEGER nid, nvarid, idim1, idim2, idim3
100      INTEGER ierr
101      INTEGER length
102      PARAMETER (length=100)
103      REAL tab_cntrl(length)
104c
105      INTEGER isoil, nsrf
106      CHARACTER*7 str7
107      CHARACTER*2 str2
108c
109      REAL,SAVE,ALLOCATABLE :: rlat_mpi(:), rlon_mpi(:)
110      REAL,SAVE,ALLOCATABLE :: tsol_mpi(:,:)
111      REAL,SAVE,ALLOCATABLE :: tsoil_mpi(:,:,:)
112      REAL,SAVE,ALLOCATABLE :: tslab_mpi(:)
113      REAL,SAVE,ALLOCATABLE :: seaice_mpi(:)
114      REAL,SAVE,ALLOCATABLE :: qsurf_mpi(:,:)
115      REAL,SAVE,ALLOCATABLE :: qsol_mpi(:)
116      REAL,SAVE,ALLOCATABLE :: snow_mpi(:,:)
117      REAL,SAVE,ALLOCATABLE :: albedo_mpi(:,:)
118      REAL,SAVE,ALLOCATABLE :: alblw_mpi(:,:)
119      REAL,SAVE,ALLOCATABLE :: evap_mpi(:,:)
120      REAL,SAVE,ALLOCATABLE :: rain_fall_mpi(:)
121      REAL,SAVE,ALLOCATABLE :: snow_fall_mpi(:)
122      REAL,SAVE,ALLOCATABLE :: solsw_mpi(:)
123      REAL,SAVE,ALLOCATABLE :: sollw_mpi(:)
124      REAL,SAVE,ALLOCATABLE :: fder_mpi(:)
125      REAL,SAVE,ALLOCATABLE :: radsol_mpi(:)
126      REAL,SAVE,ALLOCATABLE :: frugs_mpi(:,:)
127      REAL,SAVE,ALLOCATABLE :: agesno_mpi(:,:)
128      REAL,SAVE,ALLOCATABLE :: zmea_mpi(:)
129      REAL,SAVE,ALLOCATABLE :: zstd_mpi(:)
130      REAL,SAVE,ALLOCATABLE :: zsig_mpi(:)
131      REAL,SAVE,ALLOCATABLE :: zgam_mpi(:)
132      REAL,SAVE,ALLOCATABLE :: zthe_mpi(:)
133      REAL,SAVE,ALLOCATABLE :: zpic_mpi(:)
134      REAL,SAVE,ALLOCATABLE :: zval_mpi(:)
135      REAL,SAVE,ALLOCATABLE :: rugsrel_mpi(:)
136      REAL,SAVE,ALLOCATABLE :: pctsrf_mpi(:, :)
137      REAL,SAVE,ALLOCATABLE :: t_ancien_mpi(:,:), q_ancien_mpi(:,:)
138      REAL,SAVE,ALLOCATABLE :: clwcon_mpi(:,:),rnebcon_mpi(:,:)
139      REAL,SAVE,ALLOCATABLE :: ratqs_mpi(:,:)
140      REAL,SAVE,ALLOCATABLE :: run_off_lic_0_mpi(:)
141      REAL,SAVE,ALLOCATABLE :: masq_mpi(:)
142     
143c$OMP MASTER
144      ALLOCATE( rlat_mpi(klon_mpi), rlon_mpi(klon_mpi))
145      ALLOCATE( tsol_mpi(klon_mpi,nbsrf))
146      ALLOCATE( tsoil_mpi(klon_mpi,nsoilmx,nbsrf))
147      ALLOCATE( tslab_mpi(klon_mpi))
148      ALLOCATE( seaice_mpi(klon_mpi))
149      ALLOCATE( qsurf_mpi(klon_mpi,nbsrf))
150      ALLOCATE( qsol_mpi(klon_mpi))
151      ALLOCATE( snow_mpi(klon_mpi,nbsrf))
152      ALLOCATE( albedo_mpi(klon_mpi,nbsrf))
153      ALLOCATE( alblw_mpi(klon_mpi,nbsrf))
154      ALLOCATE( evap_mpi(klon_mpi,nbsrf))
155      ALLOCATE( rain_fall_mpi(klon_mpi))
156      ALLOCATE( snow_fall_mpi(klon_mpi))
157      ALLOCATE( solsw_mpi(klon_mpi))
158      ALLOCATE( sollw_mpi(klon_mpi))
159      ALLOCATE( fder_mpi(klon_mpi))
160      ALLOCATE( radsol_mpi(klon_mpi))
161      ALLOCATE( frugs_mpi(klon_mpi,nbsrf))
162      ALLOCATE( agesno_mpi(klon_mpi,nbsrf))
163      ALLOCATE( zmea_mpi(klon_mpi))
164      ALLOCATE( zstd_mpi(klon_mpi))
165      ALLOCATE( zsig_mpi(klon_mpi))
166      ALLOCATE( zgam_mpi(klon_mpi))
167      ALLOCATE( zthe_mpi(klon_mpi))
168      ALLOCATE( zpic_mpi(klon_mpi))
169      ALLOCATE( zval_mpi(klon_mpi))
170      ALLOCATE( rugsrel_mpi(klon_mpi))
171      ALLOCATE( pctsrf_mpi(klon_mpi, nbsrf))
172      ALLOCATE( t_ancien_mpi(klon_mpi,klev))
173      ALLOCATE( q_ancien_mpi(klon_mpi,klev))
174      ALLOCATE( clwcon_mpi(klon_mpi,klev),rnebcon_mpi(klon_mpi,klev))
175      ALLOCATE( ratqs_mpi(klon_mpi,klev))
176      ALLOCATE( run_off_lic_0_mpi(klon_mpi))
177      ALLOCATE( masq_mpi(klon_mpi))
178c$OMP END MASTER
179c$OMP BARRIER     
180
181      call GatherField_omp( rlat_p,rlat_mpi,1)
182      call GatherField_omp( rlon_p,rlon_mpi,1)
183      call GatherField_omp( tsol_p,tsol_mpi,nbsrf)
184      call GatherField_omp( tsoil_p,tsoil_mpi,nsoilmx*nbsrf)
185      call GatherField_omp( tslab_p,tslab_mpi,1)
186      call GatherField_omp( seaice_p,seaice_mpi,1)
187      call GatherField_omp( qsurf_p,qsurf_mpi,nbsrf)
188      call GatherField_omp( qsol_p,qsol_mpi,1)
189      call GatherField_omp( snow_p,snow_mpi,nbsrf)
190      call GatherField_omp( albedo_p,albedo_mpi,nbsrf)
191      call GatherField_omp( alblw_p,alblw_mpi,nbsrf)
192      call GatherField_omp( evap_p,evap_mpi,nbsrf)
193      call GatherField_omp( radsol_p,radsol_mpi,1)
194      call GatherField_omp( rain_fall_p,rain_fall_mpi,1)
195      call GatherField_omp( snow_fall_p,snow_fall_mpi,1)
196      call GatherField_omp( sollw_p,sollw_mpi,1)
197      call GatherField_omp( solsw_p,solsw_mpi,1)
198      call GatherField_omp( fder_p,fder_mpi,1)
199      call GatherField_omp( frugs_p,frugs_mpi,nbsrf)
200      call GatherField_omp( agesno_p,agesno_mpi,nbsrf)
201      call GatherField_omp( zmea_p,zmea_mpi,1)
202      call GatherField_omp( zstd_p,zstd_mpi,1)
203      call GatherField_omp( zsig_p,zsig_mpi,1)
204      call GatherField_omp( zgam_p,zgam_mpi,1)
205      call GatherField_omp( zthe_p,zthe_mpi,1)
206      call GatherField_omp( zpic_p,zpic_mpi,1)
207      call GatherField_omp( zval_p,zval_mpi,1)
208      call GatherField_omp( rugsrel_p,rugsrel_mpi,1)
209      call GatherField_omp( pctsrf_p,pctsrf_mpi, nbsrf)
210      call GatherField_omp( run_off_lic_0_p,run_off_lic_0_mpi,1)
211      call GatherField_omp( t_ancien_p,t_ancien_mpi,klev)
212      call GatherField_omp( q_ancien_p,q_ancien_mpi,klev)
213      call GatherField_omp( rnebcon_p,rnebcon_mpi,klev)
214      call GatherField_omp( clwcon_p,clwcon_mpi,klev)
215      call GatherField_omp( ratqs_p,ratqs_mpi,klev)
216      call GatherField_omp( zmasq,masq_mpi,1)
217
218c$OMP MASTER     
219      call GatherField( rlat_mpi,rlat,1)
220      call GatherField( rlon_mpi,rlon,1)
221      call GatherField( tsol_mpi,tsol,nbsrf)
222      call GatherField( tsoil_mpi,tsoil,nsoilmx*nbsrf)
223      call GatherField( tslab_mpi,tslab,1)
224      call GatherField( seaice_mpi,seaice,1)
225      call GatherField( qsurf_mpi,qsurf,nbsrf)
226      call GatherField( qsol_mpi,qsol,1)
227      call GatherField( snow_mpi,snow,nbsrf)
228      call GatherField( albedo_mpi,albedo,nbsrf)
229      call GatherField( alblw_mpi,alblw,nbsrf)
230      call GatherField( evap_mpi,evap,nbsrf)
231      call GatherField( radsol_mpi,radsol,1)
232      call GatherField( rain_fall_mpi,rain_fall,1)
233      call GatherField( snow_fall_mpi,snow_fall,1)
234      call GatherField( sollw_mpi,sollw,1)
235      call GatherField( solsw_mpi,solsw,1)
236      call GatherField( fder_mpi,fder,1)
237      call GatherField( frugs_mpi,frugs,nbsrf)
238      call GatherField( agesno_mpi,agesno,nbsrf)
239      call GatherField( zmea_mpi,zmea,1)
240      call GatherField( zstd_mpi,zstd,1)
241      call GatherField( zsig_mpi,zsig,1)
242      call GatherField( zgam_mpi,zgam,1)
243      call GatherField( zthe_mpi,zthe,1)
244      call GatherField( zpic_mpi,zpic,1)
245      call GatherField( zval_mpi,zval,1)
246      call GatherField( rugsrel_mpi,rugsrel,1)
247      call GatherField( pctsrf_mpi,pctsrf, nbsrf)
248      call GatherField( run_off_lic_0_mpi,run_off_lic_0,1)
249      call GatherField( t_ancien_mpi,t_ancien,klev)
250      call GatherField( q_ancien_mpi,q_ancien,klev)
251      call GatherField( rnebcon_mpi,rnebcon,klev)
252      call GatherField( clwcon_mpi,clwcon,klev)
253      call GatherField( ratqs_mpi,ratqs,klev)
254      call GatherField( masq_mpi,masq,1)
255     
256      DEALLOCATE( rlat_mpi, rlon_mpi)
257      DEALLOCATE( tsol_mpi)
258      DEALLOCATE( tsoil_mpi)
259      DEALLOCATE( tslab_mpi)
260      DEALLOCATE( seaice_mpi)
261      DEALLOCATE( qsurf_mpi)
262      DEALLOCATE( qsol_mpi)
263      DEALLOCATE( snow_mpi)
264      DEALLOCATE( albedo_mpi)
265      DEALLOCATE( alblw_mpi)
266      DEALLOCATE( evap_mpi)
267      DEALLOCATE( rain_fall_mpi)
268      DEALLOCATE( snow_fall_mpi)
269      DEALLOCATE( solsw_mpi)
270      DEALLOCATE( sollw_mpi)
271      DEALLOCATE( fder_mpi)
272      DEALLOCATE( radsol_mpi)
273      DEALLOCATE( frugs_mpi)
274      DEALLOCATE( agesno_mpi)
275      DEALLOCATE( zmea_mpi)
276      DEALLOCATE( zstd_mpi)
277      DEALLOCATE( zsig_mpi)
278      DEALLOCATE( zgam_mpi)
279      DEALLOCATE( zthe_mpi)
280      DEALLOCATE( zpic_mpi)
281      DEALLOCATE( zval_mpi)
282      DEALLOCATE( rugsrel_mpi)
283      DEALLOCATE( pctsrf_mpi)
284      DEALLOCATE( t_ancien_mpi)
285      DEALLOCATE( q_ancien_mpi)
286      DEALLOCATE( clwcon_mpi,rnebcon_mpi)
287      DEALLOCATE( ratqs_mpi)
288      DEALLOCATE( run_off_lic_0_mpi)
289      DEALLOCATE( masq_mpi)
290
291      if (phy_rank==0) then
292     
293      ierr = NF_CREATE(fichnom, NF_CLOBBER, nid)
294      IF (ierr.NE.NF_NOERR) THEN
295        write(6,*)' Pb d''ouverture du fichier '//fichnom
296        write(6,*)' ierr = ', ierr
297        CALL ABORT
298      ENDIF
299c
300      ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 28,
301     .                       "Fichier redemmarage physique")
302c
303      ierr = NF_DEF_DIM (nid, "index", length, idim1)
304      ierr = NF_DEF_DIM (nid, "points_physiques", klon2, idim2)
305      ierr = NF_DEF_DIM (nid, "horizon_vertical", klon2*klev, idim3)
306c
307      ierr = NF_ENDDEF(nid)
308c
309      DO ierr = 1, length
310         tab_cntrl(ierr) = 0.0
311      ENDDO
312      tab_cntrl(1) = dtime
313      tab_cntrl(2) = radpas
314      tab_cntrl(3) = co2_ppm
315      tab_cntrl(4) = solaire
316      tab_cntrl(5) = iflag_con
317      tab_cntrl(6) = nbapp_rad
318
319      IF( cycle_diurne ) tab_cntrl( 7 ) = 1.
320      IF(   soil_model ) tab_cntrl( 8 ) = 1.
321      IF(     new_oliq ) tab_cntrl( 9 ) = 1.
322      IF(     ok_orodr ) tab_cntrl(10 ) = 1.
323      IF(     ok_orolf ) tab_cntrl(11 ) = 1.
324
325      tab_cntrl(13) = day_end
326      tab_cntrl(14) = annee_ref
327      tab_cntrl(15) = itau_phy
328c
329      ierr = NF_REDEF (nid)
330#ifdef NC_DOUBLE
331      ierr = NF_DEF_VAR (nid, "controle", NF_DOUBLE, 1, idim1,nvarid)
332#else
333      ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1, idim1,nvarid)
334#endif
335      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22,
336     .                        "Parametres de controle")
337      ierr = NF_ENDDEF(nid)
338#ifdef NC_DOUBLE
339      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)
340#else
341      ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)
342#endif
343c
344      ierr = NF_REDEF (nid)
345#ifdef NC_DOUBLE
346      ierr = NF_DEF_VAR (nid, "longitude", NF_DOUBLE, 1, idim2,nvarid)
347#else
348      ierr = NF_DEF_VAR (nid, "longitude", NF_FLOAT, 1, idim2,nvarid)
349#endif
350      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32,
351     .                        "Longitudes de la grille physique")
352      ierr = NF_ENDDEF(nid)
353#ifdef NC_DOUBLE
354      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlon)
355#else
356      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlon)
357#endif
358c
359      ierr = NF_REDEF (nid)
360#ifdef NC_DOUBLE
361      ierr = NF_DEF_VAR (nid, "latitude", NF_DOUBLE, 1, idim2,nvarid)
362#else
363      ierr = NF_DEF_VAR (nid, "latitude", NF_FLOAT, 1, idim2,nvarid)
364#endif
365      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 31,
366     .                        "Latitudes de la grille physique")
367      ierr = NF_ENDDEF(nid)
368#ifdef NC_DOUBLE
369      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlat)
370#else
371      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlat)
372#endif
373c
374C PB ajout du masque terre/mer
375C
376      ierr = NF_REDEF (nid)
377#ifdef NC_DOUBLE
378      ierr = NF_DEF_VAR (nid, "masque", NF_DOUBLE, 1, idim2,nvarid)
379#else
380      ierr = NF_DEF_VAR (nid, "masque", NF_FLOAT, 1, idim2,nvarid)
381#endif
382      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 16,
383     .                        "masque terre mer")
384      ierr = NF_ENDDEF(nid)
385#ifdef NC_DOUBLE
386      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,masq)
387#else
388      ierr = NF_PUT_VAR_REAL (nid,nvarid,masq)
389#endif     
390c BP ajout des fraction de chaque sous-surface
391C
392C 1. fraction de terre
393C
394      ierr = NF_REDEF (nid)
395#ifdef NC_DOUBLE
396      ierr = NF_DEF_VAR (nid, "FTER", NF_DOUBLE, 1, idim2,nvarid)
397#else
398      ierr = NF_DEF_VAR (nid, "FTER", NF_FLOAT, 1, idim2,nvarid)
399#endif
400      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
401     .                        "fraction de continent")
402      ierr = NF_ENDDEF(nid)
403#ifdef NC_DOUBLE
404      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon2, is_ter))
405#else
406      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon2, is_ter))
407#endif
408C
409C 2. Fraction de glace de terre
410C
411      ierr = NF_REDEF (nid)
412#ifdef NC_DOUBLE
413      ierr = NF_DEF_VAR (nid, "FLIC", NF_DOUBLE, 1, idim2,nvarid)
414#else
415      ierr = NF_DEF_VAR (nid, "FLIC", NF_FLOAT, 1, idim2,nvarid)
416#endif
417      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 24,
418     .                        "fraction glace de terre")
419      ierr = NF_ENDDEF(nid)
420#ifdef NC_DOUBLE
421      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon2,is_lic))
422#else
423      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon2, is_lic))
424#endif
425C
426C 3. fraction ocean
427C
428      ierr = NF_REDEF (nid)
429#ifdef NC_DOUBLE
430      ierr = NF_DEF_VAR (nid, "FOCE", NF_DOUBLE, 1, idim2,nvarid)
431#else
432      ierr = NF_DEF_VAR (nid, "FOCE", NF_FLOAT, 1, idim2,nvarid)
433#endif
434      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14,
435     .                        "fraction ocean")
436      ierr = NF_ENDDEF(nid)
437#ifdef NC_DOUBLE
438      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon2, is_oce))
439#else
440      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon2, is_oce))
441#endif
442C
443C 4. Fraction glace de mer
444C
445      ierr = NF_REDEF (nid)
446#ifdef NC_DOUBLE
447      ierr = NF_DEF_VAR (nid, "FSIC", NF_DOUBLE, 1, idim2,nvarid)
448#else
449      ierr = NF_DEF_VAR (nid, "FSIC", NF_FLOAT, 1, idim2,nvarid)
450#endif
451      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 18,
452     .                        "fraction glace mer")
453      ierr = NF_ENDDEF(nid)
454#ifdef NC_DOUBLE
455      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pctsrf(1 : klon2, is_sic))
456#else
457      ierr = NF_PUT_VAR_REAL (nid,nvarid,pctsrf(1 : klon2, is_sic))
458#endif
459C
460C
461c
462      DO nsrf = 1, nbsrf
463        IF (nsrf.LE.99) THEN
464        WRITE(str2,'(i2.2)') nsrf
465        ierr = NF_REDEF (nid)
466#ifdef NC_DOUBLE
467        ierr = NF_DEF_VAR (nid, "TS"//str2, NF_DOUBLE, 1, idim2,nvarid)
468#else
469        ierr = NF_DEF_VAR (nid, "TS"//str2, NF_FLOAT, 1, idim2,nvarid)
470#endif
471        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
472     .                        "Temperature de surface No."//str2)
473        ierr = NF_ENDDEF(nid)
474        ELSE
475        PRINT*, "Trop de sous-mailles"
476        CALL abort
477        ENDIF
478#ifdef NC_DOUBLE
479        ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tsol(1,nsrf))
480#else
481        ierr = NF_PUT_VAR_REAL (nid,nvarid,tsol(1,nsrf))
482#endif
483      ENDDO
484c
485      DO nsrf = 1, nbsrf
486      DO isoil=1, nsoilmx
487        IF (isoil.LE.99 .AND. nsrf.LE.99) THEN
488        WRITE(str7,'(i2.2,"srf",i2.2)') isoil,nsrf
489        ierr = NF_REDEF (nid)
490#ifdef NC_DOUBLE
491        ierr = NF_DEF_VAR (nid, "Tsoil"//str7,NF_DOUBLE,1,idim2,nvarid)
492#else
493        ierr = NF_DEF_VAR (nid, "Tsoil"//str7,NF_FLOAT,1,idim2,nvarid)
494#endif
495        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 29,
496     .                        "Temperature du sol No."//str7)
497        ierr = NF_ENDDEF(nid)
498        ELSE
499        PRINT*, "Trop de couches"
500        CALL abort
501        ENDIF
502#ifdef NC_DOUBLE
503        ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tsoil(1,isoil,nsrf))
504#else
505        ierr = NF_PUT_VAR_REAL (nid,nvarid,tsoil(1,isoil,nsrf))
506#endif
507      ENDDO
508      ENDDO
509c
510cIM "slab" ocean
511      ierr = NF_REDEF (nid)
512#ifdef NC_DOUBLE
513      ierr = NF_DEF_VAR (nid, "TSLAB", NF_DOUBLE, 1, idim2,nvarid)
514#else
515      ierr = NF_DEF_VAR (nid, "TSLAB", NF_FLOAT, 1, idim2,nvarid)
516#endif
517      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33,
518     .                        "Ecart de la SST (pour slab-ocean)")
519      ierr = NF_ENDDEF(nid)
520#ifdef NC_DOUBLE
521      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tslab)
522#else
523      ierr = NF_PUT_VAR_REAL (nid,nvarid,tslab)
524#endif
525c
526      ierr = NF_REDEF (nid)
527#ifdef NC_DOUBLE
528      ierr = NF_DEF_VAR (nid, "SEAICE", NF_DOUBLE, 1, idim2,nvarid)
529#else
530      ierr = NF_DEF_VAR (nid, "SEAICE", NF_FLOAT, 1, idim2,nvarid)
531#endif
532      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33,
533     .                        "Glace de mer kg/m2 (pour slab-ocean)")
534      ierr = NF_ENDDEF(nid)
535#ifdef NC_DOUBLE
536      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,seaice)
537#else
538      ierr = NF_PUT_VAR_REAL (nid,nvarid,seaice)
539#endif
540c
541      DO nsrf = 1, nbsrf
542        IF (nsrf.LE.99) THEN
543        WRITE(str2,'(i2.2)') nsrf
544        ierr = NF_REDEF (nid)
545#ifdef NC_DOUBLE
546        ierr = NF_DEF_VAR (nid,"QS"//str2,NF_DOUBLE,1,idim2,nvarid)
547#else
548        ierr = NF_DEF_VAR (nid,"QS"//str2,NF_FLOAT,1,idim2,nvarid)
549#endif
550        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 25,
551     .                        "Humidite de surface No."//str2)
552        ierr = NF_ENDDEF(nid)
553        ELSE
554        PRINT*, "Trop de sous-mailles"
555        CALL abort
556        ENDIF
557#ifdef NC_DOUBLE
558      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,qsurf(1,nsrf))
559#else
560      ierr = NF_PUT_VAR_REAL (nid,nvarid,qsurf(1,nsrf))
561#endif
562      END DO
563C
564      ierr = NF_REDEF (nid)
565#ifdef NC_DOUBLE
566      ierr = NF_DEF_VAR (nid,"QSOL",NF_DOUBLE,1,idim2,nvarid)
567#else
568      ierr = NF_DEF_VAR (nid,"QSOL",NF_FLOAT,1,idim2,nvarid)
569#endif
570      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
571     .    "Eau dans le sol (mm)")
572      ierr = NF_ENDDEF(nid)
573#ifdef NC_DOUBLE
574      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,qsol)
575#else
576      ierr = NF_PUT_VAR_REAL (nid,nvarid,qsol)
577#endif
578c
579      DO nsrf = 1, nbsrf
580        IF (nsrf.LE.99) THEN
581        WRITE(str2,'(i2.2)') nsrf
582        ierr = NF_REDEF (nid)
583#ifdef NC_DOUBLE
584        ierr = NF_DEF_VAR (nid,"ALBE"//str2,NF_DOUBLE,1,idim2,nvarid)
585#else
586        ierr = NF_DEF_VAR (nid,"ALBE"//str2,NF_FLOAT,1,idim2,nvarid)
587#endif
588        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23,
589     .                        "albedo de surface No."//str2)
590        ierr = NF_ENDDEF(nid)
591        ELSE
592        PRINT*, "Trop de sous-mailles"
593        CALL abort
594        ENDIF
595#ifdef NC_DOUBLE
596      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,albedo(1,nsrf))
597#else
598      ierr = NF_PUT_VAR_REAL (nid,nvarid,albedo(1,nsrf))
599#endif
600      ENDDO
601
602cIM BEG albedo LW
603        DO nsrf = 1, nbsrf
604        IF (nsrf.LE.99) THEN
605        WRITE(str2,'(i2.2)') nsrf
606        ierr = NF_REDEF (nid)
607#ifdef NC_DOUBLE
608        ierr = NF_DEF_VAR (nid,"ALBLW"//str2,NF_DOUBLE,1,idim2,nvarid)
609#else
610        ierr = NF_DEF_VAR (nid,"ALBLW"//str2,NF_FLOAT,1,idim2,nvarid)
611#endif
612        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23,
613     .                        "albedo LW de surface No."//str2)
614        ierr = NF_ENDDEF(nid)
615        ELSE
616        PRINT*, "Trop de sous-mailles"
617        CALL abort
618        ENDIF
619#ifdef NC_DOUBLE
620      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,alblw(1,nsrf))
621#else
622      ierr = NF_PUT_VAR_REAL (nid,nvarid,alblw(1,nsrf))
623#endif
624      ENDDO
625cIM END albedo LW
626c
627      DO nsrf = 1, nbsrf
628        IF (nsrf.LE.99) THEN
629        WRITE(str2,'(i2.2)') nsrf
630        ierr = NF_REDEF (nid)
631#ifdef NC_DOUBLE
632        ierr = NF_DEF_VAR (nid,"EVAP"//str2,NF_DOUBLE,1,idim2,nvarid)
633#else
634        ierr = NF_DEF_VAR (nid,"EVAP"//str2,NF_FLOAT,1,idim2,nvarid)
635#endif
636        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
637     .                        "Evaporation de surface No."//str2)
638        ierr = NF_ENDDEF(nid)
639        ELSE
640        PRINT*, "Trop de sous-mailles"
641        CALL abort
642        ENDIF
643#ifdef NC_DOUBLE
644      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,evap(1,nsrf))
645#else
646      ierr = NF_PUT_VAR_REAL (nid,nvarid,evap(1,nsrf))
647#endif
648      ENDDO
649
650c
651      DO nsrf = 1, nbsrf
652        IF (nsrf.LE.99) THEN
653        WRITE(str2,'(i2.2)') nsrf
654        ierr = NF_REDEF (nid)
655#ifdef NC_DOUBLE
656        ierr = NF_DEF_VAR (nid,"SNOW"//str2,NF_DOUBLE,1,idim2,nvarid)
657#else
658        ierr = NF_DEF_VAR (nid,"SNOW"//str2,NF_FLOAT,1,idim2,nvarid)
659#endif
660        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22,
661     .                        "Neige de surface No."//str2)
662        ierr = NF_ENDDEF(nid)
663        ELSE
664        PRINT*, "Trop de sous-mailles"
665        CALL abort
666        ENDIF
667#ifdef NC_DOUBLE
668      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,snow(1,nsrf))
669#else
670      ierr = NF_PUT_VAR_REAL (nid,nvarid,snow(1,nsrf))
671#endif
672      ENDDO
673
674c
675      ierr = NF_REDEF (nid)
676#ifdef NC_DOUBLE
677      ierr = NF_DEF_VAR (nid, "RADS", NF_DOUBLE, 1, idim2,nvarid)
678#else
679      ierr = NF_DEF_VAR (nid, "RADS", NF_FLOAT, 1, idim2,nvarid)
680#endif
681      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
682     .                        "Rayonnement net a la surface")
683      ierr = NF_ENDDEF(nid)
684#ifdef NC_DOUBLE
685      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,radsol)
686#else
687      ierr = NF_PUT_VAR_REAL (nid,nvarid,radsol)
688#endif
689c
690      ierr = NF_REDEF (nid)
691#ifdef NC_DOUBLE
692      ierr = NF_DEF_VAR (nid, "solsw", NF_DOUBLE, 1, idim2,nvarid)
693#else
694      ierr = NF_DEF_VAR (nid, "solsw", NF_FLOAT, 1, idim2,nvarid)
695#endif
696      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32,
697     .                        "Rayonnement solaire a la surface")
698      ierr = NF_ENDDEF(nid)
699#ifdef NC_DOUBLE
700      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,solsw)
701#else
702      ierr = NF_PUT_VAR_REAL (nid,nvarid,solsw)
703#endif
704c
705      ierr = NF_REDEF (nid)
706#ifdef NC_DOUBLE
707      ierr = NF_DEF_VAR (nid, "sollw", NF_DOUBLE, 1, idim2,nvarid)
708#else
709      ierr = NF_DEF_VAR (nid, "sollw", NF_FLOAT, 1, idim2,nvarid)
710#endif
711      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 27,
712     .                        "Rayonnement IF a la surface")
713      ierr = NF_ENDDEF(nid)
714#ifdef NC_DOUBLE
715      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,sollw)
716#else
717      ierr = NF_PUT_VAR_REAL (nid,nvarid,sollw)
718#endif
719c
720      ierr = NF_REDEF (nid)
721#ifdef NC_DOUBLE
722      ierr = NF_DEF_VAR (nid, "fder", NF_DOUBLE, 1, idim2,nvarid)
723#else
724      ierr = NF_DEF_VAR (nid, "fder", NF_FLOAT, 1, idim2,nvarid)
725#endif
726      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14,
727     .                        "Derive de flux")
728      ierr = NF_ENDDEF(nid)
729#ifdef NC_DOUBLE
730      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,fder)
731#else
732      ierr = NF_PUT_VAR_REAL (nid,nvarid,fder)
733#endif
734c
735      ierr = NF_REDEF (nid)
736#ifdef NC_DOUBLE
737      ierr = NF_DEF_VAR (nid, "rain_f", NF_DOUBLE, 1, idim2,nvarid)
738#else
739      ierr = NF_DEF_VAR (nid, "rain_f", NF_FLOAT, 1, idim2,nvarid)
740#endif
741      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
742     .                        "precipitation liquide")
743      ierr = NF_ENDDEF(nid)
744#ifdef NC_DOUBLE
745      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rain_fall)
746#else
747      ierr = NF_PUT_VAR_REAL (nid,nvarid,rain_fall)
748#endif
749c
750      ierr = NF_REDEF (nid)
751#ifdef NC_DOUBLE
752      ierr = NF_DEF_VAR (nid, "snow_f", NF_DOUBLE, 1, idim2,nvarid)
753#else
754      ierr = NF_DEF_VAR (nid, "snow_f", NF_FLOAT, 1, idim2,nvarid)
755#endif
756      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
757     .                        "precipitation solide")
758      ierr = NF_ENDDEF(nid)
759#ifdef NC_DOUBLE
760      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,snow_fall)
761#else
762      ierr = NF_PUT_VAR_REAL (nid,nvarid,snow_fall)
763#endif
764c
765       endif
766c$OMP END MASTER
767cc ----> necessaire pour eviter bug openMP sur SX6
768c$OMP MASTER
769      if (phy_rank==0) then
770      DO nsrf = 1, nbsrf
771        IF (nsrf.LE.99) THEN
772        WRITE(str2,'(i2.2)') nsrf
773        ierr = NF_REDEF (nid)
774#ifdef NC_DOUBLE
775        ierr = NF_DEF_VAR (nid,"RUG"//str2,NF_DOUBLE,1,idim2,nvarid)
776#else
777        ierr = NF_DEF_VAR (nid,"RUG"//str2,NF_FLOAT,1,idim2,nvarid)
778#endif
779        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23,
780     .                        "rugosite de surface No."//str2)
781        ierr = NF_ENDDEF(nid)
782        ELSE
783        PRINT*, "Trop de sous-mailles"
784        CALL abort
785        ENDIF
786#ifdef NC_DOUBLE
787      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,frugs(1,nsrf))
788#else
789      ierr = NF_PUT_VAR_REAL (nid,nvarid,frugs(1,nsrf))
790#endif
791      ENDDO
792c
793      DO nsrf = 1, nbsrf
794        IF (nsrf.LE.99) THEN
795            WRITE(str2,'(i2.2)') nsrf
796            ierr = NF_REDEF (nid)
797#ifdef NC_DOUBLE
798            ierr = NF_DEF_VAR (nid,"AGESNO"//str2,NF_DOUBLE,1,idim2
799     $          ,nvarid)
800#else
801            ierr = NF_DEF_VAR (nid,"AGESNO"//str2,NF_FLOAT,1,idim2
802     $          ,nvarid)
803#endif
804            ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 15,
805     .                        "Age de la neige surface No."//str2)
806            ierr = NF_ENDDEF(nid)
807        ELSE
808            PRINT*, "Trop de sous-mailles"
809            CALL abort
810        ENDIF
811#ifdef NC_DOUBLE
812        ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,agesno(1,nsrf))
813#else
814      ierr = NF_PUT_VAR_REAL (nid,nvarid,agesno(1,nsrf))
815#endif
816      ENDDO
817c
818      ierr = NF_REDEF (nid)
819#ifdef NC_DOUBLE
820      ierr = NF_DEF_VAR (nid, "ZMEA", NF_DOUBLE, 1, idim2,nvarid)
821#else
822      ierr = NF_DEF_VAR (nid, "ZMEA", NF_FLOAT, 1, idim2,nvarid)
823#endif
824      ierr = NF_ENDDEF(nid)
825#ifdef NC_DOUBLE
826      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zmea)
827#else
828      ierr = NF_PUT_VAR_REAL (nid,nvarid,zmea)
829#endif
830c
831      ierr = NF_REDEF (nid)
832#ifdef NC_DOUBLE
833      ierr = NF_DEF_VAR (nid, "ZSTD", NF_DOUBLE, 1, idim2,nvarid)
834#else
835      ierr = NF_DEF_VAR (nid, "ZSTD", NF_FLOAT, 1, idim2,nvarid)
836#endif
837      ierr = NF_ENDDEF(nid)
838#ifdef NC_DOUBLE
839      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zstd)
840#else
841      ierr = NF_PUT_VAR_REAL (nid,nvarid,zstd)
842#endif
843      ierr = NF_REDEF (nid)
844#ifdef NC_DOUBLE
845      ierr = NF_DEF_VAR (nid, "ZSIG", NF_DOUBLE, 1, idim2,nvarid)
846#else
847      ierr = NF_DEF_VAR (nid, "ZSIG", NF_FLOAT, 1, idim2,nvarid)
848#endif
849      ierr = NF_ENDDEF(nid)
850#ifdef NC_DOUBLE
851      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zsig)
852#else
853      ierr = NF_PUT_VAR_REAL (nid,nvarid,zsig)
854#endif
855      ierr = NF_REDEF (nid)
856#ifdef NC_DOUBLE
857      ierr = NF_DEF_VAR (nid, "ZGAM", NF_DOUBLE, 1, idim2,nvarid)
858#else
859      ierr = NF_DEF_VAR (nid, "ZGAM", NF_FLOAT, 1, idim2,nvarid)
860#endif
861      ierr = NF_ENDDEF(nid)
862#ifdef NC_DOUBLE
863      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zgam)
864#else
865      ierr = NF_PUT_VAR_REAL (nid,nvarid,zgam)
866#endif
867      ierr = NF_REDEF (nid)
868#ifdef NC_DOUBLE
869      ierr = NF_DEF_VAR (nid, "ZTHE", NF_DOUBLE, 1, idim2,nvarid)
870#else
871      ierr = NF_DEF_VAR (nid, "ZTHE", NF_FLOAT, 1, idim2,nvarid)
872#endif
873      ierr = NF_ENDDEF(nid)
874#ifdef NC_DOUBLE
875      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zthe)
876#else
877      ierr = NF_PUT_VAR_REAL (nid,nvarid,zthe)
878#endif
879      ierr = NF_REDEF (nid)
880#ifdef NC_DOUBLE
881      ierr = NF_DEF_VAR (nid, "ZPIC", NF_DOUBLE, 1, idim2,nvarid)
882#else
883      ierr = NF_DEF_VAR (nid, "ZPIC", NF_FLOAT, 1, idim2,nvarid)
884#endif
885      ierr = NF_ENDDEF(nid)
886#ifdef NC_DOUBLE
887      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zpic)
888#else
889      ierr = NF_PUT_VAR_REAL (nid,nvarid,zpic)
890#endif
891      ierr = NF_REDEF (nid)
892#ifdef NC_DOUBLE
893      ierr = NF_DEF_VAR (nid, "ZVAL", NF_DOUBLE, 1, idim2,nvarid)
894#else
895      ierr = NF_DEF_VAR (nid, "ZVAL", NF_FLOAT, 1, idim2,nvarid)
896#endif
897      ierr = NF_ENDDEF(nid)
898#ifdef NC_DOUBLE
899      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,zval)
900#else
901      ierr = NF_PUT_VAR_REAL (nid,nvarid,zval)
902#endif
903      ierr = NF_REDEF (nid)
904#ifdef NC_DOUBLE
905      ierr = NF_DEF_VAR (nid, "RUGSREL", NF_DOUBLE, 1, idim2,nvarid)
906#else
907      ierr = NF_DEF_VAR (nid, "RUGSREL", NF_FLOAT, 1, idim2,nvarid)
908#endif
909      ierr = NF_ENDDEF(nid)
910#ifdef NC_DOUBLE
911      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rugsrel)
912#else
913      ierr = NF_PUT_VAR_REAL (nid,nvarid,rugsrel)
914#endif
915c
916      ierr = NF_REDEF (nid)
917#ifdef NC_DOUBLE
918      ierr = NF_DEF_VAR (nid, "TANCIEN", NF_DOUBLE, 1, idim3,nvarid)
919#else
920      ierr = NF_DEF_VAR (nid, "TANCIEN", NF_FLOAT, 1, idim3,nvarid)
921#endif
922      ierr = NF_ENDDEF(nid)
923#ifdef NC_DOUBLE
924      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,t_ancien)
925#else
926      ierr = NF_PUT_VAR_REAL (nid,nvarid,t_ancien)
927#endif
928c
929      ierr = NF_REDEF (nid)
930#ifdef NC_DOUBLE
931      ierr = NF_DEF_VAR (nid, "QANCIEN", NF_DOUBLE, 1, idim3,nvarid)
932#else
933      ierr = NF_DEF_VAR (nid, "QANCIEN", NF_FLOAT, 1, idim3,nvarid)
934#endif
935      ierr = NF_ENDDEF(nid)
936#ifdef NC_DOUBLE
937      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q_ancien)
938#else
939      ierr = NF_PUT_VAR_REAL (nid,nvarid,q_ancien)
940#endif
941c
942      ierr = NF_REDEF (nid)
943#ifdef NC_DOUBLE
944      ierr = NF_DEF_VAR (nid, "RUGMER", NF_DOUBLE, 1, idim2,nvarid)
945#else
946      ierr = NF_DEF_VAR (nid, "RUGMER", NF_FLOAT, 1, idim2,nvarid)
947#endif
948      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
949     .                        "Longueur de rugosite sur mer")
950      ierr = NF_ENDDEF(nid)
951#ifdef NC_DOUBLE
952      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,frugs(1,is_oce))
953#else
954      ierr = NF_PUT_VAR_REAL (nid,nvarid,frugs(1,is_oce))
955#endif
956c
957      ierr = NF_REDEF (nid)
958#ifdef NC_DOUBLE
959      ierr = NF_DEF_VAR (nid, "CLWCON", NF_DOUBLE, 1, idim2,nvarid)
960#else
961      ierr = NF_DEF_VAR (nid, "CLWCON", NF_FLOAT, 1, idim2,nvarid)
962#endif
963      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
964     .                        "Eau liquide convective")
965      ierr = NF_ENDDEF(nid)
966#ifdef NC_DOUBLE
967      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,clwcon)
968#else
969      ierr = NF_PUT_VAR_REAL (nid,nvarid,clwcon)
970#endif
971c
972      ierr = NF_REDEF (nid)
973#ifdef NC_DOUBLE
974      ierr = NF_DEF_VAR (nid, "RNEBCON", NF_DOUBLE, 1, idim2,nvarid)
975#else
976      ierr = NF_DEF_VAR (nid, "RNEBCON", NF_FLOAT, 1, idim2,nvarid)
977#endif
978      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
979     .                        "Nebulosite convective")
980      ierr = NF_ENDDEF(nid)
981#ifdef NC_DOUBLE
982      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rnebcon)
983#else
984      ierr = NF_PUT_VAR_REAL (nid,nvarid,rnebcon)
985#endif
986c
987      ierr = NF_REDEF (nid)
988#ifdef NC_DOUBLE
989      ierr = NF_DEF_VAR (nid, "RATQS", NF_DOUBLE, 1, idim2,nvarid)
990#else
991      ierr = NF_DEF_VAR (nid, "RATQS", NF_FLOAT, 1, idim2,nvarid)
992#endif
993      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
994     .                        "Ratqs")
995      ierr = NF_ENDDEF(nid)
996#ifdef NC_DOUBLE
997      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ratqs)
998#else
999      ierr = NF_PUT_VAR_REAL (nid,nvarid,ratqs)
1000#endif
1001c
1002c run_off_lic_0
1003c
1004      ierr = NF_REDEF (nid)
1005#ifdef NC_DOUBLE
1006      ierr=NF_DEF_VAR(nid,"RUNOFFLIC0",NF_DOUBLE,1,idim2,nvarid)
1007#else
1008      ierr=NF_DEF_VAR(nid,"RUNOFFLIC0",NF_FLOAT, 1,idim2,nvarid)
1009#endif
1010      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
1011     .                        "Runofflic0")
1012      ierr = NF_ENDDEF(nid)
1013#ifdef NC_DOUBLE
1014      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,run_off_lic_0)
1015#else
1016      ierr = NF_PUT_VAR_REAL (nid,nvarid,run_off_lic_0)
1017#endif
1018c
1019c
1020      ierr = NF_CLOSE(nid)
1021c
1022      endif   ! phy_rank==0
1023c$OMP END MASTER
1024      RETURN
1025      END
Note: See TracBrowser for help on using the repository browser.