--- trunk/libf/dyn3d/dynredem1.f90 2008/08/07 12:29:13 18 +++ trunk/dyn3d/dynredem1.f90 2013/11/15 18:45:49 76 @@ -1,118 +1,83 @@ -SUBROUTINE dynredem1(fichnom, time, vcov, ucov, teta, q, masse, ps) +module dynredem1_m - ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30 + IMPLICIT NONE - ! Ecriture du fichier de redémarrage au format NetCDF +contains - USE dimens_m, ONLY : llm, nqmx - USE paramet_m, ONLY : ip1jm, ip1jmp1 - USE temps, ONLY : itaufin, itau_dyn - USE abort_gcm_m, ONLY : abort_gcm - USE iniadvtrac_m, ONLY : tname - use netcdf, only: nf90_open, nf90_write, nf90_noerr, nf90_inq_varid + SUBROUTINE dynredem1(fichnom, vcov, ucov, teta, q, masse, ps, itau) - IMPLICIT NONE + ! From dyn3d/dynredem.F, version 1.2 2004/06/22 11:45:30 + ! Ecriture du fichier de redémarrage au format NetCDF + + USE dimens_m, ONLY: iim, jjm, llm, nqmx + USE iniadvtrac_m, ONLY: tname + use netcdf, only: nf90_write + use netcdf95, only: nf95_close, nf95_inq_varid, nf95_open, nf95_put_var + use nr_util, only: assert + + CHARACTER(len=*), INTENT (IN):: fichnom + REAL, INTENT(IN):: vcov(:, :, :) ! (iim + 1, jjm, llm) + REAL, INTENT(IN):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm) + REAL, INTENT(IN):: teta(:, :, :) ! (iim + 1, jjm + 1, llm) + REAL, INTENT(IN):: q(:, :, :, :) ! (iim + 1, jjm + 1, llm, nqmx) + REAL, INTENT(IN):: masse(:, :, :) ! (iim + 1, jjm + 1, llm) + REAL, INTENT(IN):: ps(:, :) ! (iim + 1, jjm + 1) + INTEGER, INTENT(IN):: itau + + ! Variables local to the procedure: + INTEGER ncid, varid + INTEGER iq + INTEGER:: nb = 0 + + !--------------------------------------------------------- + + PRINT *, 'Call sequence information: dynredem1' + + call assert((/size(vcov, 1), size(ucov, 1), size(teta, 1), size(q, 1), & + size(masse, 1), size(ps, 1)/) == iim + 1, "dynredem1 iim") + call assert((/size(vcov, 2) + 1, size(ucov, 2), size(teta, 2), size(q, 2), & + size(masse, 2), size(ps, 2)/) == jjm + 1, "dynredem1 jjm") + call assert((/size(vcov, 3), size(ucov, 3), size(teta, 3), size(q, 3), & + size(masse, 3)/) == llm, "dynredem1 llm") + call assert(size(q, 4) == nqmx, "dynredem1 nqmx") + + call nf95_open(fichnom, nf90_write, ncid) + + ! Écriture/extension de la coordonnée temps + nb = nb + 1 + call nf95_inq_varid(ncid, 'temps', varid) + call nf95_put_var(ncid, varid, values=0., start=(/nb/)) + PRINT *, "Enregistrement pour nb = ", nb + + ! Récriture du tableau de contrôle, "itaufin" n'est pas défini quand + ! on passe dans "dynredem0" + call nf95_inq_varid(ncid, 'controle', varid) + call nf95_put_var(ncid, varid, real(itau), start=(/31/)) + + ! Écriture des champs + + call nf95_inq_varid(ncid, 'ucov', varid) + call nf95_put_var(ncid, varid, ucov) + + call nf95_inq_varid(ncid, 'vcov', varid) + call nf95_put_var(ncid, varid, vcov) + + call nf95_inq_varid(ncid, 'teta', varid) + call nf95_put_var(ncid, varid, teta) + + DO iq = 1, nqmx + call nf95_inq_varid(ncid, tname(iq), varid) + call nf95_put_var(ncid, varid, q(:, :, :, iq)) + END DO + + call nf95_inq_varid(ncid, 'masse', varid) + call nf95_put_var(ncid, varid, masse) + + call nf95_inq_varid(ncid, 'ps', varid) + call nf95_put_var(ncid, varid, ps) - REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm), teta(ip1jmp1, llm) - REAL, INTENT (IN) :: ps(ip1jmp1), masse(ip1jmp1, llm) - REAL, INTENT (IN) :: q(ip1jmp1, llm, nqmx) - CHARACTER(len=*), INTENT (IN) :: fichnom - - INCLUDE 'netcdf.inc' - - REAL :: time - INTEGER :: nid, nvarid - INTEGER :: ierr - INTEGER :: iq - INTEGER :: length - PARAMETER (length=100) - REAL :: tab_cntrl(length) ! tableau des parametres du run - CHARACTER (len=20) :: modname - CHARACTER (len=80) :: abort_message - INTEGER :: nb = 0 - - !--------------------------------------------------------- - - PRINT *, 'Call sequence information: dynredem1' - - modname = 'dynredem1' - ierr = nf90_open(fichnom, nf90_write, nid) - IF (ierr/=nf90_noerr) THEN - PRINT *, 'Pb. d ouverture ' // fichnom - STOP 1 - END IF - - ! Ecriture/extension de la coordonnee temps - - nb = nb + 1 - 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) - END IF - ierr = nf_put_var1_real(nid, nvarid, nb, time) - PRINT *, 'Enregistrement pour ', nb, time - - - ! Re-ecriture du tableau de controle, itaufin n'est plus defini quand - ! on passe dans dynredem0 - 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) - END IF - ierr = nf_get_var_real(nid, nvarid, tab_cntrl) - tab_cntrl(31) = real(itau_dyn+itaufin) - ierr = nf_put_var_real(nid, nvarid, tab_cntrl) - - ! Ecriture des champs - - 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 = 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 = 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 = 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 = 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 = nf90_inq_varid(nid, 'ps', nvarid) - IF (ierr/=nf90_noerr) THEN - PRINT *, 'Variable ps n est pas definie' - STOP 1 - END IF - ierr = nf_put_var_real(nid, nvarid, ps) + call nf95_close(ncid) - ierr = nf_close(nid) + END SUBROUTINE dynredem1 -END SUBROUTINE dynredem1 +end module dynredem1_m