source: LMDZ4/trunk/libf/phylmd/write_histins.h @ 719

Last change on this file since 719 was 684, checked in by lmdzadmin, 19 years ago

Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def +
Lecture freq. ecriture en nombre de jours
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 9.7 KB
Line 
1!
2! $Header$
3!
4      IF (ok_instan) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8c
9c Champs 2D:
10c
11      itau_w = itau_phy + itap
12
13      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
14      CALL histwrite(nid_ins,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
15c
16      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
17      CALL histwrite(nid_ins,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
18
19      DO i = 1, klon
20         zx_tmp_fi2d(i) = paprs(i,1)
21      ENDDO
22      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
23      CALL histwrite(nid_ins,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
24c
25      DO i = 1, klon
26         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
27      ENDDO
28      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
29      CALL histwrite(nid_ins,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
30c
31      DO i = 1, klon
32         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
33      ENDDO
34      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
35      CALL histwrite(nid_ins,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
36c
37      DO i = 1, klon
38         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
39      ENDDO
40      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
41      CALL histwrite(nid_ins,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
42
43      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
44      CALL histwrite(nid_ins,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
45cccIM
46      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m, zx_tmp_2d)
47      CALL histwrite(nid_ins,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
48c
49      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m, zx_tmp_2d)
50      CALL histwrite(nid_ins,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
51c
52      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m, zx_tmp_2d)
53      CALL histwrite(nid_ins,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
54C
55      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m, zx_tmp_2d)
56      CALL histwrite(nid_ins,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
57c
58      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
59      CALL histwrite(nid_ins,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
60
61      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
62      CALL histwrite(nid_ins,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
63c
64      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
65      CALL histwrite(nid_ins,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
66c
67      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
68      CALL histwrite(nid_ins,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
69c
70      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
71      CALL histwrite(nid_ins,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
72c
73      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
74      CALL histwrite(nid_ins,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
75c
76      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
77      CALL histwrite(nid_ins,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
78c
79      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
80      CALL histwrite(nid_ins,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
81     .                ndex2d)
82c
83      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
84      CALL histwrite(nid_ins,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
85c
86      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
87c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
88      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
89      CALL histwrite(nid_ins,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
90c
91      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
92      CALL histwrite(nid_ins,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
93c
94      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_oce),zx_tmp_2d)
95      CALL histwrite(nid_ins,"dtsvdfo",itau_w,zx_tmp_2d,iim*jjmp1,
96     .               ndex2d)
97c
98      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_ter),zx_tmp_2d)
99      CALL histwrite(nid_ins,"dtsvdft",itau_w,zx_tmp_2d,iim*jjmp1,
100     .               ndex2d)
101c
102      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_lic),zx_tmp_2d)
103      CALL histwrite(nid_ins,"dtsvdfg",itau_w,zx_tmp_2d,iim*jjmp1,
104     .               ndex2d)
105c
106      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_sic),zx_tmp_2d)
107      CALL histwrite(nid_ins,"dtsvdfi",itau_w,zx_tmp_2d,iim*jjmp1,
108     .               ndex2d)
109
110      DO nsrf = 1, nbsrf
111CXXX
112        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
113        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
114        CALL histwrite(nid_ins,"pourc_"//clnsurf(nsrf),itau_w,
115     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
116C
117        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
118        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
119        CALL histwrite(nid_ins,"fract_"//clnsurf(nsrf),itau_w,
120     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
121C
122        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
123        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
124        CALL histwrite(nid_ins,"sens_"//clnsurf(nsrf),itau_w,
125     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
126C
127        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
128        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
129        CALL histwrite(nid_ins,"lat_"//clnsurf(nsrf),itau_w,
130     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
131C
132        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
133        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
134        CALL histwrite(nid_ins,"tsol_"//clnsurf(nsrf),itau_w,
135     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
136C
137        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
138        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
139        CALL histwrite(nid_ins,"taux_"//clnsurf(nsrf),itau_w,
140     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
141C     
142        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
143        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
144        CALL histwrite(nid_ins,"tauy_"//clnsurf(nsrf),itau_w,
145     $      zx_tmp_2d,iim*jjmp1,ndex2d)
146C
147        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
148        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
149        CALL histwrite(nid_ins,"rugs_"//clnsurf(nsrf),itau_w,
150     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
151C
152        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
153        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
154        CALL histwrite(nid_ins,"albe_"//clnsurf(nsrf),itau_w,
155     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
156C
157      END DO 
158      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
159      CALL histwrite(nid_ins,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
160      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
161      CALL histwrite(nid_ins,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
162c
163cIM: 140404   CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
164c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d)
165c     CALL histwrite(nid_ins,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
166c    .               ndex2d)
167c
168      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxrugs,zx_tmp_2d)
169      CALL histwrite(nid_ins,"rugs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
170c
171cIM cf. AM 081204 BEG
172c
173cHBTM2
174c
175c
176      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d)
177      CALL histwrite(nid_ins,"s_pblh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
178c
179      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d)
180      CALL histwrite(nid_ins,"s_pblt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
181c
182      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d)
183      CALL histwrite(nid_ins,"s_lcl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
184c
185      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d)
186      CALL histwrite(nid_ins,"s_capCL",itau_w,zx_tmp_2d,iim*jjmp1,
187     .               ndex2d)
188c
189      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d)
190      CALL histwrite(nid_ins,"s_oliqCL",itau_w,zx_tmp_2d,iim*jjmp1,
191     .               ndex2d)
192c
193      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d)
194      CALL histwrite(nid_ins,"s_cteiCL",itau_w,zx_tmp_2d,iim*jjmp1,
195     .               ndex2d)
196c
197      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d)
198      CALL histwrite(nid_ins,"s_therm",itau_w,zx_tmp_2d,iim*jjmp1,
199     .               ndex2d)
200c
201      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d)
202      CALL histwrite(nid_ins,"s_trmb1",itau_w,zx_tmp_2d,iim*jjmp1,
203     .               ndex2d)
204c
205      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d)
206      CALL histwrite(nid_ins,"s_trmb2",itau_w,zx_tmp_2d,iim*jjmp1,
207     .               ndex2d)
208c
209      CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d)
210      CALL histwrite(nid_ins,"s_trmb3",itau_w,zx_tmp_2d,iim*jjmp1,
211     .               ndex2d)
212c
213cIM cf. AM 081204 END
214c
215c Champs 3D:
216c
217      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
218      CALL histwrite(nid_ins,"temp",itau_w,zx_tmp_3d,
219     .                                   iim*jjmp1*klev,ndex3d)
220c
221      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
222      CALL histwrite(nid_ins,"vitu",itau_w,zx_tmp_3d,
223     .                                   iim*jjmp1*klev,ndex3d)
224c
225      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
226      CALL histwrite(nid_ins,"vitv",itau_w,zx_tmp_3d,
227     .                                   iim*jjmp1*klev,ndex3d)
228c
229      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
230      CALL histwrite(nid_ins,"geop",itau_w,zx_tmp_3d,
231     .                                   iim*jjmp1*klev,ndex3d)
232c
233      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
234      CALL histwrite(nid_ins,"pres",itau_w,zx_tmp_3d,
235     .                                   iim*jjmp1*klev,ndex3d)
236c
237      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
238      CALL histwrite(nid_ins,"dtvdf",itau_w,zx_tmp_3d,
239     .                                   iim*jjmp1*klev,ndex3d)
240c
241      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
242      CALL histwrite(nid_ins,"dqvdf",itau_w,zx_tmp_3d,
243     .                                   iim*jjmp1*klev,ndex3d)
244
245c
246      if (ok_sync) then
247        call histsync(nid_ins)
248      endif
249      ENDIF
Note: See TracBrowser for help on using the repository browser.