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