Context Navigation

source: LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h @ 675

Last change on this file since 675 was 644, checked in by Laurent Fairhead, 19 years ago
Synchronisation avec tous les diagnostiques de Ionela IM Inclusion du slab ocean IM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 6.2 KB

Line
1	!
2	! $Header$
3	!
4	IF (ok_mensuel) THEN
5	c
6	zstophy = dtime
7	zstohf = dtime * FLOAT(ecrit_hf)
8	zstomth = dtime * FLOAT(ecrit_mth)
9	zout = dtime * FLOAT(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: