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

Contents of /trunk/dyn3d/conf_gcm.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 263 - (show annotations)
Wed Mar 7 14:41:46 2018 UTC (6 years, 2 months ago) by guez
File size: 2666 byte(s)
Allow output of dynamics at frequency greater than one per day. Change
meaning of iecri.

Remove arguments tstep, nq of procedure inithist: use directly
module variables dtvr and nqmx instead.

Change meaning of last dummy argument of writehist. Simpler and
clearer, because time suggested a dimensional quantity.

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

  ViewVC Help
Powered by ViewVC 1.1.21