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
RevLine 
[1351]1!
2! $Header$
3!
4c$OMP MASTER
[1374]5      IF (ok_histNMC(2)) THEN
6c
[1351]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
[1398]32         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
33     .                 nlevSTD, rlevSTD/100., nvert,"down")
[1351]34        ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
[1398]35         CALL histvert(nid_daynmc, "plev", "pressure", "mb",
36     .                 nlevSTD8, rlevSTD8/100., nvert,"down")
[1351]37        ENDIF
[1398]38c
39cIM Astuce MAF: remplacer inst par ave pour les variables NMC pour avoir
40cIM             le time_counter et les bounds
41c
[1351]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,
[1398]56     .                "ave(X)", zout,zout)
[1351]57c
58          CALL histdef(nid_daynmc, "ta",
59     .                 "Air temperature","K",
60     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]61     .                "ave(X)", zout,zout)
[1351]62c
63         CALL histdef(nid_daynmc, "zg",
64     .                "Geopotential height", "m",
65     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]66     .                "ave(X)", zout,zout)
[1351]67c
68          CALL histdef(nid_daynmc, "hus",
69     .                 "Specific humidity","1",
70     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]71     .                "ave(X)", zout,zout)
[1351]72c
73         CALL histdef(nid_daynmc, "hur",
74     .                 "Relative humidity", "%",
75     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]76     .                "ave(X)", zout,zout)
[1351]77c
78          CALL histdef(nid_daynmc, "ua",
79     .                 "Eastward wind","m s-1",
80     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]81     .                "ave(X)", zout,zout)
[1351]82c
83          CALL histdef(nid_daynmc, "va",
84     .                 "Northward wind","m s-1",
85     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[1398]86     .                "ave(X)", zout,zout)
[1351]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,
[1398]91     .                "ave(X)", zout,zout)
[1351]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,
[1398]98     .                "ave(X)", zout,zout)
[1351]99c
100          CALL histdef(nid_daynmc, "ta",
101     .                 "Air temperature","K",
102     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]103     .                "ave(X)", zout,zout)
[1351]104c
105         CALL histdef(nid_daynmc, "zg",
106     .                "Geopotential height", "m",
107     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]108     .                "ave(X)", zout,zout)
[1351]109c
110          CALL histdef(nid_daynmc, "hus",
111     .                 "Specific humidity","1",
112     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]113     .                "ave(X)", zout,zout)
[1351]114c
115         CALL histdef(nid_daynmc, "hur",
116     .                 "Relative humidity", "%",
117     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]118     .                "ave(X)", zout,zout)
[1351]119c
120          CALL histdef(nid_daynmc, "ua",
121     .                 "Eastward wind","m s-1",
122     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]123     .                "ave(X)", zout,zout)
[1351]124c
125          CALL histdef(nid_daynmc, "va",
126     .                 "Northward wind","m s-1",
127     .                iim,jj_nb,nhori, nlevSTD8,1,nlevSTD8, nvert, 32,
[1398]128     .                "ave(X)", zout,zout)
[1351]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,
[1398]133     .                "ave(X)", zout,zout)
[1351]134c
135        ENDIF
136c
137         CALL histend(nid_daynmc)
138c
[1374]139      ENDIF !(ok_histNMC(2)) THEN
[1351]140c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.