source: trunk/libf/phytitan/write_histins.h @ 21

Last change on this file since 21 was 3, checked in by slebonnois, 15 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 8.6 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/write_histins.h,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
3!
4      IF (ok_instan) THEN
5c
6      ndex2d = 0
7      ndex3d = 0
8      zx_tmp_2d = 0.
9      z4_tmp_2d = 0.
10      zx_tmp_3d = 0.
11      z4_tmp_3d = 0.
12      zx_tmp_fi2d=0.
13      zx_tmp_fi3d=0.
14c
15          zsto = dtime * FLOAT(ecrit_ins)
16          zout = dtime * FLOAT(ecrit_ins)
17         itau_w = itau_phy + itap
18
19c
20c-------------------------------------------------------
21      IF(lev_histday.GE.1) THEN
22c
23ccccccccccccc 2D fields, invariables
24c
25      CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
26      z4_tmp_2d=zx_tmp_2d
27      CALL histwrite(nid_ins,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
28C
29      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
30      z4_tmp_2d=zx_tmp_2d
31      CALL histwrite(nid_ins,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
32c
33ccccccc axe Ls
34      do j=1,jjmp1
35       do i=1,iim
36        z4_tmp_2d(i,j)=zls*180./RPI      ! zls est en radians !!
37       enddo
38      enddo
39      CALL histwrite(nid_ins,"ls",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
40c
41ccccccccccccc 2D fields, variables
42c
43      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d)
44      z4_tmp_2d=zx_tmp_2d
45      CALL histwrite(nid_ins,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
46c
47      DO i = 1, klon
48         zx_tmp_fi2d(i) = paprs(i,1)
49      ENDDO
50      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
51      z4_tmp_2d=zx_tmp_2d
52      CALL histwrite(nid_ins,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
53c
54c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d)
55c     z4_tmp_2d=zx_tmp_2d
56c     CALL histwrite(nid_ins,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
57c
58c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d)
59c     z4_tmp_2d=zx_tmp_2d
60c     CALL histwrite(nid_ins,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
61c
62      ENDIF !lev_histday.GE.1
63c
64c-------------------------------------------------------
65      IF(lev_histday.GE.2) THEN
66c
67ccccccccccccc 3D fields, basics
68c
69      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
70      z4_tmp_3d=zx_tmp_3d
71      CALL histwrite(nid_ins,"temp",itau_w,z4_tmp_3d,
72     .                                   iim*jjmp1*klev,ndex3d)
73c
74      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
75      z4_tmp_3d=zx_tmp_3d
76      CALL histwrite(nid_ins,"pres",itau_w,z4_tmp_3d,
77     .                                   iim*jjmp1*klev,ndex3d)
78c
79      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
80      z4_tmp_3d=zx_tmp_3d
81      CALL histwrite(nid_ins,"geop",itau_w,z4_tmp_3d,
82     .                                   iim*jjmp1*klev,ndex3d)
83c
84      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
85      z4_tmp_3d=zx_tmp_3d
86      CALL histwrite(nid_ins,"vitu",itau_w,z4_tmp_3d,
87     .                                   iim*jjmp1*klev,ndex3d)
88c
89      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
90      z4_tmp_3d=zx_tmp_3d
91      CALL histwrite(nid_ins,"vitv",itau_w,z4_tmp_3d,
92     .                                   iim*jjmp1*klev,ndex3d)
93c
94      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
95      z4_tmp_3d=zx_tmp_3d
96      CALL histwrite(nid_ins,"vitw",itau_w,z4_tmp_3d,
97     .                                   iim*jjmp1*klev,ndex3d)
98c
99      CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
100      z4_tmp_2d=zx_tmp_2d
101      CALL histwrite(nid_ins,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
102c
103      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d)
104      z4_tmp_3d=zx_tmp_3d
105      CALL histwrite(nid_ins,"duvdf",itau_w,z4_tmp_3d,
106     .                                   iim*jjmp1*klev,ndex3d)
107c
108      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d)
109      z4_tmp_3d=zx_tmp_3d
110      CALL histwrite(nid_ins,"dudyn",itau_w,z4_tmp_3d,
111     .                                   iim*jjmp1*klev,ndex3d)
112c
113      ENDIF !lev_histday.GE.2
114c
115c-------------------------------------------------------
116      IF(lev_histday.GE.3) THEN
117c
118cccccccccccccccccc  Tracers
119c
120         if (iflag_trac.eq.1) THEN
121          if (microfi.eq.1) then
122           DO iq=1,nmicro
123       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d)
124       z4_tmp_3d=zx_tmp_3d
125       CALL histwrite(nid_ins,tname(iq),itau_w,z4_tmp_3d,
126     .                                   iim*jjmp1*klev,ndex3d)
127           ENDDO
128          endif
129          if (nmicro.lt.nqmax) then
130           DO iq=nmicro+1,nqmax
131       CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d)
132       z4_tmp_3d=zx_tmp_3d
133       CALL histwrite(nid_ins,tname(iq),itau_w,z4_tmp_3d,
134     .                                   iim*jjmp1*klev,ndex3d)
135           ENDDO
136          endif
137         endif
138c
139cccccccccccccccccc  Radiative transfer
140c
141c 2D
142c
143      CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
144      z4_tmp_2d=zx_tmp_2d
145      CALL histwrite(nid_ins,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
146c
147      CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
148      z4_tmp_2d=zx_tmp_2d
149      CALL histwrite(nid_ins,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
150c
151      CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
152      z4_tmp_2d=zx_tmp_2d
153      CALL histwrite(nid_ins,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
154c
155c 3D
156c
157      zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev)
158      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
159      z4_tmp_3d=zx_tmp_3d
160      CALL histwrite(nid_ins,"SWnet",itau_w,z4_tmp_3d,
161     .                                   iim*jjmp1*klev,ndex3d)
162c
163      zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev)
164      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
165      z4_tmp_3d=zx_tmp_3d
166      CALL histwrite(nid_ins,"LWnet",itau_w,z4_tmp_3d,
167     .                                   iim*jjmp1*klev,ndex3d)
168c
169      ENDIF !lev_histday.GE.3
170c
171c-------------------------------------------------------
172      IF(lev_histday.GE.4) THEN
173c
174      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d)
175      z4_tmp_3d=zx_tmp_3d
176      CALL histwrite(nid_ins,"dtdyn",itau_w,z4_tmp_3d,
177     .                                   iim*jjmp1*klev,ndex3d)
178c
179      CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d)
180      z4_tmp_3d=zx_tmp_3d
181      CALL histwrite(nid_ins,"dtphy",itau_w,z4_tmp_3d,
182     .                                   iim*jjmp1*klev,ndex3d)
183c K/s
184      zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)
185      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
186      z4_tmp_3d=zx_tmp_3d
187      CALL histwrite(nid_ins,"dtvdf",itau_w,z4_tmp_3d,
188     .                                   iim*jjmp1*klev,ndex3d)
189c
190c K/s
191      zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)
192      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
193      z4_tmp_3d=zx_tmp_3d
194      CALL histwrite(nid_ins,"dtajs",itau_w,z4_tmp_3d,
195     .                                   iim*jjmp1*klev,ndex3d)
196c
197c K/s
198      zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)
199      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
200      z4_tmp_3d=zx_tmp_3d
201      CALL histwrite(nid_ins,"dtswr",itau_w,z4_tmp_3d,
202     .                                   iim*jjmp1*klev,ndex3d)
203c
204c K/s     
205      zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)
206      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
207      z4_tmp_3d=zx_tmp_3d
208      CALL histwrite(nid_ins,"dtlwr",itau_w,z4_tmp_3d,
209     .                                   iim*jjmp1*klev,ndex3d)
210c K/s     
211      zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)
212      CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
213      z4_tmp_3d=zx_tmp_3d
214      CALL histwrite(nid_ins,"dtec",itau_w,z4_tmp_3d,
215     .                                   iim*jjmp1*klev,ndex3d)
216c
217c     CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d)
218c     z4_tmp_3d=zx_tmp_3d
219c     CALL histwrite(nid_ins,"dvvdf",itau_w,z4_tmp_3d,
220c    .                                   iim*jjmp1*klev,ndex3d)
221c
222      ENDIF !lev_histday.GE.4
223c
224c-------------------------------------------------------
225      IF(lev_histday.GE.5) THEN
226c
227c
228c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d)
229c      z4_tmp_2d=zx_tmp_2d
230c      CALL histwrite(nid_ins,"taux_",itau_w,
231c    $      z4_tmp_2d,iim*jjmp1,ndex2d) 
232c     
233c      CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d)
234c      z4_tmp_2d=zx_tmp_2d
235c      CALL histwrite(nid_ins,"tauy_",itau_w,
236c    $      z4_tmp_2d,iim*jjmp1,ndex2d)
237c
238c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d)
239c      z4_tmp_2d=zx_tmp_2d
240c     CALL histwrite(nid_ins,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
241c
242c     CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d)
243c      z4_tmp_2d=zx_tmp_2d
244c     CALL histwrite(nid_ins,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d)
245c
246      ENDIF !lev_histday.GE.5
247c-------------------------------------------------------
248c
249      if (ok_sync) then
250        call histsync(nid_ins)
251      endif
252      ENDIF
Note: See TracBrowser for help on using the repository browser.