source: LMDZ.3.3/trunk/libf/phylmd/ini_histhf.h @ 4944

Last change on this file since 4944 was 416, checked in by lmdzadmin, 22 years ago

Inclusion initiale

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.1 KB
Line 
1
2      IF (ok_hf) THEN
3c
4         PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf
5
6cccIM    CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
7         CALL ymds2ju(annee_ref, 1, 1, 0.0, zjulian)
8         zjulian = zjulian + day_ini
9c
10         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
11         DO i = 1, iim
12            zx_lon(i,1) = rlon(i+1)
13            zx_lon(i,jjmp1) = rlon(i+1)
14         ENDDO
15
16         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
17
18cccIM      CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat,
19         CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
20     .                 1,iim,1,jjmp1, 0, zjulian, dtime, 
21     .                 nhori, nid_hf)
22
23         CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb",
24     .                 klev, presnivs, nvert)
25         zsto = dtime
26
27c   pour les champs instantannes, il faut mettre la meme valeur pour
28c   zout et tsto.
29c   dtime est passe par ailleurs a histbeg
30
31         zout = dtime * FLOAT(NINT(86400./dtime*ecrit_hf))
32         zsto = zout
33         print*,'zout,zsto=',zout,zsto
34
35c
36c        CALL histdef(nid_hf, "phis", "Surface geop. height", "-",
37c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
38c    .                "once", zsto,zout)
39c
40c        CALL histdef(nid_hf, "aire", "Grid area", "-",
41c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
42c    .                "once", zsto,zout)
43c
44c Champs 2D:
45c
46         CALL histdef(nid_hf, "tsol", "Surface Temperature", "K",
47     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
48     .                "inst(X)", zsto,zout)
49c
50         CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa",
51     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
52     .                "inst(X)", zsto,zout)
53c
54         print*,'ATTENTION METTRE AVE(X) POUR LES PRECIPS'
55
56         CALL histdef(nid_hf, "rain", "Precipitation", "mm/d",
57     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
58     .                "inst(X)", zsto,zout)
59c
60         CALL histdef(nid_hf, "u850", "Zonal wind 850mb", "m/s",
61     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
62     .                "inst(X)", zsto,zout)
63
64         CALL histdef(nid_hf, "v850", "Meridional wind 850mb", "m/s",
65     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
66     .                "inst(X)", zsto,zout)
67c
68         CALL histdef(nid_hf, "u500", "Zonal wind 500mb", "m/s",
69     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
70     .                "inst(X)", zsto,zout)
71
72         CALL histdef(nid_hf, "v500", "Meridional wind 500mb", "m/s",
73     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
74     .                "inst(X)", zsto,zout)
75
76         CALL histdef(nid_hf, "u200", "Zonal wind 200mb", "m/s",
77     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
78     .                "inst(X)", zsto,zout)
79
80         CALL histdef(nid_hf, "v200", "Meridional wind 200mb", "m/s",
81     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
82     .                "inst(X)", zsto,zout)
83
84         CALL histdef(nid_hf, "phi500", "Geopotentiel à 500mb", "m2/s2",
85     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
86     .                "inst(X)", zsto,zout)
87
88c
89         CALL histend(nid_hf)
90
91      endif ! ok_hf
Note: See TracBrowser for help on using the repository browser.