Sources/dyn3d/iniadvtrac.f

module iniadvtrac_m

  ! From advtrac.h, version 1.1.1.1 2004/05/19 12:53:06

  ! iq = 1 pour l'eau vapeur
  ! iq = 2 pour l'eau liquide
  ! et \'eventuellement iq = 3, ..., nqmx pour les autres traceurs

  use dimens_m, only: nqmx

  implicit none

  private nqmx

  INTEGER iadv(nqmx) ! indice du sch\'ema d'advection pour l'eau et les traceurs

  character(len=10) tname(nqmx)
  ! nom du traceur pour fichiers restart et historiques

  character(len=13) ttext(nqmx) ! nom long du traceur pour sorties

contains

  subroutine iniadvtrac

    ! From dyn3d/iniadvtrac.F, version 1.3 2005/04/13 08:58:34

    ! Authors: P. Le Van, L. Fairhead, F. Hourdin, F. Codron,
    ! F. Forget, M.-A. Filiberti

    ! Initialisation des traceurs
    ! Choix du sch\'ema pour l'advection dans le fichier "traceur.def"

    use nr_util, only: assert
    use jumble, only: new_unit

    ! Local:
    character(len=3) descrq(0:14)
    integer iq, iostat, nq_local, unit

    integer, parameter:: allowed_adv(5) = (/0, 10, 12, 13, 14/)
    ! Allowed values for iadv:
    ! 10: schema Van-leer (retenu pour l'eau vapeur et liquide)
    ! 12: schema Frederic Hourdin I
    ! 13: schema Frederic Hourdin II
    ! 14: schema Van-leer + humidite specifique 

    !-----------------------------------------------------------------------

    print *, "Call sequence information: iniadvtrac"

    ! Initializations:
    descrq(0)=''
    descrq(10)='VL1'
    descrq(12)='FH1'
    descrq(13)='FH2'
    descrq(14)='VLH'

    ! Choix du sch\'ema pour l'advection dans fichier "traceur.def"
    call new_unit(unit)
    open(unit, file='traceur.def', status='old', action="read", &
         position="rewind", iostat=iostat)
    if (iostat == 0) then
       print *, 'Ouverture de "traceur.def" ok'
       read(unit, fmt=*) nq_local
       print *, 'nombre de traceurs ', nq_local
       call assert(nq_local == nqmx, "iniadvtrac nq_local")

       do iq=1, nqmx
          read(unit, fmt=*) iadv(iq), tname(iq)
          if (.not. any(iadv(iq) == allowed_adv)) then
             print *, "bad number for advection scheme"
             stop 1
          end if
       end do
       close(unit) 
    else
       print *, 'Could not open "traceur.def".'
       print *, 'Using default values.'
       call assert(nqmx == 4, "iniadvtrac nqmx")
       iadv(:4) = (/14, 10, 10, 10/)
       tname(1) = 'H2Ov'
       tname(2) = 'H2Ol'
       tname(3) = 'RN'
       tname(4) = 'PB'
       do iq = 1, nqmx
          print *, iadv(iq), tname(iq)
       end do
    ENDIF

    ! \`A partir du nom court du traceur et du sch\'ema d'advection, on
    ! d\'etermine le nom long :
    do iq = 1, nqmx
       ttext(iq)=trim(tname(iq)) // descrq(iadv(iq))
    end do

  END subroutine iniadvtrac

end module iniadvtrac_m
1	module iniadvtrac_m
2
3	! From advtrac.h, version 1.1.1.1 2004/05/19 12:53:06
4
5	! iq = 1 pour l'eau vapeur
6	! iq = 2 pour l'eau liquide
7	! et \'eventuellement iq = 3, ..., nqmx pour les autres traceurs
8
9	use dimens_m, only: nqmx
10
11	implicit none
12
13	private nqmx
14
15	INTEGER iadv(nqmx) ! indice du sch\'ema d'advection pour l'eau et les traceurs
16
17	character(len=10) tname(nqmx)
18	! nom du traceur pour fichiers restart et historiques
19
20	character(len=13) ttext(nqmx) ! nom long du traceur pour sorties
21
22	contains
23
24	subroutine iniadvtrac
25
26	! From dyn3d/iniadvtrac.F, version 1.3 2005/04/13 08:58:34
27
28	! Authors: P. Le Van, L. Fairhead, F. Hourdin, F. Codron,
29	! F. Forget, M.-A. Filiberti
30
31	! Initialisation des traceurs
32	! Choix du sch\'ema pour l'advection dans le fichier "traceur.def"
33
34	use nr_util, only: assert
35	use jumble, only: new_unit
36
37	! Local:
38	character(len=3) descrq(0:14)
39	integer iq, iostat, nq_local, unit
40
41	integer, parameter:: allowed_adv(5) = (/0, 10, 12, 13, 14/)
42	! Allowed values for iadv:
43	! 10: schema Van-leer (retenu pour l'eau vapeur et liquide)
44	! 12: schema Frederic Hourdin I
45	! 13: schema Frederic Hourdin II
46	! 14: schema Van-leer + humidite specifique
47
48	!-----------------------------------------------------------------------
49
50	print *, "Call sequence information: iniadvtrac"
51
52	! Initializations:
53	descrq(0)=''
54	descrq(10)='VL1'
55	descrq(12)='FH1'
56	descrq(13)='FH2'
57	descrq(14)='VLH'
58
59	! Choix du sch\'ema pour l'advection dans fichier "traceur.def"
60	call new_unit(unit)
61	open(unit, file='traceur.def', status='old', action="read", &
62	position="rewind", iostat=iostat)
63	if (iostat == 0) then
64	print *, 'Ouverture de "traceur.def" ok'
65	read(unit, fmt=*) nq_local
66	print *, 'nombre de traceurs ', nq_local
67	call assert(nq_local == nqmx, "iniadvtrac nq_local")
68
69	do iq=1, nqmx
70	read(unit, fmt=*) iadv(iq), tname(iq)
71	if (.not. any(iadv(iq) == allowed_adv)) then
72	print *, "bad number for advection scheme"
73	stop 1
74	end if
75	end do
76	close(unit)
77	else
78	print *, 'Could not open "traceur.def".'
79	print *, 'Using default values.'
80	call assert(nqmx == 4, "iniadvtrac nqmx")
81	iadv(:4) = (/14, 10, 10, 10/)
82	tname(1) = 'H2Ov'
83	tname(2) = 'H2Ol'
84	tname(3) = 'RN'
85	tname(4) = 'PB'
86	do iq = 1, nqmx
87	print *, iadv(iq), tname(iq)
88	end do
89	ENDIF
90
91	! \`A partir du nom court du traceur et du sch\'ema d'advection, on
92	! d\'etermine le nom long :
93	do iq = 1, nqmx
94	ttext(iq)=trim(tname(iq)) // descrq(iadv(iq))
95	end do
96
97	END subroutine iniadvtrac
98
99	end module iniadvtrac_m