source: lmdz_wrf/trunk/WRFV3/lmdz/ini_histdayNMC.h @ 1672

Last change on this file since 1672 was 1, checked in by lfita, 11 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: 6.8 KB
Line 
1!
2! $Header$
3!
4!$OMP MASTER
5      IF (ok_histNMC(2)) THEN
6!c
7       zstophy = dtime
8       zstohf = ecrit_hf
9       zstomth = ecrit_mth
10       zout = freq_outNMC(2)
11!c
12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14!c
15!cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
16!cym         DO i = 1, iim
17!cym            zx_lon(i,1) = rlon(i+1)
18!cym            zx_lon(i,jjmp1) = rlon(i+1)
19!cym         ENDDO
20         DO ll=1,klev
21            znivsig(ll)=float(ll)
22         ENDDO
23!cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24!cym         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
25!cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
26!cym     .                 nhori, nid_daynmc)
27
28         CALL histbeg_phy("histdayNMC",itau_phy, zjulian, dtime,                     &
29       &                 nhori, nid_daynmc)
30!c
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
38!c
39!IM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
40!IM             le time_counter et les bounds
41!c
42!ccc Champs 3D interpolles sur des niveaux de pression du NMC
43!ccc
44!c
45!c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
46!c             sur les niveaux du NMC et on somme & moyenne
47!c             toutes les freq_moyNMC secondes par des routines undefSTD et
48!c             moy_undefSTD pour eliminer les valeurs "undef"
49!c             de la moyenne mensuelle
50!c ======> le "inst(X)" ci-dessous est par consequence factice !
51!c
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)
57!c
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)
62!c
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)
67!c
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)
72!c
73         CALL histdef(nid_daynmc, "hur",                                             &
74       &                 "Relative humidity", "%",                                   &
75       &                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,               &
76       &                "ave(X)", zout,zout)
77!c
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)
82!c
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)
87!c
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)
92!c
93        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
94!c
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)
99!c
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)
104!c
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)
109!c
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)
114!c
115         CALL histdef(nid_daynmc, "hur",                                             &
116       &                 "Relative humidity", "%",                                   &
117       &                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,             &
118       &                "ave(X)", zout,zout)
119!c
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)
124!c
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)
129!c
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)
134!c
135        ENDIF
136!c
137         CALL histend(nid_daynmc)
138!c
139      ENDIF !(ok_histNMC(2)) THEN
140!$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.