/[lmdze]/trunk/Sources/phylmd/phyetat0.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/phyetat0.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 43 by guez, Fri Apr 8 12:43:31 2011 UTC
# Line 1  Line 1 
1  module phyetat0_m  module phyetat0_m
2    
3    use dimphy, only: klon, klev, zmasq    use dimphy, only: klon
4      
5    IMPLICIT none    IMPLICIT none
6    
7    REAL, save:: rlat(klon), rlon(klon)    REAL, save:: rlat(klon), rlon(klon) ! latitude and longitude, in degrees
   ! latitude et longitude pour chaque point, in degrees  
8    
9    private klon, klev, zmasq    private klon
10    
11  contains  contains
12    
13    SUBROUTINE phyetat0(fichnom,dtime,co2_ppm_etat0,solaire_etat0, &    SUBROUTINE phyetat0(fichnom, pctsrf, tsol,tsoil, ocean, tslab,seaice, &
        pctsrf, tsol,tsoil, &  
        ocean, tslab,seaice, & !IM "slab" ocean  
14         qsurf,qsol,snow, &         qsurf,qsol,snow, &
15         albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, &         albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, &
16         fder,radsol,frugs,agesno,clesphy0, &         fder,radsol,frugs,agesno, &
17         zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0, &         zmea,zstd,zsig,zgam,zthe,zpic,zval, &
18         t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon, &         t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon, &
19         run_off_lic_0)         run_off_lic_0)
20    
21      ! From phylmd/phyetat0.F,v 1.4 2005/06/03 10:03:07      ! From phylmd/phyetat0.F,v 1.4 2005/06/03 10:03:07
   
     use dimens_m  
     use indicesol  
     use dimsoil  
     use temps  
     use clesphys  
   
22      ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818      ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
23      ! Objet: Lecture de l'etat initial pour la physique      ! Objet: Lecture de l'etat initial pour la physique
24    
25        USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf
26        USE dimsoil, ONLY : nsoilmx
27        USE temps, ONLY : itau_phy
28        use netcdf, only: nf90_get_att, nf90_global, nf90_inq_varid, NF90_NOERR
29        use netcdf95, only: handle_err, nf95_get_var
30        use dimphy, only: zmasq, klev
31    
32      include "netcdf.inc"      include "netcdf.inc"
33    
34      CHARACTER*(*) fichnom      CHARACTER(len=*) fichnom
     REAL, intent(out):: dtime  
     INTEGER radpas  
     REAL, intent(out):: co2_ppm_etat0  
     REAL, intent(out):: solaire_etat0  
35      REAL tsol(klon,nbsrf)      REAL tsol(klon,nbsrf)
36      REAL tsoil(klon,nsoilmx,nbsrf)      REAL tsoil(klon,nsoilmx,nbsrf)
37      !IM "slab" ocean      !IM "slab" ocean
# Line 58  contains Line 51  contains
51      REAL frugs(klon,nbsrf)      REAL frugs(klon,nbsrf)
52      REAL agesno(klon,nbsrf)      REAL agesno(klon,nbsrf)
53      REAL zmea(klon)      REAL zmea(klon)
54      REAL zstd(klon)      REAL, intent(out):: zstd(klon)
55      REAL zsig(klon)      REAL, intent(out):: zsig(klon)
56      REAL zgam(klon)      REAL zgam(klon)
57      REAL zthe(klon)      REAL zthe(klon)
58      REAL zpic(klon)      REAL zpic(klon)
59      REAL zval(klon)      REAL zval(klon)
     REAL rugsrel(klon)  
60      REAL pctsrf(klon, nbsrf)      REAL pctsrf(klon, nbsrf)
61      REAL fractint(klon)      REAL fractint(klon)
62      REAL run_off_lic_0(klon)      REAL run_off_lic_0(klon)
# Line 73  contains Line 65  contains
65      real rnebcon(klon,klev),clwcon(klon,klev),ratqs(klon,klev)      real rnebcon(klon,klev),clwcon(klon,klev),ratqs(klon,klev)
66      LOGICAL ancien_ok      LOGICAL ancien_ok
67    
68      CHARACTER*6 ocean      CHARACTER(len=*), intent(in):: ocean
   
     INTEGER        longcles  
     PARAMETER    ( longcles = 20 )  
     REAL clesphy0( longcles )  
69    
70      REAL xmin, xmax      REAL xmin, xmax
71    
72      INTEGER nid, nvarid      INTEGER nid, nvarid
73      INTEGER ierr, i, nsrf, isoil      INTEGER ierr, i, nsrf, isoil
     INTEGER length  
     PARAMETER (length=100)  
     REAL tab_cntrl(length), tabcntr0(length)  
74      CHARACTER*7 str7      CHARACTER*7 str7
75      CHARACTER*2 str2      CHARACTER*2 str2
76    
# Line 97  contains Line 82  contains
82    
83      print *, 'fichnom = ', fichnom      print *, 'fichnom = ', fichnom
84      ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)      ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)
85      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
86         write(6,*)' Pb d''ouverture du fichier '//fichnom         write(6,*)' Pb d''ouverture du fichier '//fichnom
87         write(6,*)' ierr = ', ierr         write(6,*)' ierr = ', ierr
88         STOP 1         STOP 1
89      ENDIF      ENDIF
90    
91      ! Lecture des parametres de controle:      ierr = nf90_get_att(nid, nf90_global, "itau_phy", itau_phy)
92        call handle_err("phyetat0 itau_phy", ierr, nid, nf90_global)
     ierr = NF_INQ_VARID (nid, "controle", nvarid)  
     IF (ierr.NE.NF_NOERR) THEN  
        PRINT*, 'phyetat0: Le champ <controle> est absent'  
        stop 1  
     ENDIF  
     ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)  
     IF (ierr.NE.NF_NOERR) THEN  
        PRINT*, 'phyetat0: Lecture echouee pour <controle>'  
        stop 1  
     ELSE  
        DO i = 1, length  
           tabcntr0( i ) = tab_cntrl( i )  
        ENDDO  
   
        cycle_diurne   = .FALSE.  
        soil_model     = .FALSE.  
        new_oliq       = .FALSE.  
        ok_orodr       = .FALSE.  
        ok_orolf       = .FALSE.  
        ok_limitvrai   = .FALSE.  
   
   
        IF( clesphy0(1).NE.tab_cntrl( 5 ) )  THEN  
           tab_cntrl( 5 ) = clesphy0(1)  
        ENDIF  
   
        IF( clesphy0(2).NE.tab_cntrl( 6 ) )  THEN  
           tab_cntrl( 6 ) = clesphy0(2)  
        ENDIF  
   
        IF( clesphy0(3).NE.tab_cntrl( 7 ) )  THEN  
           tab_cntrl( 7 ) = clesphy0(3)  
        ENDIF  
   
        IF( clesphy0(4).NE.tab_cntrl( 8 ) )  THEN  
           tab_cntrl( 8 ) = clesphy0(4)  
        ENDIF  
   
        IF( clesphy0(5).NE.tab_cntrl( 9 ) )  THEN  
           tab_cntrl( 9 ) = clesphy0( 5 )  
        ENDIF  
   
        IF( clesphy0(6).NE.tab_cntrl( 10 ) )  THEN  
           tab_cntrl( 10 ) = clesphy0( 6 )  
        ENDIF  
   
        IF( clesphy0(7).NE.tab_cntrl( 11 ) )  THEN  
           tab_cntrl( 11 ) = clesphy0( 7 )  
        ENDIF  
   
        IF( clesphy0(8).NE.tab_cntrl( 12 ) )  THEN  
           tab_cntrl( 12 ) = clesphy0( 8 )  
        ENDIF  
   
   
        dtime        = tab_cntrl(1)  
        radpas       = int(tab_cntrl(2))  
        co2_ppm_etat0      = tab_cntrl(3)  
        solaire_etat0      = tab_cntrl(4)  
        iflag_con    = tab_cntrl(5)  
        nbapp_rad    = tab_cntrl(6)  
   
   
        cycle_diurne    = .FALSE.  
        soil_model      = .FALSE.  
        new_oliq        = .FALSE.  
        ok_orodr        = .FALSE.  
        ok_orolf        = .FALSE.  
        ok_limitvrai    = .FALSE.  
   
        IF( tab_cntrl( 7) .EQ. 1. )    cycle_diurne  = .TRUE.  
        IF( tab_cntrl( 8) .EQ. 1. )       soil_model = .TRUE.  
        IF( tab_cntrl( 9) .EQ. 1. )         new_oliq = .TRUE.  
        IF( tab_cntrl(10) .EQ. 1. )         ok_orodr = .TRUE.  
        IF( tab_cntrl(11) .EQ. 1. )         ok_orolf = .TRUE.  
        IF( tab_cntrl(12) .EQ. 1. )     ok_limitvrai = .TRUE.  
     ENDIF  
   
     itau_phy = tab_cntrl(15)  
   
93    
94      ! Lecture des latitudes (coordonnees):      ! Lecture des latitudes (coordonnees):
95    
96      ierr = NF_INQ_VARID (nid, "latitude", nvarid)      ierr = NF90_INQ_VARID (nid, "latitude", nvarid)
97      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
98         PRINT*, 'phyetat0: Le champ <latitude> est absent'         PRINT *, 'phyetat0: Le champ <latitude> est absent'
99         stop 1         stop 1
100      ENDIF      ENDIF
101      ierr = NF_GET_VAR_REAL(nid, nvarid, rlat)      ierr = NF_GET_VAR_REAL(nid, nvarid, rlat)
102      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
103         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'         PRINT *, 'phyetat0: Lecture echouee pour <latitude>'
104         stop 1         stop 1
105      ENDIF      ENDIF
106    
107      ! Lecture des longitudes (coordonnees):      ! Lecture des longitudes (coordonnees):
108    
109      ierr = NF_INQ_VARID (nid, "longitude", nvarid)      ierr = NF90_INQ_VARID (nid, "longitude", nvarid)
110      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
111         PRINT*, 'phyetat0: Le champ <longitude> est absent'         PRINT *, 'phyetat0: Le champ <longitude> est absent'
112         stop 1         stop 1
113      ENDIF      ENDIF
114      ierr = NF_GET_VAR_REAL(nid, nvarid, rlon)      ierr = NF_GET_VAR_REAL(nid, nvarid, rlon)
115      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
116         PRINT*, 'phyetat0: Lecture echouee pour <latitude>'         PRINT *, 'phyetat0: Lecture echouee pour <latitude>'
117         stop 1         stop 1
118      ENDIF      ENDIF
119    
120    
121      ! Lecture du masque terre mer      ! Lecture du masque terre mer
122    
123      ierr = NF_INQ_VARID (nid, "masque", nvarid)      ierr = NF90_INQ_VARID (nid, "masque", nvarid)
124      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
125         ierr = NF_GET_VAR_REAL(nid, nvarid, zmasq)         call nf95_get_var(nid, nvarid, zmasq)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <masque>'  
           stop 1  
        ENDIF  
126      else      else
127         PRINT*, 'phyetat0: Le champ <masque> est absent'         PRINT *, 'phyetat0: Le champ <masque> est absent'
128         PRINT*, 'fichier startphy non compatible avec phyetat0'         PRINT *, 'fichier startphy non compatible avec phyetat0'
129         !      stop 1         !      stop 1
130      ENDIF      ENDIF
131      ! Lecture des fractions pour chaque sous-surface      ! Lecture des fractions pour chaque sous-surface
# Line 235  contains Line 136  contains
136    
137      ! fraction de terre      ! fraction de terre
138    
139      ierr = NF_INQ_VARID (nid, "FTER", nvarid)      ierr = NF90_INQ_VARID (nid, "FTER", nvarid)
140      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
141         ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_ter))         call nf95_get_var(nid, nvarid, pctsrf(1 : klon,is_ter))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <FTER>'  
           stop 1  
        ENDIF  
142      else      else
143         PRINT*, 'phyetat0: Le champ <FTER> est absent'         PRINT *, 'phyetat0: Le champ <FTER> est absent'
144         !$$$         stop 1         !$$$         stop 1
145      ENDIF      ENDIF
146    
147      ! fraction de glace de terre      ! fraction de glace de terre
148    
149      ierr = NF_INQ_VARID (nid, "FLIC", nvarid)      ierr = NF90_INQ_VARID (nid, "FLIC", nvarid)
150      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
151         ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_lic))         call nf95_get_var(nid, nvarid, pctsrf(1 : klon,is_lic))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <FLIC>'  
           stop 1  
        ENDIF  
152      else      else
153         PRINT*, 'phyetat0: Le champ <FLIC> est absent'         PRINT *, 'phyetat0: Le champ <FLIC> est absent'
154         !$$$         stop 1         !$$$         stop 1
155      ENDIF      ENDIF
156    
157      ! fraction d'ocean      ! fraction d'ocean
158    
159      ierr = NF_INQ_VARID (nid, "FOCE", nvarid)      ierr = NF90_INQ_VARID (nid, "FOCE", nvarid)
160      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
161         ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_oce))         call nf95_get_var(nid, nvarid, pctsrf(1 : klon,is_oce))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <FOCE>'  
           stop 1  
        ENDIF  
162      else      else
163         PRINT*, 'phyetat0: Le champ <FOCE> est absent'         PRINT *, 'phyetat0: Le champ <FOCE> est absent'
164         !$$$         stop 1         !$$$         stop 1
165      ENDIF      ENDIF
166    
167      ! fraction glace de mer      ! fraction glace de mer
168    
169      ierr = NF_INQ_VARID (nid, "FSIC", nvarid)      ierr = NF90_INQ_VARID (nid, "FSIC", nvarid)
170      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
171         ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon, is_sic))         call nf95_get_var(nid, nvarid, pctsrf(1 : klon, is_sic))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <FSIC>'  
           stop 1  
        ENDIF  
172      else      else
173         PRINT*, 'phyetat0: Le champ <FSIC> est absent'         PRINT *, 'phyetat0: Le champ <FSIC> est absent'
174         !$$$         stop 1         !$$$         stop 1
175      ENDIF      ENDIF
176    
# Line 312  contains Line 197  contains
197    
198      ! Lecture des temperatures du sol:      ! Lecture des temperatures du sol:
199    
200      ierr = NF_INQ_VARID (nid, "TS", nvarid)      ierr = NF90_INQ_VARID (nid, "TS", nvarid)
201      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
202         PRINT*, 'phyetat0: Le champ <TS> est absent'         PRINT *, 'phyetat0: Le champ <TS> est absent'
203         PRINT*, '          Mais je vais essayer de lire TS**'         PRINT *, '          Mais je vais essayer de lire TS**'
204         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
205            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
206               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
207               stop 1               stop 1
208            ENDIF            ENDIF
209            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
210            ierr = NF_INQ_VARID (nid, "TS"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "TS"//str2, nvarid)
211            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
212               PRINT*, "phyetat0: Le champ <TS"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <TS"//str2//"> est absent"
213               stop 1               stop 1
214            ENDIF            ENDIF
215            ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,nsrf))
216            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
217               PRINT*, "phyetat0: Lecture echouee pour <TS"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <TS"//str2//">"
218               stop 1               stop 1
219            ENDIF            ENDIF
220            xmin = 1.0E+20            xmin = 1.0E+20
# Line 338  contains Line 223  contains
223               xmin = MIN(tsol(i,nsrf),xmin)               xmin = MIN(tsol(i,nsrf),xmin)
224               xmax = MAX(tsol(i,nsrf),xmax)               xmax = MAX(tsol(i,nsrf),xmax)
225            ENDDO            ENDDO
226            PRINT*,'Temperature du sol TS**:', nsrf, xmin, xmax            PRINT *,'Temperature du sol TS**:', nsrf, xmin, xmax
227         ENDDO         ENDDO
228      ELSE      ELSE
229         PRINT*, 'phyetat0: Le champ <TS> est present'         PRINT *, 'phyetat0: Le champ <TS> est present'
230         PRINT*, '          J ignore donc les autres temperatures TS**'         PRINT *, '          J ignore donc les autres temperatures TS**'
231         ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,1))         call nf95_get_var(nid, nvarid, tsol(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <TS>"  
           stop 1  
        ENDIF  
232         xmin = 1.0E+20         xmin = 1.0E+20
233         xmax = -1.0E+20         xmax = -1.0E+20
234         DO i = 1, klon         DO i = 1, klon
235            xmin = MIN(tsol(i,1),xmin)            xmin = MIN(tsol(i,1),xmin)
236            xmax = MAX(tsol(i,1),xmax)            xmax = MAX(tsol(i,1),xmax)
237         ENDDO         ENDDO
238         PRINT*,'Temperature du sol <TS>', xmin, xmax         PRINT *,'Temperature du sol <TS>', xmin, xmax
239         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
240            DO i = 1, klon            DO i = 1, klon
241               tsol(i,nsrf) = tsol(i,1)               tsol(i,nsrf) = tsol(i,1)
# Line 367  contains Line 248  contains
248      DO nsrf = 1, nbsrf      DO nsrf = 1, nbsrf
249         DO isoil=1, nsoilmx         DO isoil=1, nsoilmx
250            IF (isoil.GT.99 .AND. nsrf.GT.99) THEN            IF (isoil.GT.99 .AND. nsrf.GT.99) THEN
251               PRINT*, "Trop de couches ou sous-mailles"               PRINT *, "Trop de couches ou sous-mailles"
252               stop 1               stop 1
253            ENDIF            ENDIF
254            WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf            WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf
255            ierr = NF_INQ_VARID (nid, 'Tsoil'//str7, nvarid)            ierr = NF90_INQ_VARID (nid, 'Tsoil'//str7, nvarid)
256            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
257               PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"               PRINT *, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
258               PRINT*, "          Il prend donc la valeur de surface"               PRINT *, "          Il prend donc la valeur de surface"
259               DO i=1, klon               DO i=1, klon
260                  tsoil(i,isoil,nsrf)=tsol(i,nsrf)                  tsoil(i,isoil,nsrf)=tsol(i,nsrf)
261               ENDDO               ENDDO
262            ELSE            ELSE
263               ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf))               ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf))
264               IF (ierr.NE.NF_NOERR) THEN               IF (ierr.NE.NF90_NOERR) THEN
265                  PRINT*, "Lecture echouee pour <Tsoil"//str7//">"                  PRINT *, "Lecture echouee pour <Tsoil"//str7//">"
266                  stop 1                  stop 1
267               ENDIF               ENDIF
268            ENDIF            ENDIF
# Line 393  contains Line 274  contains
274      ! Lecture de tslab (pour slab ocean seulement):            ! Lecture de tslab (pour slab ocean seulement):      
275    
276      IF (ocean .eq. 'slab  ') then      IF (ocean .eq. 'slab  ') then
277         ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)         ierr = NF90_INQ_VARID (nid, "TSLAB", nvarid)
278         IF (ierr.NE.NF_NOERR) THEN         IF (ierr.NE.NF90_NOERR) THEN
279            PRINT*, "phyetat0: Le champ <TSLAB> est absent"            PRINT *, "phyetat0: Le champ <TSLAB> est absent"
           stop 1  
        ENDIF  
        ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)  
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"  
280            stop 1            stop 1
281         ENDIF         ENDIF
282           call nf95_get_var(nid, nvarid, tslab)
283         xmin = 1.0E+20         xmin = 1.0E+20
284         xmax = -1.0E+20         xmax = -1.0E+20
285         DO i = 1, klon         DO i = 1, klon
286            xmin = MIN(tslab(i),xmin)            xmin = MIN(tslab(i),xmin)
287            xmax = MAX(tslab(i),xmax)            xmax = MAX(tslab(i),xmax)
288         ENDDO         ENDDO
289         PRINT*,'Ecart de la SST tslab:', xmin, xmax         PRINT *,'Ecart de la SST tslab:', xmin, xmax
290    
291         ! Lecture de seaice (pour slab ocean seulement):         ! Lecture de seaice (pour slab ocean seulement):
292    
293         ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)         ierr = NF90_INQ_VARID (nid, "SEAICE", nvarid)
294         IF (ierr.NE.NF_NOERR) THEN         IF (ierr.NE.NF90_NOERR) THEN
295            PRINT*, "phyetat0: Le champ <SEAICE> est absent"            PRINT *, "phyetat0: Le champ <SEAICE> est absent"
           stop 1  
        ENDIF  
        ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)  
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"  
296            stop 1            stop 1
297         ENDIF         ENDIF
298           call nf95_get_var(nid, nvarid, seaice)
299         xmin = 1.0E+20         xmin = 1.0E+20
300         xmax = -1.0E+20         xmax = -1.0E+20
301         DO i = 1, klon         DO i = 1, klon
302            xmin = MIN(seaice(i),xmin)            xmin = MIN(seaice(i),xmin)
303            xmax = MAX(seaice(i),xmax)            xmax = MAX(seaice(i),xmax)
304         ENDDO         ENDDO
305         PRINT*,'Masse de la glace de mer seaice:', xmin, xmax         PRINT *,'Masse de la glace de mer seaice:', xmin, xmax
306      ELSE      ELSE
307         tslab = 0.         tslab = 0.
308         seaice = 0.         seaice = 0.
# Line 437  contains Line 310  contains
310    
311      ! Lecture de l'humidite de l'air juste au dessus du sol:      ! Lecture de l'humidite de l'air juste au dessus du sol:
312    
313      ierr = NF_INQ_VARID (nid, "QS", nvarid)      ierr = NF90_INQ_VARID (nid, "QS", nvarid)
314      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
315         PRINT*, 'phyetat0: Le champ <QS> est absent'         PRINT *, 'phyetat0: Le champ <QS> est absent'
316         PRINT*, '          Mais je vais essayer de lire QS**'         PRINT *, '          Mais je vais essayer de lire QS**'
317         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
318            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
319               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
320               stop 1               stop 1
321            ENDIF            ENDIF
322            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
323            ierr = NF_INQ_VARID (nid, "QS"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "QS"//str2, nvarid)
324            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
325               PRINT*, "phyetat0: Le champ <QS"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <QS"//str2//"> est absent"
326               stop 1               stop 1
327            ENDIF            ENDIF
328            ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))
329            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
330               PRINT*, "phyetat0: Lecture echouee pour <QS"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <QS"//str2//">"
331               stop 1               stop 1
332            ENDIF            ENDIF
333            xmin = 1.0E+20            xmin = 1.0E+20
# Line 463  contains Line 336  contains
336               xmin = MIN(qsurf(i,nsrf),xmin)               xmin = MIN(qsurf(i,nsrf),xmin)
337               xmax = MAX(qsurf(i,nsrf),xmax)               xmax = MAX(qsurf(i,nsrf),xmax)
338            ENDDO            ENDDO
339            PRINT*,'Humidite pres du sol QS**:', nsrf, xmin, xmax            PRINT *,'Humidite pres du sol QS**:', nsrf, xmin, xmax
340         ENDDO         ENDDO
341      ELSE      ELSE
342         PRINT*, 'phyetat0: Le champ <QS> est present'         PRINT *, 'phyetat0: Le champ <QS> est present'
343         PRINT*, '          J ignore donc les autres humidites QS**'         PRINT *, '          J ignore donc les autres humidites QS**'
344         ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1))         call nf95_get_var(nid, nvarid, qsurf(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <QS>"  
           stop 1  
        ENDIF  
345         xmin = 1.0E+20         xmin = 1.0E+20
346         xmax = -1.0E+20         xmax = -1.0E+20
347         DO i = 1, klon         DO i = 1, klon
348            xmin = MIN(qsurf(i,1),xmin)            xmin = MIN(qsurf(i,1),xmin)
349            xmax = MAX(qsurf(i,1),xmax)            xmax = MAX(qsurf(i,1),xmax)
350         ENDDO         ENDDO
351         PRINT*,'Humidite pres du sol <QS>', xmin, xmax         PRINT *,'Humidite pres du sol <QS>', xmin, xmax
352         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
353            DO i = 1, klon            DO i = 1, klon
354               qsurf(i,nsrf) = qsurf(i,1)               qsurf(i,nsrf) = qsurf(i,1)
# Line 489  contains Line 358  contains
358    
359      ! Eau dans le sol (pour le modele de sol "bucket")      ! Eau dans le sol (pour le modele de sol "bucket")
360    
361      ierr = NF_INQ_VARID (nid, "QSOL", nvarid)      ierr = NF90_INQ_VARID(nid, "QSOL", nvarid)
362      IF (ierr .EQ.  NF_NOERR) THEN      IF (ierr ==  NF90_NOERR) THEN
363         ierr = NF_GET_VAR_REAL(nid, nvarid, qsol)         call nf95_get_var(nid, nvarid, qsol)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <QSOL>'  
           stop 1  
        ENDIF  
364      else      else
365         PRINT*, 'phyetat0: Le champ <QSOL> est absent'         PRINT *, 'phyetat0: Le champ <QSOL> est absent'
366         PRINT*, '          Valeur par defaut nulle'         PRINT *, '          Valeur par defaut nulle'
367         qsol(:)=0.         qsol = 0.
        !$$$         stop 1  
368      ENDIF      ENDIF
369      xmin = 1.0E+20      xmin = 1.0E+20
370      xmax = -1.0E+20      xmax = -1.0E+20
# Line 508  contains Line 372  contains
372         xmin = MIN(qsol(i),xmin)         xmin = MIN(qsol(i),xmin)
373         xmax = MAX(qsol(i),xmax)         xmax = MAX(qsol(i),xmax)
374      ENDDO      ENDDO
375      PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax      PRINT *,'Eau dans le sol (mm) <QSOL>', xmin, xmax
376    
377      ! Lecture de neige au sol:      ! Lecture de neige au sol:
378    
379      ierr = NF_INQ_VARID (nid, "SNOW", nvarid)      ierr = NF90_INQ_VARID (nid, "SNOW", nvarid)
380      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
381         PRINT*, 'phyetat0: Le champ <SNOW> est absent'         PRINT *, 'phyetat0: Le champ <SNOW> est absent'
382         PRINT*, '          Mais je vais essayer de lire SNOW**'         PRINT *, '          Mais je vais essayer de lire SNOW**'
383         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
384            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
385               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
386               stop 1               stop 1
387            ENDIF            ENDIF
388            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
389            ierr = NF_INQ_VARID (nid, "SNOW"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "SNOW"//str2, nvarid)
390            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
391               PRINT*, "phyetat0: Le champ <SNOW"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <SNOW"//str2//"> est absent"
392               stop 1               stop 1
393            ENDIF            ENDIF
394            ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf))
395            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
396               PRINT*, "phyetat0: Lecture echouee pour <SNOW"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <SNOW"//str2//">"
397               stop 1               stop 1
398            ENDIF            ENDIF
399            xmin = 1.0E+20            xmin = 1.0E+20
# Line 538  contains Line 402  contains
402               xmin = MIN(snow(i,nsrf),xmin)               xmin = MIN(snow(i,nsrf),xmin)
403               xmax = MAX(snow(i,nsrf),xmax)               xmax = MAX(snow(i,nsrf),xmax)
404            ENDDO            ENDDO
405            PRINT*,'Neige du sol SNOW**:', nsrf, xmin, xmax            PRINT *,'Neige du sol SNOW**:', nsrf, xmin, xmax
406         ENDDO         ENDDO
407      ELSE      ELSE
408         PRINT*, 'phyetat0: Le champ <SNOW> est present'         PRINT *, 'phyetat0: Le champ <SNOW> est present'
409         PRINT*, '          J ignore donc les autres neiges SNOW**'         PRINT *, '          J ignore donc les autres neiges SNOW**'
410         ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,1))         call nf95_get_var(nid, nvarid, snow(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <SNOW>"  
           stop 1  
        ENDIF  
411         xmin = 1.0E+20         xmin = 1.0E+20
412         xmax = -1.0E+20         xmax = -1.0E+20
413         DO i = 1, klon         DO i = 1, klon
414            xmin = MIN(snow(i,1),xmin)            xmin = MIN(snow(i,1),xmin)
415            xmax = MAX(snow(i,1),xmax)            xmax = MAX(snow(i,1),xmax)
416         ENDDO         ENDDO
417         PRINT*,'Neige du sol <SNOW>', xmin, xmax         PRINT *,'Neige du sol <SNOW>', xmin, xmax
418         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
419            DO i = 1, klon            DO i = 1, klon
420               snow(i,nsrf) = snow(i,1)               snow(i,nsrf) = snow(i,1)
# Line 564  contains Line 424  contains
424    
425      ! Lecture de albedo au sol:      ! Lecture de albedo au sol:
426    
427      ierr = NF_INQ_VARID (nid, "ALBE", nvarid)      ierr = NF90_INQ_VARID (nid, "ALBE", nvarid)
428      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
429         PRINT*, 'phyetat0: Le champ <ALBE> est absent'         PRINT *, 'phyetat0: Le champ <ALBE> est absent'
430         PRINT*, '          Mais je vais essayer de lire ALBE**'         PRINT *, '          Mais je vais essayer de lire ALBE**'
431         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
432            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
433               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
434               stop 1               stop 1
435            ENDIF            ENDIF
436            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
437            ierr = NF_INQ_VARID (nid, "ALBE"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "ALBE"//str2, nvarid)
438            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
439               PRINT*, "phyetat0: Le champ <ALBE"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <ALBE"//str2//"> est absent"
440               stop 1               stop 1
441            ENDIF            ENDIF
442            ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,nsrf))
443            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
444               PRINT*, "phyetat0: Lecture echouee pour <ALBE"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <ALBE"//str2//">"
445               stop 1               stop 1
446            ENDIF            ENDIF
447            xmin = 1.0E+20            xmin = 1.0E+20
# Line 590  contains Line 450  contains
450               xmin = MIN(albe(i,nsrf),xmin)               xmin = MIN(albe(i,nsrf),xmin)
451               xmax = MAX(albe(i,nsrf),xmax)               xmax = MAX(albe(i,nsrf),xmax)
452            ENDDO            ENDDO
453            PRINT*,'Albedo du sol ALBE**:', nsrf, xmin, xmax            PRINT *,'Albedo du sol ALBE**:', nsrf, xmin, xmax
454         ENDDO         ENDDO
455      ELSE      ELSE
456         PRINT*, 'phyetat0: Le champ <ALBE> est present'         PRINT *, 'phyetat0: Le champ <ALBE> est present'
457         PRINT*, '          J ignore donc les autres ALBE**'         PRINT *, '          J ignore donc les autres ALBE**'
458         ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,1))         call nf95_get_var(nid, nvarid, albe(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <ALBE>"  
           stop 1  
        ENDIF  
459         xmin = 1.0E+20         xmin = 1.0E+20
460         xmax = -1.0E+20         xmax = -1.0E+20
461         DO i = 1, klon         DO i = 1, klon
462            xmin = MIN(albe(i,1),xmin)            xmin = MIN(albe(i,1),xmin)
463            xmax = MAX(albe(i,1),xmax)            xmax = MAX(albe(i,1),xmax)
464         ENDDO         ENDDO
465         PRINT*,'Neige du sol <ALBE>', xmin, xmax         PRINT *,'Neige du sol <ALBE>', xmin, xmax
466         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
467            DO i = 1, klon            DO i = 1, klon
468               albe(i,nsrf) = albe(i,1)               albe(i,nsrf) = albe(i,1)
# Line 617  contains Line 473  contains
473    
474      ! Lecture de albedo au sol LW:      ! Lecture de albedo au sol LW:
475    
476      ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)      ierr = NF90_INQ_VARID (nid, "ALBLW", nvarid)
477      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
478         PRINT*, 'phyetat0: Le champ <ALBLW> est absent'         PRINT *, 'phyetat0: Le champ <ALBLW> est absent'
479         !        PRINT*, '          Mais je vais essayer de lire ALBLW**'         !        PRINT *, '          Mais je vais essayer de lire ALBLW**'
480         PRINT*, '          Mais je vais prendre ALBE**'         PRINT *, '          Mais je vais prendre ALBE**'
481         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
482            DO i = 1, klon            DO i = 1, klon
483               alblw(i,nsrf) = albe(i,nsrf)               alblw(i,nsrf) = albe(i,nsrf)
484            ENDDO            ENDDO
485         ENDDO         ENDDO
486      ELSE      ELSE
487         PRINT*, 'phyetat0: Le champ <ALBLW> est present'         PRINT *, 'phyetat0: Le champ <ALBLW> est present'
488         PRINT*, '          J ignore donc les autres ALBLW**'         PRINT *, '          J ignore donc les autres ALBLW**'
489         ierr = NF_GET_VAR_REAL(nid, nvarid, alblw(1,1))         call nf95_get_var(nid, nvarid, alblw(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <ALBLW>"  
           stop 1  
        ENDIF  
490         xmin = 1.0E+20         xmin = 1.0E+20
491         xmax = -1.0E+20         xmax = -1.0E+20
492         DO i = 1, klon         DO i = 1, klon
493            xmin = MIN(alblw(i,1),xmin)            xmin = MIN(alblw(i,1),xmin)
494            xmax = MAX(alblw(i,1),xmax)            xmax = MAX(alblw(i,1),xmax)
495         ENDDO         ENDDO
496         PRINT*,'Neige du sol <ALBLW>', xmin, xmax         PRINT *,'Neige du sol <ALBLW>', xmin, xmax
497         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
498            DO i = 1, klon            DO i = 1, klon
499               alblw(i,nsrf) = alblw(i,1)               alblw(i,nsrf) = alblw(i,1)
# Line 651  contains Line 503  contains
503    
504      ! Lecture de evaporation:        ! Lecture de evaporation:  
505    
506      ierr = NF_INQ_VARID (nid, "EVAP", nvarid)      ierr = NF90_INQ_VARID (nid, "EVAP", nvarid)
507      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
508         PRINT*, 'phyetat0: Le champ <EVAP> est absent'         PRINT *, 'phyetat0: Le champ <EVAP> est absent'
509         PRINT*, '          Mais je vais essayer de lire EVAP**'         PRINT *, '          Mais je vais essayer de lire EVAP**'
510         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
511            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
512               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
513               stop 1               stop 1
514            ENDIF            ENDIF
515            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
516            ierr = NF_INQ_VARID (nid, "EVAP"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "EVAP"//str2, nvarid)
517            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
518               PRINT*, "phyetat0: Le champ <EVAP"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <EVAP"//str2//"> est absent"
519               stop 1               stop 1
520            ENDIF            ENDIF
521            ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf))
522            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
523               PRINT*, "phyetat0: Lecture echouee pour <EVAP"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <EVAP"//str2//">"
524               stop 1               stop 1
525            ENDIF            ENDIF
526            xmin = 1.0E+20            xmin = 1.0E+20
# Line 677  contains Line 529  contains
529               xmin = MIN(evap(i,nsrf),xmin)               xmin = MIN(evap(i,nsrf),xmin)
530               xmax = MAX(evap(i,nsrf),xmax)               xmax = MAX(evap(i,nsrf),xmax)
531            ENDDO            ENDDO
532            PRINT*,'evap du sol EVAP**:', nsrf, xmin, xmax            PRINT *,'evap du sol EVAP**:', nsrf, xmin, xmax
533         ENDDO         ENDDO
534      ELSE      ELSE
535         PRINT*, 'phyetat0: Le champ <EVAP> est present'         PRINT *, 'phyetat0: Le champ <EVAP> est present'
536         PRINT*, '          J ignore donc les autres EVAP**'         PRINT *, '          J ignore donc les autres EVAP**'
537         ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,1))         call nf95_get_var(nid, nvarid, evap(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <EVAP>"  
           stop 1  
        ENDIF  
538         xmin = 1.0E+20         xmin = 1.0E+20
539         xmax = -1.0E+20         xmax = -1.0E+20
540         DO i = 1, klon         DO i = 1, klon
541            xmin = MIN(evap(i,1),xmin)            xmin = MIN(evap(i,1),xmin)
542            xmax = MAX(evap(i,1),xmax)            xmax = MAX(evap(i,1),xmax)
543         ENDDO         ENDDO
544         PRINT*,'Evap du sol <EVAP>', xmin, xmax         PRINT *,'Evap du sol <EVAP>', xmin, xmax
545         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
546            DO i = 1, klon            DO i = 1, klon
547               evap(i,nsrf) = evap(i,1)               evap(i,nsrf) = evap(i,1)
# Line 703  contains Line 551  contains
551    
552      ! Lecture precipitation liquide:      ! Lecture precipitation liquide:
553    
554      ierr = NF_INQ_VARID (nid, "rain_f", nvarid)      ierr = NF90_INQ_VARID (nid, "rain_f", nvarid)
555      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
556         PRINT*, 'phyetat0: Le champ <rain_f> est absent'         PRINT *, 'phyetat0: Le champ <rain_f> est absent'
557         stop 1         stop 1
558      ENDIF      ENDIF
559      ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall)      ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall)
560      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
561         PRINT*, 'phyetat0: Lecture echouee pour <rain_f>'         PRINT *, 'phyetat0: Lecture echouee pour <rain_f>'
562         stop 1         stop 1
563      ENDIF      ENDIF
564      xmin = 1.0E+20      xmin = 1.0E+20
# Line 719  contains Line 567  contains
567         xmin = MIN(rain_fall(i),xmin)         xmin = MIN(rain_fall(i),xmin)
568         xmax = MAX(rain_fall(i),xmax)         xmax = MAX(rain_fall(i),xmax)
569      ENDDO      ENDDO
570      PRINT*,'Precipitation liquide rain_f:', xmin, xmax      PRINT *,'Precipitation liquide rain_f:', xmin, xmax
571    
572      ! Lecture precipitation solide:      ! Lecture precipitation solide:
573    
574      ierr = NF_INQ_VARID (nid, "snow_f", nvarid)      ierr = NF90_INQ_VARID (nid, "snow_f", nvarid)
575      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
576         PRINT*, 'phyetat0: Le champ <snow_f> est absent'         PRINT *, 'phyetat0: Le champ <snow_f> est absent'
577         stop 1         stop 1
578      ENDIF      ENDIF
579      ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall)      ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall)
580      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
581         PRINT*, 'phyetat0: Lecture echouee pour <snow_f>'         PRINT *, 'phyetat0: Lecture echouee pour <snow_f>'
582         stop 1         stop 1
583      ENDIF      ENDIF
584      xmin = 1.0E+20      xmin = 1.0E+20
# Line 739  contains Line 587  contains
587         xmin = MIN(snow_fall(i),xmin)         xmin = MIN(snow_fall(i),xmin)
588         xmax = MAX(snow_fall(i),xmax)         xmax = MAX(snow_fall(i),xmax)
589      ENDDO      ENDDO
590      PRINT*,'Precipitation solide snow_f:', xmin, xmax      PRINT *,'Precipitation solide snow_f:', xmin, xmax
591    
592      ! Lecture rayonnement solaire au sol:      ! Lecture rayonnement solaire au sol:
593    
594      ierr = NF_INQ_VARID (nid, "solsw", nvarid)      ierr = NF90_INQ_VARID (nid, "solsw", nvarid)
595      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
596         PRINT*, 'phyetat0: Le champ <solsw> est absent'         PRINT *, 'phyetat0: Le champ <solsw> est absent'
597         PRINT*, 'mis a zero'         PRINT *, 'mis a zero'
598         solsw = 0.         solsw = 0.
599      ELSE      ELSE
600         ierr = NF_GET_VAR_REAL(nid, nvarid, solsw)         call nf95_get_var(nid, nvarid, solsw)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <solsw>'  
           stop 1  
        ENDIF  
601      ENDIF      ENDIF
602      xmin = 1.0E+20      xmin = 1.0E+20
603      xmax = -1.0E+20      xmax = -1.0E+20
# Line 761  contains Line 605  contains
605         xmin = MIN(solsw(i),xmin)         xmin = MIN(solsw(i),xmin)
606         xmax = MAX(solsw(i),xmax)         xmax = MAX(solsw(i),xmax)
607      ENDDO      ENDDO
608      PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax      PRINT *,'Rayonnement solaire au sol solsw:', xmin, xmax
609    
610      ! Lecture rayonnement IF au sol:      ! Lecture rayonnement IF au sol:
611    
612      ierr = NF_INQ_VARID (nid, "sollw", nvarid)      ierr = NF90_INQ_VARID (nid, "sollw", nvarid)
613      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
614         PRINT*, 'phyetat0: Le champ <sollw> est absent'         PRINT *, 'phyetat0: Le champ <sollw> est absent'
615         PRINT*, 'mis a zero'         PRINT *, 'mis a zero'
616         sollw = 0.         sollw = 0.
617      ELSE      ELSE
618         ierr = NF_GET_VAR_REAL(nid, nvarid, sollw)         call nf95_get_var(nid, nvarid, sollw)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <sollw>'  
           stop 1  
        ENDIF  
619      ENDIF      ENDIF
620      xmin = 1.0E+20      xmin = 1.0E+20
621      xmax = -1.0E+20      xmax = -1.0E+20
# Line 783  contains Line 623  contains
623         xmin = MIN(sollw(i),xmin)         xmin = MIN(sollw(i),xmin)
624         xmax = MAX(sollw(i),xmax)         xmax = MAX(sollw(i),xmax)
625      ENDDO      ENDDO
626      PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax      PRINT *,'Rayonnement IF au sol sollw:', xmin, xmax
627    
628    
629      ! Lecture derive des flux:      ! Lecture derive des flux:
630    
631      ierr = NF_INQ_VARID (nid, "fder", nvarid)      ierr = NF90_INQ_VARID (nid, "fder", nvarid)
632      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
633         PRINT*, 'phyetat0: Le champ <fder> est absent'         PRINT *, 'phyetat0: Le champ <fder> est absent'
634         PRINT*, 'mis a zero'         PRINT *, 'mis a zero'
635         fder = 0.         fder = 0.
636      ELSE      ELSE
637         ierr = NF_GET_VAR_REAL(nid, nvarid, fder)         call nf95_get_var(nid, nvarid, fder)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, 'phyetat0: Lecture echouee pour <fder>'  
           stop 1  
        ENDIF  
638      ENDIF      ENDIF
639      xmin = 1.0E+20      xmin = 1.0E+20
640      xmax = -1.0E+20      xmax = -1.0E+20
# Line 806  contains Line 642  contains
642         xmin = MIN(fder(i),xmin)         xmin = MIN(fder(i),xmin)
643         xmax = MAX(fder(i),xmax)         xmax = MAX(fder(i),xmax)
644      ENDDO      ENDDO
645      PRINT*,'Derive des flux fder:', xmin, xmax      PRINT *,'Derive des flux fder:', xmin, xmax
646    
647    
648      ! Lecture du rayonnement net au sol:      ! Lecture du rayonnement net au sol:
649    
650      ierr = NF_INQ_VARID (nid, "RADS", nvarid)      ierr = NF90_INQ_VARID (nid, "RADS", nvarid)
651      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
652         PRINT*, 'phyetat0: Le champ <RADS> est absent'         PRINT *, 'phyetat0: Le champ <RADS> est absent'
653         stop 1         stop 1
654      ENDIF      ENDIF
655      ierr = NF_GET_VAR_REAL(nid, nvarid, radsol)      ierr = NF_GET_VAR_REAL(nid, nvarid, radsol)
656      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
657         PRINT*, 'phyetat0: Lecture echouee pour <RADS>'         PRINT *, 'phyetat0: Lecture echouee pour <RADS>'
658         stop 1         stop 1
659      ENDIF      ENDIF
660      xmin = 1.0E+20      xmin = 1.0E+20
# Line 827  contains Line 663  contains
663         xmin = MIN(radsol(i),xmin)         xmin = MIN(radsol(i),xmin)
664         xmax = MAX(radsol(i),xmax)         xmax = MAX(radsol(i),xmax)
665      ENDDO      ENDDO
666      PRINT*,'Rayonnement net au sol radsol:', xmin, xmax      PRINT *,'Rayonnement net au sol radsol:', xmin, xmax
667    
668      ! Lecture de la longueur de rugosite      ! Lecture de la longueur de rugosite
669    
670    
671      ierr = NF_INQ_VARID (nid, "RUG", nvarid)      ierr = NF90_INQ_VARID (nid, "RUG", nvarid)
672      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
673         PRINT*, 'phyetat0: Le champ <RUG> est absent'         PRINT *, 'phyetat0: Le champ <RUG> est absent'
674         PRINT*, '          Mais je vais essayer de lire RUG**'         PRINT *, '          Mais je vais essayer de lire RUG**'
675         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
676            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
677               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
678               stop 1               stop 1
679            ENDIF            ENDIF
680            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
681            ierr = NF_INQ_VARID (nid, "RUG"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "RUG"//str2, nvarid)
682            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
683               PRINT*, "phyetat0: Le champ <RUG"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <RUG"//str2//"> est absent"
684               stop 1               stop 1
685            ENDIF            ENDIF
686            ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf))
687            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
688               PRINT*, "phyetat0: Lecture echouee pour <RUG"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <RUG"//str2//">"
689               stop 1               stop 1
690            ENDIF            ENDIF
691            xmin = 1.0E+20            xmin = 1.0E+20
# Line 858  contains Line 694  contains
694               xmin = MIN(frugs(i,nsrf),xmin)               xmin = MIN(frugs(i,nsrf),xmin)
695               xmax = MAX(frugs(i,nsrf),xmax)               xmax = MAX(frugs(i,nsrf),xmax)
696            ENDDO            ENDDO
697            PRINT*,'rugosite du sol RUG**:', nsrf, xmin, xmax            PRINT *,'rugosite du sol RUG**:', nsrf, xmin, xmax
698         ENDDO         ENDDO
699      ELSE      ELSE
700         PRINT*, 'phyetat0: Le champ <RUG> est present'         PRINT *, 'phyetat0: Le champ <RUG> est present'
701         PRINT*, '          J ignore donc les autres RUG**'         PRINT *, '          J ignore donc les autres RUG**'
702         ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,1))         call nf95_get_var(nid, nvarid, frugs(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <RUG>"  
           stop 1  
        ENDIF  
703         xmin = 1.0E+20         xmin = 1.0E+20
704         xmax = -1.0E+20         xmax = -1.0E+20
705         DO i = 1, klon         DO i = 1, klon
706            xmin = MIN(frugs(i,1),xmin)            xmin = MIN(frugs(i,1),xmin)
707            xmax = MAX(frugs(i,1),xmax)            xmax = MAX(frugs(i,1),xmax)
708         ENDDO         ENDDO
709         PRINT*,'rugosite <RUG>', xmin, xmax         PRINT *,'rugosite <RUG>', xmin, xmax
710         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
711            DO i = 1, klon            DO i = 1, klon
712               frugs(i,nsrf) = frugs(i,1)               frugs(i,nsrf) = frugs(i,1)
# Line 885  contains Line 717  contains
717    
718      ! Lecture de l'age de la neige:      ! Lecture de l'age de la neige:
719    
720      ierr = NF_INQ_VARID (nid, "AGESNO", nvarid)      ierr = NF90_INQ_VARID (nid, "AGESNO", nvarid)
721      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
722         PRINT*, 'phyetat0: Le champ <AGESNO> est absent'         PRINT *, 'phyetat0: Le champ <AGESNO> est absent'
723         PRINT*, '          Mais je vais essayer de lire AGESNO**'         PRINT *, '          Mais je vais essayer de lire AGESNO**'
724         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
725            IF (nsrf.GT.99) THEN            IF (nsrf.GT.99) THEN
726               PRINT*, "Trop de sous-mailles"               PRINT *, "Trop de sous-mailles"
727               stop 1               stop 1
728            ENDIF            ENDIF
729            WRITE(str2,'(i2.2)') nsrf            WRITE(str2,'(i2.2)') nsrf
730            ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid)            ierr = NF90_INQ_VARID (nid, "AGESNO"//str2, nvarid)
731            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
732               PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent"               PRINT *, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
733               agesno = 50.0               agesno = 50.0
734            ENDIF            ENDIF
735            ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))            ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))
736            IF (ierr.NE.NF_NOERR) THEN            IF (ierr.NE.NF90_NOERR) THEN
737               PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"               PRINT *, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
738               stop 1               stop 1
739            ENDIF            ENDIF
740            xmin = 1.0E+20            xmin = 1.0E+20
# Line 911  contains Line 743  contains
743               xmin = MIN(agesno(i,nsrf),xmin)               xmin = MIN(agesno(i,nsrf),xmin)
744               xmax = MAX(agesno(i,nsrf),xmax)               xmax = MAX(agesno(i,nsrf),xmax)
745            ENDDO            ENDDO
746            PRINT*,'Age de la neige AGESNO**:', nsrf, xmin, xmax            PRINT *,'Age de la neige AGESNO**:', nsrf, xmin, xmax
747         ENDDO         ENDDO
748      ELSE      ELSE
749         PRINT*, 'phyetat0: Le champ <AGESNO> est present'         PRINT *, 'phyetat0: Le champ <AGESNO> est present'
750         PRINT*, '          J ignore donc les autres AGESNO**'         PRINT *, '          J ignore donc les autres AGESNO**'
751         ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1))         call nf95_get_var(nid, nvarid, agesno(:,1))
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <AGESNO>"  
           stop 1  
        ENDIF  
752         xmin = 1.0E+20         xmin = 1.0E+20
753         xmax = -1.0E+20         xmax = -1.0E+20
754         DO i = 1, klon         DO i = 1, klon
755            xmin = MIN(agesno(i,1),xmin)            xmin = MIN(agesno(i,1),xmin)
756            xmax = MAX(agesno(i,1),xmax)            xmax = MAX(agesno(i,1),xmax)
757         ENDDO         ENDDO
758         PRINT*,'Age de la neige <AGESNO>', xmin, xmax         PRINT *,'Age de la neige <AGESNO>', xmin, xmax
759         DO nsrf = 2, nbsrf         DO nsrf = 2, nbsrf
760            DO i = 1, klon            DO i = 1, klon
761               agesno(i,nsrf) = agesno(i,1)               agesno(i,nsrf) = agesno(i,1)
# Line 936  contains Line 764  contains
764      ENDIF      ENDIF
765    
766    
767      ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)      ierr = NF90_INQ_VARID (nid, "ZMEA", nvarid)
768      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
769         PRINT*, 'phyetat0: Le champ <ZMEA> est absent'         PRINT *, 'phyetat0: Le champ <ZMEA> est absent'
770         stop 1         stop 1
771      ENDIF      ENDIF
772      ierr = NF_GET_VAR_REAL(nid, nvarid, zmea)      ierr = NF_GET_VAR_REAL(nid, nvarid, zmea)
773      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
774         PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>'         PRINT *, 'phyetat0: Lecture echouee pour <ZMEA>'
775         stop 1         stop 1
776      ENDIF      ENDIF
777      xmin = 1.0E+20      xmin = 1.0E+20
# Line 952  contains Line 780  contains
780         xmin = MIN(zmea(i),xmin)         xmin = MIN(zmea(i),xmin)
781         xmax = MAX(zmea(i),xmax)         xmax = MAX(zmea(i),xmax)
782      ENDDO      ENDDO
783      PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
784    
785    
786      ierr = NF_INQ_VARID (nid, "ZSTD", nvarid)      ierr = NF90_INQ_VARID (nid, "ZSTD", nvarid)
787      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
788         PRINT*, 'phyetat0: Le champ <ZSTD> est absent'         PRINT *, 'phyetat0: Le champ <ZSTD> est absent'
789         stop 1         stop 1
790      ENDIF      ENDIF
791      ierr = NF_GET_VAR_REAL(nid, nvarid, zstd)      ierr = NF_GET_VAR_REAL(nid, nvarid, zstd)
792      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
793         PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>'         PRINT *, 'phyetat0: Lecture echouee pour <ZSTD>'
794         stop 1         stop 1
795      ENDIF      ENDIF
796      xmin = 1.0E+20      xmin = 1.0E+20
# Line 971  contains Line 799  contains
799         xmin = MIN(zstd(i),xmin)         xmin = MIN(zstd(i),xmin)
800         xmax = MAX(zstd(i),xmax)         xmax = MAX(zstd(i),xmax)
801      ENDDO      ENDDO
802      PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
803    
804    
805      ierr = NF_INQ_VARID (nid, "ZSIG", nvarid)      ierr = NF90_INQ_VARID (nid, "ZSIG", nvarid)
806      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
807         PRINT*, 'phyetat0: Le champ <ZSIG> est absent'         PRINT *, 'phyetat0: Le champ <ZSIG> est absent'
808         stop 1         stop 1
809      ENDIF      ENDIF
810      ierr = NF_GET_VAR_REAL(nid, nvarid, zsig)      ierr = NF_GET_VAR_REAL(nid, nvarid, zsig)
811      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
812         PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>'         PRINT *, 'phyetat0: Lecture echouee pour <ZSIG>'
813         stop 1         stop 1
814      ENDIF      ENDIF
815      xmin = 1.0E+20      xmin = 1.0E+20
# Line 990  contains Line 818  contains
818         xmin = MIN(zsig(i),xmin)         xmin = MIN(zsig(i),xmin)
819         xmax = MAX(zsig(i),xmax)         xmax = MAX(zsig(i),xmax)
820      ENDDO      ENDDO
821      PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
822    
823    
824      ierr = NF_INQ_VARID (nid, "ZGAM", nvarid)      ierr = NF90_INQ_VARID (nid, "ZGAM", nvarid)
825      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
826         PRINT*, 'phyetat0: Le champ <ZGAM> est absent'         PRINT *, 'phyetat0: Le champ <ZGAM> est absent'
827         stop 1         stop 1
828      ENDIF      ENDIF
829      ierr = NF_GET_VAR_REAL(nid, nvarid, zgam)      ierr = NF_GET_VAR_REAL(nid, nvarid, zgam)
830      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
831         PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>'         PRINT *, 'phyetat0: Lecture echouee pour <ZGAM>'
832         stop 1         stop 1
833      ENDIF      ENDIF
834      xmin = 1.0E+20      xmin = 1.0E+20
# Line 1009  contains Line 837  contains
837         xmin = MIN(zgam(i),xmin)         xmin = MIN(zgam(i),xmin)
838         xmax = MAX(zgam(i),xmax)         xmax = MAX(zgam(i),xmax)
839      ENDDO      ENDDO
840      PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
841    
842    
843      ierr = NF_INQ_VARID (nid, "ZTHE", nvarid)      ierr = NF90_INQ_VARID (nid, "ZTHE", nvarid)
844      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
845         PRINT*, 'phyetat0: Le champ <ZTHE> est absent'         PRINT *, 'phyetat0: Le champ <ZTHE> est absent'
846         stop 1         stop 1
847      ENDIF      ENDIF
848      ierr = NF_GET_VAR_REAL(nid, nvarid, zthe)      ierr = NF_GET_VAR_REAL(nid, nvarid, zthe)
849      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
850         PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>'         PRINT *, 'phyetat0: Lecture echouee pour <ZTHE>'
851         stop 1         stop 1
852      ENDIF      ENDIF
853      xmin = 1.0E+20      xmin = 1.0E+20
# Line 1028  contains Line 856  contains
856         xmin = MIN(zthe(i),xmin)         xmin = MIN(zthe(i),xmin)
857         xmax = MAX(zthe(i),xmax)         xmax = MAX(zthe(i),xmax)
858      ENDDO      ENDDO
859      PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
860    
861    
862      ierr = NF_INQ_VARID (nid, "ZPIC", nvarid)      ierr = NF90_INQ_VARID (nid, "ZPIC", nvarid)
863      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
864         PRINT*, 'phyetat0: Le champ <ZPIC> est absent'         PRINT *, 'phyetat0: Le champ <ZPIC> est absent'
865         stop 1         stop 1
866      ENDIF      ENDIF
867      ierr = NF_GET_VAR_REAL(nid, nvarid, zpic)      ierr = NF_GET_VAR_REAL(nid, nvarid, zpic)
868      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
869         PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>'         PRINT *, 'phyetat0: Lecture echouee pour <ZPIC>'
870         stop 1         stop 1
871      ENDIF      ENDIF
872      xmin = 1.0E+20      xmin = 1.0E+20
# Line 1047  contains Line 875  contains
875         xmin = MIN(zpic(i),xmin)         xmin = MIN(zpic(i),xmin)
876         xmax = MAX(zpic(i),xmax)         xmax = MAX(zpic(i),xmax)
877      ENDDO      ENDDO
878      PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
879    
880      ierr = NF_INQ_VARID (nid, "ZVAL", nvarid)      ierr = NF90_INQ_VARID (nid, "ZVAL", nvarid)
881      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
882         PRINT*, 'phyetat0: Le champ <ZVAL> est absent'         PRINT *, 'phyetat0: Le champ <ZVAL> est absent'
883         stop 1         stop 1
884      ENDIF      ENDIF
885      ierr = NF_GET_VAR_REAL(nid, nvarid, zval)      ierr = NF_GET_VAR_REAL(nid, nvarid, zval)
886      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
887         PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>'         PRINT *, 'phyetat0: Lecture echouee pour <ZVAL>'
888         stop 1         stop 1
889      ENDIF      ENDIF
890      xmin = 1.0E+20      xmin = 1.0E+20
# Line 1065  contains Line 893  contains
893         xmin = MIN(zval(i),xmin)         xmin = MIN(zval(i),xmin)
894         xmax = MAX(zval(i),xmax)         xmax = MAX(zval(i),xmax)
895      ENDDO      ENDDO
896      PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax      PRINT *,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
   
   
     ierr = NF_INQ_VARID (nid, "RUGSREL", nvarid)  
     IF (ierr.NE.NF_NOERR) THEN  
        PRINT*, 'phyetat0: Le champ <RUGSREL> est absent'  
        stop 1  
     ENDIF  
     ierr = NF_GET_VAR_REAL(nid, nvarid, rugsrel)  
     IF (ierr.NE.NF_NOERR) THEN  
        PRINT*, 'phyetat0: Lecture echouee pour <RUGSREL>'  
        stop 1  
     ENDIF  
     xmin = 1.0E+20  
     xmax = -1.0E+20  
     DO i = 1, klon  
        xmin = MIN(rugsrel(i),xmin)  
        xmax = MAX(rugsrel(i),xmax)  
     ENDDO  
     PRINT*,'Rugosite relief (ecart-type) rugsrel:', xmin, xmax  
   
897    
898      ancien_ok = .TRUE.      ancien_ok = .TRUE.
899    
900      ierr = NF_INQ_VARID (nid, "TANCIEN", nvarid)      ierr = NF90_INQ_VARID (nid, "TANCIEN", nvarid)
901      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
902         PRINT*, "phyetat0: Le champ <TANCIEN> est absent"         PRINT *, "phyetat0: Le champ <TANCIEN> est absent"
903         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
904         ancien_ok = .FALSE.         ancien_ok = .FALSE.
905      ELSE      ELSE
906         ierr = NF_GET_VAR_REAL(nid, nvarid, t_ancien)         call nf95_get_var(nid, nvarid, t_ancien)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <TANCIEN>"  
           stop 1  
        ENDIF  
907      ENDIF      ENDIF
908    
909      ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)      ierr = NF90_INQ_VARID (nid, "QANCIEN", nvarid)
910      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
911         PRINT*, "phyetat0: Le champ <QANCIEN> est absent"         PRINT *, "phyetat0: Le champ <QANCIEN> est absent"
912         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
913         ancien_ok = .FALSE.         ancien_ok = .FALSE.
914      ELSE      ELSE
915         ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)         call nf95_get_var(nid, nvarid, q_ancien)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"  
           stop 1  
        ENDIF  
916      ENDIF      ENDIF
917    
918      ierr = NF_INQ_VARID (nid, "CLWCON", nvarid)      ierr = NF90_INQ_VARID (nid, "CLWCON", nvarid)
919      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
920         PRINT*, "phyetat0: Le champ CLWCON est absent"         PRINT *, "phyetat0: Le champ CLWCON est absent"
921         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
922         clwcon = 0.         clwcon = 0.
923      ELSE      ELSE
924         ierr = NF_GET_VAR_REAL(nid, nvarid, clwcon)         call nf95_get_var(nid, nvarid, clwcon)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <CLWCON>"  
           stop 1  
        ENDIF  
925      ENDIF      ENDIF
926      xmin = 1.0E+20      xmin = 1.0E+20
927      xmax = -1.0E+20      xmax = -1.0E+20
928      xmin = MINval(clwcon)      xmin = MINval(clwcon)
929      xmax = MAXval(clwcon)      xmax = MAXval(clwcon)
930      PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax      PRINT *,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
931    
932      ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid)      ierr = NF90_INQ_VARID (nid, "RNEBCON", nvarid)
933      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
934         PRINT*, "phyetat0: Le champ RNEBCON est absent"         PRINT *, "phyetat0: Le champ RNEBCON est absent"
935         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
936         rnebcon = 0.         rnebcon = 0.
937      ELSE      ELSE
938         ierr = NF_GET_VAR_REAL(nid, nvarid, rnebcon)         call nf95_get_var(nid, nvarid, rnebcon)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <RNEBCON>"  
           stop 1  
        ENDIF  
939      ENDIF      ENDIF
940      xmin = 1.0E+20      xmin = 1.0E+20
941      xmax = -1.0E+20      xmax = -1.0E+20
942      xmin = MINval(rnebcon)      xmin = MINval(rnebcon)
943      xmax = MAXval(rnebcon)      xmax = MAXval(rnebcon)
944      PRINT*,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax      PRINT *,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax
945    
946    
947      ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)      ierr = NF90_INQ_VARID (nid, "QANCIEN", nvarid)
948      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
949         PRINT*, "phyetat0: Le champ <QANCIEN> est absent"         PRINT *, "phyetat0: Le champ <QANCIEN> est absent"
950         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
951         ancien_ok = .FALSE.         ancien_ok = .FALSE.
952      ELSE      ELSE
953         ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)         call nf95_get_var(nid, nvarid, q_ancien)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"  
           stop 1  
        ENDIF  
954      ENDIF      ENDIF
955    
956      ! Lecture ratqs      ! Lecture ratqs
957    
958      ierr = NF_INQ_VARID (nid, "RATQS", nvarid)      ierr = NF90_INQ_VARID (nid, "RATQS", nvarid)
959      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
960         PRINT*, "phyetat0: Le champ <RATQS> est absent"         PRINT *, "phyetat0: Le champ <RATQS> est absent"
961         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
962         ratqs = 0.         ratqs = 0.
963      ELSE      ELSE
964         ierr = NF_GET_VAR_REAL(nid, nvarid, ratqs)         call nf95_get_var(nid, nvarid, ratqs)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <RATQS>"  
           stop 1  
        ENDIF  
965      ENDIF      ENDIF
966      xmin = 1.0E+20      xmin = 1.0E+20
967      xmax = -1.0E+20      xmax = -1.0E+20
968      xmin = MINval(ratqs)      xmin = MINval(ratqs)
969      xmax = MAXval(ratqs)      xmax = MAXval(ratqs)
970      PRINT*,'(ecart-type) ratqs:', xmin, xmax      PRINT *,'(ecart-type) ratqs:', xmin, xmax
971    
972      ! Lecture run_off_lic_0      ! Lecture run_off_lic_0
973    
974      ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid)      ierr = NF90_INQ_VARID (nid, "RUNOFFLIC0", nvarid)
975      IF (ierr.NE.NF_NOERR) THEN      IF (ierr.NE.NF90_NOERR) THEN
976         PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent"         PRINT *, "phyetat0: Le champ <RUNOFFLIC0> est absent"
977         PRINT*, "Depart legerement fausse. Mais je continue"         PRINT *, "Depart legerement fausse. Mais je continue"
978         run_off_lic_0 = 0.         run_off_lic_0 = 0.
979      ELSE      ELSE
980         ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0)         call nf95_get_var(nid, nvarid, run_off_lic_0)
        IF (ierr.NE.NF_NOERR) THEN  
           PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>"  
           stop 1  
        ENDIF  
981      ENDIF      ENDIF
982      xmin = 1.0E+20      xmin = 1.0E+20
983      xmax = -1.0E+20      xmax = -1.0E+20
984      xmin = MINval(run_off_lic_0)      xmin = MINval(run_off_lic_0)
985      xmax = MAXval(run_off_lic_0)      xmax = MAXval(run_off_lic_0)
986      PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax      PRINT *,'(ecart-type) run_off_lic_0:', xmin, xmax
987    
988      ! Fermer le fichier:      ! Fermer le fichier:
989    
Legend:
Removed from v.3  
changed lines
  Added in v.43
  ViewVC Help
Powered by ViewVC 1.1.21