source: LMDZ4/trunk/libf/phy_IPCC_AR4/phyetat0.F @ 902

Last change on this file since 902 was 868, checked in by Laurent Fairhead, 17 years ago

Preparation du remplacement de la physique utilisee pour l'exercice IPCC_AR4
par la version de la physique avec thermique. On garde le repertoire phylmd
pour un petit moment pour que les utilisateurs ne soient pas trop perdus ...
phy_IPCC_AR4 = phylmd
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 45.7 KB
RevLine 
[868]1!
2! $Header$
3!
4c
5c
6      SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
7     .           rlat_p, rlon_p, pctsrf_p, tsol_p,
8     .           ocean_in, ok_veget_in,
9     .           albe_p, alblw_p,
10     .           rain_fall_p, snow_fall_p,solsw_p, sollw_p,
11     .           radsol_p,clesphy0,
12     .           zmea_p,zstd_p,zsig_p,zgam_p,zthe_p,zpic_p,zval_p,
13     .           rugsrel_p,tabcntr0,
14     .           t_ancien_p,q_ancien_p,ancien_ok_p, rnebcon_p, ratqs_p,
15     .           clwcon_p)
16
17      USE dimphy
18      USE mod_grid_phy_lmdz
19      USE mod_phys_lmdz_para
20      USE iophy
21      USE ocean_slab_mod,   ONLY : ocean_slab_init
22      USE ocean_cpl_mod,    ONLY : ocean_cpl_init
23      USE ocean_forced_mod, ONLY : ocean_forced_init
24      USE fonte_neige_mod,  ONLY : fonte_neige_init
25      USE pbl_surface_mod,  ONLY : pbl_surface_init
26      USE surface_data,     ONLY : ocean, ok_veget
27
28      IMPLICIT none
29c======================================================================
30c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
31c Objet: Lecture de l'etat initial pour la physique
32c======================================================================
33#include "dimensions.h"
34#include "netcdf.inc"
35#include "indicesol.h"
36#include "dimsoil.h"
37#include "clesphys.h"
38#include "temps.h"
39c======================================================================
40      CHARACTER*(*) fichnom
41      REAL dtime
42      INTEGER radpas
43      REAL rlat_p(klon), rlon_p(klon)
44      REAL co2_ppm_etat0
45      REAL solaire_etat0
46      REAL tsol_p(klon,nbsrf)
47      REAL tsoil_p(klon,nsoilmx,nbsrf)
48      REAL tslab_p(klon), seaice_p(klon)
49      REAL qsurf_p(klon,nbsrf)
50      REAL qsol_p(klon)
51      REAL snow_p(klon,nbsrf)
52      REAL albe_p(klon,nbsrf)
53cIM BEG alblw
54      REAL alblw_p(klon,nbsrf)
55cIM END alblw
56      REAL evap_p(klon,nbsrf)
57      REAL radsol_p(klon)
58      REAL rain_fall_p(klon)
59      REAL snow_fall_p(klon)
60      REAL sollw_p(klon)
61      real solsw_p(klon)
62      real fder_p(klon)
63      REAL frugs_p(klon,nbsrf)
64      REAL agesno_p(klon,nbsrf)
65      REAL zmea_p(klon)
66      REAL zstd_p(klon)
67      REAL zsig_p(klon)
68      REAL zgam_p(klon)
69      REAL zthe_p(klon)
70      REAL zpic_p(klon)
71      REAL zval_p(klon)
72      REAL rugsrel_p(klon)
73      REAL pctsrf_p(klon, nbsrf)
74      REAL run_off_lic_0_p(klon)
75     
76      REAL t_ancien_p(klon,klev), q_ancien_p(klon,klev)
77      real rnebcon_p(klon,klev),clwcon_p(klon,klev)
78      real ratqs_p(klon,klev)
79      LOGICAL,SAVE ::  ancien_ok
80      LOGICAL ::  ancien_ok_p
81     
82      REAL zmasq_glo(klon_glo)
83      REAL rlat(klon_glo), rlon(klon_glo)
84      REAL tsol(klon_glo,nbsrf)
85      REAL tsoil(klon_glo,nsoilmx,nbsrf)
86cIM "slab" ocean
87      REAL tslab(klon_glo), seaice(klon_glo)
88      REAL qsurf(klon_glo,nbsrf)
89      REAL qsol(klon_glo)
90      REAL snow(klon_glo,nbsrf)
91      REAL albe(klon_glo,nbsrf)
92      REAL alblw(klon_glo,nbsrf)
93      REAL evap(klon_glo,nbsrf)
94      REAL radsol(klon_glo)
95      REAL rain_fall(klon_glo)
96      REAL snow_fall(klon_glo)
97      REAL sollw(klon_glo)
98      real solsw(klon_glo)
99      real fder(klon_glo)
100      REAL frugs(klon_glo,nbsrf)
101      REAL agesno(klon_glo,nbsrf)
102      REAL zmea(klon_glo)
103      REAL zstd(klon_glo)
104      REAL zsig(klon_glo)
105      REAL zgam(klon_glo)
106      REAL zthe(klon_glo)
107      REAL zpic(klon_glo)
108      REAL zval(klon_glo)
109      REAL rugsrel(klon_glo)
110      REAL pctsrf(klon_glo, nbsrf)
111      REAL fractint(klon_glo)
112      REAL run_off_lic_0(klon_glo)
113      REAL t_ancien(klon_glo,klev)
114      REAL q_ancien(klon_glo,klev)
115      real rnebcon(klon_glo,klev)
116      real clwcon(klon_glo,klev)
117      real ratqs(klon_glo,klev)
118
119      CHARACTER*6 ocean_in
120      LOGICAL ok_veget_in
121
122      INTEGER        longcles
123      PARAMETER    ( longcles = 20 )
124      REAL clesphy0( longcles )
125c
126      REAL xmin, xmax
127c
128      INTEGER nid, nvarid
129      INTEGER ierr, i, nsrf, isoil
130      INTEGER length
131      PARAMETER (length=100)
132      REAL tab_cntrl(length), tabcntr0(length)
133      REAL,SAVE :: tab_cntrl_omp(length)
134      CHARACTER*7 str7
135      CHARACTER*2 str2
136      real iolat(jjm+1)
137c
138c Ouvrir le fichier contenant l'etat initial:
139c
140
141c$OMP MASTER
142      print *,'MASTER -x , omp_rank=',omp_rank
143c$OMP END MASTER
144
145c$OMP MASTER
146      IF (is_mpi_root) THEN
147        print*,'fichnom ',fichnom
148        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)
149        IF (ierr.NE.NF_NOERR) THEN
150          write(6,*)' Pb d''ouverture du fichier '//fichnom
151          write(6,*)' ierr = ', ierr
152          CALL ABORT
153        ENDIF
154      ENDIF
155c$OMP END MASTER
156c
157c Lecture des parametres de controle:
158c
159c$OMP MASTER
160      IF (is_mpi_root) THEN
161     
162      ierr = NF_INQ_VARID (nid, "controle", nvarid)
163      IF (ierr.NE.NF_NOERR) THEN
164         PRINT*, 'phyetat0: Le champ <controle> est absent'
165         CALL abort
166      ENDIF
167#ifdef NC_DOUBLE
168      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
169#else
170      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
171#endif
172      IF (ierr.NE.NF_NOERR) THEN
173         PRINT*, 'phyetat0: Lecture echouee pour <controle>'
174         CALL abort
175      ENDIF
176      ENDIF
177
178c$OMP END MASTER
179       
180       CALL bcast(tab_cntrl)
181       
182c
183       
184         DO i = 1, length
185           tabcntr0( i ) = tab_cntrl( i )
186         ENDDO
187c
188         cycle_diurne   = .FALSE.
189         soil_model     = .FALSE.
190         new_oliq       = .FALSE.
191         ok_orodr       = .FALSE.
192         ok_orolf       = .FALSE.
193         ok_limitvrai   = .FALSE.
194
195
196         IF( clesphy0(1).NE.tab_cntrl( 5 ) )  THEN
197             tab_cntrl( 5 ) = clesphy0(1)
198         ENDIF
199
200         IF( clesphy0(2).NE.tab_cntrl( 6 ) )  THEN
201             tab_cntrl( 6 ) = clesphy0(2)
202         ENDIF
203
204         IF( clesphy0(3).NE.tab_cntrl( 7 ) )  THEN
205             tab_cntrl( 7 ) = clesphy0(3)
206         ENDIF
207
208         IF( clesphy0(4).NE.tab_cntrl( 8 ) )  THEN
209             tab_cntrl( 8 ) = clesphy0(4)
210         ENDIF
211
212         IF( clesphy0(5).NE.tab_cntrl( 9 ) )  THEN
213             tab_cntrl( 9 ) = clesphy0( 5 )
214         ENDIF
215
216         IF( clesphy0(6).NE.tab_cntrl( 10 ) )  THEN
217             tab_cntrl( 10 ) = clesphy0( 6 )
218         ENDIF
219
220         IF( clesphy0(7).NE.tab_cntrl( 11 ) )  THEN
221             tab_cntrl( 11 ) = clesphy0( 7 )
222         ENDIF
223
224         IF( clesphy0(8).NE.tab_cntrl( 12 ) )  THEN
225             tab_cntrl( 12 ) = clesphy0( 8 )
226         ENDIF
227
228
229         dtime        = tab_cntrl(1)
230         radpas       = tab_cntrl(2)
231         co2_ppm_etat0      = tab_cntrl(3)
232         solaire_etat0      = tab_cntrl(4)
233         iflag_con    = tab_cntrl(5)
234         nbapp_rad    = tab_cntrl(6)
235
236
237         cycle_diurne    = .FALSE.
238         soil_model      = .FALSE.
239         new_oliq        = .FALSE.
240         ok_orodr        = .FALSE.
241         ok_orolf        = .FALSE.
242         ok_limitvrai    = .FALSE.
243
244         IF( tab_cntrl( 7) .EQ. 1. )    cycle_diurne  = .TRUE.
245         IF( tab_cntrl( 8) .EQ. 1. )       soil_model = .TRUE.
246         IF( tab_cntrl( 9) .EQ. 1. )         new_oliq = .TRUE.
247         IF( tab_cntrl(10) .EQ. 1. )         ok_orodr = .TRUE.
248         IF( tab_cntrl(11) .EQ. 1. )         ok_orolf = .TRUE.
249         IF( tab_cntrl(12) .EQ. 1. )     ok_limitvrai = .TRUE.
250
251
252      itau_phy = tab_cntrl(15)
253
254c
255c Lecture des latitudes (coordonnees):
256c
257c$OMP MASTER
258      IF (is_mpi_root) THEN
259     
260      ierr = NF_INQ_VARID (nid, "latitude", nvarid)
261      IF (ierr.NE.NF_NOERR) THEN
262         PRINT*, 'phyetat0: Le champ <latitude> est absent'
263         CALL abort
264      ENDIF
265#ifdef NC_DOUBLE
266      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlat)
267#else
268      ierr = NF_GET_VAR_REAL(nid, nvarid, rlat)
269#endif
270      IF (ierr.NE.NF_NOERR) THEN
271         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
272         CALL abort
273      ENDIF
274
275c
276c Lecture des longitudes (coordonnees):
277c
278      ierr = NF_INQ_VARID (nid, "longitude", nvarid)
279      IF (ierr.NE.NF_NOERR) THEN
280         PRINT*, 'phyetat0: Le champ <longitude> est absent'
281         CALL abort
282      ENDIF
283#ifdef NC_DOUBLE
284      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlon)
285#else
286      ierr = NF_GET_VAR_REAL(nid, nvarid, rlon)
287#endif
288      IF (ierr.NE.NF_NOERR) THEN
289         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
290         CALL abort
291      ENDIF
292C
293C
294C Lecture du masque terre mer
295C
296
297      ierr = NF_INQ_VARID (nid, "masque", nvarid)
298      IF (ierr .EQ.  NF_NOERR) THEN
299#ifdef NC_DOUBLE
300          ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmasq_glo)
301#else
302          ierr = NF_GET_VAR_REAL(nid, nvarid, zmasq_glo)
303#endif
304          IF (ierr.NE.NF_NOERR) THEN
305              PRINT*, 'phyetat0: Lecture echouee pour <masque>'
306              CALL abort
307          ENDIF
308      else
309          PRINT*, 'phyetat0: Le champ <masque> est absent'
310          PRINT*, 'fichier startphy non compatible avec phyetat0'
311C      CALL abort
312      ENDIF
313
314       
315C Lecture des fractions pour chaque sous-surface
316C
317C initialisation des sous-surfaces
318C
319      pctsrf = 0.
320C
321C fraction de terre
322C
323
324      ierr = NF_INQ_VARID (nid, "FTER", nvarid)
325      IF (ierr .EQ.  NF_NOERR) THEN
326#ifdef NC_DOUBLE
327          ierr = NF_GET_VAR_DOUBLE(nid, nvarid,       
328     .                             pctsrf(1 : klon_glo,is_ter))
329#else
330          ierr = NF_GET_VAR_REAL(nid, nvarid,
331     .                           pctsrf(1 : klon_glo,is_ter))
332#endif
333          IF (ierr.NE.NF_NOERR) THEN
334              PRINT*, 'phyetat0: Lecture echouee pour <FTER>'
335              CALL abort
336          ENDIF
337      else
338          PRINT*, 'phyetat0: Le champ <FTER> est absent'
339c@$$         CALL abort
340      ENDIF
341
342C
343C fraction de glace de terre
344C
345      ierr = NF_INQ_VARID (nid, "FLIC", nvarid)
346      IF (ierr .EQ.  NF_NOERR) THEN
347#ifdef NC_DOUBLE
348          ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
349     .                             pctsrf(1 : klon_glo,is_lic))
350#else
351          ierr = NF_GET_VAR_REAL(nid, nvarid,
352     .                           pctsrf(1 : klon_glo,is_lic))
353#endif
354          IF (ierr.NE.NF_NOERR) THEN
355              PRINT*, 'phyetat0: Lecture echouee pour <FLIC>'
356              CALL abort
357          ENDIF
358      else
359          PRINT*, 'phyetat0: Le champ <FLIC> est absent'
360c@$$         CALL abort
361      ENDIF
362C
363C fraction d'ocean
364C
365      ierr = NF_INQ_VARID (nid, "FOCE", nvarid)
366      IF (ierr .EQ.  NF_NOERR) THEN
367#ifdef NC_DOUBLE
368          ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
369     .                             pctsrf(1 : klon_glo,is_oce))
370#else
371          ierr = NF_GET_VAR_REAL(nid, nvarid,
372     .                           pctsrf(1 : klon_glo,is_oce))
373#endif
374          IF (ierr.NE.NF_NOERR) THEN
375              PRINT*, 'phyetat0: Lecture echouee pour <FOCE>'
376              CALL abort
377          ENDIF
378      else
379          PRINT*, 'phyetat0: Le champ <FOCE> est absent'
380c@$$         CALL abort
381      ENDIF
382
383C
384C fraction glace de mer
385C
386      ierr = NF_INQ_VARID (nid, "FSIC", nvarid)
387      IF (ierr .EQ.  NF_NOERR) THEN
388#ifdef NC_DOUBLE
389          ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
390     .                             pctsrf(1 : klon_glo,is_sic))
391#else
392          ierr = NF_GET_VAR_REAL(nid, nvarid,
393     .                           pctsrf(1 : klon_glo, is_sic))
394#endif
395          IF (ierr.NE.NF_NOERR) THEN
396              PRINT*, 'phyetat0: Lecture echouee pour <FSIC>'
397              CALL abort
398          ENDIF
399      else
400          PRINT*, 'phyetat0: Le champ <FSIC> est absent'
401c@$$         CALL abort
402      ENDIF
403
404C
405C  Verification de l'adequation entre le masque et les sous-surfaces
406C
407      fractint( 1 : klon_glo) = pctsrf(1 : klon_glo, is_ter)
408     $    + pctsrf(1 : klon_glo, is_lic)
409      DO i = 1 , klon_glo
410        IF ( abs(fractint(i) - zmasq_glo(i) ) .GT. EPSFRA ) THEN
411            WRITE(*,*) 'phyetat0: attention fraction terre pas ',
412     $          'coherente ', i, zmasq_glo(i), pctsrf(i, is_ter)
413     $          ,pctsrf(i, is_lic)
414        ENDIF
415      END DO
416      fractint (1 : klon_glo) =  pctsrf(1 : klon_glo, is_oce)
417     $    + pctsrf(1 : klon_glo, is_sic)
418      DO i = 1 , klon_glo
419        IF ( abs( fractint(i) - (1. - zmasq_glo(i))) .GT. EPSFRA ) THEN
420            WRITE(*,*) 'phyetat0 attention fraction ocean pas ',
421     $          'coherente ', i, zmasq_glo(i) , pctsrf(i, is_oce)
422     $          ,pctsrf(i, is_sic)
423        ENDIF
424      END DO
425
426C
427c Lecture des temperatures du sol:
428c
429
430      ierr = NF_INQ_VARID (nid, "TS", nvarid)
431      IF (ierr.NE.NF_NOERR) THEN
432         PRINT*, 'phyetat0: Le champ <TS> est absent'
433         PRINT*, '          Mais je vais essayer de lire TS**'
434         DO nsrf = 1, nbsrf
435           IF (nsrf.GT.99) THEN
436             PRINT*, "Trop de sous-mailles"
437             CALL abort
438           ENDIF
439           WRITE(str2,'(i2.2)') nsrf
440           ierr = NF_INQ_VARID (nid, "TS"//str2, nvarid)
441           IF (ierr.NE.NF_NOERR) THEN
442              PRINT*, "phyetat0: Le champ <TS"//str2//"> est absent"
443              CALL abort
444           ENDIF
445#ifdef NC_DOUBLE
446           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol(1,nsrf))
447#else
448           ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,nsrf))
449#endif
450           IF (ierr.NE.NF_NOERR) THEN
451             PRINT*, "phyetat0: Lecture echouee pour <TS"//str2//">"
452             CALL abort
453           ENDIF
454
455           xmin = 1.0E+20
456           xmax = -1.0E+20
457           DO i = 1, klon_glo
458              xmin = MIN(tsol(i,nsrf),xmin)
459              xmax = MAX(tsol(i,nsrf),xmax)
460           ENDDO
461           PRINT*,'Temperature du sol TS**:', nsrf, xmin, xmax
462         ENDDO
463      ELSE
464         PRINT*, 'phyetat0: Le champ <TS> est present'
465         PRINT*, '          J ignore donc les autres temperatures TS**'
466#ifdef NC_DOUBLE
467         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol(1,1))
468#else
469         ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,1))
470#endif
471         IF (ierr.NE.NF_NOERR) THEN
472            PRINT*, "phyetat0: Lecture echouee pour <TS>"
473            CALL abort
474         ENDIF
475         xmin = 1.0E+20
476         xmax = -1.0E+20
477         DO i = 1, klon_glo
478            xmin = MIN(tsol(i,1),xmin)
479            xmax = MAX(tsol(i,1),xmax)
480         ENDDO
481         PRINT*,'Temperature du sol <TS>', xmin, xmax
482         DO nsrf = 2, nbsrf
483         DO i = 1, klon_glo
484            tsol(i,nsrf) = tsol(i,1)
485         ENDDO
486         ENDDO
487      ENDIF
488
489c
490c Lecture des temperatures du sol profond:
491c
492      DO nsrf = 1, nbsrf
493      DO isoil=1, nsoilmx
494      IF (isoil.GT.99 .AND. nsrf.GT.99) THEN
495         PRINT*, "Trop de couches ou sous-mailles"
496         CALL abort
497      ENDIF
498      WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf
499      ierr = NF_INQ_VARID (nid, 'Tsoil'//str7, nvarid)
500      IF (ierr.NE.NF_NOERR) THEN
501         PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
502         PRINT*, "          Il prend donc la valeur de surface"
503         DO i=1, klon_glo
504             tsoil(i,isoil,nsrf)=tsol(i,nsrf)
505         ENDDO
506      ELSE
507#ifdef NC_DOUBLE
508         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsoil(1,isoil,nsrf))
509#else
510         ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf))
511#endif
512         IF (ierr.NE.NF_NOERR) THEN
513            PRINT*, "Lecture echouee pour <Tsoil"//str7//">"
514            CALL abort
515         ENDIF
516      ENDIF
517      ENDDO
518      ENDDO
519c
520cIM "slab" ocean
521c
522c Lecture de tslab (pour slab ocean seulement):     
523c
524      IF (ocean_in .eq. 'slab  ') then
525        ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
526        IF (ierr.NE.NF_NOERR) THEN
527          PRINT*, "phyetat0: Le champ <TSLAB> est absent"
528          CALL abort
529        ENDIF
530#ifdef NC_DOUBLE
531        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
532#else
533        ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
534#endif
535        IF (ierr.NE.NF_NOERR) THEN
536          PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"
537          CALL abort
538        ENDIF
539        xmin = 1.0E+20
540        xmax = -1.0E+20
541        DO i = 1, klon_glo
542          xmin = MIN(tslab(i),xmin)
543          xmax = MAX(tslab(i),xmax)
544        ENDDO
545        PRINT*,'Min, Max tslab (utilise si OCEAN=slab )', xmin, xmax
546c
547c Lecture de seaice (pour slab ocean seulement):
548c
549        ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)
550        IF (ierr.NE.NF_NOERR) THEN
551          PRINT*, "phyetat0: Le champ <SEAICE> est absent"
552          CALL abort
553        ENDIF
554#ifdef NC_DOUBLE
555        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
556#else
557        ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
558#endif
559        IF (ierr.NE.NF_NOERR) THEN
560          PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"
561          CALL abort
562        ENDIF
563        xmin = 1.0E+20
564        xmax = -1.0E+20
565        DO i = 1, klon_glo
566          xmin = MIN(seaice(i),xmin)
567          xmax = MAX(seaice(i),xmax)
568        ENDDO
569        PRINT*,'Masse de la glace de mer (utilise si OCEAN=slab)',
570     $  xmin, xmax
571      ELSE
572        tslab = 0.
573        seaice = 0.
574      ENDIF
575c
576c Lecture de l'humidite de l'air juste au dessus du sol:
577c
578      ierr = NF_INQ_VARID (nid, "QS", nvarid)
579      IF (ierr.NE.NF_NOERR) THEN
580         PRINT*, 'phyetat0: Le champ <QS> est absent'
581         PRINT*, '          Mais je vais essayer de lire QS**'
582         DO nsrf = 1, nbsrf
583           IF (nsrf.GT.99) THEN
584             PRINT*, "Trop de sous-mailles"
585             CALL abort
586           ENDIF
587           WRITE(str2,'(i2.2)') nsrf
588           ierr = NF_INQ_VARID (nid, "QS"//str2, nvarid)
589           IF (ierr.NE.NF_NOERR) THEN
590              PRINT*, "phyetat0: Le champ <QS"//str2//"> est absent"
591              CALL abort
592           ENDIF
593#ifdef NC_DOUBLE
594           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf))
595#else
596           ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))
597#endif
598           IF (ierr.NE.NF_NOERR) THEN
599             PRINT*, "phyetat0: Lecture echouee pour <QS"//str2//">"
600             CALL abort
601           ENDIF
602           xmin = 1.0E+20
603           xmax = -1.0E+20
604           DO i = 1, klon_glo
605              xmin = MIN(qsurf(i,nsrf),xmin)
606              xmax = MAX(qsurf(i,nsrf),xmax)
607           ENDDO
608           PRINT*,'Humidite pres du sol QS**:', nsrf, xmin, xmax
609         ENDDO
610      ELSE
611         PRINT*, 'phyetat0: Le champ <QS> est present'
612         PRINT*, '          J ignore donc les autres humidites QS**'
613#ifdef NC_DOUBLE
614         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1))
615#else
616         ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1))
617#endif
618         IF (ierr.NE.NF_NOERR) THEN
619            PRINT*, "phyetat0: Lecture echouee pour <QS>"
620            CALL abort
621         ENDIF
622         xmin = 1.0E+20
623         xmax = -1.0E+20
624         DO i = 1, klon_glo
625            xmin = MIN(qsurf(i,1),xmin)
626            xmax = MAX(qsurf(i,1),xmax)
627         ENDDO
628         PRINT*,'Humidite pres du sol <QS>', xmin, xmax
629         DO nsrf = 2, nbsrf
630         DO i = 1, klon_glo
631            qsurf(i,nsrf) = qsurf(i,1)
632         ENDDO
633         ENDDO
634      ENDIF
635C
636C Eau dans le sol (pour le modele de sol "bucket")
637C
638      ierr = NF_INQ_VARID (nid, "QSOL", nvarid)
639      IF (ierr .EQ.  NF_NOERR) THEN
640#ifdef NC_DOUBLE
641          ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol)
642#else
643          ierr = NF_GET_VAR_REAL(nid, nvarid, qsol)
644#endif
645          IF (ierr.NE.NF_NOERR) THEN
646              PRINT*, 'phyetat0: Lecture echouee pour <QSOL>'
647              CALL abort
648          ENDIF
649      else
650          PRINT*, 'phyetat0: Le champ <QSOL> est absent'
651          PRINT*, '          Valeur par defaut nulle'
652          qsol(:)=0.
653c@$$         CALL abort
654      ENDIF
655      xmin = 1.0E+20
656      xmax = -1.0E+20
657      DO i = 1, klon_glo
658        xmin = MIN(qsol(i),xmin)
659        xmax = MAX(qsol(i),xmax)
660      ENDDO
661      PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax
662c
663c Lecture de neige au sol:
664c
665      ierr = NF_INQ_VARID (nid, "SNOW", nvarid)
666      IF (ierr.NE.NF_NOERR) THEN
667         PRINT*, 'phyetat0: Le champ <SNOW> est absent'
668         PRINT*, '          Mais je vais essayer de lire SNOW**'
669         DO nsrf = 1, nbsrf
670           IF (nsrf.GT.99) THEN
671             PRINT*, "Trop de sous-mailles"
672             CALL abort
673           ENDIF
674           WRITE(str2,'(i2.2)') nsrf
675           ierr = NF_INQ_VARID (nid, "SNOW"//str2, nvarid)
676           IF (ierr.NE.NF_NOERR) THEN
677              PRINT*, "phyetat0: Le champ <SNOW"//str2//"> est absent"
678              CALL abort
679           ENDIF
680#ifdef NC_DOUBLE
681           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,nsrf))
682#else
683           ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf))
684#endif
685           IF (ierr.NE.NF_NOERR) THEN
686             PRINT*, "phyetat0: Lecture echouee pour <SNOW"//str2//">"
687             CALL abort
688           ENDIF
689           xmin = 1.0E+20
690           xmax = -1.0E+20
691           DO i = 1, klon_glo
692              xmin = MIN(snow(i,nsrf),xmin)
693              xmax = MAX(snow(i,nsrf),xmax)
694           ENDDO
695           PRINT*,'Neige du sol SNOW**:', nsrf, xmin, xmax
696         ENDDO
697      ELSE
698         PRINT*, 'phyetat0: Le champ <SNOW> est present'
699         PRINT*, '          J ignore donc les autres neiges SNOW**'
700#ifdef NC_DOUBLE
701         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,1))
702#else
703         ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,1))
704#endif
705         IF (ierr.NE.NF_NOERR) THEN
706            PRINT*, "phyetat0: Lecture echouee pour <SNOW>"
707            CALL abort
708         ENDIF
709         xmin = 1.0E+20
710         xmax = -1.0E+20
711         DO i = 1, klon_glo
712            xmin = MIN(snow(i,1),xmin)
713            xmax = MAX(snow(i,1),xmax)
714         ENDDO
715         PRINT*,'Neige du sol <SNOW>', xmin, xmax
716         DO nsrf = 2, nbsrf
717         DO i = 1, klon_glo
718            snow(i,nsrf) = snow(i,1)
719         ENDDO
720         ENDDO
721      ENDIF
722c
723c Lecture de albedo au sol:
724c
725      ierr = NF_INQ_VARID (nid, "ALBE", nvarid)
726      IF (ierr.NE.NF_NOERR) THEN
727         PRINT*, 'phyetat0: Le champ <ALBE> est absent'
728         PRINT*, '          Mais je vais essayer de lire ALBE**'
729         DO nsrf = 1, nbsrf
730           IF (nsrf.GT.99) THEN
731             PRINT*, "Trop de sous-mailles"
732             CALL abort
733           ENDIF
734           WRITE(str2,'(i2.2)') nsrf
735           ierr = NF_INQ_VARID (nid, "ALBE"//str2, nvarid)
736           IF (ierr.NE.NF_NOERR) THEN
737              PRINT*, "phyetat0: Le champ <ALBE"//str2//"> est absent"
738              CALL abort
739           ENDIF
740#ifdef NC_DOUBLE
741           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,nsrf))
742#else
743           ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,nsrf))
744#endif
745           IF (ierr.NE.NF_NOERR) THEN
746             PRINT*, "phyetat0: Lecture echouee pour <ALBE"//str2//">"
747             CALL abort
748           ENDIF
749           xmin = 1.0E+20
750           xmax = -1.0E+20
751           DO i = 1, klon_glo
752              xmin = MIN(albe(i,nsrf),xmin)
753              xmax = MAX(albe(i,nsrf),xmax)
754           ENDDO
755           PRINT*,'Albedo du sol ALBE**:', nsrf, xmin, xmax
756         ENDDO
757      ELSE
758         PRINT*, 'phyetat0: Le champ <ALBE> est present'
759         PRINT*, '          J ignore donc les autres ALBE**'
760#ifdef NC_DOUBLE
761         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,1))
762#else
763         ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,1))
764#endif
765         IF (ierr.NE.NF_NOERR) THEN
766            PRINT*, "phyetat0: Lecture echouee pour <ALBE>"
767            CALL abort
768         ENDIF
769         xmin = 1.0E+20
770         xmax = -1.0E+20
771         DO i = 1, klon_glo
772            xmin = MIN(albe(i,1),xmin)
773            xmax = MAX(albe(i,1),xmax)
774         ENDDO
775         PRINT*,'Neige du sol <ALBE>', xmin, xmax
776         DO nsrf = 2, nbsrf
777         DO i = 1, klon_glo
778            albe(i,nsrf) = albe(i,1)
779         ENDDO
780         ENDDO
781      ENDIF
782
783c
784c Lecture de albedo au sol LW:
785c
786      ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)
787      IF (ierr.NE.NF_NOERR) THEN
788         PRINT*, 'phyetat0: Le champ <ALBLW> est absent'
789c        PRINT*, '          Mais je vais essayer de lire ALBLW**'
790         PRINT*, '          Mais je vais prendre ALBE**'
791         DO nsrf = 1, nbsrf
792           DO i = 1, klon_glo
793             alblw(i,nsrf) = albe(i,nsrf)
794           ENDDO
795         ENDDO
796      ELSE
797         PRINT*, 'phyetat0: Le champ <ALBLW> est present'
798         PRINT*, '          J ignore donc les autres ALBLW**'
799#ifdef NC_DOUBLE
800         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alblw(1,1))
801#else
802         ierr = NF_GET_VAR_REAL(nid, nvarid, alblw(1,1))
803#endif
804         IF (ierr.NE.NF_NOERR) THEN
805            PRINT*, "phyetat0: Lecture echouee pour <ALBLW>"
806            CALL abort
807         ENDIF
808         xmin = 1.0E+20
809         xmax = -1.0E+20
810         DO i = 1, klon_glo
811            xmin = MIN(alblw(i,1),xmin)
812            xmax = MAX(alblw(i,1),xmax)
813         ENDDO
814         PRINT*,'Neige du sol <ALBLW>', xmin, xmax
815         DO nsrf = 2, nbsrf
816         DO i = 1, klon_glo
817            alblw(i,nsrf) = alblw(i,1)
818         ENDDO
819         ENDDO
820      ENDIF
821c
822c Lecture de evaporation: 
823c
824      ierr = NF_INQ_VARID (nid, "EVAP", nvarid)
825      IF (ierr.NE.NF_NOERR) THEN
826         PRINT*, 'phyetat0: Le champ <EVAP> est absent'
827         PRINT*, '          Mais je vais essayer de lire EVAP**'
828         DO nsrf = 1, nbsrf
829           IF (nsrf.GT.99) THEN
830             PRINT*, "Trop de sous-mailles"
831             CALL abort
832           ENDIF
833           WRITE(str2,'(i2.2)') nsrf
834           ierr = NF_INQ_VARID (nid, "EVAP"//str2, nvarid)
835           IF (ierr.NE.NF_NOERR) THEN
836              PRINT*, "phyetat0: Le champ <EVAP"//str2//"> est absent"
837              CALL abort
838           ENDIF
839#ifdef NC_DOUBLE
840           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,nsrf))
841#else
842           ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf))
843#endif
844           IF (ierr.NE.NF_NOERR) THEN
845             PRINT*, "phyetat0: Lecture echouee pour <EVAP"//str2//">"
846             CALL abort
847           ENDIF
848           xmin = 1.0E+20
849           xmax = -1.0E+20
850           DO i = 1, klon_glo
851              xmin = MIN(evap(i,nsrf),xmin)
852              xmax = MAX(evap(i,nsrf),xmax)
853           ENDDO
854           PRINT*,'evap du sol EVAP**:', nsrf, xmin, xmax
855         ENDDO
856      ELSE
857         PRINT*, 'phyetat0: Le champ <EVAP> est present'
858         PRINT*, '          J ignore donc les autres EVAP**'
859#ifdef NC_DOUBLE
860         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,1))
861#else
862         ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,1))
863#endif
864         IF (ierr.NE.NF_NOERR) THEN
865            PRINT*, "phyetat0: Lecture echouee pour <EVAP>"
866            CALL abort
867         ENDIF
868         xmin = 1.0E+20
869         xmax = -1.0E+20
870         DO i = 1, klon_glo
871            xmin = MIN(evap(i,1),xmin)
872            xmax = MAX(evap(i,1),xmax)
873         ENDDO
874         PRINT*,'Evap du sol <EVAP>', xmin, xmax
875         DO nsrf = 2, nbsrf
876         DO i = 1, klon_glo
877            evap(i,nsrf) = evap(i,1)
878         ENDDO
879         ENDDO
880      ENDIF
881c
882c Lecture precipitation liquide:
883c
884      ierr = NF_INQ_VARID (nid, "rain_f", nvarid)
885      IF (ierr.NE.NF_NOERR) THEN
886         PRINT*, 'phyetat0: Le champ <rain_f> est absent'
887         CALL abort
888      ENDIF
889#ifdef NC_DOUBLE
890      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rain_fall)
891#else
892      ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall)
893#endif
894      IF (ierr.NE.NF_NOERR) THEN
895         PRINT*, 'phyetat0: Lecture echouee pour <rain_f>'
896         CALL abort
897      ENDIF
898      xmin = 1.0E+20
899      xmax = -1.0E+20
900      DO i = 1, klon_glo
901         xmin = MIN(rain_fall(i),xmin)
902         xmax = MAX(rain_fall(i),xmax)
903      ENDDO
904      PRINT*,'Precipitation liquide rain_f:', xmin, xmax
905c
906c Lecture precipitation solide:
907c
908      ierr = NF_INQ_VARID (nid, "snow_f", nvarid)
909      IF (ierr.NE.NF_NOERR) THEN
910         PRINT*, 'phyetat0: Le champ <snow_f> est absent'
911         CALL abort
912      ENDIF
913#ifdef NC_DOUBLE
914      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow_fall)
915#else
916      ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall)
917#endif
918      IF (ierr.NE.NF_NOERR) THEN
919         PRINT*, 'phyetat0: Lecture echouee pour <snow_f>'
920         CALL abort
921      ENDIF
922      xmin = 1.0E+20
923      xmax = -1.0E+20
924      DO i = 1, klon_glo
925         xmin = MIN(snow_fall(i),xmin)
926         xmax = MAX(snow_fall(i),xmax)
927      ENDDO
928      PRINT*,'Precipitation solide snow_f:', xmin, xmax
929c
930c Lecture rayonnement solaire au sol:
931c
932      ierr = NF_INQ_VARID (nid, "solsw", nvarid)
933      IF (ierr.NE.NF_NOERR) THEN
934         PRINT*, 'phyetat0: Le champ <solsw> est absent'
935         PRINT*, 'mis a zero'
936         solsw = 0.
937      ELSE
938#ifdef NC_DOUBLE
939        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, solsw)
940#else
941        ierr = NF_GET_VAR_REAL(nid, nvarid, solsw)
942#endif
943        IF (ierr.NE.NF_NOERR) THEN
944          PRINT*, 'phyetat0: Lecture echouee pour <solsw>'
945          CALL abort
946        ENDIF
947      ENDIF
948      xmin = 1.0E+20
949      xmax = -1.0E+20
950      DO i = 1, klon_glo
951         xmin = MIN(solsw(i),xmin)
952         xmax = MAX(solsw(i),xmax)
953      ENDDO
954      PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax
955c
956c Lecture rayonnement IF au sol:
957c
958      ierr = NF_INQ_VARID (nid, "sollw", nvarid)
959      IF (ierr.NE.NF_NOERR) THEN
960         PRINT*, 'phyetat0: Le champ <sollw> est absent'
961         PRINT*, 'mis a zero'
962         sollw = 0.
963      ELSE
964#ifdef NC_DOUBLE
965        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, sollw)
966#else
967        ierr = NF_GET_VAR_REAL(nid, nvarid, sollw)
968#endif
969        IF (ierr.NE.NF_NOERR) THEN
970          PRINT*, 'phyetat0: Lecture echouee pour <sollw>'
971          CALL abort
972        ENDIF
973      ENDIF
974      xmin = 1.0E+20
975      xmax = -1.0E+20
976      DO i = 1, klon_glo
977         xmin = MIN(sollw(i),xmin)
978         xmax = MAX(sollw(i),xmax)
979      ENDDO
980      PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax
981     
982      ENDIF  ! is_mpi_root
983c$OMP END MASTER
984
985
986c$OMP MASTER
987      IF (is_mpi_root) THEN
988c
989c Lecture derive des flux:
990c
991      ierr = NF_INQ_VARID (nid, "fder", nvarid)
992      IF (ierr.NE.NF_NOERR) THEN
993         PRINT*, 'phyetat0: Le champ <fder> est absent'
994         PRINT*, 'mis a zero'
995         fder = 0.
996      ELSE
997#ifdef NC_DOUBLE
998        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder)
999#else
1000        ierr = NF_GET_VAR_REAL(nid, nvarid, fder)
1001#endif
1002        IF (ierr.NE.NF_NOERR) THEN
1003          PRINT*, 'phyetat0: Lecture echouee pour <fder>'
1004          CALL abort
1005        ENDIF
1006      ENDIF
1007      xmin = 1.0E+20
1008      xmax = -1.0E+20
1009      DO i = 1, klon_glo
1010         xmin = MIN(fder(i),xmin)
1011         xmax = MAX(fder(i),xmax)
1012      ENDDO
1013      PRINT*,'Derive des flux fder:', xmin, xmax
1014
1015c
1016c Lecture du rayonnement net au sol:
1017c
1018      ierr = NF_INQ_VARID (nid, "RADS", nvarid)
1019      IF (ierr.NE.NF_NOERR) THEN
1020         PRINT*, 'phyetat0: Le champ <RADS> est absent'
1021         CALL abort
1022      ENDIF
1023#ifdef NC_DOUBLE
1024      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, radsol)
1025#else
1026      ierr = NF_GET_VAR_REAL(nid, nvarid, radsol)
1027#endif
1028      IF (ierr.NE.NF_NOERR) THEN
1029         PRINT*, 'phyetat0: Lecture echouee pour <RADS>'
1030         CALL abort
1031      ENDIF
1032      xmin = 1.0E+20
1033      xmax = -1.0E+20
1034      DO i = 1, klon_glo
1035         xmin = MIN(radsol(i),xmin)
1036         xmax = MAX(radsol(i),xmax)
1037      ENDDO
1038      PRINT*,'Rayonnement net au sol radsol:', xmin, xmax
1039c
1040c Lecture de la longueur de rugosite
1041c
1042c
1043      ierr = NF_INQ_VARID (nid, "RUG", nvarid)
1044      IF (ierr.NE.NF_NOERR) THEN
1045         PRINT*, 'phyetat0: Le champ <RUG> est absent'
1046         PRINT*, '          Mais je vais essayer de lire RUG**'
1047         DO nsrf = 1, nbsrf
1048           IF (nsrf.GT.99) THEN
1049             PRINT*, "Trop de sous-mailles"
1050             CALL abort
1051           ENDIF
1052           WRITE(str2,'(i2.2)') nsrf
1053           ierr = NF_INQ_VARID (nid, "RUG"//str2, nvarid)
1054           IF (ierr.NE.NF_NOERR) THEN
1055              PRINT*, "phyetat0: Le champ <RUG"//str2//"> est absent"
1056              CALL abort
1057           ENDIF
1058#ifdef NC_DOUBLE
1059           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,nsrf))
1060#else
1061           ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf))
1062#endif
1063           IF (ierr.NE.NF_NOERR) THEN
1064             PRINT*, "phyetat0: Lecture echouee pour <RUG"//str2//">"
1065             CALL abort
1066           ENDIF
1067           xmin = 1.0E+20
1068           xmax = -1.0E+20
1069           DO i = 1, klon_glo
1070              xmin = MIN(frugs(i,nsrf),xmin)
1071              xmax = MAX(frugs(i,nsrf),xmax)
1072           ENDDO
1073           PRINT*,'rugosite du sol RUG**:', nsrf, xmin, xmax
1074         ENDDO
1075      ELSE
1076         PRINT*, 'phyetat0: Le champ <RUG> est present'
1077         PRINT*, '          J ignore donc les autres RUG**'
1078#ifdef NC_DOUBLE
1079         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,1))
1080#else
1081         ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,1))
1082#endif
1083         IF (ierr.NE.NF_NOERR) THEN
1084            PRINT*, "phyetat0: Lecture echouee pour <RUG>"
1085            CALL abort
1086         ENDIF
1087         xmin = 1.0E+20
1088         xmax = -1.0E+20
1089         DO i = 1, klon_glo
1090            xmin = MIN(frugs(i,1),xmin)
1091            xmax = MAX(frugs(i,1),xmax)
1092         ENDDO
1093         PRINT*,'rugosite <RUG>', xmin, xmax
1094         DO nsrf = 2, nbsrf
1095         DO i = 1, klon_glo
1096            frugs(i,nsrf) = frugs(i,1)
1097         ENDDO
1098         ENDDO
1099      ENDIF
1100
1101c
1102c Lecture de l'age de la neige:
1103c
1104      ierr = NF_INQ_VARID (nid, "AGESNO", nvarid)
1105      IF (ierr.NE.NF_NOERR) THEN
1106         PRINT*, 'phyetat0: Le champ <AGESNO> est absent'
1107         PRINT*, '          Mais je vais essayer de lire AGESNO**'
1108         DO nsrf = 1, nbsrf
1109           IF (nsrf.GT.99) THEN
1110             PRINT*, "Trop de sous-mailles"
1111             CALL abort
1112           ENDIF
1113           WRITE(str2,'(i2.2)') nsrf
1114           ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid)
1115           IF (ierr.NE.NF_NOERR) THEN
1116              PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
1117              agesno = 50.0
1118           ENDIF
1119#ifdef NC_DOUBLE
1120           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,nsrf))
1121#else
1122           ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))
1123#endif
1124           IF (ierr.NE.NF_NOERR) THEN
1125             PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
1126             CALL abort
1127           ENDIF
1128           xmin = 1.0E+20
1129           xmax = -1.0E+20
1130           DO i = 1, klon_glo
1131              xmin = MIN(agesno(i,nsrf),xmin)
1132              xmax = MAX(agesno(i,nsrf),xmax)
1133           ENDDO
1134           PRINT*,'Age de la neige AGESNO**:', nsrf, xmin, xmax
1135         ENDDO
1136      ELSE
1137         PRINT*, 'phyetat0: Le champ <AGESNO> est present'
1138         PRINT*, '          J ignore donc les autres AGESNO**'
1139#ifdef NC_DOUBLE
1140         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,1))
1141#else
1142         ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1))
1143#endif
1144         IF (ierr.NE.NF_NOERR) THEN
1145            PRINT*, "phyetat0: Lecture echouee pour <AGESNO>"
1146            CALL abort
1147         ENDIF
1148         xmin = 1.0E+20
1149         xmax = -1.0E+20
1150         DO i = 1, klon_glo
1151            xmin = MIN(agesno(i,1),xmin)
1152            xmax = MAX(agesno(i,1),xmax)
1153         ENDDO
1154         PRINT*,'Age de la neige <AGESNO>', xmin, xmax
1155         DO nsrf = 2, nbsrf
1156         DO i = 1, klon_glo
1157            agesno(i,nsrf) = agesno(i,1)
1158         ENDDO
1159         ENDDO
1160      ENDIF
1161
1162c
1163      ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)
1164      IF (ierr.NE.NF_NOERR) THEN
1165         PRINT*, 'phyetat0: Le champ <ZMEA> est absent'
1166         CALL abort
1167      ENDIF
1168#ifdef NC_DOUBLE
1169      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmea)
1170#else
1171      ierr = NF_GET_VAR_REAL(nid, nvarid, zmea)
1172#endif
1173      IF (ierr.NE.NF_NOERR) THEN
1174         PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>'
1175         CALL abort
1176      ENDIF
1177      xmin = 1.0E+20
1178      xmax = -1.0E+20
1179      DO i = 1, klon_glo
1180         xmin = MIN(zmea(i),xmin)
1181         xmax = MAX(zmea(i),xmax)
1182      ENDDO
1183      PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
1184c
1185c
1186      ierr = NF_INQ_VARID (nid, "ZSTD", nvarid)
1187      IF (ierr.NE.NF_NOERR) THEN
1188         PRINT*, 'phyetat0: Le champ <ZSTD> est absent'
1189         CALL abort
1190      ENDIF
1191#ifdef NC_DOUBLE
1192      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zstd)
1193#else
1194      ierr = NF_GET_VAR_REAL(nid, nvarid, zstd)
1195#endif
1196      IF (ierr.NE.NF_NOERR) THEN
1197         PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>'
1198         CALL abort
1199      ENDIF
1200      xmin = 1.0E+20
1201      xmax = -1.0E+20
1202      DO i = 1, klon_glo
1203         xmin = MIN(zstd(i),xmin)
1204         xmax = MAX(zstd(i),xmax)
1205      ENDDO
1206      PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
1207c
1208c
1209      ierr = NF_INQ_VARID (nid, "ZSIG", nvarid)
1210      IF (ierr.NE.NF_NOERR) THEN
1211         PRINT*, 'phyetat0: Le champ <ZSIG> est absent'
1212         CALL abort
1213      ENDIF
1214#ifdef NC_DOUBLE
1215      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zsig)
1216#else
1217      ierr = NF_GET_VAR_REAL(nid, nvarid, zsig)
1218#endif
1219      IF (ierr.NE.NF_NOERR) THEN
1220         PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>'
1221         CALL abort
1222      ENDIF
1223      xmin = 1.0E+20
1224      xmax = -1.0E+20
1225      DO i = 1, klon_glo
1226         xmin = MIN(zsig(i),xmin)
1227         xmax = MAX(zsig(i),xmax)
1228      ENDDO
1229      PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
1230c
1231c
1232      ierr = NF_INQ_VARID (nid, "ZGAM", nvarid)
1233      IF (ierr.NE.NF_NOERR) THEN
1234         PRINT*, 'phyetat0: Le champ <ZGAM> est absent'
1235         CALL abort
1236      ENDIF
1237#ifdef NC_DOUBLE
1238      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zgam)
1239#else
1240      ierr = NF_GET_VAR_REAL(nid, nvarid, zgam)
1241#endif
1242      IF (ierr.NE.NF_NOERR) THEN
1243         PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>'
1244         CALL abort
1245      ENDIF
1246      xmin = 1.0E+20
1247      xmax = -1.0E+20
1248      DO i = 1, klon_glo
1249         xmin = MIN(zgam(i),xmin)
1250         xmax = MAX(zgam(i),xmax)
1251      ENDDO
1252      PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
1253c
1254c
1255      ierr = NF_INQ_VARID (nid, "ZTHE", nvarid)
1256      IF (ierr.NE.NF_NOERR) THEN
1257         PRINT*, 'phyetat0: Le champ <ZTHE> est absent'
1258         CALL abort
1259      ENDIF
1260#ifdef NC_DOUBLE
1261      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zthe)
1262#else
1263      ierr = NF_GET_VAR_REAL(nid, nvarid, zthe)
1264#endif
1265      IF (ierr.NE.NF_NOERR) THEN
1266         PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>'
1267         CALL abort
1268      ENDIF
1269      xmin = 1.0E+20
1270      xmax = -1.0E+20
1271      DO i = 1, klon_glo
1272         xmin = MIN(zthe(i),xmin)
1273         xmax = MAX(zthe(i),xmax)
1274      ENDDO
1275      PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
1276c
1277c
1278      ierr = NF_INQ_VARID (nid, "ZPIC", nvarid)
1279      IF (ierr.NE.NF_NOERR) THEN
1280         PRINT*, 'phyetat0: Le champ <ZPIC> est absent'
1281         CALL abort
1282      ENDIF
1283#ifdef NC_DOUBLE
1284      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zpic)
1285#else
1286      ierr = NF_GET_VAR_REAL(nid, nvarid, zpic)
1287#endif
1288      IF (ierr.NE.NF_NOERR) THEN
1289         PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>'
1290         CALL abort
1291      ENDIF
1292      xmin = 1.0E+20
1293      xmax = -1.0E+20
1294      DO i = 1, klon_glo
1295         xmin = MIN(zpic(i),xmin)
1296         xmax = MAX(zpic(i),xmax)
1297      ENDDO
1298      PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
1299c
1300      ierr = NF_INQ_VARID (nid, "ZVAL", nvarid)
1301      IF (ierr.NE.NF_NOERR) THEN
1302         PRINT*, 'phyetat0: Le champ <ZVAL> est absent'
1303         CALL abort
1304      ENDIF
1305#ifdef NC_DOUBLE
1306      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zval)
1307#else
1308      ierr = NF_GET_VAR_REAL(nid, nvarid, zval)
1309#endif
1310      IF (ierr.NE.NF_NOERR) THEN
1311         PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>'
1312         CALL abort
1313      ENDIF
1314      xmin = 1.0E+20
1315      xmax = -1.0E+20
1316      DO i = 1, klon_glo
1317         xmin = MIN(zval(i),xmin)
1318         xmax = MAX(zval(i),xmax)
1319      ENDDO
1320      PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
1321c
1322c
1323      ierr = NF_INQ_VARID (nid, "RUGSREL", nvarid)
1324      IF (ierr.NE.NF_NOERR) THEN
1325         PRINT*, 'phyetat0: Le champ <RUGSREL> est absent'
1326         CALL abort
1327      ENDIF
1328#ifdef NC_DOUBLE
1329      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rugsrel)
1330#else
1331      ierr = NF_GET_VAR_REAL(nid, nvarid, rugsrel)
1332#endif
1333      IF (ierr.NE.NF_NOERR) THEN
1334         PRINT*, 'phyetat0: Lecture echouee pour <RUGSREL>'
1335         CALL abort
1336      ENDIF
1337      xmin = 1.0E+20
1338      xmax = -1.0E+20
1339      DO i = 1, klon_glo
1340         xmin = MIN(rugsrel(i),xmin)
1341         xmax = MAX(rugsrel(i),xmax)
1342      ENDDO
1343      PRINT*,'Rugosite relief (ecart-type) rugsrel:', xmin, xmax
1344c
1345c
1346      ancien_ok = .TRUE.
1347c
1348      ierr = NF_INQ_VARID (nid, "TANCIEN", nvarid)
1349      IF (ierr.NE.NF_NOERR) THEN
1350         PRINT*, "phyetat0: Le champ <TANCIEN> est absent"
1351         PRINT*, "Depart legerement fausse. Mais je continue"
1352         ancien_ok = .FALSE.
1353      ELSE
1354#ifdef NC_DOUBLE
1355         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, t_ancien)
1356#else
1357         ierr = NF_GET_VAR_REAL(nid, nvarid, t_ancien)
1358#endif
1359         IF (ierr.NE.NF_NOERR) THEN
1360            PRINT*, "phyetat0: Lecture echouee pour <TANCIEN>"
1361            CALL abort
1362         ENDIF
1363      ENDIF
1364c
1365      ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
1366      IF (ierr.NE.NF_NOERR) THEN
1367         PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
1368         PRINT*, "Depart legerement fausse. Mais je continue"
1369         ancien_ok = .FALSE.
1370      ELSE
1371#ifdef NC_DOUBLE
1372         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien)
1373#else
1374         ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)
1375#endif
1376         IF (ierr.NE.NF_NOERR) THEN
1377            PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
1378            CALL abort
1379         ENDIF
1380      ENDIF
1381c
1382      clwcon=0.
1383      ierr = NF_INQ_VARID (nid, "CLWCON", nvarid)
1384      IF (ierr.NE.NF_NOERR) THEN
1385         PRINT*, "phyetat0: Le champ CLWCON est absent"
1386         PRINT*, "Depart legerement fausse. Mais je continue"
1387         clwcon = 0.
1388      ELSE
1389#ifdef NC_DOUBLE
1390         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, clwcon)
1391#else
1392         ierr = NF_GET_VAR_REAL(nid, nvarid, clwcon)
1393#endif
1394         IF (ierr.NE.NF_NOERR) THEN
1395            PRINT*, "phyetat0: Lecture echouee pour <CLWCON>"
1396            CALL abort
1397         ENDIF
1398      ENDIF
1399      xmin = 1.0E+20
1400      xmax = -1.0E+20
1401      xmin = MINval(clwcon)
1402      xmax = MAXval(clwcon)
1403      PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
1404c
1405      rnebcon=0.
1406      ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid)
1407      IF (ierr.NE.NF_NOERR) THEN
1408         PRINT*, "phyetat0: Le champ RNEBCON est absent"
1409         PRINT*, "Depart legerement fausse. Mais je continue"
1410         rnebcon = 0.
1411      ELSE
1412#ifdef NC_DOUBLE
1413         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rnebcon)
1414#else
1415         ierr = NF_GET_VAR_REAL(nid, nvarid, rnebcon)
1416#endif
1417         IF (ierr.NE.NF_NOERR) THEN
1418            PRINT*, "phyetat0: Lecture echouee pour <RNEBCON>"
1419            CALL abort
1420         ENDIF
1421      ENDIF
1422      xmin = 1.0E+20
1423      xmax = -1.0E+20
1424      xmin = MINval(rnebcon)
1425      xmax = MAXval(rnebcon)
1426      PRINT*,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax
1427
1428c
1429      ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
1430      IF (ierr.NE.NF_NOERR) THEN
1431         PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
1432         PRINT*, "Depart legerement fausse. Mais je continue"
1433         ancien_ok = .FALSE.
1434      ELSE
1435#ifdef NC_DOUBLE
1436         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien)
1437#else
1438         ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)
1439#endif
1440         IF (ierr.NE.NF_NOERR) THEN
1441            PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
1442            CALL abort
1443         ENDIF
1444      ENDIF
1445c
1446c Lecture ratqs
1447c
1448       ratqs=0.
1449      ierr = NF_INQ_VARID (nid, "RATQS", nvarid)
1450      IF (ierr.NE.NF_NOERR) THEN
1451         PRINT*, "phyetat0: Le champ <RATQS> est absent"
1452         PRINT*, "Depart legerement fausse. Mais je continue"
1453         ratqs = 0.
1454      ELSE
1455#ifdef NC_DOUBLE
1456         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ratqs)
1457#else
1458         ierr = NF_GET_VAR_REAL(nid, nvarid, ratqs)
1459#endif
1460         IF (ierr.NE.NF_NOERR) THEN
1461            PRINT*, "phyetat0: Lecture echouee pour <RATQS>"
1462            CALL abort
1463         ENDIF
1464      ENDIF
1465      xmin = 1.0E+20
1466      xmax = -1.0E+20
1467      xmin = MINval(ratqs)
1468      xmax = MAXval(ratqs)
1469      PRINT*,'(ecart-type) ratqs:', xmin, xmax
1470c
1471c Lecture run_off_lic_0
1472c
1473      ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid)
1474      IF (ierr.NE.NF_NOERR) THEN
1475         PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent"
1476         PRINT*, "Depart legerement fausse. Mais je continue"
1477         run_off_lic_0 = 0.
1478      ELSE
1479#ifdef NC_DOUBLE
1480         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, run_off_lic_0)
1481#else
1482         ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0)
1483#endif
1484         IF (ierr.NE.NF_NOERR) THEN
1485            PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>"
1486            CALL abort
1487         ENDIF
1488      ENDIF
1489      xmin = 1.0E+20
1490      xmax = -1.0E+20
1491      xmin = MINval(run_off_lic_0)
1492      xmax = MAXval(run_off_lic_0)
1493      PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax
1494c
1495c Fermer le fichier:
1496c
1497      ierr = NF_CLOSE(nid)
1498      ENDIF ! is_mpi_root
1499c
1500c$OMP END MASTER
1501
1502c$OMP MASTER
1503cym  en attendant mieux
1504        iolat(1)=rlat(1)
1505        iolat(jjm+1)=rlat(klon_glo)
1506        do i=2,jjm
1507          iolat(i)=rlat(2+(i-2)*iim)
1508        enddo
1509        CALL bcast_mpi(iolat)
1510        CALL bcast_mpi(rlon)
1511        call init_iophy(iolat,rlon(2:iim+1))
1512       
1513c$OMP END MASTER
1514       
1515      call Scatter( rlat,rlat_p)
1516      call Scatter( rlon,rlon_p)
1517      call Scatter( tsol,tsol_p)
1518      call Scatter( tsoil,tsoil_p)
1519      call Scatter( tslab,tslab_p)
1520      call Scatter( seaice,seaice_p)
1521      call Scatter( qsurf,qsurf_p)
1522      call Scatter( qsol,qsol_p)
1523      call Scatter( snow,snow_p)
1524      call Scatter( albe,albe_p)
1525      call Scatter( alblw,alblw_p)
1526      call Scatter( evap,evap_p)
1527      call Scatter( radsol,radsol_p)
1528      call Scatter( rain_fall,rain_fall_p)
1529      call Scatter( snow_fall,snow_fall_p)
1530      call Scatter( sollw,sollw_p)
1531      call Scatter( solsw,solsw_p)
1532      call Scatter( fder,fder_p)
1533      call Scatter( frugs,frugs_p)
1534      call Scatter( agesno,agesno_p)
1535      call Scatter( zmea,zmea_p)
1536      call Scatter( zstd,zstd_p)
1537      call Scatter( zsig,zsig_p)
1538      call Scatter( zgam,zgam_p)
1539      call Scatter( zthe,zthe_p)
1540      call Scatter( zpic,zpic_p)
1541      call Scatter( zval,zval_p)
1542      call Scatter( rugsrel,rugsrel_p)
1543      call Scatter( pctsrf,pctsrf_p)
1544      call Scatter( run_off_lic_0,run_off_lic_0_p)
1545      call Scatter( t_ancien,t_ancien_p)
1546      call Scatter( q_ancien,q_ancien_p)
1547      call Scatter( rnebcon,rnebcon_p)
1548      call Scatter( clwcon,clwcon_p)
1549      call Scatter( ratqs,ratqs_p)
1550      call Scatter( zmasq_glo,zmasq)
1551
1552c
1553c Initilalize variables in module surface_data
1554c
1555      ok_veget = ok_veget_in
1556      ocean    = ocean_in
1557c
1558c Initialize module pbl_surface_mod
1559c
1560      CALL pbl_surface_init(qsol_p, fder_p, snow_p, qsurf_p,
1561     $     evap_p, frugs_p, agesno_p, tsoil_p)
1562
1563c Initialize ocean module according to ocean type
1564      IF ( ocean == 'slab' ) THEN
1565c        initilalize module ocean_slab_init
1566         CALL ocean_slab_init(dtime, tslab_p, seaice_p, pctsrf_p)
1567      ELSEIF ( ocean == 'couple' ) THEN
1568c        initilalize module ocean_cpl_init
1569         CALL ocean_cpl_init(dtime, rlon_p, rlat_p)
1570      ELSE
1571c        initilalize module ocean_forced_init
1572         CALL ocean_forced_init
1573      ENDIF
1574c
1575c Initilialize module fonte_neige_mod     
1576c
1577      CALL fonte_neige_init(run_off_lic_0_p)
1578
1579
1580      RETURN
1581      END
Note: See TracBrowser for help on using the repository browser.