source: LMDZ4/tags/LMDZ4_V2/libf/phylmd/write_histins.h @ 741

Last change on this file since 741 was 741, checked in by (none), 17 years ago

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