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

Annotation of /trunk/dyn3d/gcm.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 115 - (hide annotations)
Fri Sep 19 17:36:20 2014 UTC (9 years, 7 months ago) by guez
File size: 5010 byte(s) 
Extracted code from tau2alpha for first call into new procedure
init_tau2alpha. dxdys, dxdyu, dxdyv are now local variables if guide
computed by init_tau2alpha. This allows us to remove terrible argument
type of tau2alpha: we just give to tau2alpha the right dxdy and
rlat.

In module conf_guide_m, changed default values of tau_min_*, because
0.02 is too small for the default daystep = 240, iperiod = 5. Changed
default values of guide_[uv] to false. Moved variable ok_guide from
conf_gcm_m to conf_guide_m, ok_guide is no longer an input parameter,
it is computed from guide_*. Had then to move test on ok_guide and
day_step from conf_gcm_m to conf_guide_m. Added checks on input
nudging parameters in procedure conf_guide. Upgraded variable factt to
module conf_guide_m in order to check nudging parameters. Bug fix:
variable guide_q was not in namelist conf_guide_nml.

Removed unused variables aire_min, aire_max of MODULE guide_m.

Moved the call to conf_guide from guide to gcm. This was needed to
define ok_guide before getting into guide (since ok_guide is no longer
in conf_gcm_m). Moved test on grossismx and grossismy from tau2alpha
to guide. It is clearer now that only tau_max is used for a regular
grid, and we do not have to repeat this test in each call to
tau2alpha. In guide, we only call writefield when alpha is not a
constant.
1 guez 3 PROGRAM gcm
2    
3 guez 36 ! Authors: P. Le Van, L. Fairhead, F. Hourdin
4     ! From "gcm.F", version 1.4, 2006/04/04 15:05:16
5 guez 3
6 guez 38 ! General circulation model of LMD. Avec coordonnée verticale
7     ! hybride, avec nouveaux opérateurs de dissipation "*" ("gradiv2",
8     ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour
9     ! l'advection de "q", en modifiant "iadv" dans "traceur.def".
10    
11 guez 79 use comconst, only: daysec, dtvr, iniconst
12 guez 78 use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom
13 guez 37 use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
14 guez 3 use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
15 guez 57 nday, raz_date, periodav, conf_gcm, iflag_phys
16 guez 115 use conf_guide_m, only: conf_guide
17 guez 26 use dimens_m, only: iim, jjm, llm, nqmx
18     use dimphy, only: klon
19 guez 79 USE disvert_m, ONLY : disvert
20 guez 25 use dynetat0_m, only: dynetat0, day_ini
21 guez 26 use dynredem0_m, only: dynredem0
22 guez 3 use grid_change, only: dyn_phy, init_dyn_phy
23 guez 61 use histclo_m, only: histclo
24 guez 18 use iniadvtrac_m, only: iniadvtrac
25 guez 26 use inidissip_m, only: inidissip
26 guez 54 use inifilr_m, only: inifilr
27 guez 26 use initdynav_m, only: initdynav
28     use inithist_m, only: inithist
29 guez 57 use init_dynzon_m, only: init_dynzon
30 guez 92 USE ioconf_calendar_m, only: ioconf_calendar
31 guez 57 use jumble, only: new_unit
32 guez 3 use leapfrog_m, only: leapfrog
33 guez 108 use netcdf95, only: nf95_close
34 guez 37 use suphec_m, only: suphec
35 guez 26 use temps, only: day_ref, annee_ref, day_end, itau_dyn
36     use tracstoke, only: istdyn, istphy
37 guez 57 use unit_nml_m, only: unit_nml
38 guez 38 use yoethf_m, only: yoethf
39 guez 110 use createnewfield_m, only: NbField, Ncid
40 guez 3
41     IMPLICIT NONE
42    
43     ! Variables dynamiques :
44 guez 55 REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm) ! vent covariant
45 guez 43 REAL teta(iim + 1, jjm + 1, llm) ! température potentielle
46 guez 40 REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés
47 guez 39 REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
48 guez 55 REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
49 guez 7 REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol
50 guez 3
51     ! Variables pour le fichier histoire :
52 guez 25 REAL time_0 ! time in day, as a fraction of day, in [0, 1[
53 guez 3
54     ! Calendrier :
55     LOGICAL:: true_calendar = .false. ! default value
56    
57     logical mask_v(iim + 1, jjm)
58     ! (mask for points in the "v" grid, first index is for longitude,
59     ! second index is for latitude)
60    
61 guez 108 integer i
62    
63 guez 3 namelist /main_nml/true_calendar
64    
65     !------------------------------------------------------------
66    
67 guez 57 call new_unit(unit_nml)
68 guez 79 open(unit_nml, file="used_namelists.txt", status="replace", action="write")
69 guez 57
70     CALL conf_gcm
71    
72 guez 3 print *, "Enter namelist 'main_nml'."
73     read (unit=*, nml=main_nml)
74 guez 57 write(unit_nml, nml=main_nml)
75 guez 3
76     ! Choix du calendrier :
77     if (true_calendar) then
78     call ioconf_calendar('gregorian')
79     else
80     call ioconf_calendar('360d')
81     endif
82    
83     ! Initialisation des traceurs
84     ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
85 guez 23 call iniadvtrac
86 guez 3
87 guez 79 CALL iniconst
88    
89 guez 3 ! Lecture du fichier "start.nc" :
90     CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
91    
92 guez 98 ! On remet le calendrier à zéro si demandé :
93     if (raz_date) then
94     print *, 'On réinitialise à la date lue dans la namelist.'
95     annee_ref = anneeref
96     day_ref = dayref
97     day_ini = dayref
98     itau_dyn = 0
99     time_0 = 0.
100     else
101     print *, 'On garde les dates du fichier "start".'
102 guez 3 endif
103    
104 guez 79 CALL disvert
105 guez 36 CALL inigeom ! initialisation de la géometrie
106     CALL inifilr ! initialisation du filtre
107 guez 27 CALL inidissip
108 guez 3 call init_dyn_phy
109    
110     ! Initialisation de la physique :
111     IF (iflag_phys == 1) THEN
112 guez 37 rlatd(1)=rlatu(1)
113     rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
114     rlatd(klon)= rlatu(jjm + 1)
115 guez 3
116 guez 37 rlond(1)=0.
117     rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
118     rlond(klon)= 0.
119 guez 3
120 guez 37 cuphy = pack(cu_2d, dyn_phy)
121 guez 3
122     ! Construct a mask for points in the "v" grid:
123     mask_v = .true.
124     mask_v(2:, 1) = .false.
125     mask_v(iim + 1, 2:) = .false.
126    
127 guez 37 cvphy(:klon - 1) = pack(cv_2d, mask_v)
128     cvphy(klon) = cv_2d(1, jjm)
129     ! (that value of "cv_2d" is used twice in "cvphy")
130 guez 3
131 guez 37 airephy = pack(aire_2d, dyn_phy)
132     CALL suphec
133 guez 38 call yoethf
134 guez 3 ENDIF
135    
136     ! Initialisation des entrées-sorties :
137     day_end = day_ini + nday
138     print *, "day_ini = ", day_ini
139     print *, "day_end = ", day_end
140    
141 guez 5 CALL dynredem0("restart.nc", day_end, phis)
142 guez 79 CALL inithist(day_ref, annee_ref, dtvr, nqmx, t_ops = iecri * daysec, &
143 guez 56 t_wrt = iecri * daysec)
144 guez 79 CALL initdynav(day_ref, annee_ref, dtvr, nqmx, t_ops = iperiod * dtvr, &
145 guez 56 t_wrt = periodav * daysec)
146 guez 57 call init_dynzon(dt_app = dtvr * iperiod)
147 guez 3
148     ! Choix des fréquences de stockage pour le hors-ligne :
149 guez 36 istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
150 guez 3 istphy = istdyn / iphysiq
151    
152 guez 115 CALL conf_guide
153    
154 guez 3 ! Intégration temporelle du modèle :
155 guez 23 CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
156 guez 3
157 guez 68 close(unit_nml)
158 guez 10 call histclo
159 guez 108
160     do i = 1, nbfield
161 guez 109 call nf95_close(Ncid(i))
162 guez 108 end do
163    
164 guez 10 print *, 'Simulation finished'
165     print *, 'Everything is cool'
166    
167 guez 3 END PROGRAM gcm
  ViewVC Help
Powered by ViewVC 1.1.21