My Project
 All Classes Files Functions Variables Macros
ini_histday_seri.h
Go to the documentation of this file.
1 c
2 c $Id: ini_histday_seri.h 1403 2010-07-01 09:02:53Z fairhead $
3 c
4 cym Ne fonctionnera pas en mode parallele
5  IF (is_sequential) THEN
6 
7  IF (type_run.EQ."AMIP") THEN
8 c
9  zstophy = dtime
10  zout = ecrit_day
11 c
12  idayref = day_ref
13  CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14 c
15  CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
16  DO i = 1, iim
17  zx_lon(i,1) = rlon(i+1)
18  zx_lon(i,jjmp1) = rlon(i+1)
19  ENDDO
20  DO ll=1,klev
21  znivsig(ll)=REAL(ll)
22  ENDDO
23  CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
24 c
25  imin_debut=1
26  nbpti=1
27  jmin_debut=1
28  nbptj=1
29 c
30  CALL histbeg("histday_seri.nc",
31  . iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
32  . imin_debut,nbpti,jmin_debut,nbptj,
33  . itau_phy, zjulian, dtime,
34  . nhori, nid_day_seri)
35 c
36  CALL histvert(nid_day_seri, "presnivs",
37  . "Vertical levels","mb",
38  . klev, presnivs/100., nvert)
39 c
40  CALL histdef(nid_day_seri, "bilTOA",
41  . "Net radiation at model top", "W/m2",
42  . iim,jjmp1,nhori, 1,1,1, -99, 32,
43  . "ave(X)", zstophy,zout)
44 c
45  CALL histdef(nid_day_seri, "bils",
46  . "Net downward energy flux at surface","W/m2",
47  . iim,jjmp1,nhori, 1,1,1, -99, 32,
48  . "ave(X)", zstophy,zout)
49 c
50  CALL histdef(nid_day_seri, "ecin",
51  . "Total kinetic energy (per unit area)","J/m2",
52  . iim,jjmp1,nhori, 1,1,1, -99, 32,
53  . "ave(X)", zstophy,zout)
54 c
55 cIM 151004 BEG
56  IF(1.EQ.0) THEN
57 c
58  CALL histdef(nid_day_seri, "momang",
59  . "Total relative angular momentum (per unit area)",
60  . "kg/s",
61  . iim,jjmp1,nhori, 1,1,1, -99, 32,
62  . "ave(X)", zstophy,zout)
63 c
64  CALL histdef(nid_day_seri, "frictor",
65  . "Friction torque (per unit area)", "N/m",
66  . iim,jjmp1,nhori, 1,1,1, -99, 32,
67  . "ave(X)", zstophy,zout)
68 c
69  CALL histdef(nid_day_seri, "mountor",
70  . "Mountain torque (per unit area)", "N/m",
71  . iim,jjmp1,nhori, 1,1,1, -99, 32,
72  . "ave(X)", zstophy,zout)
73 c
74  ENDIF !(1.EQ.0) THEN
75 c
76  CALL histdef(nid_day_seri, "momang",
77  . "Axial angular momentum (per unit area)",
78  . "kg/s",
79  . iim,jjmp1,nhori, 1,1,1, -99, 32,
80  . "ave(X)", zstophy,zout)
81 c
82  CALL histdef(nid_day_seri, "torsfc",
83  . "Total surface torque (including mountain torque)", "N/m",
84  . iim,jjmp1,nhori, 1,1,1, -99, 32,
85  . "ave(X)", zstophy,zout)
86 c
87 cIM 151004 END
88 c
89  CALL histdef(nid_day_seri, "tamv",
90  . "Temperature (mass-weighted vert. ave)", "K",
91  . iim,jjmp1,nhori, 1,1,1, -99, 32,
92  . "ave(X)", zstophy,zout)
93 c
94  CALL histdef(nid_day_seri, "psol",
95  . "Surface pressure", "Pa",
96  . iim,jjmp1,nhori, 1,1,1, -99, 32,
97  . "ave(X)", zstophy,zout)
98 c
99  CALL histdef(nid_day_seri, "evap",
100  . "Evaporation and sublimation (per unit area)",
101  . "kg/(m2*s)",
102  . iim,jjmp1,nhori, 1,1,1, -99, 32,
103  . "ave(X)", zstophy,zout)
104 c
105 c call histdef(nid_day_seri,
106 c . "SnowFrac",
107 c . "Snow-covered area ", "%",
108 c . iim,jjmp1,nhori, 1,1,1, -99, 32,
109 c . "ave(X)", zstophy,zout)
110 c
111 c CALL histdef(nid_day_seri, "snow_depth",
112 cIM 080904 . "Snow Depth (water equivalent)", "m",
113 cIM 191104 . "Snow Depth (water equivalent)", "kg/m2",
114 c . "Snow Mass", "kg/m2",
115 c . iim,jjmp1,nhori, 1,1,1, -99, 32,
116 c . "ave(X)", zstophy,zout)
117 c
118  call histdef(nid_day_seri,
119  . "tsol_"//clnsurf(is_oce),
120  . "SST over open (ice-free) ocean ", "K",
121  . iim,jjmp1,nhori, 1,1,1, -99, 32,
122  . "ave(X)", zstophy,zout)
123 c
124 c=================================================================
125 c
126  CALL histend(nid_day_seri)
127 c
128 c=================================================================
129  ENDIF ! fin de test sur type_run.EQ.AMIP
130 
131  ENDIF ! is_sequential