My Project
 All Classes Files Functions Variables Macros
ini_histhfNMC.h
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4 c$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 cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
37 cIM 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 c$OMP END MASTER