source: LMDZ4/trunk/libf/phylmd/ini_histrac.h @ 1097

Last change on this file since 1097 was 1030, checked in by lmdzadmin, 16 years ago

Nettoyage de commentaires + correction d'un bug dans le cas sans convection
FH/IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.0 KB
RevLine 
[524]1!
2! $Header$
3!
[959]4      IF (config_inca == 'none') THEN
[766]5c$OMP MASTER
[524]6         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
7c
[766]8         CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,
[524]9     .                 nhori, nid_tra)
10         CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
11     .                 klev, presnivs, nvert)
12
13
14
15         zsto = pdtphys
[956]16         zout = ecrit_tra
[524]17c
18         CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
[776]19     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[524]20     .                "once",  zsto,zout)
21c
22         CALL histdef(nid_tra, "aire", "Grid area", "-",
[776]23     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[524]24     .                "once",  zsto,zout)
25         DO it=1,nqmax
26C champ 2D
27         iq=it+2
28         iiq=niadv(iq)
29         CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga",
[776]30     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]31     .                "ave(X)", zsto,zout)
32         if (lessivage) THEN
33         CALL histdef(nid_tra, "fl"//tnom(iq),"Flux "//ttext(iiq),
[776]34     .              "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]35     .              "ave(X)", zsto,zout)
36         endif
[549]37
38c---Ajout Olivia
39         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq),
40     .                "tendance thermique"// ttext(iiq), "?",
[776]41     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]42     .                "ave(X)", zsto,zout)
[973]43c
44         if(iflag_con.GE.2) then
[549]45         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq),
46     .                "tendance convection"// ttext(iiq), "?",
[776]47     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]48     .                "ave(X)", zsto,zout)
[973]49             endif !(iflag_con.GE.2) then
[549]50         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq),
51     .                "tendance couche limite"// ttext(iiq), "?",
[776]52     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[549]53     .                "ave(X)", zsto,zout)
54c---fin Olivia   
55
[524]56         ENDDO
[549]57
[524]58         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
[776]59     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[766]60     .                "inst(X)",  zout,zout)
[524]61
62         CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
[776]63     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]64     .                "inst(X)",  zout,zout)
[524]65         CALL histdef(nid_tra, "psrf1", "nature sol", "-",
[776]66     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]67     .                "inst(X)",  zout,zout)
[524]68         CALL histdef(nid_tra, "psrf2", "nature sol", "-",
[776]69     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]70     .                "inst(X)",  zout,zout)
[524]71         CALL histdef(nid_tra, "psrf3", "nature sol", "-",
[776]72     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]73     .                "inst(X)",  zout,zout)
[524]74         CALL histdef(nid_tra, "psrf4", "nature sol", "-",
[776]75     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]76     .                "inst(X)",  zout,zout)
[524]77         CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
[776]78     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]79     .                "inst(X)",  zout,zout)
[524]80         CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
[776]81     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]82     .                "inst(X)",  zout,zout)
[524]83         CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
[776]84     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[766]85     .                "inst",  zout,zout)
[524]86         CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
[776]87     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
[657]88     .                "inst(X)",  zout,zout)
[524]89         CALL histdef(nid_tra, "pplay", "flux u mont","-",
[776]90     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[657]91     .                "inst(X)", zout,zout)
[524]92         CALL histdef(nid_tra, "t", "flux u mont","-",
[776]93     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[657]94     .                "inst(X)", zout,zout)
[524]95         CALL histdef(nid_tra, "mfu", "flux u mont","-",
[776]96     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]97     .                "ave(X)", zsto,zout)
98         CALL histdef(nid_tra, "mfd", "flux u decen","-",
[776]99     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]100     .                "ave(X)", zsto,zout)
101         CALL histdef(nid_tra, "en_u", "flux u mont","-",
[776]102     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]103     .                "ave(X)", zsto,zout)
104         CALL histdef(nid_tra, "en_d", "flux u mont","-",
[776]105     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]106     .                "ave(X)", zsto,zout)
107         CALL histdef(nid_tra, "de_d", "flux u mont","-",
[776]108     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]109     .                "ave(X)", zsto,zout)
110         CALL histdef(nid_tra, "de_u", "flux u decen","-",
[776]111     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]112     .                "ave(X)", zsto,zout)
113         CALL histdef(nid_tra, "coefh", "turbulent coef","-",
[776]114     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
[524]115     .                "ave(X)", zsto,zout)
116
117c
118         CALL histend(nid_tra)
119         ndex2d = 0
120         ndex3d = 0
121         ndex = 0
[766]122c$OMP END MASTER
[959]123      END IF
Note: See TracBrowser for help on using the repository browser.