source: LMDZ6/trunk/libf/phylmd/ini_histday_seri.h @ 5020

Last change on this file since 5020 was 2399, checked in by Ehouarn Millour, 9 years ago

Follow-up from commit 2395: get rid of rlon and rlat, longitude_deg and latitude_deg (from module geometry_mod) should be used instead. Longitudes and latitudes are no longer loaded from startphy.nc but inherited from dynamics (and compatibility with values in startphy.nc is checked). This will change bench results because of roundoffs differences between the two.
EM
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.7 KB
Line 
1!
2! $Id: ini_histday_seri.h 2399 2015-11-20 16:23:28Z acozic $
3!
4!ym Ne fonctionnera pas en mode parallele
5      IF (is_sequential) THEN
6     
7      IF (type_run.EQ."AMIP") THEN
8!
9       zstophy = pdtphys
10       zout = ecrit_day
11!
12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
14!
15         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,longitude_deg,zx_lon)
16         DO i = 1, nbp_lon
17            zx_lon(i,1) = longitude_deg(i+1)
18            zx_lon(i,nbp_lat) = longitude_deg(i+1)
19         ENDDO
20         DO ll=1,klev
21            znivsig(ll)=REAL(ll)
22         ENDDO
23         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,latitude_deg,zx_lat)
24!
25         imin_debut=1 
26         nbpti=1
27         jmin_debut=1 
28         nbptj=1
29!
30         CALL histbeg("histday_seri.nc",  &
31                       nbp_lon,zx_lon(:,1),nbp_lat,zx_lat(1,:), &
32                       imin_debut,nbpti,jmin_debut,nbptj, &
33                       itau_phy, zjulian, dtime, &
34                       nhori, nid_day_seri)
35!
36         CALL histvert(nid_day_seri, "presnivs",  &
37                      "Vertical levels","mb", &
38                       klev, presnivs/100., nvert)
39!
40         CALL histdef(nid_day_seri, "bilTOA",  &
41                      "Net radiation at model top", "W/m2", &
42                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
43                      "ave(X)", zstophy,zout)
44!
45         CALL histdef(nid_day_seri, "bils",  &
46                      "Net downward energy flux at surface","W/m2", &
47                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
48                      "ave(X)", zstophy,zout)
49!
50         CALL histdef(nid_day_seri, "ecin",  &
51                      "Total kinetic energy (per unit area)","J/m2", &
52                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
53                      "ave(X)", zstophy,zout)
54!
55!IM 151004 BEG
56         IF(1.EQ.0) THEN
57!
58         CALL histdef(nid_day_seri, "momang",  &
59                     "Total relative angular momentum (per unit area)", &
60                     "kg/s", &
61                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
62                      "ave(X)", zstophy,zout)
63!
64         CALL histdef(nid_day_seri, "frictor",  &
65                     "Friction torque (per unit area)", "N/m", &
66                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
67                      "ave(X)", zstophy,zout)
68!
69         CALL histdef(nid_day_seri, "mountor",  &
70                     "Mountain torque (per unit area)", "N/m", &
71                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
72                      "ave(X)", zstophy,zout)
73!
74         ENDIF !(1.EQ.0) THEN
75!
76         CALL histdef(nid_day_seri, "momang",  &
77                     "Axial angular momentum (per unit area)", &
78                     "kg/s", &
79                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
80                      "ave(X)", zstophy,zout)
81!
82         CALL histdef(nid_day_seri, "torsfc",  &
83              "Total surface torque (including mountain torque)", "N/m", &
84                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
85                      "ave(X)", zstophy,zout)
86!
87!IM 151004 END       
88!
89         CALL histdef(nid_day_seri, "tamv",  &
90                      "Temperature (mass-weighted vert. ave)", "K", &
91                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
92                      "ave(X)", zstophy,zout)
93!
94         CALL histdef(nid_day_seri, "psol",  &
95                      "Surface pressure", "Pa", &
96                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
97                      "ave(X)", zstophy,zout)
98!
99         CALL histdef(nid_day_seri, "evap",  &
100                      "Evaporation and sublimation (per unit area)",  &
101                      "kg/(m2*s)", &
102                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
103                      "ave(X)", zstophy,zout)
104!
105!          call histdef(nid_day_seri, 
106!    .         "SnowFrac", 
107!    .         "Snow-covered area ", "%", 
108!    .         nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,
109!    .         "ave(X)", zstophy,zout)
110!
111!        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!    .                "Snow Mass", "kg/m2",
115!    .                nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, 
116!    .               "ave(X)", zstophy,zout)
117!
118           call histdef(nid_day_seri,  &
119               "tsol_"//clnsurf(is_oce),  &
120               "SST over open (ice-free) ocean ", "K",   &
121               nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
122               "ave(X)", zstophy,zout)
123!
124!=================================================================
125!
126         CALL histend(nid_day_seri)
127!
128!=================================================================
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.