source: LMDZ4/trunk/libf/phylmd/ini_histmthNMC.h @ 1374

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

Add logical flags ok_histNMC(3) to manage output of the
standard pressure files hist*NMC.nc
IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.9 KB
RevLine 
[524]1!
2! $Header$
3!
[766]4c$OMP MASTER
[1374]5c
6      IF (ok_histNMC(1)) THEN
7c
[1352]8       zout = freq_outNMC(1)
[524]9c
[1352]10       idayref = day_ref
11       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
[524]12c
[766]13cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
14cym         DO i = 1, iim
15cym            zx_lon(i,1) = rlon(i+1)
16cym            zx_lon(i,jjmp1) = rlon(i+1)
17cym         ENDDO
[524]18         DO ll=1,klev
19            znivsig(ll)=float(ll)
20         ENDDO
[766]21cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
22cym         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
23cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
[1352]24cym     .                 nhori, nid_mthnmc)
[766]25
[1352]26         CALL histbeg_phy("histmthNMC",itau_phy, zjulian, dtime, 
27     .                 nhori, nid_mthnmc)
[644]28c
[1352]29         CALL histvert(nid_mthnmc, "plev", "pressure", "Pa",
30     .                 nlevSTD, rlevSTD, nvert)
[644]31ccc
32ccc Champs 3D interpolles sur des niveaux de pression du NMC
33ccc
[524]34c
[644]35c ATTENTION : pour AMIP2 on interpole t,u,v,wphi,q,rh
36c             sur les niveaux du NMC et on somme & moyenne
[1352]37c             toutes les freq_moyNMC secondes par des routines undefSTD et
[644]38c             moy_undefSTD pour eliminer les valeurs "undef"
39c             de la moyenne mensuelle
40c ======> le "inst(X)" ci-dessous est par consequence factice !
[524]41c
[1352]42c
43          CALL histdef(nid_mthnmc, "tnondef",
44     .                 "Valeurs non-definies","-",
[776]45     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]46     .                "inst(X)", zout,zout)
[524]47c
[1352]48          CALL histdef(nid_mthnmc, "ta",
49     .                 "Air temperature","K",
[776]50     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]51     .                "inst(X)", zout,zout)
[524]52c
[1352]53         CALL histdef(nid_mthnmc, "zg",
54     .                "Geopotential height", "m",
[776]55     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]56     .                "inst(X)", zout,zout)
57c
[1352]58          CALL histdef(nid_mthnmc, "hus",
59     .                 "Specific humidity","1",
60     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
61     .                "inst(X)", zout,zout)
62c
63         CALL histdef(nid_mthnmc, "hur",
[644]64     .                 "Relative humidity", "%",
[776]65     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]66     .                "inst(X)", zout,zout)
67c
[1352]68          CALL histdef(nid_mthnmc, "ua",
69     .                 "Eastward wind","m s-1",
[776]70     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]71     .                "inst(X)", zout,zout)
72c
[1352]73          CALL histdef(nid_mthnmc, "va",
74     .                 "Northward wind","m s-1",
[776]75     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]76     .                "inst(X)", zout,zout)
77c
[1352]78          CALL histdef(nid_mthnmc, "wap",
79     .                 "Lagrangian tendency of air pressure","Pa s-1",
[776]80     .                iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]81     .                "inst(X)", zout,zout)
82c
[1352]83          CALL histdef(nid_mthnmc, "psbg",
[644]84     .         "Pressure sfce below ground","%",
[776]85     .         iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]86     .         "inst(X)", zout,zout)
87c
[1352]88          CALL histdef(nid_mthnmc, "uv",
[644]89     .         "uv ",
[776]90     .         "m2/s2",iim,jj_nb,nhori, nlevSTD,1,nlevSTD, nvert, 32,
[644]91     .         "inst(X)", zout,zout)
92c
[1352]93          CALL histdef(nid_mthnmc, "vq",
[644]94     .         "vq ",
[776]95     .         "m/s * (kg/kg)",iim,jj_nb,nhori, 
[644]96     .          nlevSTD,1,nlevSTD, nvert, 32,
97     .         "inst(X)", zout,zout)
98c
[1352]99          CALL histdef(nid_mthnmc, "vT",
[644]100     .         "vT ", 
[776]101     .         "mK/s",iim,jj_nb,nhori, 
[644]102     .          nlevSTD,1,nlevSTD, nvert, 32,
103     .         "inst(X)", zout,zout)
104c
[1352]105          CALL histdef(nid_mthnmc, "wq",
[644]106     .         "wq ", 
[776]107     .         "(Pa/s)*(kg/kg)",iim,jj_nb,nhori,
[644]108     .          nlevSTD,1,nlevSTD, nvert, 32,
109     .         "inst(X)", zout,zout)
110c
[1352]111          CALL histdef(nid_mthnmc, "vphi",
[644]112     .         "vphi ", 
[776]113     .         "m2/s",iim,jj_nb,nhori, 
[644]114     .          nlevSTD,1,nlevSTD, nvert, 32,
115     .         "inst(X)", zout,zout)
116c
[1352]117          CALL histdef(nid_mthnmc, "wT",
[644]118     .         "wT ", 
[776]119     .         "K*Pa/s",iim,jj_nb,nhori,
[644]120     .          nlevSTD,1,nlevSTD, nvert, 32,
121     .         "inst(X)", zout,zout)
122c
[1352]123          CALL histdef(nid_mthnmc, "uxu",
[644]124     .         "u2 ", 
[776]125     .         "m2/s2",iim,jj_nb,nhori,
[644]126     .          nlevSTD,1,nlevSTD, nvert, 32,
127     .         "inst(X)", zout,zout)
128c
[1352]129          CALL histdef(nid_mthnmc, "vxv",
[644]130     .         "v2 ", 
[776]131     .         "m2/s2",iim,jj_nb,nhori,
[644]132     .          nlevSTD,1,nlevSTD, nvert, 32,
133     .         "inst(X)", zout,zout)
134c
[1352]135          CALL histdef(nid_mthnmc, "TxT",
[644]136     .         "T2 ", 
[776]137     .         "K2",iim,jj_nb,nhori,
[644]138     .          nlevSTD,1,nlevSTD, nvert, 32,
139     .         "inst(X)", zout,zout)
140c
[524]141
[1352]142         CALL histend(nid_mthnmc)
[524]143c
[1374]144      ENDIF !(ok_histNMC(1)) THEN
145c
[766]146c$OMP END MASTER
Note: See TracBrowser for help on using the repository browser.