source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday_seri.h @ 728

Last change on this file since 728 was 723, checked in by lmdzadmin, 18 years ago

On passe a des ecrit_ins, ecrit_day, etc en nombre de jours (REAL)
On lit frequence ecriture traceurs ecrit_trac dans physiq.def
Correction petits pbs ini_histrac.h
IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.4 KB
Line 
1c
2c $Header$
3c
4       IF (type_run.EQ."AMIP") THEN
5c
6       zstophy = dtime
7       zout = dtime * ecrit_day
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)
21c
22         imin_debut=1 
23         nbpti=1
24         jmin_debut=1 
25         nbptj=1
26c
27         CALL histbeg("histday_seri.nc", 
28     .                 iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
29     .                 imin_debut,nbpti,jmin_debut,nbptj,
30     .                 itau_phy, zjulian, dtime,
31     .                 nhori, nid_day_seri)
32c
33         CALL histvert(nid_day_seri, "presnivs", 
34     .                "Vertical levels","mb",
35     .                 klev, presnivs/100., nvert)
36c
37         CALL histdef(nid_day_seri, "bilTOA", 
38     .                "Net radiation at model top", "W/m2",
39     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
40     .                "ave(X)", zstophy,zout)
41c
42         CALL histdef(nid_day_seri, "bils", 
43     .                "Net downward energy flux at surface","W/m2",
44     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
45     .                "ave(X)", zstophy,zout)
46c
47         CALL histdef(nid_day_seri, "ecin", 
48     .                "Total kinetic energy (per unit area)","J/m2",
49     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
50     .                "ave(X)", zstophy,zout)
51c
52cIM 151004 BEG
53         IF(1.EQ.0) THEN
54c
55         CALL histdef(nid_day_seri, "momang", 
56     .               "Total relative angular momentum (per unit area)",
57     .               "kg/s",
58     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
59     .                "ave(X)", zstophy,zout)
60c
61         CALL histdef(nid_day_seri, "frictor", 
62     .               "Friction torque (per unit area)", "N/m",
63     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
64     .                "ave(X)", zstophy,zout)
65c
66         CALL histdef(nid_day_seri, "mountor", 
67     .               "Mountain torque (per unit area)", "N/m",
68     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
69     .                "ave(X)", zstophy,zout)
70c
71         ENDIF !(1.EQ.0) THEN
72c
73         CALL histdef(nid_day_seri, "momang", 
74     .               "Axial angular momentum (per unit area)",
75     .               "kg/s",
76     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
77     .                "ave(X)", zstophy,zout)
78c
79         CALL histdef(nid_day_seri, "torsfc", 
80     .        "Total surface torque (including mountain torque)", "N/m",
81     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
82     .                "ave(X)", zstophy,zout)
83c
84cIM 151004 END       
85c
86         CALL histdef(nid_day_seri, "tamv", 
87     .                "Temperature (mass-weighted vert. ave)", "K",
88     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
89     .                "ave(X)", zstophy,zout)
90c
91         CALL histdef(nid_day_seri, "psol", 
92     .                "Surface pressure", "Pa",
93     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
94     .                "ave(X)", zstophy,zout)
95c
96         CALL histdef(nid_day_seri, "evap", 
97     .                "Evaporation and sublimation (per unit area)", 
98     .                "kg/(m2*s)",
99     .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
100     .                "ave(X)", zstophy,zout)
101c
102c          call histdef(nid_day_seri, 
103c    .         "SnowFrac", 
104c    .         "Snow-covered area ", "%", 
105c    .         iim,jjmp1,nhori, 1,1,1, -99, 32,
106c    .         "ave(X)", zstophy,zout)
107c
108c        CALL histdef(nid_day_seri, "snow_depth", 
109cIM 080904  .                "Snow Depth (water equivalent)", "m",
110cIM 191104  .                "Snow Depth (water equivalent)", "kg/m2",
111c    .                "Snow Mass", "kg/m2",
112c    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
113c    .               "ave(X)", zstophy,zout)
114c
115           call histdef(nid_day_seri, 
116     .         "tsol_"//clnsurf(is_oce),
117     .         "SST over open (ice-free) ocean ", "K", 
118     .         iim,jjmp1,nhori, 1,1,1, -99, 32,
119     .         "ave(X)", zstophy,zout)
120c
121c=================================================================
122c
123         CALL histend(nid_day_seri)
124c
125c=================================================================
126      ENDIF ! fin de test sur type_run.EQ.AMIP
Note: See TracBrowser for help on using the repository browser.