source: LMDZ5/branches/testing/libf/phylmd/ini_histday_seri.h @ 5455

Last change on this file since 5455 was 2435, checked in by Laurent Fairhead, 9 years ago

Merged trunk changes r2396:2434 into testing branch

  • 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
RevLine 
[1910]1!
2! $Id: ini_histday_seri.h 2435 2016-01-28 16:02:13Z aborella $
3!
4!ym Ne fonctionnera pas en mode parallele
[776]5      IF (is_sequential) THEN
[766]6     
[644]7      IF (type_run.EQ."AMIP") THEN
[1910]8!
[2408]9       zstophy = pdtphys
[684]10       zout = ecrit_day
[1910]11!
[644]12         idayref = day_ref
13         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
[1910]14!
[2435]15         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,longitude_deg,zx_lon)
[2408]16         DO i = 1, nbp_lon
[2435]17            zx_lon(i,1) = longitude_deg(i+1)
18            zx_lon(i,nbp_lat) = longitude_deg(i+1)
[644]19         ENDDO
20         DO ll=1,klev
[1403]21            znivsig(ll)=REAL(ll)
[644]22         ENDDO
[2435]23         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,latitude_deg,zx_lat)
[1910]24!
[644]25         imin_debut=1 
26         nbpti=1
27         jmin_debut=1 
28         nbptj=1
[1910]29!
30         CALL histbeg("histday_seri.nc",  &
[2408]31                       nbp_lon,zx_lon(:,1),nbp_lat,zx_lat(1,:), &
[1910]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", &
[2408]42                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
[1910]43                      "ave(X)", zstophy,zout)
44!
45         CALL histdef(nid_day_seri, "bils",  &
46                      "Net downward energy flux at surface","W/m2", &
[2408]47                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
[1910]48                      "ave(X)", zstophy,zout)
49!
50         CALL histdef(nid_day_seri, "ecin",  &
51                      "Total kinetic energy (per unit area)","J/m2", &
[2408]52                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]53                      "ave(X)", zstophy,zout)
54!
55!IM 151004 BEG
[644]56         IF(1.EQ.0) THEN
[1910]57!
58         CALL histdef(nid_day_seri, "momang",  &
59                     "Total relative angular momentum (per unit area)", &
60                     "kg/s", &
[2408]61                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]62                      "ave(X)", zstophy,zout)
63!
64         CALL histdef(nid_day_seri, "frictor",  &
65                     "Friction torque (per unit area)", "N/m", &
[2408]66                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]67                      "ave(X)", zstophy,zout)
68!
69         CALL histdef(nid_day_seri, "mountor",  &
70                     "Mountain torque (per unit area)", "N/m", &
[2408]71                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]72                      "ave(X)", zstophy,zout)
73!
[644]74         ENDIF !(1.EQ.0) THEN
[1910]75!
76         CALL histdef(nid_day_seri, "momang",  &
77                     "Axial angular momentum (per unit area)", &
78                     "kg/s", &
[2408]79                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]80                      "ave(X)", zstophy,zout)
81!
82         CALL histdef(nid_day_seri, "torsfc",  &
83              "Total surface torque (including mountain torque)", "N/m", &
[2408]84                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]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", &
[2408]91                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]92                      "ave(X)", zstophy,zout)
93!
94         CALL histdef(nid_day_seri, "psol",  &
95                      "Surface pressure", "Pa", &
[2408]96                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
[1910]97                      "ave(X)", zstophy,zout)
98!
99         CALL histdef(nid_day_seri, "evap",  &
100                      "Evaporation and sublimation (per unit area)",  &
101                      "kg/(m2*s)", &
[2408]102                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
[1910]103                      "ave(X)", zstophy,zout)
104!
105!          call histdef(nid_day_seri, 
106!    .         "SnowFrac", 
107!    .         "Snow-covered area ", "%", 
[2408]108!    .         nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,
[1910]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",
[2408]115!    .                nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, 
[1910]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",   &
[2408]121               nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
[1910]122               "ave(X)", zstophy,zout)
123!
124!=================================================================
125!
[644]126         CALL histend(nid_day_seri)
[1910]127!
128!=================================================================
[644]129      ENDIF ! fin de test sur type_run.EQ.AMIP
[766]130     
[776]131      ENDIF ! is_sequential
Note: See TracBrowser for help on using the repository browser.