/[lmdze]/trunk/Sources/phylmd/phyredem.f
ViewVC logotype

Contents of /trunk/Sources/phylmd/phyredem.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 12 - (show annotations)
Mon Jul 21 16:05:07 2008 UTC (15 years, 10 months ago) by guez
Original Path: trunk/libf/phylmd/phyredem.f90
File size: 13041 byte(s) 
-- Minor modification of input/output:

Created procedure "read_logic". Variables of module "logic" are read
by "read_logic" instead of "conf_gcm". Variable "offline" of module
"conf_gcm" is read from namelist instead of "*.def".

Deleted arguments "dtime", "co2_ppm_etat0", "solaire_etat0",
"tabcntr0" and local variables "radpas", "tab_cntrl" of
"phyetat0". "phyetat0" does not read "controle" in "startphy.nc" any
longer. "phyetat0" now reads global attribute "itau_phy" from
"startphy.nc". "phyredem" does not create variable "controle" in
"startphy.nc" any longer. "phyredem" now writes global attribute
"itau_phy" of "startphy.nc". Deleted argument "tabcntr0" of
"printflag". Removed diagnostic messages written by "printflag" for
comparison of the variable "controle" of "startphy.nc" and the
variables read from "*.def" or namelist input.

-- Removing unwanted functionality:

Removed variable "lunout" from module "iniprint", replaced everywhere
by standard output.

Removed case "ocean == 'couple'" in "clmain", "interfsurf_hq" and
"physiq". Removed procedure "interfoce_cpl".

-- Should not change anything at run time:

Automated creation of graphs in documentation. More documentation on
input files.

Converted Fortran files to free format: "phyredem.f90", "printflag.f90".

Split module "clesphy" into "clesphys" and "clesphys2".

Removed variables "conser", "leapf", "forward", "apphys", "apdiss" and
"statcl" from module "logic". Added arguments "conser" to "advect",
"leapf" to "integrd". Added local variables "forward", "leapf",
"apphys", "conser", "apdiss" in "leapfrog".

Added intent attributes.

Deleted arguments "dtime" of "phyredem", "pdtime" of "flxdtdq", "sh"
of "phytrac", "dt" of "yamada".

Deleted local variables "dtime", "co2_ppm_etat0", "solaire_etat0",
"length", "tabcntr0" in "physiq". Replaced all references to "dtime"
by references to "pdtphys".
1 SUBROUTINE phyredem(fichnom,radpas,rlat,rlon,pctsrf,tsol,tsoil, &
2 tslab,seaice,qsurf,qsol,snow,albedo,alblw,evap,rain_fall,snow_fall, &
3 solsw,sollw,fder,radsol,frugs,agesno,zmea,zstd,zsig,zgam,zthe,zpic, &
4 zval,rugsrel,t_ancien,q_ancien,rnebcon,ratqs,clwcon,run_off_lic_0)
5
6 ! From phylmd/phyredem.F,v 1.3 2005/05/25 13:10:09
7 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
8 ! Objet: Ecriture de l'etat de démarrage ou redémarrage pour la physique
9
10 USE indicesol, ONLY : is_lic, is_oce, is_sic, is_ter, nbsrf
11 USE dimphy, ONLY : klev, klon, zmasq
12 USE dimsoil, ONLY : nsoilmx
13 USE temps, ONLY : itau_phy
14 USE netcdf95, ONLY : nf95_create, nf95_put_att
15 USE netcdf, ONLY : nf90_clobber, nf90_global
16
17 IMPLICIT NONE
18
19 INCLUDE 'netcdf.inc'
20
21 CHARACTER(len=*) fichnom
22 INTEGER :: radpas
23 REAL, INTENT (IN) :: rlat(klon), rlon(klon)
24 REAL :: tsol(klon,nbsrf)
25 REAL :: tsoil(klon,nsoilmx,nbsrf)
26 !IM "slab" ocean
27 REAL :: tslab(klon), seaice(klon) !IM "slab" ocean
28 REAL :: qsurf(klon,nbsrf)
29 REAL :: qsol(klon)
30 REAL :: snow(klon,nbsrf)
31 REAL :: albedo(klon,nbsrf)
32 !IM BEG
33 REAL :: alblw(klon,nbsrf)
34 !IM END
35 REAL :: evap(klon,nbsrf)
36 REAL :: rain_fall(klon)
37 REAL :: snow_fall(klon)
38 REAL :: solsw(klon)
39 REAL :: sollw(klon)
40 REAL :: fder(klon)
41 REAL :: radsol(klon)
42 REAL :: frugs(klon,nbsrf)
43 REAL :: agesno(klon,nbsrf)
44 REAL :: zmea(klon)
45 REAL :: zstd(klon)
46 REAL :: zsig(klon)
47 REAL :: zgam(klon)
48 REAL :: zthe(klon)
49 REAL :: zpic(klon)
50 REAL :: zval(klon)
51 REAL :: rugsrel(klon)
52 REAL :: pctsrf(klon,nbsrf)
53 REAL :: t_ancien(klon,klev), q_ancien(klon,klev)
54 REAL :: clwcon(klon,klev), rnebcon(klon,klev), ratqs(klon,klev)
55 REAL :: run_off_lic_0(klon)
56
57 INTEGER :: nid, nvarid, idim2, idim3
58 INTEGER :: ierr
59
60 INTEGER :: isoil, nsrf
61 CHARACTER (7) :: str7
62 CHARACTER (2) :: str2
63
64 !------------------------------------------------------------
65
66 PRINT *, 'Call sequence information: phyredem'
67 CALL nf95_create(fichnom,nf90_clobber,nid)
68
69 call nf95_put_att(nid, nf90_global, 'title', 'Fichier redémarrage physique')
70 call nf95_put_att(nid, nf90_global, "itau_phy", itau_phy)
71
72 ierr = nf_def_dim(nid,'points_physiques',klon,idim2)
73 ierr = nf_def_dim(nid,'horizon_vertical',klon*klev,idim3)
74
75 ierr = nf_def_var(nid,'longitude',nf_float,1,idim2,nvarid)
76 ierr = nf_put_att_text(nid,nvarid,'title',32, &
77 'Longitudes de la grille physique')
78 ierr = nf_enddef(nid)
79 ierr = nf_put_var_real(nid,nvarid,rlon)
80
81 ierr = nf_redef(nid)
82 ierr = nf_def_var(nid,'latitude',nf_float,1,idim2,nvarid)
83 ierr = nf_put_att_text(nid,nvarid,'title',31, &
84 'Latitudes de la grille physique')
85 ierr = nf_enddef(nid)
86 ierr = nf_put_var_real(nid,nvarid,rlat)
87
88 ! PB ajout du masque terre/mer
89
90 ierr = nf_redef(nid)
91 ierr = nf_def_var(nid,'masque',nf_float,1,idim2,nvarid)
92 ierr = nf_put_att_text(nid,nvarid,'title',16,'masque terre mer')
93 ierr = nf_enddef(nid)
94 ierr = nf_put_var_real(nid,nvarid,zmasq)
95 ! BP ajout des fraction de chaque sous-surface
96
97 ! 1. fraction de terre
98
99 ierr = nf_redef(nid)
100 ierr = nf_def_var(nid,'FTER',nf_float,1,idim2,nvarid)
101 ierr = nf_put_att_text(nid,nvarid,'title',21,'fraction de continent')
102 ierr = nf_enddef(nid)
103 ierr = nf_put_var_real(nid,nvarid,pctsrf(1:klon,is_ter))
104
105 ! 2. Fraction de glace de terre
106
107 ierr = nf_redef(nid)
108 ierr = nf_def_var(nid,'FLIC',nf_float,1,idim2,nvarid)
109 ierr = nf_put_att_text(nid,nvarid,'title',24,'fraction glace de terre')
110 ierr = nf_enddef(nid)
111 ierr = nf_put_var_real(nid,nvarid,pctsrf(1:klon,is_lic))
112
113 ! 3. fraction ocean
114
115 ierr = nf_redef(nid)
116 ierr = nf_def_var(nid,'FOCE',nf_float,1,idim2,nvarid)
117 ierr = nf_put_att_text(nid,nvarid,'title',14,'fraction ocean')
118 ierr = nf_enddef(nid)
119 ierr = nf_put_var_real(nid,nvarid,pctsrf(1:klon,is_oce))
120
121 ! 4. Fraction glace de mer
122
123 ierr = nf_redef(nid)
124 ierr = nf_def_var(nid,'FSIC',nf_float,1,idim2,nvarid)
125 ierr = nf_put_att_text(nid,nvarid,'title',18,'fraction glace mer')
126 ierr = nf_enddef(nid)
127 ierr = nf_put_var_real(nid,nvarid,pctsrf(1:klon,is_sic))
128
129
130
131 DO nsrf = 1, nbsrf
132 IF (nsrf<=99) THEN
133 WRITE (str2,'(i2.2)') nsrf
134 ierr = nf_redef(nid)
135 ierr = nf_def_var(nid,'TS'//str2,nf_float,1,idim2,nvarid)
136 ierr = nf_put_att_text(nid,nvarid,'title',28, &
137 'Temperature de surface No.'//str2)
138 ierr = nf_enddef(nid)
139 ELSE
140 PRINT *, 'Trop de sous-mailles'
141 STOP 1
142 END IF
143 ierr = nf_put_var_real(nid,nvarid,tsol(1,nsrf))
144 END DO
145
146 DO nsrf = 1, nbsrf
147 DO isoil = 1, nsoilmx
148 IF (isoil<=99 .AND. nsrf<=99) THEN
149 WRITE (str7,'(i2.2,"srf",i2.2)') isoil, nsrf
150 ierr = nf_redef(nid)
151 ierr = nf_def_var(nid,'Tsoil'//str7,nf_float,1,idim2,nvarid)
152 ierr = nf_put_att_text(nid,nvarid,'title',29, &
153 'Temperature du sol No.'//str7)
154 ierr = nf_enddef(nid)
155 ELSE
156 PRINT *, 'Trop de couches'
157 STOP 1
158 END IF
159 ierr = nf_put_var_real(nid,nvarid,tsoil(1,isoil,nsrf))
160 END DO
161 END DO
162
163 !IM "slab" ocean
164 ierr = nf_redef(nid)
165 ierr = nf_def_var(nid,'TSLAB',nf_float,1,idim2,nvarid)
166 ierr = nf_put_att_text(nid,nvarid,'title',33, &
167 'Ecart de la SST (pour slab-ocean)')
168 ierr = nf_enddef(nid)
169 ierr = nf_put_var_real(nid,nvarid,tslab)
170
171 ierr = nf_redef(nid)
172 ierr = nf_def_var(nid,'SEAICE',nf_float,1,idim2,nvarid)
173 ierr = nf_put_att_text(nid,nvarid,'title',33, &
174 'Glace de mer kg/m2 (pour slab-ocean)')
175 ierr = nf_enddef(nid)
176 ierr = nf_put_var_real(nid,nvarid,seaice)
177
178 DO nsrf = 1, nbsrf
179 IF (nsrf<=99) THEN
180 WRITE (str2,'(i2.2)') nsrf
181 ierr = nf_redef(nid)
182 ierr = nf_def_var(nid,'QS'//str2,nf_float,1,idim2,nvarid)
183 ierr = nf_put_att_text(nid,nvarid,'title',25, &
184 'Humidite de surface No.'//str2)
185 ierr = nf_enddef(nid)
186 ELSE
187 PRINT *, 'Trop de sous-mailles'
188 STOP 1
189 END IF
190 ierr = nf_put_var_real(nid,nvarid,qsurf(1,nsrf))
191 END DO
192
193 ierr = nf_redef(nid)
194 ierr = nf_def_var(nid,'QSOL',nf_float,1,idim2,nvarid)
195 ierr = nf_put_att_text(nid,nvarid,'title',20,'Eau dans le sol (mm)')
196 ierr = nf_enddef(nid)
197 ierr = nf_put_var_real(nid,nvarid,qsol)
198
199 DO nsrf = 1, nbsrf
200 IF (nsrf<=99) THEN
201 WRITE (str2,'(i2.2)') nsrf
202 ierr = nf_redef(nid)
203 ierr = nf_def_var(nid,'ALBE'//str2,nf_float,1,idim2,nvarid)
204 ierr = nf_put_att_text(nid,nvarid,'title',23, &
205 'albedo de surface No.'//str2)
206 ierr = nf_enddef(nid)
207 ELSE
208 PRINT *, 'Trop de sous-mailles'
209 STOP 1
210 END IF
211 ierr = nf_put_var_real(nid,nvarid,albedo(1,nsrf))
212 END DO
213
214 !IM BEG albedo LW
215 DO nsrf = 1, nbsrf
216 IF (nsrf<=99) THEN
217 WRITE (str2,'(i2.2)') nsrf
218 ierr = nf_redef(nid)
219 ierr = nf_def_var(nid,'ALBLW'//str2,nf_float,1,idim2,nvarid)
220 ierr = nf_put_att_text(nid,nvarid,'title',23, &
221 'albedo LW de surface No.'//str2)
222 ierr = nf_enddef(nid)
223 ELSE
224 PRINT *, 'Trop de sous-mailles'
225 STOP 1
226 END IF
227 ierr = nf_put_var_real(nid,nvarid,alblw(1,nsrf))
228 END DO
229 !IM END albedo LW
230
231 DO nsrf = 1, nbsrf
232 IF (nsrf<=99) THEN
233 WRITE (str2,'(i2.2)') nsrf
234 ierr = nf_redef(nid)
235 ierr = nf_def_var(nid,'EVAP'//str2,nf_float,1,idim2,nvarid)
236 ierr = nf_put_att_text(nid,nvarid,'title',28, &
237 'Evaporation de surface No.'//str2)
238 ierr = nf_enddef(nid)
239 ELSE
240 PRINT *, 'Trop de sous-mailles'
241 STOP 1
242 END IF
243 ierr = nf_put_var_real(nid,nvarid,evap(1,nsrf))
244 END DO
245
246
247 DO nsrf = 1, nbsrf
248 IF (nsrf<=99) THEN
249 WRITE (str2,'(i2.2)') nsrf
250 ierr = nf_redef(nid)
251 ierr = nf_def_var(nid,'SNOW'//str2,nf_float,1,idim2,nvarid)
252 ierr = nf_put_att_text(nid,nvarid,'title',22, &
253 'Neige de surface No.'//str2)
254 ierr = nf_enddef(nid)
255 ELSE
256 PRINT *, 'Trop de sous-mailles'
257 STOP 1
258 END IF
259 ierr = nf_put_var_real(nid,nvarid,snow(1,nsrf))
260 END DO
261
262
263 ierr = nf_redef(nid)
264 ierr = nf_def_var(nid,'RADS',nf_float,1,idim2,nvarid)
265 ierr = nf_put_att_text(nid,nvarid,'title',28, &
266 'Rayonnement net a la surface')
267 ierr = nf_enddef(nid)
268 ierr = nf_put_var_real(nid,nvarid,radsol)
269
270 ierr = nf_redef(nid)
271 ierr = nf_def_var(nid,'solsw',nf_float,1,idim2,nvarid)
272 ierr = nf_put_att_text(nid,nvarid,'title',32, &
273 'Rayonnement solaire a la surface')
274 ierr = nf_enddef(nid)
275 ierr = nf_put_var_real(nid,nvarid,solsw)
276
277 ierr = nf_redef(nid)
278 ierr = nf_def_var(nid,'sollw',nf_float,1,idim2,nvarid)
279 ierr = nf_put_att_text(nid,nvarid,'title',27, &
280 'Rayonnement IF a la surface')
281 ierr = nf_enddef(nid)
282 ierr = nf_put_var_real(nid,nvarid,sollw)
283
284 ierr = nf_redef(nid)
285 ierr = nf_def_var(nid,'fder',nf_float,1,idim2,nvarid)
286 ierr = nf_put_att_text(nid,nvarid,'title',14,'Derive de flux')
287 ierr = nf_enddef(nid)
288 ierr = nf_put_var_real(nid,nvarid,fder)
289
290 ierr = nf_redef(nid)
291 ierr = nf_def_var(nid,'rain_f',nf_float,1,idim2,nvarid)
292 ierr = nf_put_att_text(nid,nvarid,'title',21,'precipitation liquide')
293 ierr = nf_enddef(nid)
294 ierr = nf_put_var_real(nid,nvarid,rain_fall)
295
296 ierr = nf_redef(nid)
297 ierr = nf_def_var(nid,'snow_f',nf_float,1,idim2,nvarid)
298 ierr = nf_put_att_text(nid,nvarid,'title',20,'precipitation solide')
299 ierr = nf_enddef(nid)
300 ierr = nf_put_var_real(nid,nvarid,snow_fall)
301
302 DO nsrf = 1, nbsrf
303 IF (nsrf<=99) THEN
304 WRITE (str2,'(i2.2)') nsrf
305 ierr = nf_redef(nid)
306 ierr = nf_def_var(nid,'RUG'//str2,nf_float,1,idim2,nvarid)
307 ierr = nf_put_att_text(nid,nvarid,'title',23, &
308 'rugosite de surface No.'//str2)
309 ierr = nf_enddef(nid)
310 ELSE
311 PRINT *, 'Trop de sous-mailles'
312 STOP 1
313 END IF
314 ierr = nf_put_var_real(nid,nvarid,frugs(1,nsrf))
315 END DO
316
317 DO nsrf = 1, nbsrf
318 IF (nsrf<=99) THEN
319 WRITE (str2,'(i2.2)') nsrf
320 ierr = nf_redef(nid)
321 ierr = nf_def_var(nid,'AGESNO'//str2,nf_float,1,idim2,nvarid)
322 ierr = nf_put_att_text(nid,nvarid,'title',15, &
323 'Age de la neige surface No.'//str2)
324 ierr = nf_enddef(nid)
325 ELSE
326 PRINT *, 'Trop de sous-mailles'
327 STOP 1
328 END IF
329 ierr = nf_put_var_real(nid,nvarid,agesno(1,nsrf))
330 END DO
331
332 ierr = nf_redef(nid)
333 ierr = nf_def_var(nid,'ZMEA',nf_float,1,idim2,nvarid)
334 ierr = nf_enddef(nid)
335 ierr = nf_put_var_real(nid,nvarid,zmea)
336
337 ierr = nf_redef(nid)
338 ierr = nf_def_var(nid,'ZSTD',nf_float,1,idim2,nvarid)
339 ierr = nf_enddef(nid)
340 ierr = nf_put_var_real(nid,nvarid,zstd)
341 ierr = nf_redef(nid)
342 ierr = nf_def_var(nid,'ZSIG',nf_float,1,idim2,nvarid)
343 ierr = nf_enddef(nid)
344 ierr = nf_put_var_real(nid,nvarid,zsig)
345 ierr = nf_redef(nid)
346 ierr = nf_def_var(nid,'ZGAM',nf_float,1,idim2,nvarid)
347 ierr = nf_enddef(nid)
348 ierr = nf_put_var_real(nid,nvarid,zgam)
349 ierr = nf_redef(nid)
350 ierr = nf_def_var(nid,'ZTHE',nf_float,1,idim2,nvarid)
351 ierr = nf_enddef(nid)
352 ierr = nf_put_var_real(nid,nvarid,zthe)
353 ierr = nf_redef(nid)
354 ierr = nf_def_var(nid,'ZPIC',nf_float,1,idim2,nvarid)
355 ierr = nf_enddef(nid)
356 ierr = nf_put_var_real(nid,nvarid,zpic)
357 ierr = nf_redef(nid)
358 ierr = nf_def_var(nid,'ZVAL',nf_float,1,idim2,nvarid)
359 ierr = nf_enddef(nid)
360 ierr = nf_put_var_real(nid,nvarid,zval)
361 ierr = nf_redef(nid)
362 ierr = nf_def_var(nid,'RUGSREL',nf_float,1,idim2,nvarid)
363 ierr = nf_enddef(nid)
364 ierr = nf_put_var_real(nid,nvarid,rugsrel)
365
366 ierr = nf_redef(nid)
367 ierr = nf_def_var(nid,'TANCIEN',nf_float,1,idim3,nvarid)
368 ierr = nf_enddef(nid)
369 ierr = nf_put_var_real(nid,nvarid,t_ancien)
370
371 ierr = nf_redef(nid)
372 ierr = nf_def_var(nid,'QANCIEN',nf_float,1,idim3,nvarid)
373 ierr = nf_enddef(nid)
374 ierr = nf_put_var_real(nid,nvarid,q_ancien)
375
376 ierr = nf_redef(nid)
377 ierr = nf_def_var(nid,'RUGMER',nf_float,1,idim2,nvarid)
378 ierr = nf_put_att_text(nid,nvarid,'title',28, &
379 'Longueur de rugosite sur mer')
380 ierr = nf_enddef(nid)
381 ierr = nf_put_var_real(nid,nvarid,frugs(1,is_oce))
382
383 ierr = nf_redef(nid)
384 ierr = nf_def_var(nid,'CLWCON',nf_float,1,idim2,nvarid)
385 ierr = nf_put_att_text(nid,nvarid,'title',28,'Eau liquide convective')
386 ierr = nf_enddef(nid)
387 ierr = nf_put_var_real(nid,nvarid,clwcon)
388
389 ierr = nf_redef(nid)
390 ierr = nf_def_var(nid,'RNEBCON',nf_float,1,idim2,nvarid)
391 ierr = nf_put_att_text(nid,nvarid,'title',28,'Nebulosite convective')
392 ierr = nf_enddef(nid)
393 ierr = nf_put_var_real(nid,nvarid,rnebcon)
394
395 ierr = nf_redef(nid)
396 ierr = nf_def_var(nid,'RATQS',nf_float,1,idim2,nvarid)
397 ierr = nf_put_att_text(nid,nvarid,'title',5,'Ratqs')
398 ierr = nf_enddef(nid)
399 ierr = nf_put_var_real(nid,nvarid,ratqs)
400
401 ! run_off_lic_0
402
403 ierr = nf_redef(nid)
404 ierr = nf_def_var(nid,'RUNOFFLIC0',nf_float,1,idim2,nvarid)
405 ierr = nf_put_att_text(nid,nvarid,'title',10,'Runofflic0')
406 ierr = nf_enddef(nid)
407 ierr = nf_put_var_real(nid,nvarid,run_off_lic_0)
408
409
410 ierr = nf_close(nid)
411
412 END SUBROUTINE phyredem
  ViewVC Help
Powered by ViewVC 1.1.21