source: LMDZ4/trunk/libf/phylmd/ini_histhfNMC.h @ 1398

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

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

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