source: LMDZ.3.3/trunk/libf/phylmd/ini_histmthNMC.h @ 5322

Last change on this file since 5322 was 485, checked in by lmdzadmin, 21 years ago

Initial import
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 7.7 KB
Line 
1c     IF (ok_mensuelNMC) THEN
2      IF (ok_mensuel) THEN
3c
4       IF(region.EQ.'3d') THEN
5         zsto = dtime
6         zout = dtime * ecrit_mth
7       ELSEIF(region.EQ.'1d') THEN
8         zsto = dtime
9c        zout = dtime
10         zout = dtime *4.
11       ENDIF
12c
13         idayref = day_ref
14         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
15c
16         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
17         DO i = 1, iim
18            zx_lon(i,1) = rlon(i+1)
19            zx_lon(i,jjmp1) = rlon(i+1)
20         ENDDO
21         DO ll=1,klev
22            znivsig(ll)=float(ll)
23         ENDDO
24         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
25         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
26c071103
27c    .                 1,iim,1,jjmp1, 
28     .                 imin_debut,nbpti,
29     .                 jmin_debut,nbptj,
30c071103
31     .                 itau_phy, zjulian, dtime, 
32     .                 nhori, nid_nmc)
33         write(*,*)'Mensuel NMC ', itau_phy, zjulian
34         CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb",
35     .                 klev, presnivs, nvert)
36c        call histvert(nid_nmc, 'sig_s', 'Niveaux sigma','-',
37c    .              klev, znivsig, nvert)
38c
39c Champs 2D:
40
41c Champs interpolles sur des niveaux de pression du NMC
42         CALL histdef(nid_nmc, "u1000",
43     . "Zonal wind 1000mb", 
44     .                "m/s",
45     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
46     .                "ave(X)", zsto,zout)
47
48         CALL histdef(nid_nmc, "v1000",
49     .                "Meridional wind 1000mb", "m/s",
50     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
51     .                "ave(X)", zsto,zout)
52
53         CALL histdef(nid_nmc, "u925", "Zonal wind 925mb", 
54     .                "m/s",
55     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
56     .                "ave(X)", zsto,zout)
57
58         CALL histdef(nid_nmc, "v925", "Meridional wind 925mb", 
59     .                "m/s",
60     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
61     .                "ave(X)", zsto,zout)
62
63         CALL histdef(nid_nmc, "u850", "Zonal wind 850mb", 
64     .                "m/s",
65     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
66     .                "ave(X)", zsto,zout)
67
68         CALL histdef(nid_nmc, "v850", "Meridional wind 850mb", 
69     .                "m/s",
70     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
71     .                "ave(X)", zsto,zout)
72c
73         CALL histdef(nid_nmc, "u700", "Zonal wind 700mb", 
74     .                "m/s",
75     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
76     .                "ave(X)", zsto,zout)
77
78         CALL histdef(nid_nmc, "v700", "Meridional wind 700mb", 
79     .                "m/s",
80     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
81     .                "ave(X)", zsto,zout)
82c
83         CALL histdef(nid_nmc, "u600", "Zonal wind 600mb", 
84     .                "m/s",
85     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
86     .                "ave(X)", zsto,zout)
87
88         CALL histdef(nid_nmc, "v600", "Meridional wind 600mb", 
89     .                "m/s",
90     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
91     .                "ave(X)", zsto,zout)
92c
93         CALL histdef(nid_nmc, "u500", "Zonal wind 500mb", 
94     .                "m/s",
95     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
96     .                "ave(X)", zsto,zout)
97
98         CALL histdef(nid_nmc, "v500", "Meridional wind 500mb", 
99     .                "m/s",
100     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
101     .                "ave(X)", zsto,zout)
102c
103         CALL histdef(nid_nmc, "u400", "Zonal wind 400mb", 
104     .                "m/s",
105     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
106     .                "ave(X)", zsto,zout)
107
108         CALL histdef(nid_nmc, "v400", "Meridional wind 400mb", 
109     .                "m/s",
110     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
111     .                "ave(X)", zsto,zout)
112c
113         CALL histdef(nid_nmc, "u300", "Zonal wind 300mb", 
114     .                "m/s",
115     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
116     .                "ave(X)", zsto,zout)
117
118         CALL histdef(nid_nmc, "v300", "Meridional wind 300mb", 
119     .                "m/s",
120     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
121     .                "ave(X)", zsto,zout)
122c
123         CALL histdef(nid_nmc, "u250", "Zonal wind 250mb", 
124     .                "m/s",
125     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
126     .                "ave(X)", zsto,zout)
127
128         CALL histdef(nid_nmc, "v250", "Meridional wind 250mb", 
129     .                "m/s",
130     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
131     .                "ave(X)", zsto,zout)
132c
133         CALL histdef(nid_nmc, "u200", "Zonal wind 200mb", 
134     .                "m/s",
135     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
136     .                "ave(X)", zsto,zout)
137
138         CALL histdef(nid_nmc, "v200", "Meridional wind 200mb", 
139     .                "m/s",
140     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
141     .                "ave(X)", zsto,zout)
142c
143         CALL histdef(nid_nmc, "u150", "Zonal wind 150mb", 
144     .                "m/s",
145     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
146     .                "ave(X)", zsto,zout)
147
148         CALL histdef(nid_nmc, "v150", "Meridional wind 150mb", 
149     .                "m/s",
150     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
151     .                "ave(X)", zsto,zout)
152c
153         CALL histdef(nid_nmc, "u100", "Zonal wind 100mb", 
154     .                "m/s",
155     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
156     .                "ave(X)", zsto,zout)
157
158         CALL histdef(nid_nmc, "v100", "Meridional wind 100mb", 
159     .                "m/s",
160     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
161     .                "ave(X)", zsto,zout)
162c
163         CALL histdef(nid_nmc, "u70", "Zonal wind 70mb", 
164     .                "m/s",
165     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
166     .                "ave(X)", zsto,zout)
167
168         CALL histdef(nid_nmc, "v70", "Meridional wind 70mb", 
169     .                "m/s",
170     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
171     .                "ave(X)", zsto,zout)
172c
173         CALL histdef(nid_nmc, "u50", "Zonal wind 50mb", 
174     .                "m/s",
175     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
176     .                "ave(X)", zsto,zout)
177
178         CALL histdef(nid_nmc, "v50", "Meridional wind 50mb", 
179     .                "m/s",
180     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
181     .                "ave(X)", zsto,zout)
182c
183         CALL histdef(nid_nmc, "u30", "Zonal wind 30mb", 
184     .                "m/s",
185     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
186     .                "ave(X)", zsto,zout)
187
188         CALL histdef(nid_nmc, "v30", "Meridional wind 30mb", 
189     .                "m/s",
190     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
191     .                "ave(X)", zsto,zout)
192c
193         CALL histdef(nid_nmc, "u20", "Zonal wind 20mb", 
194     .                "m/s",
195     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
196     .                "ave(X)", zsto,zout)
197
198         CALL histdef(nid_nmc, "v20", "Meridional wind 20mb", 
199     .                "m/s",
200     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
201     .                "ave(X)", zsto,zout)
202c
203         CALL histdef(nid_nmc, "u10", "Zonal wind 10mb", 
204     .                "m/s",
205     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
206     .                "ave(X)", zsto,zout)
207
208         CALL histdef(nid_nmc, "v10", "Meridional wind 10mb", 
209     .                "m/s",
210     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
211     .                "ave(X)", zsto,zout)
212c
213         CALL histend(nid_nmc)
214c
215c        ndex2d = 0
216c
217      ENDIF ! fin de test sur ok_mensuel
218
Note: See TracBrowser for help on using the repository browser.