Context Navigation

source: LMDZ4/branches/pre_V3/libf/phylmd/ini_histmthNMC.h @ 4096

Last change on this file since 4096 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 6.1 KB

Line
1	!
2	! $Header$
3	!
4	IF (ok_mensuel) THEN
5	c
6	zstophy = dtime
7	zstohf = ecrit_hf
8	zstomth = ecrit_mth
9	zout = ecrit_mth
10	c
11	idayref = day_ref
12	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
13	c
14	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
15	DO i = 1, iim
16	zx_lon(i,1) = rlon(i+1)
17	zx_lon(i,jjmp1) = rlon(i+1)
18	ENDDO
19	DO ll=1,klev
20	znivsig(ll)=float(ll)
21	ENDDO
22	CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
23	CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
24	. 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
25	. nhori, nid_nmc)
26	c
27	CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb",
28	. nlevSTD, rlevSTD/100., nvert)
29	ccc
30	ccc Champs 3D interpolles sur des niveaux de pression du NMC
31	ccc
32	IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
33	c
34	CALL histdef(nid_nmc, "temp",
35	. "Temperature","K",
36	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
37	. "inst(X)", zout,zout)
38	c
39	CALL histdef(nid_nmc, "phi",
40	. "Geopotential", "m",
41	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
42	. "inst(X)", zout,zout)
43	c
44	CALL histdef(nid_nmc, "q",
45	. "Specific humidity","kg/kg",
46	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
47	. "inst(X)", zout,zout)
48	c
49	CALL histdef(nid_nmc, "rh",
50	. "Relative humidity", "%",
51	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
52	. "inst(X)", zout,zout)
53	c
54	CALL histdef(nid_nmc, "u",
55	. "Zonal wind","m/s",
56	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
57	. "inst(X)", zout,zout)
58	c
59	CALL histdef(nid_nmc, "v",
60	. "Meridional wind","m/s",
61	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
62	. "inst(X)", zout,zout)
63	c
64	ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
65	c
66	c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
67	c sur les niveaux du NMC et on somme & moyenne
68	c toutes les 6 heures par des routines undefSTD et
69	c moy_undefSTD pour eliminer les valeurs "undef"
70	c de la moyenne mensuelle
71	c ======> le "inst(X)" ci-dessous est par consequence factice !
72	c
73	CALL histdef(nid_nmc, "temp",
74	. "Temperature","K",
75	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
76	. "inst(X)", zout,zout)
77	c
78	CALL histdef(nid_nmc, "phi",
79	. "Geopotential ", "m",
80	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
81	. "inst(X)", zout,zout)
82	c
83	CALL histdef(nid_nmc, "q",
84	. "Specific humidity","kg/kg",
85	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
86	. "inst(X)", zout,zout)
87	c
88	CALL histdef(nid_nmc, "rh",
89	. "Relative humidity", "%",
90	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
91	. "inst(X)", zout,zout)
92	c
93	CALL histdef(nid_nmc, "u",
94	. "Zonal wind","m/s",
95	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
96	. "inst(X)", zout,zout)
97	c
98	CALL histdef(nid_nmc, "v",
99	. "Meridional wind","m/s",
100	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
101	. "inst(X)", zout,zout)
102	c
103	CALL histdef(nid_nmc, "w",
104	. "Vertical motion","Pa/s",
105	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
106	. "inst(X)", zout,zout)
107	c
108	c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
109	c sur les niveaux du NMC et on somme & moyenne
110	c toutes les 6 heures par des routines undefSTD et
111	c moy_undefSTD pour eliminer les valeurs "undef"
112	c de la moyenne mensuelle
113	c ======> le "inst(X)" ci-dessus est par consequence factice !
114	c
115	c
116	CALL histdef(nid_nmc, "psbg",
117	. "Pressure sfce below ground","%",
118	. iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
119	. "inst(X)", zout,zout)
120	c
121	CALL histdef(nid_nmc, "uv",
122	. "uv ",
123	. "m2/s2",iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
124	. "inst(X)", zout,zout)
125	c
126	CALL histdef(nid_nmc, "vq",
127	. "vq ",
128	. "m/s * (kg/kg)",iim,jjmp1,nhori,
129	. nlevSTD,1,nlevSTD, nvert, 32,
130	. "inst(X)", zout,zout)
131	c
132	CALL histdef(nid_nmc, "vT",
133	. "vT ",
134	. "mK/s",iim,jjmp1,nhori,
135	. nlevSTD,1,nlevSTD, nvert, 32,
136	. "inst(X)", zout,zout)
137	c
138	CALL histdef(nid_nmc, "wq",
139	. "wq ",
140	. "(Pa/s)*(kg/kg)",iim,jjmp1,nhori,
141	. nlevSTD,1,nlevSTD, nvert, 32,
142	. "inst(X)", zout,zout)
143	c
144	CALL histdef(nid_nmc, "vphi",
145	. "vphi ",
146	. "m2/s",iim,jjmp1,nhori,
147	. nlevSTD,1,nlevSTD, nvert, 32,
148	. "inst(X)", zout,zout)
149	c
150	CALL histdef(nid_nmc, "wT",
151	. "wT ",
152	. "K*Pa/s",iim,jjmp1,nhori,
153	. nlevSTD,1,nlevSTD, nvert, 32,
154	. "inst(X)", zout,zout)
155	c
156	CALL histdef(nid_nmc, "uxu",
157	. "u2 ",
158	. "m2/s2",iim,jjmp1,nhori,
159	. nlevSTD,1,nlevSTD, nvert, 32,
160	. "inst(X)", zout,zout)
161	c
162	CALL histdef(nid_nmc, "vxv",
163	. "v2 ",
164	. "m2/s2",iim,jjmp1,nhori,
165	. nlevSTD,1,nlevSTD, nvert, 32,
166	. "inst(X)", zout,zout)
167	c
168	CALL histdef(nid_nmc, "TxT",
169	. "T2 ",
170	. "K2",iim,jjmp1,nhori,
171	. nlevSTD,1,nlevSTD, nvert, 32,
172	. "inst(X)", zout,zout)
173	c
174	ENDIF !(type_run.EQ."AMIP")
175
176	CALL histend(nid_nmc)
177	c
178	ENDIF ! fin de test sur ok_mensuel

Note: See TracBrowser for help on using the repository browser.

Download in other formats: