source: dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/ini_paramLMDZ_phy.h @ 3825

Last change on this file since 3825 was 3825, checked in by ymipsl, 10 years ago

Reorganize geometry and grid modules. Prepare physics for unstructutured grid support. Simplify initialization of physics from dynamic.
Compiled only with dynd3dmem, but not tested for moment.

YM
  • Property svn:eol-style set to native
  • Property svn:executable set to *
File size: 4.8 KB
Line 
1!IM    Implemente en modes sequentiel et parallele
2
3       CALL gather(lat_degrees,rlat_glo)
4       CALL bcast(rlat_glo)
5       CALL gather(lon_degrees,rlon_glo)
6       CALL bcast(rlon_glo)
7
8!$OMP MASTER
9      if (is_mpi_root) then
10!
11       zstophy = dtime
12       zout = mth_len*un_jour
13!
14       idayref = day_ref
15       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
16!
17       CALL grid1dTo2d_glo(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 grid1dTo2d_glo(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.