source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histdayNMC.h @ 5360

Last change on this file since 5360 was 1349, checked in by musat, 15 years ago

Add new output files for variables on standard pressure levels at daily and high
frequency (3hr)
IM
  • Property svn:executable set to *
File size: 5.0 KB
RevLine 
[1349]1!
2! $Header$
3!
4c$OMP MASTER
5
6       zstophy = dtime
7       zstohf = ecrit_hf
8       zstomth = ecrit_mth
9       zout = freq_outNMC(2)
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_daynmc)
26
27         CALL histbeg_phy("histdayNMC",itau_phy, zjulian, dtime, 
28     .                 nhori, nid_daynmc)
29c
30        IF(lev_histdayNMC.EQ.nlevSTD) THEN
31         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
32     .                 nlevSTD, rlevSTD, nvert)
33        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
34         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
35     .                 nlevSTD8, rlevSTD8, nvert)
36        ENDIF
37ccc
38ccc Champs 3D interpolles sur des niveaux de pression du NMC
39ccc
40c
41c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
42c             sur les niveaux du NMC et on somme & moyenne
43c             toutes les freq_moyNMC secondes par des routines undefSTD et
44c             moy_undefSTD pour eliminer les valeurs "undef"
45c             de la moyenne mensuelle
46c ======> le "inst(X)" ci-dessous est par consequence factice !
47c
48        IF(lev_histdayNMC.EQ.nlevSTD) THEN
49          CALL histdef(nid_daynmc, "tnondef",
50     .                 "Valeurs non-definies","-",
51     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
52     .                "inst(X)", zout,zout)
53c
54          CALL histdef(nid_daynmc, "ta",
55     .                 "Air temperature","K",
56     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
57     .                "inst(X)", zout,zout)
58c
59         CALL histdef(nid_daynmc, "zg",
60     .                "Geopotential height", "m",
61     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
62     .                "inst(X)", zout,zout)
63c
64          CALL histdef(nid_daynmc, "hus",
65     .                 "Specific humidity","1",
66     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
67     .                "inst(X)", zout,zout)
68c
69         CALL histdef(nid_daynmc, "hur",
70     .                 "Relative humidity", "%",
71     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
72     .                "inst(X)", zout,zout)
73c
74          CALL histdef(nid_daynmc, "ua",
75     .                 "Eastward wind","m s-1",
76     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
77     .                "inst(X)", zout,zout)
78c
79          CALL histdef(nid_daynmc, "va",
80     .                 "Northward wind","m s-1",
81     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
82     .                "inst(X)", zout,zout)
83c
84          CALL histdef(nid_daynmc, "wap",
85     .                 "Lagrangian tendency of air pressure","Pa s-1",
86     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
87     .                "inst(X)", zout,zout)
88c
89        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
90c
91          CALL histdef(nid_daynmc, "tnondef",
92     .                 "Valeurs non-definies","-",
93     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
94     .                "inst(X)", zout,zout)
95c
96          CALL histdef(nid_daynmc, "ta",
97     .                 "Air temperature","K",
98     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
99     .                "inst(X)", zout,zout)
100c
101         CALL histdef(nid_daynmc, "zg",
102     .                "Geopotential height", "m",
103     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
104     .                "inst(X)", zout,zout)
105c
106          CALL histdef(nid_daynmc, "hus",
107     .                 "Specific humidity","1",
108     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
109     .                "inst(X)", zout,zout)
110c
111         CALL histdef(nid_daynmc, "hur",
112     .                 "Relative humidity", "%",
113     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
114     .                "inst(X)", zout,zout)
115c
116          CALL histdef(nid_daynmc, "ua",
117     .                 "Eastward wind","m s-1",
118     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
119     .                "inst(X)", zout,zout)
120c
121          CALL histdef(nid_daynmc, "va",
122     .                 "Northward wind","m s-1",
123     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
124     .                "inst(X)", zout,zout)
125c
126          CALL histdef(nid_daynmc, "wap",
127     .                 "Lagrangian tendency of air pressure","Pa s-1",
128     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
129     .                "inst(X)", zout,zout)
130c
131        ENDIF
132c
133         CALL histend(nid_daynmc)
134c
135c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.