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 |
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) |
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 ! ' |
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 |
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 |