trunk/dyn3d/conf_gcm.f

module conf_gcm_m

  IMPLICIT NONE

  INTEGER:: nday = 1 ! 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. ! output variables for offline transport

  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

  logical:: iflag_phys = .true. ! call parameterizations of physics
  INTEGER, SAVE:: lmt_pas ! number of time steps of "physics" per day

contains

  SUBROUTINE conf_gcm

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

    use abort_gcm_m, only: abort_gcm
    use nr_util, only: assert
    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)

    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("conf_gcm", &
         'Il faut choisir un nombre de pas par jour multiple de "iperiod".')

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

    call assert(mod(iphysiq, iperiod) == 0, &
         "conf_gcm -- iphysiq must be multiple of iperiod")

    lmt_pas = day_step / iphysiq
    print *, 'Number of time steps of "physics" per day: ', lmt_pas

  END SUBROUTINE conf_gcm

end module conf_gcm_m
1	guez	3	module conf_gcm_m
2
3			IMPLICIT NONE
4
5	guez	212	INTEGER:: nday = 1 ! nombre de jours d'int\'egration
6	guez	154	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	guez	3
9	guez	39	integer:: iapp_tracvl = 5
10	guez	3	! Should normally be equal to "iperiod"
11			! frequence du groupement des flux (en pas de temps)
12
13	guez	39	integer:: iconser = 240
14			! number of time steps between output of control variables
15	guez	3
16	guez	39	integer:: iecri = 1 ! time interval between outputs to "dyn_hist.nc" (in days)
17	guez	3
18	guez	40	integer:: iphysiq = 5
19			! number of time steps of dynamics between calls to physics
20
21	guez	15	logical:: raz_date = .false.
22	guez	154	! 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	guez	98	! fichier restart)
25	guez	3
26	guez	56	integer:: periodav = 1
27			! time interval between outputs in the dynamical part, in days
28	guez	3
29	guez	223	logical:: offline = .FALSE. ! output variables for offline transport
30	guez	3
31	guez	154	integer:: prt_level = 0 ! niveau d'impression souhait\'e (0 = minimum)
32	guez	57
33			LOGICAL:: purmats= .FALSE.
34			! Help = Choix du schema d'integration temporel.
35			! y = pure Matsuno sinon c'est du Matsuno-leapfrog
36
37	guez	208	logical:: iflag_phys = .true. ! call parameterizations of physics
38	guez	202	INTEGER, SAVE:: lmt_pas ! number of time steps of "physics" per day
39
40	guez	3	contains
41
42	guez	13	SUBROUTINE conf_gcm
43	guez	3
44			! Auteurs : L. Fairhead, P. Le Van
45			! Version du 29/04/97
46
47	guez	33	use abort_gcm_m, only: abort_gcm
48	guez	212	use nr_util, only: assert
49	guez	57	use unit_nml_m, only: unit_nml
50	guez	3
51	guez	129	namelist /conf_gcm_nml/ raz_date, nday, day_step, iperiod, iapp_tracvl, &
52			iconser, iecri, periodav, iphysiq, offline
53	guez	3
54	guez	129	namelist /iniprint_nml/ prt_level
55	guez	57
56	guez	115	namelist /logic_nml/ purmats, iflag_phys
57	guez	57
58	guez	3	!------------------------------------
59
60			print *, "Call sequence information: conf_gcm"
61	guez	12
62	guez	57	print *, "Enter namelist 'iniprint_nml'."
63			read(unit=*, nml=iniprint_nml)
64			write(unit_nml, nml=iniprint_nml)
65
66			print *, "Enter namelist 'logic_nml'."
67			read(unit=*, nml=logic_nml)
68			write(unit_nml, nml=logic_nml)
69
70	guez	3	print *, "Enter namelist 'conf_gcm_nml'."
71			read(unit=*, nml=conf_gcm_nml)
72	guez	57	write(unit_nml, nml=conf_gcm_nml)
73	guez	3
74	guez	171	IF (MOD(day_step, iperiod) /= 0) call abort_gcm("conf_gcm", &
75			'Il faut choisir un nombre de pas par jour multiple de "iperiod".')
76	guez	33
77	guez	171	IF (MOD(day_step, iphysiq)/= 0) call abort_gcm("conf_gcm", &
78			'Il faut choisir un nombre de pas par jour multiple de "iphysiq".')
79	guez	33
80	guez	212	call assert(mod(iphysiq, iperiod) == 0, &
81			"conf_gcm -- iphysiq must be multiple of iperiod")
82
83	guez	202	lmt_pas = day_step / iphysiq
84			print *, 'Number of time steps of "physics" per day: ', lmt_pas
85
86	guez	3	END SUBROUTINE conf_gcm
87
88			end module conf_gcm_m