source: LMDZ4/trunk/libf/phytherm/write_histins.h @ 1080

Last change on this file since 1080 was 814, checked in by Laurent Fairhead, 17 years ago

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