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

Diff of /trunk/dyn3d/iniadvtrac.f

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

revision 130 by guez, Fri Feb 13 18:22:38 2015 UTC revision 131 by guez, Fri Feb 27 16:44:07 2015 UTC
# Line 30  module iniadvtrac_m Line 30  module iniadvtrac_m
30    
31    INTEGER hadv(nqmx) ! indice schéma transport horizontal    INTEGER hadv(nqmx) ! indice schéma transport horizontal
32    INTEGER vadv(nqmx) ! indice schéma transport vertical    INTEGER vadv(nqmx) ! indice schéma transport vertical
33    character(len=8) tnom(nqmx) ! nom court du traceur  
34    character(len=10) tname(nqmx) ! nom du traceur pour restart    character(len=10) tname(nqmx)
35      ! nom du traceur pour fichiers restart et historiques
36    
37    character(len=13) ttext(nqmx) ! nom long du traceur pour sorties    character(len=13) ttext(nqmx) ! nom long du traceur pour sorties
38    
39  contains  contains
# Line 79  contains Line 81  contains
81         call assert(nq_local == nqmx, "iniadvtrac nq_local")         call assert(nq_local == nqmx, "iniadvtrac nq_local")
82    
83         do iq=1, nqmx         do iq=1, nqmx
84            read(unit, fmt=*) hadv(iq), vadv(iq), tnom(iq)            read(unit, fmt=*) hadv(iq), vadv(iq), tname(iq)
85            if (.not. any(hadv(iq) == allowed_adv) &            if (.not. any(hadv(iq) == allowed_adv) &
86                 .or. .not. any(vadv(iq) == allowed_adv)) then                 .or. .not. any(vadv(iq) == allowed_adv)) then
87               print *, "bad number for advection scheme"               print *, "bad number for advection scheme"
# Line 93  contains Line 95  contains
95         call assert(nqmx == 4, "iniadvtrac nqmx")         call assert(nqmx == 4, "iniadvtrac nqmx")
96         hadv(:4) = (/14, 10, 10, 10/)         hadv(:4) = (/14, 10, 10, 10/)
97         vadv(:4) = hadv(:4)         vadv(:4) = hadv(:4)
98         tnom(1) = 'H2Ov'         tname(1) = 'H2Ov'
99         tnom(2) = 'H2Ol'         tname(2) = 'H2Ol'
100         tnom(3) = 'RN'         tname(3) = 'RN'
101         tnom(4) = 'PB'         tname(4) = 'PB'
102         do iq = 1, nqmx         do iq = 1, nqmx
103            print *, hadv(iq), vadv(iq), tnom(iq)            print *, hadv(iq), vadv(iq), tname(iq)
104         end do         end do
105      ENDIF      ENDIF
106    
     tname = tnom  
   
107      ! À partir du nom court du traceur et du schéma d'advection, on      ! À partir du nom court du traceur et du schéma d'advection, on
108      ! détermine le nom long :      ! détermine le nom long :
109      do iq = 1, nqmx      do iq = 1, nqmx
# Line 119  contains Line 119  contains
119         endif         endif
120    
121         IF (iadv(iq) == 0) THEN         IF (iadv(iq) == 0) THEN
122            ttext(iq) = tnom(iq)            ttext(iq) = tname(iq)
123         ELSE         ELSE
124            ttext(iq)=trim(tnom(iq)) // descrq(iadv(iq))            ttext(iq)=trim(tname(iq)) // descrq(iadv(iq))
125         endif         endif
126      end do      end do
127    
Legend:
Removed from v.130  
changed lines
  Added in v.131
  ViewVC Help
Powered by ViewVC 1.1.21