source: LMDZ5/branches/AI-cosp/libf/phylmd/ini_paramLMDZ_phy.h @ 3717

Last change on this file since 3717 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:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 4.8 KB
Line 
1!IM    Implemente en modes sequentiel et parallele
2
3       CALL gather(latitude_deg,rlat_glo)
4       CALL bcast(rlat_glo)
5       CALL gather(longitude_deg,rlon_glo)
6       CALL bcast(rlon_glo)
7
8!$OMP MASTER
9      if (is_mpi_root) then
10!
11       zstophy = pdtphys
12       zout = mth_len*un_jour
13!
14       idayref = day_ref
15       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
16!
17       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon_glo,zx_lon)
18       if (nbp_lon.gt.1) then
19       DO i = 1, nbp_lon
20         zx_lon(i,1) = rlon_glo(i+1)
21         zx_lon(i,nbp_lat) = rlon_glo(i+1)
22       ENDDO
23       endif
24       CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat_glo,zx_lat)
25!
26       CALL histbeg("paramLMDZ_phy.nc",  &
27                       np,zx_lon(np:np,1), np,zx_lat(1,np:np), &
28                       1,1,1,1, &
29                       itau_phy, zjulian, dtime, &
30                       nhori, nid_ctesGCM)
31!
32       CALL histdef(nid_ctesGCM, "R_ecc",  &
33                      "Excentricite","-", &
34                      1,1,nhori, 1,1,1, -99, 32, &
35                      "ave", zstophy,zout)
36!
37       CALL histdef(nid_ctesGCM, "R_peri",  &
38                      "Equinoxe","-", &
39                      1,1,nhori, 1,1,1, -99, 32, &
40                      "ave", zstophy,zout)
41!
42       CALL histdef(nid_ctesGCM, "R_incl",  &
43                      "Inclinaison","deg", &
44                      1,1,nhori, 1,1,1, -99, 32, &
45                      "ave", zstophy,zout)
46!
47       CALL histdef(nid_ctesGCM, "solaire",  &
48                      "Constante solaire","W/m2", &
49                      1,1,nhori, 1,1,1, -99, 32, &
50                      "ave", zstophy,zout)
51!
52       CALL histdef(nid_ctesGCM, "co2_ppm",  &
53                      "Concentration du CO2", "ppm", &
54                      1,1,nhori, 1,1,1, -99, 32,  &
55                      "ave(X)", zstophy,zout)
56!
57       CALL histdef(nid_ctesGCM, "CH4_ppb",  &
58                      "Concentration du CH4", "ppb", &
59                      1,1,nhori, 1,1,1, -99, 32,  &
60                      "ave(X)", zstophy,zout)
61!
62       CALL histdef(nid_ctesGCM, "N2O_ppb", &
63                      "Concentration du N2O", "ppb", &
64                      1,1,nhori, 1,1,1, -99, 32, &
65                      "ave(X)", zstophy,zout)
66!
67       CALL histdef(nid_ctesGCM, "CFC11_ppt", &
68                      "Concentration du CFC11", "ppt", &
69                      1,1,nhori, 1,1,1, -99, 32, &
70                      "ave(X)", zstophy,zout)
71!
72       CALL histdef(nid_ctesGCM, "CFC12_ppt", &
73                      "Concentration du CFC12", "ppt", &
74                      1,1,nhori, 1,1,1, -99, 32, &
75                      "ave(X)", zstophy,zout)
76!
77       CALL histdef(nid_ctesGCM, "bils", &
78                      "Surface total heat flux", "W m-2", &
79                      1,1,nhori, 1,1,1, -99, 32, &
80                      "ave", zstophy,zout)
81!
82       CALL histdef(nid_ctesGCM, "evap", &
83                      "Evaporation", "kg m-2 s-1", &
84                      1,1,nhori, 1,1,1, -99, 32, &
85                      "ave", zstophy,zout)
86!
87       CALL histdef(nid_ctesGCM, "evap_land", &
88                      "Land evaporation", "kg m-2 s-1", &
89                      1,1,nhori, 1,1,1, -99, 32, &
90                      "ave", zstophy,zout)
91!
92       CALL histdef(nid_ctesGCM, "flat", &
93                      "Latent heat flux", "W m-2", &
94                      1,1,nhori, 1,1,1, -99, 32, &
95                      "ave", zstophy,zout)
96!
97       CALL histdef(nid_ctesGCM, "nettop0", &
98                      "Clear sky net downward radiatif flux at TOA",  &
99                      "W m-2", &
100                      1,1,nhori, 1,1,1, -99, 32, &
101                      "ave", zstophy,zout)
102!
103       CALL histdef(nid_ctesGCM, "nettop", &
104                      "Net downward radiatif flux at TOA", "W m-2", &
105                      1,1,nhori, 1,1,1, -99, 32, &
106                      "ave", zstophy,zout)
107!
108       CALL histdef(nid_ctesGCM, "precip", &
109                      "Total precipitation (liq+sol)", "kg m-2 s-1", &
110                      1,1,nhori, 1,1,1, -99, 32, &
111                      "ave", zstophy,zout)
112!
113       CALL histdef(nid_ctesGCM, "tsol", &
114                      "Surface temperature", "K", &
115                      1,1,nhori, 1,1,1, -99, 32, &
116                      "ave", zstophy,zout)
117!
118       CALL histdef(nid_ctesGCM, "t2m", &
119                      "Temperature at 2m", "K", &
120                      1,1,nhori, 1,1,1, -99, 32, &
121                      "ave", zstophy,zout)
122!
123       CALL histdef(nid_ctesGCM, "prw", &
124                      "Precipitable water", "kg m-2", &
125                      1,1,nhori, 1,1,1, -99, 32, &
126                      "ave", zstophy,zout) 
127!=================================================================
128!
129       CALL histend(nid_ctesGCM)
130       
131       endif !(is_mpi_root)
132!$OMP END MASTER
133!=================================================================
Note: See TracBrowser for help on using the repository browser.