source: LMDZ5/trunk/libf/phylmd/ini_histdayNMC.h @ 1641

Last change on this file since 1641 was 1400, checked in by musat, 14 years ago

Change units for vertical axes from "mb" to (SI) "Pa"
Change titles for variables at standard pressure levels from "mb" to "hPa"
IM

  • Property svn:executable set to *
File size: 5.2 KB
Line 
1!
2! $Header$
3!
4c$OMP MASTER
5      IF (ok_histNMC(2)) THEN
6c
7       zstophy = dtime
8       zstohf = ecrit_hf
9       zstomth = ecrit_mth
10       zout = freq_outNMC(2)
11c
12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14c
15cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
16cym         DO i = 1, iim
17cym            zx_lon(i,1) = rlon(i+1)
18cym            zx_lon(i,jjmp1) = rlon(i+1)
19cym         ENDDO
20         DO ll=1,klev
21            znivsig(ll)=float(ll)
22         ENDDO
23cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24cym         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
25cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
26cym     .                 nhori, nid_daynmc)
27
28         CALL histbeg_phy("histdayNMC",itau_phy, zjulian, dtime, 
29     .                 nhori, nid_daynmc)
30c
31        IF(lev_histdayNMC.EQ.nlevSTD) THEN
32         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
33     .                 nlevSTD, rlevSTD, nvert,"down")
34        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
35         CALL histvert(nid_daynmc, "plev", "pressure", "Pa",
36     .                 nlevSTD8, rlevSTD8, nvert,"down")
37        ENDIF
38c
39cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
40cIM             le time_counter et les bounds
41c
42ccc Champs 3D interpolles sur des niveaux de pression du NMC
43ccc
44c
45c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
46c             sur les niveaux du NMC et on somme & moyenne
47c             toutes les freq_moyNMC secondes par des routines undefSTD et
48c             moy_undefSTD pour eliminer les valeurs "undef"
49c             de la moyenne mensuelle
50c ======> le "inst(X)" ci-dessous est par consequence factice !
51c
52        IF(lev_histdayNMC.EQ.nlevSTD) THEN
53          CALL histdef(nid_daynmc, "tnondef",
54     .                 "Valeurs non-definies","-",
55     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
56     .                "ave(X)", zout,zout)
57c
58          CALL histdef(nid_daynmc, "ta",
59     .                 "Air temperature","K",
60     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
61     .                "ave(X)", zout,zout)
62c
63         CALL histdef(nid_daynmc, "zg",
64     .                "Geopotential height", "m",
65     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
66     .                "ave(X)", zout,zout)
67c
68          CALL histdef(nid_daynmc, "hus",
69     .                 "Specific humidity","1",
70     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
71     .                "ave(X)", zout,zout)
72c
73         CALL histdef(nid_daynmc, "hur",
74     .                 "Relative humidity", "%",
75     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
76     .                "ave(X)", zout,zout)
77c
78          CALL histdef(nid_daynmc, "ua",
79     .                 "Eastward wind","m s-1",
80     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
81     .                "ave(X)", zout,zout)
82c
83          CALL histdef(nid_daynmc, "va",
84     .                 "Northward wind","m s-1",
85     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
86     .                "ave(X)", zout,zout)
87c
88          CALL histdef(nid_daynmc, "wap",
89     .                 "Lagrangian tendency of air pressure","Pa s-1",
90     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
91     .                "ave(X)", zout,zout)
92c
93        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
94c
95          CALL histdef(nid_daynmc, "tnondef",
96     .                 "Valeurs non-definies","-",
97     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
98     .                "ave(X)", zout,zout)
99c
100          CALL histdef(nid_daynmc, "ta",
101     .                 "Air temperature","K",
102     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
103     .                "ave(X)", zout,zout)
104c
105         CALL histdef(nid_daynmc, "zg",
106     .                "Geopotential height", "m",
107     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
108     .                "ave(X)", zout,zout)
109c
110          CALL histdef(nid_daynmc, "hus",
111     .                 "Specific humidity","1",
112     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
113     .                "ave(X)", zout,zout)
114c
115         CALL histdef(nid_daynmc, "hur",
116     .                 "Relative humidity", "%",
117     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
118     .                "ave(X)", zout,zout)
119c
120          CALL histdef(nid_daynmc, "ua",
121     .                 "Eastward wind","m s-1",
122     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
123     .                "ave(X)", zout,zout)
124c
125          CALL histdef(nid_daynmc, "va",
126     .                 "Northward wind","m s-1",
127     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
128     .                "ave(X)", zout,zout)
129c
130          CALL histdef(nid_daynmc, "wap",
131     .                 "Lagrangian tendency of air pressure","Pa s-1",
132     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
133     .                "ave(X)", zout,zout)
134c
135        ENDIF
136c
137         CALL histend(nid_daynmc)
138c
139      ENDIF !(ok_histNMC(2)) THEN
140c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.