--- trunk/libf/dyn3d/dynredem1.f90 2008/03/03 16:32:04 5 +++ trunk/libf/dyn3d/dynredem1.f90 2008/08/07 12:29:13 18 @@ -8,7 +8,8 @@ USE paramet_m, ONLY : ip1jm, ip1jmp1 USE temps, ONLY : itaufin, itau_dyn USE abort_gcm_m, ONLY : abort_gcm - USE advtrac_m, ONLY : tname + USE iniadvtrac_m, ONLY : tname + use netcdf, only: nf90_open, nf90_write, nf90_noerr, nf90_inq_varid IMPLICIT NONE @@ -35,8 +36,8 @@ PRINT *, 'Call sequence information: dynredem1' modname = 'dynredem1' - ierr = nf_open(fichnom, nf_write, nid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_open(fichnom, nf90_write, nid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Pb. d ouverture ' // fichnom STOP 1 END IF @@ -44,8 +45,8 @@ ! Ecriture/extension de la coordonnee temps nb = nb + 1 - ierr = nf_inq_varid(nid, 'temps', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'temps', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, nf_strerror(ierr) abort_message = 'Variable temps n est pas definie' CALL abort_gcm(modname, abort_message, ierr) @@ -56,8 +57,8 @@ ! Re-ecriture du tableau de controle, itaufin n'est plus defini quand ! on passe dans dynredem0 - ierr = nf_inq_varid(nid, 'controle', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'controle', nvarid) + IF (ierr/=nf90_noerr) THEN abort_message = 'dynredem1: Le champ est absent' ierr = 1 CALL abort_gcm(modname, abort_message, ierr) @@ -68,45 +69,45 @@ ! Ecriture des champs - ierr = nf_inq_varid(nid, 'ucov', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'ucov', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable ucov n est pas definie' STOP 1 END IF ierr = nf_put_var_real(nid, nvarid, ucov) - ierr = nf_inq_varid(nid, 'vcov', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'vcov', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable vcov n est pas definie' STOP 1 END IF ierr = nf_put_var_real(nid, nvarid, vcov) - ierr = nf_inq_varid(nid, 'teta', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'teta', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable teta n est pas definie' STOP 1 END IF ierr = nf_put_var_real(nid, nvarid, teta) DO iq = 1, nqmx - ierr = nf_inq_varid(nid, tname(iq), nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, tname(iq), nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable ', tname(iq), 'n''est pas définie' STOP 1 END IF ierr = nf_put_var_real(nid, nvarid, q(1, 1, iq)) END DO - ierr = nf_inq_varid(nid, 'masse', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'masse', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable masse n est pas definie' STOP 1 END IF ierr = nf_put_var_real(nid, nvarid, masse) - ierr = nf_inq_varid(nid, 'ps', nvarid) - IF (ierr/=nf_noerr) THEN + ierr = nf90_inq_varid(nid, 'ps', nvarid) + IF (ierr/=nf90_noerr) THEN PRINT *, 'Variable ps n est pas definie' STOP 1 END IF