source: LMDZ4/branches/LMDZ4_par_0/libf/phylmd/ini_histmthNMC.h @ 3400

Last change on this file since 3400 was 634, checked in by Laurent Fairhead, 20 years ago

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

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