source: LMDZ4/branches/LMDZ4-dev-20091210/libf/phylmd/ini_histrac.h @ 4104

Last change on this file since 4104 was 1192, checked in by jghattas, 15 years ago

Deplace #ifdef CPP_IOIPSL de ini_histrac et write_histrac vers phytrac.
Sur VARGAS ca ne compile pas avec des #ifdef dans des fichiers include.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.9 KB
Line 
1!
2! $Id $
3!
4  IF (ecrit_tra>0. .AND. config_inca == 'none') THEN
5!$OMP MASTER
6     CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
7     CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
8     CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",klev, presnivs, nvert)
9
10     zsto = pdtphys
11     zout = ecrit_tra
12     CALL histdef(nid_tra, "phis", "Surface geop. height", "-",   &
13          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
14     CALL histdef(nid_tra, "aire", "Grid area", "-",              &
15          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
16
17!TRACEURS
18!----------------
19     DO it = 1,nbtr
20        iiq = niadv(it+2)
21
22! CONCENTRATIONS
23        CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",    &
24             iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
25
26! TD LESSIVAGE
27        IF (lessivage .AND. aerosol(it)) THEN
28           CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
29                "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
30                "ave(X)", zsto,zout)
31        END IF
32
33! TD THERMIQUES
34        IF (iflag_thermals.gt.0) THEN
35           CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),      &
36                "tendance thermique"// ttext(iiq), "?",       &
37                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
38                "ave(X)", zsto,zout)
39        ENDIF
40
41! TD CONVECTION
42        IF (iflag_con.GE.2) THEN
43           CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),   &
44                "tendance convection"// ttext(iiq), "?",   &
45                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
46                "ave(X)", zsto,zout)
47        ENDIF
48
49! TD COUCHE-LIMITE
50        CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),      &
51             "tendance couche limite"// ttext(iiq), "?",   &
52             iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
53             "ave(X)", zsto,zout)
54     ENDDO
55!---------------   
56!
57! VENT (niveau 1)
58     CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",      &
59          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
60          "inst(X)",  zout,zout)     
61     CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",      &
62          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
63          "inst(X)",  zout,zout)
64
65! TEMPERATURE DU SOL
66     CALL histdef(nid_tra, "ftsol1", "temper sol", "-",    &
67          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
68          "inst(X)",  zout,zout)
69     CALL histdef(nid_tra, "ftsol2", "temper sol", "-",    &
70          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
71          "inst(X)",  zout,zout)
72     CALL histdef(nid_tra, "ftsol3", "temper sol", "-",    &
73          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
74          "inst",  zout,zout)
75     CALL histdef(nid_tra, "ftsol4", "temper sol", "-",    &
76          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
77          "inst(X)",  zout,zout)
78
79! NATURE DU SOL
80     CALL histdef(nid_tra, "psrf1", "nature sol", "-",     &
81          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
82          "inst(X)",  zout,zout)
83     CALL histdef(nid_tra, "psrf2", "nature sol", "-",     &
84          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
85          "inst(X)",  zout,zout)
86     CALL histdef(nid_tra, "psrf3", "nature sol", "-",     &
87          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
88          "inst(X)",  zout,zout)
89     CALL histdef(nid_tra, "psrf4", "nature sol", "-",     &
90          iim,jj_nb,nhori, 1,1,1, -99, 32,                 & 
91          "inst(X)",  zout,zout)
92! DIVERS
93     CALL histdef(nid_tra, "pplay", "flux u mont","-",     &
94          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
95          "inst(X)", zout,zout)
96     CALL histdef(nid_tra, "t", "flux u mont","-",         &
97          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
98          "inst(X)", zout,zout)
99     CALL histdef(nid_tra, "mfu", "flux u mont","-",       &
100          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
101          "ave(X)", zsto,zout)
102     CALL histdef(nid_tra, "mfd", "flux u decen","-",      &
103          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
104          "ave(X)", zsto,zout)
105     CALL histdef(nid_tra, "en_u", "flux u mont","-",      &
106          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
107          "ave(X)", zsto,zout)
108     CALL histdef(nid_tra, "en_d", "flux u mont","-",      &
109          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
110          "ave(X)", zsto,zout)
111     CALL histdef(nid_tra, "de_d", "flux u mont","-",      &
112          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
113          "ave(X)", zsto,zout)
114     CALL histdef(nid_tra, "de_u", "flux u decen","-",     &
115          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
116          "ave(X)", zsto,zout)
117     CALL histdef(nid_tra, "coefh", "turbulent coef","-",  &
118          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
119          "ave(X)", zsto,zout)   
120     
121     CALL histend(nid_tra)
122!$OMP END MASTER
123  END IF ! ecrit_tra>0. .AND. config_inca == 'none'
124 
Note: See TracBrowser for help on using the repository browser.