source: lmdz_wrf/WRFV3/lmdz/ini_histhfNMC.h @ 1

Last change on this file since 1 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:

  • Property svn:executable set to *
File size: 7.7 KB
Line 
1!
2! $Header$
3!
4!$OMP MASTER
5!c
6      IF (ok_histNMC(3)) THEN
7!c
8       zstophy = dtime
9       zstohf = ecrit_hf
10       zstomth = ecrit_mth
11!c      zout = 6 * 3600.
12       zout = freq_outNMC(3)
13!c
14         idayref = day_ref
15         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
16!c
17!cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
18!cym         DO i = 1, iim
19!cym            zx_lon(i,1) = rlon(i+1)
20!cym            zx_lon(i,jjmp1) = rlon(i+1)
21!cym         ENDDO
22         DO ll=1,klev
23            znivsig(ll)=float(ll)
24         ENDDO
25!cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
26!cym         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
27!cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
28!cym     .                 nhori, nid_hfnmc)
29
30         CALL histbeg_phy("histhfNMC",itau_phy, zjulian, dtime,                      &
31       &                 nhori, nid_hfnmc)
32!c
33         CALL histvert(nid_hfnmc, "plev", "pressure", "Pa",                          &
34       &                 nlevSTD3, rlevSTD3, nvert,"down")
35!c
36!IM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
37!IM             le time_counter et les bounds
38!c
39!ccc
40!ccc Champs 3D interpolles sur des niveaux de pression du NMC
41!ccc
42!c
43!c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
44!c             sur les niveaux du NMC et on somme & moyenne
45!c             toutes les freq_moyNMC secondes par des routines undefSTD et
46!c             moy_undefSTD pour eliminer les valeurs "undef"
47!c             de la moyenne mensuelle
48!c ======> le "inst(X)" ci-dessous est par consequence factice !
49!c
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)
54!c
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)
59!c
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)
64!c
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)
69!c
70         CALL histdef(nid_hfnmc, "hur",                                              &
71       &                 "Relative humidity", "%",                                   &
72       &                iim,jj_nb,nhori, nlevSTD3,1,nlevSTD3, nvert, 32,             &
73       &                "ave(X)", zout,zout)
74!c
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)
79!c
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)
84!c
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)
89!c
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)
94!c
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)
99!c
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)
105!c
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)
111!c
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)
117!c
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)
123!c
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)
129!c
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)
135!c
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)
141!c
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)
147!c
148         CALL histend(nid_hfnmc)
149!c
150      ENDIF !(ok_histNMC(2)) THEN
151!c
152!$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.