source: LMDZ4/branches/unlabeled-1.1.1/libf/phylmd/ini_histmthNMC.h @ 4667

Last change on this file since 4667 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.4 KB
Line 
1!
2! $Header$
3!
4      IF (ok_mensuel) THEN
5c
6         zsto = dtime
7         zout = dtime * ecrit_mth
8c
9         idayref = day_ref
10         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
11c
12         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
13         DO i = 1, iim
14            zx_lon(i,1) = rlon(i+1)
15            zx_lon(i,jjmp1) = rlon(i+1)
16         ENDDO
17         DO ll=1,klev
18            znivsig(ll)=float(ll)
19         ENDDO
20         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
21         CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
22     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 
23     .                 nhori, nid_nmc)
24         write(*,*)'Mensuel NMC ', itau_phy, zjulian
25         CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb",
26     .                 klev, presnivs, nvert)
27c        call histvert(nid_nmc, 'sig_s', 'Niveaux sigma','-',
28c    .              klev, znivsig, nvert)
29c
30c Champs 2D:
31
32c Champs interpolles sur des niveaux de pression du NMC
33c IMIMIM 110304 BEG
34c
35         DO k=1, nlevSTD
36c
37         bb=clevSTD(k)
38c
39         IF(k.GE.2) THEN
40          aa=clevSTD(k)
41          bb=aa(1:lnblnk1(aa))
42         ENDIF
43c
44          CALL histdef(nid_nmc, "t"//bb,
45     .                 "Temperature"//bb//"mb","K",
46     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
47     .                "ave(X)", zsto,zout)
48c
49         CALL histdef(nid_nmc, "phi"//bb,
50     .                "Geopotential"//bb//"mb", "m2/s2",
51     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
52     .                "ave(X)", zsto,zout)
53c
54          CALL histdef(nid_nmc, "q"//bb,
55     .                 "Specific humidity"//bb//"mb","kg/kg",
56     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
57     .                "ave(X)", zsto,zout)
58c
59         CALL histdef(nid_nmc, "rh"//bb,
60     .                 "Relative humidity"//bb//"mb", "%",
61     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
62     .                "ave(X)", zsto,zout)
63c
64          CALL histdef(nid_nmc, "u"//bb,
65     .                 "Zonal wind"//bb//"mb","K",
66     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
67     .                "ave(X)", zsto,zout)
68c
69          CALL histdef(nid_nmc, "v"//bb,
70     .                 "Meridional wind"//bb//"mb","K",
71     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
72     .                "ave(X)", zsto,zout)
73c
74         ENDDO !nlevSTD
75c IMIMIM 110304 END
76
77         CALL histend(nid_nmc)
78c
79c        ndex2d = 0
80c
81      ENDIF ! fin de test sur ok_mensuel
82
Note: See TracBrowser for help on using the repository browser.