source: lmdz_wrf/branches/LMDZ_WRFmeas/WRFV3/lmdz/ini_histmthNMC.h @ 17

Last change on this file since 17 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

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