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

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

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