source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histins.h @ 506

Last change on this file since 506 was 506, checked in by lmdzadmin, 20 years ago

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 7.9 KB
Line 
1      IF (ok_instan) THEN
2c
3      ndex2d = 0
4      ndex3d = 0
5c
6c Champs 2D:
7c
8         zsto = dtime * ecrit_ins
9         zout = dtime * ecrit_ins
10         itau_w = itau_phy + itap
11
12         i = NINT(zout/zsto)
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         i = NINT(zout/zsto)
17      CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
18      CALL histwrite(nid_ins,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
19
20      DO i = 1, klon
21         zx_tmp_fi2d(i) = paprs(i,1)
22      ENDDO
23      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
24      CALL histwrite(nid_ins,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
25c
26      DO i = 1, klon
27         zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
28      ENDDO
29      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
30      CALL histwrite(nid_ins,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
31c
32      DO i = 1, klon
33         zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
34      ENDDO
35      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
36      CALL histwrite(nid_ins,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
37c
38      DO i = 1, klon
39         zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
40      ENDDO
41      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
42      CALL histwrite(nid_ins,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
43
44      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
45      CALL histwrite(nid_ins,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
46cccIM
47      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m, zx_tmp_2d)
48      CALL histwrite(nid_ins,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
49c
50      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m, zx_tmp_2d)
51      CALL histwrite(nid_ins,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
52c
53      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m, zx_tmp_2d)
54      CALL histwrite(nid_ins,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
55C
56      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m, zx_tmp_2d)
57      CALL histwrite(nid_ins,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
58c
59      CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
60      CALL histwrite(nid_ins,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
61
62      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
63      CALL histwrite(nid_ins,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
64c
65      CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
66      CALL histwrite(nid_ins,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
67c
68      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
69      CALL histwrite(nid_ins,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
70c
71      CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
72      CALL histwrite(nid_ins,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
73c
74      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
75      CALL histwrite(nid_ins,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
76c
77      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
78      CALL histwrite(nid_ins,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
79c
80      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
81      CALL histwrite(nid_ins,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
82     .                ndex2d)
83c
84      CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
85      CALL histwrite(nid_ins,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
86c
87      zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
88c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d)
89      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
90      CALL histwrite(nid_ins,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
91c
92      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
93      CALL histwrite(nid_ins,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
94c
95      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_oce),zx_tmp_2d)
96      CALL histwrite(nid_ins,"dtsvdfo",itau_w,zx_tmp_2d,iim*jjmp1,
97     .               ndex2d)
98c
99      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_ter),zx_tmp_2d)
100      CALL histwrite(nid_ins,"dtsvdft",itau_w,zx_tmp_2d,iim*jjmp1,
101     .               ndex2d)
102c
103      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_lic),zx_tmp_2d)
104      CALL histwrite(nid_ins,"dtsvdfg",itau_w,zx_tmp_2d,iim*jjmp1,
105     .               ndex2d)
106c
107      CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_sic),zx_tmp_2d)
108      CALL histwrite(nid_ins,"dtsvdfi",itau_w,zx_tmp_2d,iim*jjmp1,
109     .               ndex2d)
110
111      DO nsrf = 1, nbsrf
112CXXX
113        zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
114        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
115        CALL histwrite(nid_ins,"pourc_"//clnsurf(nsrf),itau_w,
116     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
117C
118        zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
119        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
120        CALL histwrite(nid_ins,"sens_"//clnsurf(nsrf),itau_w,
121     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
122C
123        zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
124        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
125        CALL histwrite(nid_ins,"lat_"//clnsurf(nsrf),itau_w,
126     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
127C
128        zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
129        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
130        CALL histwrite(nid_ins,"tsol_"//clnsurf(nsrf),itau_w,
131     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
132C
133        zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
134        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
135        CALL histwrite(nid_ins,"taux_"//clnsurf(nsrf),itau_w,
136     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
137C     
138        zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
139        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
140        CALL histwrite(nid_ins,"tauy_"//clnsurf(nsrf),itau_w,
141     $      zx_tmp_2d,iim*jjmp1,ndex2d)
142C
143        zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
144        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
145        CALL histwrite(nid_ins,"rugs_"//clnsurf(nsrf),itau_w,
146     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
147C
148        zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
149        CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
150        CALL histwrite(nid_ins,"albe_"//clnsurf(nsrf),itau_w,
151     $      zx_tmp_2d,iim*jjmp1,ndex2d) 
152C
153      END DO 
154      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d)
155      CALL histwrite(nid_ins,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
156      CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d)
157      CALL histwrite(nid_ins,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
158c
159      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
160      CALL histwrite(nid_ins,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
161     .               ndex2d)
162c
163      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxrugs,zx_tmp_2d)
164      CALL histwrite(nid_ins,"rugs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
165c
166c Champs 3D:
167c
168      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
169      CALL histwrite(nid_ins,"temp",itau_w,zx_tmp_3d,
170     .                                   iim*jjmp1*klev,ndex3d)
171c
172      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
173      CALL histwrite(nid_ins,"vitu",itau_w,zx_tmp_3d,
174     .                                   iim*jjmp1*klev,ndex3d)
175c
176      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
177      CALL histwrite(nid_ins,"vitv",itau_w,zx_tmp_3d,
178     .                                   iim*jjmp1*klev,ndex3d)
179c
180      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
181      CALL histwrite(nid_ins,"geop",itau_w,zx_tmp_3d,
182     .                                   iim*jjmp1*klev,ndex3d)
183c
184      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
185      CALL histwrite(nid_ins,"pres",itau_w,zx_tmp_3d,
186     .                                   iim*jjmp1*klev,ndex3d)
187c
188      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d)
189      CALL histwrite(nid_ins,"dtvdf",itau_w,zx_tmp_3d,
190     .                                   iim*jjmp1*klev,ndex3d)
191c
192      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d)
193      CALL histwrite(nid_ins,"dqvdf",itau_w,zx_tmp_3d,
194     .                                   iim*jjmp1*klev,ndex3d)
195
196c
197      if (ok_sync) then
198        call histsync(nid_ins)
199      endif
200      ENDIF
Note: See TracBrowser for help on using the repository browser.