/[lmdze]/trunk/libf/phylmd/condsurf.f
ViewVC logotype

Diff of /trunk/libf/phylmd/condsurf.f

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

revision 60 by guez, Mon Jul 21 16:05:07 2008 UTC revision 61 by guez, Fri Apr 20 14:58:43 2012 UTC
# Line 6  c Line 6  c
6        use dimphy        use dimphy
7        use temps        use temps
8        use clesphys2, only: ok_limitvrai        use clesphys2, only: ok_limitvrai
9                use netcdf
10        IMPLICIT none        IMPLICIT none
11  c  c
12  c I. Musat 05.2005  c I. Musat 05.2005
# Line 19  c jourvrai : input  , vrai jour de la si Line 20  c jourvrai : input  , vrai jour de la si
20  c  c
21  c lmt_bils: bilan chaleur au sol (a utiliser pour "slab-ocean")  c lmt_bils: bilan chaleur au sol (a utiliser pour "slab-ocean")
22  c  c
       include "netcdf.inc"  
23        INTEGER nid, nvarid        INTEGER nid, nvarid
24        INTEGER debut(2)        INTEGER debut(2)
25        INTEGER epais(2)        INTEGER epais(2)
# Line 80  c           PRINT *,' Fichier  Limite ', Line 80  c           PRINT *,' Fichier  Limite ',
80  c  c
81  c Ouvrir le fichier en format NetCDF:  c Ouvrir le fichier en format NetCDF:
82  c  c
83        ierr = NF_OPEN (fich, NF_NOWRITE,nid)        ierr = NF90_OPEN (fich, NF90_NOWRITE,nid)
84        IF (ierr.NE.NF_NOERR) THEN        IF (ierr.NE.NF90_NOERR) THEN
85          WRITE(6,*)' Pb d''ouverture du fichier ', fich          WRITE(6,*)' Pb d''ouverture du fichier ', fich
86          WRITE(6,*)' Le fichier limit ',fich,' (avec 4 chiffres , pour'          WRITE(6,*)' Le fichier limit ',fich,' (avec 4 chiffres , pour'
87          WRITE(6,*)'       l an 2000 )  ,  n existe  pas !  '          WRITE(6,*)'       l an 2000 )  ,  n existe  pas !  '
# Line 98  c Line 98  c
98  c  c
99  c Bilan flux de chaleur au sol:  c Bilan flux de chaleur au sol:
100  c  c
101        ierr = NF_INQ_VARID (nid, "BILS", nvarid)        ierr = NF90_INQ_VARID (nid, "BILS", nvarid)
102        IF (ierr .NE. NF_NOERR) THEN        IF (ierr .NE. NF90_NOERR) THEN
103           PRINT*, "condsurf: Le champ <BILS> est absent"           PRINT*, "condsurf: Le champ <BILS> est absent"
104           stop 1           stop 1
105        ENDIF        ENDIF
106  c     PRINT*,'debut,epais',debut,epais  c     PRINT*,'debut,epais',debut,epais
107        ierr = NF_GET_VARA_REAL(nid, nvarid,debut,epais,lmt_bils)        ierr = NF90_GET_VAR(nid, nvarid,lmt_bils,debut,epais)
108        IF (ierr .NE. NF_NOERR) THEN        IF (ierr .NE. NF90_NOERR) THEN
109           PRINT*, "condsurf: Lecture echouee pour <BILS>"           PRINT*, "condsurf: Lecture echouee pour <BILS>"
110           stop 1           stop 1
111        ENDIF        ENDIF
# Line 113  c     ENDDO !k = 1, jour Line 113  c     ENDDO !k = 1, jour
113  c  c
114  c Fermer le fichier:  c Fermer le fichier:
115  c  c
116        ierr = NF_CLOSE(nid)        ierr = NF90_CLOSE(nid)
117  c  c
118  c  c
119  c     PRINT*, 'lmt_bils est lu pour jour: ', jour  c     PRINT*, 'lmt_bils est lu pour jour: ', jour
Legend:
Removed from v.60  
changed lines
  Added in v.61
  ViewVC Help
Powered by ViewVC 1.1.21