/[lmdze]/trunk/libf/dyn3d/dynredem0.f90
ViewVC logotype

Contents of /trunk/libf/dyn3d/dynredem0.f90

Parent Directory Parent Directory | Revision Log Revision Log

Revision 7 - (show annotations)
Mon Mar 31 12:24:17 2008 UTC (16 years, 1 month ago) by guez
File size: 9146 byte(s) 
This revision is not in working order. Pending some moving of files.

Important changes. In the program "etat0_lim": ozone coefficients from
Mobidic are regridded in time instead of pressure ; consequences in
"etat0". In the program "gcm", ozone coefficients from Mobidic are
read once per day only for the current day and regridded in pressure ;
consequences in "o3_chem_m", "regr_pr_coefoz", "phytrac" and
"regr_pr_comb_coefoz_m".

NetCDF95 is a library and does not export NetCDF.

New variables "nag_gl_options", "nag_fcalls_options" and
"nag_cross_options" in "nag_tools.mk".

"check_coefoz.jnl" rewritten entirely for new version of
"coefoz_LMDZ.nc".

Target "obj_etat0_lim" moved from "GNUmakefile" to "nag_rules.mk".

Added some "intent" attributes in "calfis", "clmain", "clqh",
"cltrac", "cltracrn", "cvltr", "ini_undefSTD", "moy_undefSTD",
"nflxtr", "phystokenc", "phytrac", "readsulfate", "readsulfate_preind"
and "undefSTD".

In "dynetat0", "dynredem0" and "gcm", "phis" has rank 2 instead of
1. "phis" has assumed shape in "dynredem0".

Added module containing "dynredem0". Changed some calls with NetCDF
Fortran 77 interface to calls with NetCDF95 interface.

Replaced calls to "ssum" by calls to "sum" in "inigeom".

In "make.sh", new option "-c" to change compiler.

In "aaam_bud", argument "rjour" deleted.

In "physiq": renamed some variables; deleted variable "xjour".

In "phytrac": renamed some variables; new argument "lmt_pas".
1 module dynredem0_m
2
3 IMPLICIT NONE
4
5 contains
6
7 SUBROUTINE dynredem0(fichnom, iday_end, phis)
8
9 ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30
10
11 ! Ecriture du fichier de redémarrage au format NetCDF (initialisation)
12
13 USE IOIPSL, only: ymds2ju, ju2ymds
14 use dimens_m, only: iim, jjm, llm, nqmx
15 use paramet_m, only: ip1jmp1, iip1, jjp1, llmp1
16 use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg
17 use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig
18 use logic
19 use comgeom
20 use serre
21 use temps, only: annee_ref, day_ref, itaufin, itau_dyn
22 use ener
23 use advtrac_m, only: tname, ttext
24 use netcdf95, only: nf95_create, nf95_put_att, nf95_def_dim, &
25 nf95_def_var, NF95_ENDDEF, NF95_PUT_VAR
26 use netcdf, only: NF90_CLOBBER, NF90_GLOBAL, NF90_UNLIMITED, nf90_float
27
28 CHARACTER(len=*), intent(in):: fichnom
29 INTEGER, intent(in):: iday_end
30 REAL, intent(in):: phis(:, :)
31
32 ! Local:
33
34 include "netcdf.inc"
35
36 INTEGER iq, l
37 INTEGER length
38 PARAMETER (length = 100)
39 REAL tab_cntrl(length) ! tableau des parametres du run
40 INTEGER ierr
41
42 ! Variables locales pour NetCDF:
43
44 INTEGER dims2(2), dims3(3), dims4(4)
45 INTEGER idim_index
46 INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
47 INTEGER idim_s, idim_sig
48 INTEGER idim_tim
49 INTEGER nid, nvarid
50
51 REAL zjulian, hours
52 INTEGER yyears0, jjour0, mmois0
53 character(len=30) unites
54
55 !-----------------------------------------------------------------------
56
57 print *, "Call sequence information: dynredem0"
58
59 call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
60 call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
61
62 DO l=1, length
63 tab_cntrl(l) = 0.
64 ENDDO
65 tab_cntrl(1) = REAL(iim)
66 tab_cntrl(2) = REAL(jjm)
67 tab_cntrl(3) = REAL(llm)
68 tab_cntrl(4) = REAL(day_ref)
69 tab_cntrl(5) = REAL(annee_ref)
70 tab_cntrl(6) = rad
71 tab_cntrl(7) = omeg
72 tab_cntrl(8) = g
73 tab_cntrl(9) = cpp
74 tab_cntrl(10) = kappa
75 tab_cntrl(11) = daysec
76 tab_cntrl(12) = dtvr
77 tab_cntrl(13) = etot0
78 tab_cntrl(14) = ptot0
79 tab_cntrl(15) = ztot0
80 tab_cntrl(16) = stot0
81 tab_cntrl(17) = ang0
82 tab_cntrl(18) = pa
83 tab_cntrl(19) = preff
84
85 ! Paramètres pour le zoom :
86
87 tab_cntrl(20) = clon
88 tab_cntrl(21) = clat
89 tab_cntrl(22) = grossismx
90 tab_cntrl(23) = grossismy
91
92 IF ( fxyhypb ) THEN
93 tab_cntrl(24) = 1.
94 tab_cntrl(25) = dzoomx
95 tab_cntrl(26) = dzoomy
96 tab_cntrl(27) = 0.
97 tab_cntrl(28) = taux
98 tab_cntrl(29) = tauy
99 ELSE
100 tab_cntrl(24) = 0.
101 tab_cntrl(25) = dzoomx
102 tab_cntrl(26) = dzoomy
103 tab_cntrl(27) = 0.
104 tab_cntrl(28) = 0.
105 tab_cntrl(29) = 0.
106 IF( ysinus ) tab_cntrl(27) = 1.
107 ENDIF
108
109 tab_cntrl(30) = REAL(iday_end)
110 tab_cntrl(31) = REAL(itau_dyn + itaufin)
111
112 call nf95_create(fichnom, NF90_CLOBBER, nid)
113 call nf95_put_att(nid, NF90_GLOBAL, "title", &
114 "Fichier de démarrage dynamique")
115
116 ! Definir les dimensions du fichiers:
117
118 call nf95_def_dim(nid, "index", length, idim_index)
119 call NF95_DEF_DIM(nid, "rlonu", iip1, idim_rlonu)
120 call NF95_DEF_DIM(nid, "rlatu", jjp1, idim_rlatu)
121 call NF95_DEF_DIM(nid, "rlonv", iip1, idim_rlonv)
122 call NF95_DEF_DIM(nid, "rlatv", jjm, idim_rlatv)
123 call NF95_DEF_DIM(nid, "sigs", llm, idim_s)
124 call NF95_DEF_DIM(nid, "sig", llmp1, idim_sig)
125 call NF95_DEF_DIM(nid, "temps", NF90_UNLIMITED, idim_tim)
126
127 ! Definir et enregistrer certains champs invariants:
128
129 call nf95_def_var(nid, "controle", NF90_FLOAT, idim_index, nvarid)
130 call nf95_put_att(nid, nvarid, "title", "Parametres de controle")
131 call NF95_ENDDEF(nid)
132 call NF95_PUT_VAR(nid, nvarid, tab_cntrl)
133
134 ierr = NF_REDEF (nid)
135 call nf95_def_var(nid, "rlonu", NF90_FLOAT, idim_rlonu, nvarid)
136 call nf95_put_att(nid, nvarid, "title", "Longitudes des points U")
137 call NF95_ENDDEF(nid)
138 call NF95_PUT_VAR(nid, nvarid, rlonu)
139
140 ierr = NF_REDEF (nid)
141 call nf95_def_var(nid, "rlatu", NF90_FLOAT, idim_rlatu, nvarid)
142 call nf95_put_att(nid, nvarid, "title", "Latitudes des points U")
143 call NF95_ENDDEF(nid)
144 call NF95_PUT_VAR(nid, nvarid, rlatu)
145
146 ierr = NF_REDEF (nid)
147 call nf95_def_var(nid, "rlonv", NF90_FLOAT, idim_rlonv, nvarid)
148 call nf95_put_att(nid, nvarid, "title", "Longitudes des points V")
149 call NF95_ENDDEF(nid)
150 call NF95_PUT_VAR(nid, nvarid, rlonv)
151
152 ierr = NF_REDEF (nid)
153 call nf95_def_var(nid, "rlatv", NF90_FLOAT, idim_rlatv, nvarid)
154 call nf95_put_att(nid, nvarid, "title", "Latitudes des points V")
155 call NF95_ENDDEF(nid)
156 call NF95_PUT_VAR(nid, nvarid, rlatv)
157
158 ierr = NF_REDEF (nid)
159 call nf95_def_var(nid, "nivsigs", NF90_FLOAT, idim_s, nvarid)
160 call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches s")
161 call NF95_ENDDEF(nid)
162 call NF95_PUT_VAR(nid, nvarid, nivsigs)
163
164 ierr = NF_REDEF (nid)
165 call nf95_def_var(nid, "nivsig", NF90_FLOAT, idim_sig, nvarid)
166 call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches sigma")
167 call NF95_ENDDEF(nid)
168 call NF95_PUT_VAR(nid, nvarid, nivsig)
169
170 ierr = NF_REDEF (nid)
171 call nf95_def_var(nid, "ap", NF90_FLOAT, idim_sig, nvarid)
172 call nf95_put_att(nid, nvarid, "title", "Coefficient A pour hybride")
173 call NF95_ENDDEF(nid)
174 call NF95_PUT_VAR(nid, nvarid, ap)
175
176 ierr = NF_REDEF (nid)
177 call nf95_def_var(nid, "bp", NF90_FLOAT, idim_sig, nvarid)
178 call nf95_put_att(nid, nvarid, "title", "Coefficient B pour hybride")
179 call NF95_ENDDEF(nid)
180 call NF95_PUT_VAR(nid, nvarid, bp)
181
182 ierr = NF_REDEF (nid)
183 call nf95_def_var(nid, "presnivs", NF90_FLOAT, idim_s, nvarid)
184 call NF95_ENDDEF(nid)
185 call NF95_PUT_VAR(nid, nvarid, presnivs)
186
187 ! Coefficients de passage cov. <-> contra. <--> naturel
188
189 ierr = NF_REDEF (nid)
190 dims2(1) = idim_rlonu
191 dims2(2) = idim_rlatu
192 call nf95_def_var(nid, "cu", NF90_FLOAT, dims2, nvarid)
193 call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour U")
194 call NF95_ENDDEF(nid)
195 call NF95_PUT_VAR(nid, nvarid, cu_2d)
196
197 ierr = NF_REDEF (nid)
198 dims2(1) = idim_rlonv
199 dims2(2) = idim_rlatv
200 call nf95_def_var(nid, "cv", NF90_FLOAT, dims2, nvarid)
201 call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour V")
202 call NF95_ENDDEF(nid)
203 call NF95_PUT_VAR(nid, nvarid, cv_2d)
204
205 ! Aire de chaque maille:
206
207 ierr = NF_REDEF (nid)
208 dims2(1) = idim_rlonv
209 dims2(2) = idim_rlatu
210 call nf95_def_var(nid, "aire", NF90_FLOAT, dims2, nvarid)
211 call nf95_put_att(nid, nvarid, "title", "Aires de chaque maille")
212 call NF95_ENDDEF(nid)
213 call NF95_PUT_VAR(nid, nvarid, aire_2d)
214
215 ! Geopentiel au sol:
216
217 ierr = NF_REDEF (nid)
218 dims2(1) = idim_rlonv
219 dims2(2) = idim_rlatu
220 call nf95_def_var(nid, "phisinit", NF90_FLOAT, dims2, nvarid)
221 call nf95_put_att(nid, nvarid, "title", "Geopotentiel au sol")
222 call NF95_ENDDEF(nid)
223 call NF95_PUT_VAR(nid, nvarid, phis)
224
225 ! Definir les variables pour pouvoir les enregistrer plus tard:
226
227 ierr = NF_REDEF (nid) ! entrer dans le mode de definition
228
229 call nf95_def_var(nid, "temps", NF90_FLOAT, idim_tim, nvarid)
230 call nf95_put_att(nid, nvarid, "title", "Temps de simulation")
231 write(unites, 200)yyears0, mmois0, jjour0
232 200 format('days since ', i4, '-', i2.2, '-', i2.2, ' 00:00:00')
233 call nf95_put_att(nid, nvarid, "units", unites)
234
235
236 dims4(1) = idim_rlonu
237 dims4(2) = idim_rlatu
238 dims4(3) = idim_s
239 dims4(4) = idim_tim
240 call nf95_def_var(nid, "ucov", NF90_FLOAT, dims4, nvarid)
241 call nf95_put_att(nid, nvarid, "title", "Vitesse U")
242
243 dims4(1) = idim_rlonv
244 dims4(2) = idim_rlatv
245 dims4(3) = idim_s
246 dims4(4) = idim_tim
247 call nf95_def_var(nid, "vcov", NF90_FLOAT, dims4, nvarid)
248 call nf95_put_att(nid, nvarid, "title", "Vitesse V")
249
250 dims4(1) = idim_rlonv
251 dims4(2) = idim_rlatu
252 dims4(3) = idim_s
253 dims4(4) = idim_tim
254 call nf95_def_var(nid, "teta", NF90_FLOAT, dims4, nvarid)
255 call nf95_put_att(nid, nvarid, "title", "Temperature")
256
257 dims4(1) = idim_rlonv
258 dims4(2) = idim_rlatu
259 dims4(3) = idim_s
260 dims4(4) = idim_tim
261 DO iq=1, nqmx
262 call nf95_def_var(nid, tname(iq), NF90_FLOAT, dims4, nvarid)
263 call nf95_put_att(nid, nvarid, "title", ttext(iq))
264 ENDDO
265
266 dims4(1) = idim_rlonv
267 dims4(2) = idim_rlatu
268 dims4(3) = idim_s
269 dims4(4) = idim_tim
270 call nf95_def_var(nid, "masse", NF90_FLOAT, dims4, nvarid)
271 call nf95_put_att(nid, nvarid, "title", "C est quoi ?")
272
273 dims3(1) = idim_rlonv
274 dims3(2) = idim_rlatu
275 dims3(3) = idim_tim
276 call nf95_def_var(nid, "ps", NF90_FLOAT, dims3, nvarid)
277 call nf95_put_att(nid, nvarid, "title", "Pression au sol")
278
279 ierr = NF_ENDDEF(nid) ! sortir du mode de definition
280 ierr = NF_CLOSE(nid) ! fermer le fichier
281
282 PRINT*, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
283 PRINT*, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
284
285 END SUBROUTINE dynredem0
286
287 end module dynredem0_m
  ViewVC Help
Powered by ViewVC 1.1.21