source: lmdz_wrf/WRFV3/lmdz/ini_histday_seri.h @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 6.2 KB
Line 
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!IM 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!IM 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!IM 080904  .                "Snow Depth (water equivalent)", "m",
113!IM 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
Note: See TracBrowser for help on using the repository browser.