My Project
 All Classes Files Functions Variables Macros
ini_histmthNMC.h
Go to the documentation of this file.
1 !
2 ! $Id: ini_histmthNMC.h 1403 2010-07-01 09:02:53Z fairhead $
3 !
4 c$OMP MASTER
5 c
6  IF (ok_histNMC(1)) THEN
7 c
8  zout = freq_outNMC(1)
9 c
10  idayref = day_ref
11  CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
12 c
13 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14 cym DO i = 1, iim
15 cym zx_lon(i,1) = rlon(i+1)
16 cym zx_lon(i,jjmp1) = rlon(i+1)
17 cym ENDDO
18  DO ll=1,klev
19  znivsig(ll)=REAL(ll)
20  ENDDO
21 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22 cym CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
23 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
24 cym . nhori, nid_mthnmc)
25 
26  CALL histbeg_phy("histmthNMC",itau_phy, zjulian, dtime,
27  . nhori, nid_mthnmc)
28 c
29  CALL histvert(nid_mthnmc, "plev", "pressure", "Pa",
30  . nlevSTD, rlevSTD, nvert,"down")
31 c
32 cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
33 cIM le time_counter et les bounds
34 cIM
35 ccc Champs 3D interpolles sur des niveaux de pression du NMC
36 ccc
37 c
38 c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
39 c sur les niveaux du NMC et on somme & moyenne
40 c toutes les freq_moyNMC secondes par des routines undefSTD et
41 c moy_undefSTD pour eliminer les valeurs "undef"
42 c de la moyenne mensuelle
43 c ======> le "inst(X)" ci-dessous est par consequence factice !
44 c
45 c
46  CALL histdef(nid_mthnmc, "tnondef",
47  . "Valeurs non-definies","-",
48  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
49  . "ave(X)", zout,zout)
50 c
51  CALL histdef(nid_mthnmc, "ta",
52  . "Air temperature","K",
53  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
54  . "ave(X)", zout,zout)
55 c
56  CALL histdef(nid_mthnmc, "zg",
57  . "Geopotential height", "m",
58  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
59  . "ave(X)", zout,zout)
60 c
61  CALL histdef(nid_mthnmc, "hus",
62  . "Specific humidity","1",
63  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
64  . "ave(X)", zout,zout)
65 c
66  CALL histdef(nid_mthnmc, "hur",
67  . "Relative humidity", "%",
68  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
69  . "ave(X)", zout,zout)
70 c
71  CALL histdef(nid_mthnmc, "ua",
72  . "Eastward wind","m s-1",
73  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
74  . "ave(X)", zout,zout)
75 c
76  CALL histdef(nid_mthnmc, "va",
77  . "Northward wind","m s-1",
78  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
79  . "ave(X)", zout,zout)
80 c
81  CALL histdef(nid_mthnmc, "wap",
82  . "Lagrangian tendency of air pressure","Pa s-1",
83  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
84  . "ave(X)", zout,zout)
85 c
86  CALL histdef(nid_mthnmc, "psbg",
87  . "Pressure sfce below ground","%",
88  . iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
89  . "ave(X)", zout,zout)
90 c
91  CALL histdef(nid_mthnmc, "uv",
92  . "uv ",
93  . "m2/s2",iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
94  . "ave(X)", zout,zout)
95 c
96  CALL histdef(nid_mthnmc, "vq",
97  . "vq ",
98  . "m/s * (kg/kg)",iim,jj_nb,nhori,
99  . nlevSTD,1,nlevSTD, nvert, 32,
100  . "ave(X)", zout,zout)
101 c
102  CALL histdef(nid_mthnmc, "vT",
103  . "vT ",
104  . "mK/s",iim,jj_nb,nhori,
105  . nlevSTD,1,nlevSTD, nvert, 32,
106  . "ave(X)", zout,zout)
107 c
108  CALL histdef(nid_mthnmc, "wq",
109  . "wq ",
110  . "(Pa/s)*(kg/kg)",iim,jj_nb,nhori,
111  . nlevSTD,1,nlevSTD, nvert, 32,
112  . "ave(X)", zout,zout)
113 c
114  CALL histdef(nid_mthnmc, "vphi",
115  . "vphi ",
116  . "m2/s",iim,jj_nb,nhori,
117  . nlevSTD,1,nlevSTD, nvert, 32,
118  . "ave(X)", zout,zout)
119 c
120  CALL histdef(nid_mthnmc, "wT",
121  . "wT ",
122  . "K*Pa/s",iim,jj_nb,nhori,
123  . nlevSTD,1,nlevSTD, nvert, 32,
124  . "ave(X)", zout,zout)
125 c
126  CALL histdef(nid_mthnmc, "uxu",
127  . "u2 ",
128  . "m2/s2",iim,jj_nb,nhori,
129  . nlevSTD,1,nlevSTD, nvert, 32,
130  . "ave(X)", zout,zout)
131 c
132  CALL histdef(nid_mthnmc, "vxv",
133  . "v2 ",
134  . "m2/s2",iim,jj_nb,nhori,
135  . nlevSTD,1,nlevSTD, nvert, 32,
136  . "ave(X)", zout,zout)
137 c
138  CALL histdef(nid_mthnmc, "TxT",
139  . "T2 ",
140  . "K2",iim,jj_nb,nhori,
141  . nlevSTD,1,nlevSTD, nvert, 32,
142  . "ave(X)", zout,zout)
143 c
144  CALL histdef(nid_mthnmc, "tro3",
145  . "Ozone mole fraction",
146  . "1e-9",iim,jj_nb,nhori,
147  . nlevSTD,1,nlevSTD, nvert, 32,
148  . "ave(X)", zout,zout)
149 c
150  if (read_climoz == 2) THEN
151  CALL histdef(nid_mthnmc, "tro3_daylight",
152  . "Daylight ozone mole fraction",
153  . "1e-9",iim,jj_nb,nhori,
154  . nlevSTD,1,nlevSTD, nvert, 32,
155  . "ave(X)", zout,zout)
156  endif
157 c
158  CALL histend(nid_mthnmc)
159 c
160  ENDIF !(ok_histNMC(1)) THEN
161 c
162 c$OMP END MASTER