source: LMDZ4/trunk/libf/phy_IPCC_AR4/ini_histrac.h @ 1036

Last change on this file since 1036 was 956, checked in by lmdzadmin, 17 years ago

Nettoyage du controle des parametres physiques. FH

Les parametres cycle_diurne, soil_model, new_oliq, ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad et iflag_con
sont maintenant geres par la physique uniquement.
ecritphy est elimine.
dimphy.F90 et clesphys.h ne sont plus utilises par le code dynamique.
Le test academique obtenu en compilant avec
makegcm -p nophys gcm
fonctionne. FH
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.4 KB
Line 
1!
2! $Header$
3!
4#ifndef INCA
5c$OMP MASTER
6         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
7c
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,
18     .                 nhori, nid_tra)
19         CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
20     .                 klev, presnivs, nvert)
21
22
23
24         zsto = pdtphys
25!         zout = pdtphys * ecrit_tra
26         zout = ecrit_tra
27c
28         CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
29     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
30     .                "once",  zsto,zout)
31c
32         CALL histdef(nid_tra, "aire", "Grid area", "-",
33     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
34     .                "once",  zsto,zout)
35         DO it=1,nqmax
36C champ 2D
37         iq=it+2
38         iiq=niadv(iq)
39         CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga",
40     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
41     .                "ave(X)", zsto,zout)
42         if (lessivage) THEN
43         CALL histdef(nid_tra, "fl"//tnom(iq),"Flux "//ttext(iiq),
44     .              "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
45     .              "ave(X)", zsto,zout)
46         endif
47
48c---Ajout Olivia
49         CALL histdef(nid_tra, "d_tr_th_"//tnom(iq),
50     .                "tendance thermique"// ttext(iiq), "?",
51     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
52     .                "ave(X)", zsto,zout)
53         CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq),
54     .                "tendance convection"// ttext(iiq), "?",
55     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
56     .                "ave(X)", zsto,zout)
57         CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq),
58     .                "tendance couche limite"// ttext(iiq), "?",
59     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
60     .                "ave(X)", zsto,zout)
61c---fin Olivia   
62
63         ENDDO
64
65         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
66     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
67     .                "inst(X)",  zout,zout)
68
69         CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",
70     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
71     .                "inst(X)",  zout,zout)
72         CALL histdef(nid_tra, "psrf1", "nature sol", "-",
73     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
74     .                "inst(X)",  zout,zout)
75         CALL histdef(nid_tra, "psrf2", "nature sol", "-",
76     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
77     .                "inst(X)",  zout,zout)
78         CALL histdef(nid_tra, "psrf3", "nature sol", "-",
79     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
80     .                "inst(X)",  zout,zout)
81         CALL histdef(nid_tra, "psrf4", "nature sol", "-",
82     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
83     .                "inst(X)",  zout,zout)
84         CALL histdef(nid_tra, "ftsol1", "temper sol", "-",
85     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
86     .                "inst(X)",  zout,zout)
87         CALL histdef(nid_tra, "ftsol2", "temper sol", "-",
88     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
89     .                "inst(X)",  zout,zout)
90         CALL histdef(nid_tra, "ftsol3", "temper sol", "-",
91     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
92     .                "inst",  zout,zout)
93         CALL histdef(nid_tra, "ftsol4", "temper sol", "-",
94     .                iim,jj_nb,nhori, 1,1,1, -99, 32,
95     .                "inst(X)",  zout,zout)
96         CALL histdef(nid_tra, "pplay", "flux u mont","-",
97     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
98     .                "inst(X)", zout,zout)
99         CALL histdef(nid_tra, "t", "flux u mont","-",
100     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
101     .                "inst(X)", zout,zout)
102         CALL histdef(nid_tra, "mfu", "flux u mont","-",
103     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
104     .                "ave(X)", zsto,zout)
105         CALL histdef(nid_tra, "mfd", "flux u decen","-",
106     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
107     .                "ave(X)", zsto,zout)
108         CALL histdef(nid_tra, "en_u", "flux u mont","-",
109     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
110     .                "ave(X)", zsto,zout)
111         CALL histdef(nid_tra, "en_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_d", "flux u mont","-",
115     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
116     .                "ave(X)", zsto,zout)
117         CALL histdef(nid_tra, "de_u", "flux u decen","-",
118     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
119     .                "ave(X)", zsto,zout)
120         CALL histdef(nid_tra, "coefh", "turbulent coef","-",
121     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
122     .                "ave(X)", zsto,zout)
123
124c
125         CALL histend(nid_tra)
126         ndex2d = 0
127         ndex3d = 0
128         ndex = 0
129c$OMP END MASTER
130#endif
Note: See TracBrowser for help on using the repository browser.