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 |
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" |
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 |
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 |
|
|