source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmthNMC.h @ 5220

Last change on this file since 5220 was 740, checked in by lmdzadmin, 18 years ago

Correction bogues: les ecrit_ sont des REALs lus dans conf_phys.F90 en
nombre de jours sauf pour ecrit_ins et ecrit_tra en secondes!
Les ecrit_ sont initialises dans conf_phys.F90 et peuvent etre modifies dans
physiq.def.
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
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
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)
26c
27         CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb",
28     .                 nlevSTD, rlevSTD/100., nvert)
29ccc
30ccc Champs 3D interpolles sur des niveaux de pression du NMC
31ccc
32      IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
33c
34          CALL histdef(nid_nmc, "temp",
35     .                 "Temperature","K",
36     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
37     .                "inst(X)", zout,zout)
38c
39         CALL histdef(nid_nmc, "phi",
40     .                "Geopotential", "m",
41     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
42     .                "inst(X)", zout,zout)
43c
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)
48c
49         CALL histdef(nid_nmc, "rh",
50     .                 "Relative humidity", "%",
51     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
52     .                "inst(X)", zout,zout)
53c
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)
58c
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)
63c
64      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
65c
66c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
67c             sur les niveaux du NMC et on somme & moyenne
68c             toutes les 6 heures par des routines undefSTD et
69c             moy_undefSTD pour eliminer les valeurs "undef"
70c             de la moyenne mensuelle
71c ======> le "inst(X)" ci-dessous est par consequence factice !
72c
73          CALL histdef(nid_nmc, "temp",
74     .                 "Temperature","K",
75     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
76     .                "inst(X)", zout,zout)
77c
78         CALL histdef(nid_nmc, "phi",
79     .                "Geopotential ", "m",
80     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
81     .                "inst(X)", zout,zout)
82c
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)
87c
88         CALL histdef(nid_nmc, "rh",
89     .                 "Relative humidity", "%",
90     .                iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
91     .                "inst(X)", zout,zout)
92c
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)
97c
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)
102c
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)
107c
108c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
109c             sur les niveaux du NMC et on somme & moyenne
110c             toutes les 6 heures par des routines undefSTD et
111c             moy_undefSTD pour eliminer les valeurs "undef"
112c             de la moyenne mensuelle
113c ======> le "inst(X)" ci-dessus est par consequence factice !
114c
115c
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)
120c
121          CALL histdef(nid_nmc, "uv",
122     .         "uv ",
123     .         "m2/s2",iim,jjmp1,nhori, nlevSTD,1,nlevSTD, nvert, 32,
124     .         "inst(X)", zout,zout)
125c
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)
131c
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)
137c
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)
143c
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)
149c
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)
155c
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)
161c
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)
167c
168          CALL histdef(nid_nmc, "TxT",
169     .         "T2 ", 
170     .         "K2",iim,jjmp1,nhori,
171     .          nlevSTD,1,nlevSTD, nvert, 32,
172     .         "inst(X)", zout,zout)
173c
174      ENDIF !(type_run.EQ."AMIP")
175
176         CALL histend(nid_nmc)
177c
178      ENDIF ! fin de test sur ok_mensuel
Note: See TracBrowser for help on using the repository browser.