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