source: LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h @ 1126

Last change on this file since 1126 was 743, checked in by Laurent Fairhead, 18 years ago

Nettoyage sur la version YM
LF

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