source: LMDZ4/trunk/libf/phytherm/write_histrac.h @ 1068

Last change on this file since 1068 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: 3.7 KB
Line 
1!
2! $Header$
3!
4#ifndef INCA
5      ndex = 0
6      ndex2d = 0
7      ndex3d = 0
8c
9      itau_w = itau_phy + nstep
10
11cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
12      CALL histwrite_phy(nid_tra,"phis",itau_w,pphis)
13C
14cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,airephy,zx_tmp_2d)     
15      CALL histwrite_phy(nid_tra,"aire",itau_w,airephy)
16
17      DO it=1,nqmax
18C champs 2D
19
20
21cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
22       CALL histwrite_phy(nid_tra,tnom(it+2),itau_w,tr_seri(:,:,it))
23       if (lessivage) THEN
24cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,flestottr(1,1,it),zx_tmp_3d)
25       CALL histwrite_phy(nid_tra,"fl"//tnom(it+2),itau_w,
26     .                                   flestottr(:,:,it))
27      endif
28     
29c----Olivia
30cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_th(1,1,it),zx_tmp_3d)
31       CALL histwrite_phy(nid_tra,"d_tr_th_"//tnom(it+2),itau_w,
32     .                                           d_tr_th(:,:,it))
33cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cv(1,1,it),zx_tmp_3d)
34       CALL histwrite_phy(nid_tra,"d_tr_cv_"//tnom(it+2),itau_w,
35     .                                           d_tr_cv(:,:,it))
36cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cl(1,1,it),zx_tmp_3d)
37       CALL histwrite_phy(nid_tra,"d_tr_cl_"//tnom(it+2),itau_w,
38     .                                           d_tr_cl(:,:,it))
39c---fin Olivia     
40     
41      ENDDO
42
43
44C abder
45cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,yu1,zx_tmp_2d)
46         CALL histwrite_phy(nid_tra,"pyu1",itau_w,yu1)
47
48cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,yv1,zx_tmp_2d)
49         CALL histwrite_phy(nid_tra,"pyv1",itau_w,yv1)
50
51cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol1,zx_tmp_2d)
52         CALL histwrite_phy(nid_tra,"ftsol1",itau_w,pftsol1)
53
54cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol2,zx_tmp_2d)
55         CALL histwrite_phy(nid_tra,"ftsol2",itau_w,pftsol2)
56
57cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol3,zx_tmp_2d)
58         CALL histwrite_phy(nid_tra,"ftsol3",itau_w,pftsol3)
59
60cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol4,zx_tmp_2d)
61         CALL histwrite_phy(nid_tra,"ftsol4",itau_w,pftsol4)
62
63cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf1,zx_tmp_2d)
64         CALL histwrite_phy(nid_tra,"psrf1",itau_w,ppsrf1)
65
66cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf2,zx_tmp_2d)
67         CALL histwrite_phy(nid_tra,"psrf2",itau_w,ppsrf2)
68
69cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf3,zx_tmp_2d)
70         CALL histwrite_phy(nid_tra,"psrf3",itau_w,ppsrf3)
71
72cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf4,zx_tmp_2d)
73         CALL histwrite_phy(nid_tra,"psrf4",itau_w,ppsrf4)
74cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pplay,zx_tmp_3d)
75        CALL histwrite_phy(nid_tra,"pplay",itau_w,pplay)
76
77cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,t_seri,zx_tmp_3d)
78        CALL histwrite_phy(nid_tra,"t",itau_w,t_seri)
79cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfu,zx_tmp_3d)
80        CALL histwrite_phy(nid_tra,"mfu",itau_w,pmfu)
81cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfd,zx_tmp_3d)
82        CALL histwrite_phy(nid_tra,"mfd",itau_w,pmfd)
83cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_u,zx_tmp_3d)
84        CALL histwrite_phy(nid_tra,"en_u",itau_w,pen_u)
85cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_d,zx_tmp_3d)
86        CALL histwrite_phy(nid_tra,"en_d",itau_w,pen_d)
87cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_d,zx_tmp_3d)
88        CALL histwrite_phy(nid_tra,"de_d",itau_w,pde_d)
89cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_u,zx_tmp_3d)
90        CALL histwrite_phy(nid_tra,"de_u",itau_w,pde_u)
91cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,coefh,zx_tmp_3d)
92        CALL histwrite_phy(nid_tra,"coefh",itau_w,coefh)
93
94
95c abder
96
97      if (ok_sync) then
98c$OMP MASTER
99        call histsync(nid_tra)
100c$OMP END MASTER
101       endif
102
103#endif
104
105
106
Note: See TracBrowser for help on using the repository browser.