source: LMDZ4/trunk/libf/phylmd/write_histrac.h @ 959

Last change on this file since 959 was 959, checked in by lsce, 16 years ago
  • Ajoute du parametre config_inca dans conf_gcm.F config_inca='none'(sans INCA, par defaut) config_inca='chem'(avec INCA config chemie) config_inca='aero'(avec INCA config aerosol)
  • Menage parmis les cles CPP INCA
  • Enleve le calcul d'omega dans calfis.F et active le calcul correspondant dans advtrac.F(avant uniquement pour INCA).
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.8 KB
Line 
1!
2! $Header$
3!
4
5      IF (config_inca == 'none') THEN
6      ndex = 0
7      ndex2d = 0
8      ndex3d = 0
9c
10      itau_w = itau_phy + nstep
11
12cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
13      CALL histwrite_phy(nid_tra,"phis",itau_w,pphis)
14C
15cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,airephy,zx_tmp_2d)     
16      CALL histwrite_phy(nid_tra,"aire",itau_w,airephy)
17
18      DO it=1,nqmax
19C champs 2D
20
21
22cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
23       CALL histwrite_phy(nid_tra,tnom(it+2),itau_w,tr_seri(:,:,it))
24       if (lessivage) THEN
25cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,flestottr(1,1,it),zx_tmp_3d)
26       CALL histwrite_phy(nid_tra,"fl"//tnom(it+2),itau_w,
27     .                                   flestottr(:,:,it))
28      endif
29     
30c----Olivia
31cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_th(1,1,it),zx_tmp_3d)
32       CALL histwrite_phy(nid_tra,"d_tr_th_"//tnom(it+2),itau_w,
33     .                                           d_tr_th(:,:,it))
34cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cv(1,1,it),zx_tmp_3d)
35       CALL histwrite_phy(nid_tra,"d_tr_cv_"//tnom(it+2),itau_w,
36     .                                           d_tr_cv(:,:,it))
37cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cl(1,1,it),zx_tmp_3d)
38       CALL histwrite_phy(nid_tra,"d_tr_cl_"//tnom(it+2),itau_w,
39     .                                           d_tr_cl(:,:,it))
40c---fin Olivia     
41     
42      ENDDO
43
44
45C abder
46cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,yu1,zx_tmp_2d)
47         CALL histwrite_phy(nid_tra,"pyu1",itau_w,yu1)
48
49cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,yv1,zx_tmp_2d)
50         CALL histwrite_phy(nid_tra,"pyv1",itau_w,yv1)
51
52cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol1,zx_tmp_2d)
53         CALL histwrite_phy(nid_tra,"ftsol1",itau_w,pftsol1)
54
55cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol2,zx_tmp_2d)
56         CALL histwrite_phy(nid_tra,"ftsol2",itau_w,pftsol2)
57
58cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol3,zx_tmp_2d)
59         CALL histwrite_phy(nid_tra,"ftsol3",itau_w,pftsol3)
60
61cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol4,zx_tmp_2d)
62         CALL histwrite_phy(nid_tra,"ftsol4",itau_w,pftsol4)
63
64cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf1,zx_tmp_2d)
65         CALL histwrite_phy(nid_tra,"psrf1",itau_w,ppsrf1)
66
67cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf2,zx_tmp_2d)
68         CALL histwrite_phy(nid_tra,"psrf2",itau_w,ppsrf2)
69
70cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf3,zx_tmp_2d)
71         CALL histwrite_phy(nid_tra,"psrf3",itau_w,ppsrf3)
72
73cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf4,zx_tmp_2d)
74         CALL histwrite_phy(nid_tra,"psrf4",itau_w,ppsrf4)
75cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pplay,zx_tmp_3d)
76        CALL histwrite_phy(nid_tra,"pplay",itau_w,pplay)
77
78cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,t_seri,zx_tmp_3d)
79        CALL histwrite_phy(nid_tra,"t",itau_w,t_seri)
80cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfu,zx_tmp_3d)
81        CALL histwrite_phy(nid_tra,"mfu",itau_w,pmfu)
82cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfd,zx_tmp_3d)
83        CALL histwrite_phy(nid_tra,"mfd",itau_w,pmfd)
84cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_u,zx_tmp_3d)
85        CALL histwrite_phy(nid_tra,"en_u",itau_w,pen_u)
86cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_d,zx_tmp_3d)
87        CALL histwrite_phy(nid_tra,"en_d",itau_w,pen_d)
88cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_d,zx_tmp_3d)
89        CALL histwrite_phy(nid_tra,"de_d",itau_w,pde_d)
90cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_u,zx_tmp_3d)
91        CALL histwrite_phy(nid_tra,"de_u",itau_w,pde_u)
92cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,coefh,zx_tmp_3d)
93        CALL histwrite_phy(nid_tra,"coefh",itau_w,coefh)
94
95
96c abder
97
98      if (ok_sync) then
99c$OMP MASTER
100        call histsync(nid_tra)
101c$OMP END MASTER
102       endif
103
104       END IF
105
106
107
Note: See TracBrowser for help on using the repository browser.