trunk/dyn3d/conf_gcm.f

module conf_gcm_m

  IMPLICIT NONE

  INTEGER:: nday = 10 ! nombre de jours d'int\'egration
  integer:: day_step = 240 ! nombre de pas de temps de la dynamique par jour
  integer:: iperiod = 5 ! periode pour le pas Matsuno (en pas de temps)

  integer:: iapp_tracvl = 5
  ! Should normally be equal to "iperiod"
  ! frequence du groupement des flux (en pas de temps) 

  integer:: iconser = 240
  ! number of time steps between output of control variables

  integer:: iecri = 1 ! time interval between outputs to "dyn_hist.nc" (in days)

  integer:: iphysiq = 5
  ! number of time steps of dynamics between calls to physics

  logical:: raz_date = .false.
  ! prise en compte de la date initiale de la namelist et remise \`a
  ! z\'ero des compteurs de pas de temps (sinon on garde la date du
  ! fichier restart)

  integer:: periodav = 1 
  ! time interval between outputs in the dynamical part, in days

  logical:: offline = .FALSE.
  ! permet de mettre en route la nouvelle parametrisation de l'eau liquide

  integer:: prt_level = 0 ! niveau d'impression souhait\'e (0 = minimum)

  LOGICAL:: purmats= .FALSE.
  ! Help = Choix du schema d'integration temporel.
  ! y = pure Matsuno sinon c'est du Matsuno-leapfrog

  INTEGER:: iflag_phys = 1
  ! contr\^ole l'appel \`a la physique :
  ! 0 : pas de physique
  ! 1 : physique normale (appel \`a phylmd, phymars...) (default)
  ! 2 : rappel Newtonien pour la temp\'erature + friction au sol

contains

  SUBROUTINE conf_gcm

    ! Auteurs : L. Fairhead, P. Le Van
    ! Version du 29/04/97

    use abort_gcm_m, only: abort_gcm
    use comdissnew, only: read_comdissnew
    use unit_nml_m, only: unit_nml

    namelist /conf_gcm_nml/ raz_date, nday, day_step, iperiod, iapp_tracvl, &
         iconser, iecri, periodav, iphysiq, offline

    namelist /iniprint_nml/ prt_level

    namelist /logic_nml/ purmats, iflag_phys

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

    print *, "Call sequence information: conf_gcm"

    print *, "Enter namelist 'iniprint_nml'."
    read(unit=*, nml=iniprint_nml)
    write(unit_nml, nml=iniprint_nml)

    print *, "Enter namelist 'logic_nml'."
    read(unit=*, nml=logic_nml)
    write(unit_nml, nml=logic_nml)

    call read_comdissnew

    print *, "Enter namelist 'conf_gcm_nml'."
    read(unit=*, nml=conf_gcm_nml)
    write(unit_nml, nml=conf_gcm_nml)

    IF (MOD(day_step, iperiod) /= 0) call abort_gcm(modname = "conf_gcm", &
         message = &
         'Il faut choisir un nombre de pas par jour multiple de "iperiod".', &
         ierr = 1)

    IF (MOD(day_step, iphysiq)/= 0) call abort_gcm(modname = "conf_gcm", &
         message = &
         'Il faut choisir un nombre de pas par jour multiple de "iphysiq".', &
         ierr = 1)

  END SUBROUTINE conf_gcm

end module conf_gcm_m
1	module conf_gcm_m
2
3	IMPLICIT NONE
4
5	INTEGER:: nday = 10 ! nombre de jours d'int\'egration
6	integer:: day_step = 240 ! nombre de pas de temps de la dynamique par jour
7	integer:: iperiod = 5 ! periode pour le pas Matsuno (en pas de temps)
8
9	integer:: iapp_tracvl = 5
10	! Should normally be equal to "iperiod"
11	! frequence du groupement des flux (en pas de temps)
12
13	integer:: iconser = 240
14	! number of time steps between output of control variables
15
16	integer:: iecri = 1 ! time interval between outputs to "dyn_hist.nc" (in days)
17
18	integer:: iphysiq = 5
19	! number of time steps of dynamics between calls to physics
20
21	logical:: raz_date = .false.
22	! prise en compte de la date initiale de la namelist et remise \`a
23	! z\'ero des compteurs de pas de temps (sinon on garde la date du
24	! fichier restart)
25
26	integer:: periodav = 1
27	! time interval between outputs in the dynamical part, in days
28
29	logical:: offline = .FALSE.
30	! permet de mettre en route la nouvelle parametrisation de l'eau liquide
31
32	integer:: prt_level = 0 ! niveau d'impression souhait\'e (0 = minimum)
33
34	LOGICAL:: purmats= .FALSE.
35	! Help = Choix du schema d'integration temporel.
36	! y = pure Matsuno sinon c'est du Matsuno-leapfrog
37
38	INTEGER:: iflag_phys = 1
39	! contr\^ole l'appel \`a la physique :
40	! 0 : pas de physique
41	! 1 : physique normale (appel \`a phylmd, phymars...) (default)
42	! 2 : rappel Newtonien pour la temp\'erature + friction au sol
43
44	contains
45
46	SUBROUTINE conf_gcm
47
48	! Auteurs : L. Fairhead, P. Le Van
49	! Version du 29/04/97
50
51	use abort_gcm_m, only: abort_gcm
52	use comdissnew, only: read_comdissnew
53	use unit_nml_m, only: unit_nml
54
55	namelist /conf_gcm_nml/ raz_date, nday, day_step, iperiod, iapp_tracvl, &
56	iconser, iecri, periodav, iphysiq, offline
57
58	namelist /iniprint_nml/ prt_level
59
60	namelist /logic_nml/ purmats, iflag_phys
61
62	!------------------------------------
63
64	print *, "Call sequence information: conf_gcm"
65
66	print *, "Enter namelist 'iniprint_nml'."
67	read(unit=*, nml=iniprint_nml)
68	write(unit_nml, nml=iniprint_nml)
69
70	print *, "Enter namelist 'logic_nml'."
71	read(unit=*, nml=logic_nml)
72	write(unit_nml, nml=logic_nml)
73
74	call read_comdissnew
75
76	print *, "Enter namelist 'conf_gcm_nml'."
77	read(unit=*, nml=conf_gcm_nml)
78	write(unit_nml, nml=conf_gcm_nml)
79
80	IF (MOD(day_step, iperiod) /= 0) call abort_gcm(modname = "conf_gcm", &
81	message = &
82	'Il faut choisir un nombre de pas par jour multiple de "iperiod".', &
83	ierr = 1)
84
85	IF (MOD(day_step, iphysiq)/= 0) call abort_gcm(modname = "conf_gcm", &
86	message = &
87	'Il faut choisir un nombre de pas par jour multiple de "iphysiq".', &
88	ierr = 1)
89
90	END SUBROUTINE conf_gcm
91
92	end module conf_gcm_m