source: LMDZ4/branches/V3_test/libf/phylmd/ini_histmthNMC.h @ 718

Last change on this file since 718 was 704, checked in by Laurent Fairhead, 19 years ago

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