source: LMDZ4/trunk/libf/phylmd/ini_histdayNMC.h @ 1398

Last change on this file since 1398 was 1398, checked in by musat, 14 years ago

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

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