/[lmdze]/trunk/dyn3d/wrgrads.f
ViewVC logotype

Diff of /trunk/dyn3d/wrgrads.f

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

revision 105 by guez, Thu Sep 4 10:40:24 2014 UTC revision 107 by guez, Thu Sep 11 15:09:15 2014 UTC
# Line 18  SUBROUTINE wrgrads(i_f, nl, field, name, Line 18  SUBROUTINE wrgrads(i_f, nl, field, name,
18    INTEGER, INTENT(IN):: i_f    INTEGER, INTENT(IN):: i_f
19    integer nl    integer nl
20    REAL, INTENT(IN):: field(imx*jmx*lmx)    REAL, INTENT(IN):: field(imx*jmx*lmx)
21    CHARACTER(len=10) name, file    CHARACTER(len=*) name, titlevar
22    CHARACTER(len=10) titlevar    CHARACTER(len=10) file
23    
24    ! local    ! local
25    
26    INTEGER im, jm, lm, i, j, l, lnblnk, iv, iii, iji, iif, ijf    INTEGER im, jm, lm, i, j, l, iv, iii, iji, iif, ijf
27    
28    LOGICAL writectl    LOGICAL writectl
29    
# Line 56  SUBROUTINE wrgrads(i_f, nl, field, name, Line 56  SUBROUTINE wrgrads(i_f, nl, field, name,
56        ivar(i_f) = ivar(i_f) + 1        ivar(i_f) = ivar(i_f) + 1
57        nvar(i_f) = ivar(i_f)        nvar(i_f) = ivar(i_f)
58        var(ivar(i_f), i_f) = name        var(ivar(i_f), i_f) = name
59        tvar(ivar(i_f), i_f) = titlevar(1:lnblnk(titlevar))        tvar(ivar(i_f), i_f) = trim(titlevar)
60        nld(ivar(i_f), i_f) = nl        nld(ivar(i_f), i_f) = nl
61        PRINT *, 'initialisation ecriture de ', var(ivar(i_f), i_f)        PRINT *, 'initialisation ecriture de ', var(ivar(i_f), i_f)
62        PRINT *, 'i_f ivar(i_f) nld ', i_f, ivar(i_f), nld(ivar(i_f), i_f)        PRINT *, 'i_f ivar(i_f) nld ', i_f, ivar(i_f), nld(ivar(i_f), i_f)
# Line 96  SUBROUTINE wrgrads(i_f, nl, field, name, Line 96  SUBROUTINE wrgrads(i_f, nl, field, name,
96    
97      file = fichier(i_f)      file = fichier(i_f)
98      ! WARNING! on reecrase le fichier .ctl a chaque ecriture      ! WARNING! on reecrase le fichier .ctl a chaque ecriture
99      OPEN (unit(i_f), FILE=file(1:lnblnk(file))//'.ctl', FORM='formatted', &      OPEN (unit(i_f), FILE=trim(file)//'.ctl', FORM='formatted', &
100        STATUS='unknown')        STATUS='unknown')
101      WRITE (unit(i_f), '(a5,1x,a40)') 'DSET ', '^' // file(1:lnblnk(file)) // &      WRITE (unit(i_f), '(a5,1x,a40)') 'DSET ', '^' // trim(file) // &
102        '.dat'        '.dat'
103    
104      WRITE (unit(i_f), '(a12)') 'UNDEF 1.0E30'      WRITE (unit(i_f), '(a12)') 'UNDEF 1.0E30'
Legend:
Removed from v.105  
changed lines
  Added in v.107
  ViewVC Help
Powered by ViewVC 1.1.21