/[lmdze]/trunk/phylmd/phyredem0.f
ViewVC logotype

Contents of /trunk/phylmd/phyredem0.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 308 - (show annotations)
Tue Sep 18 15:14:40 2018 UTC (5 years, 7 months ago) by guez
File size: 7133 byte(s)
In procedure calcul_fluxs, rename coef1lay to cdragh (following
LMDZ). Use named constant min_wind_speed (following LMDZ).

Change name of NetCDF variable soll to CF standard rls.

1 module phyredem0_m
2
3 IMPLICIT NONE
4
5 INTEGER ncid_restartphy
6
7 contains
8
9 SUBROUTINE phyredem0
10
11 ! From phylmd/phyredem.F, version 1.3, 2005/05/25 13:10:09
12 ! Author: Z. X. Li (LMD/CNRS)
13 ! Date: 1993/08/18
14
15 ! Objet : \'ecriture de l'\'etat de d\'emarrage ou red\'emarrage
16 ! pour la physique
17
18 use conf_gcm_m, only: nday, lmt_pas
19 USE dimphy, ONLY: klev, klon
20 USE dimsoil, ONLY: nsoilmx
21 USE indicesol, ONLY: nbsrf
22 USE netcdf, ONLY: nf90_clobber, nf90_global, nf90_float
23 USE netcdf95, ONLY: nf95_create, nf95_put_att, nf95_def_dim, &
24 nf95_def_var, nf95_enddef, nf95_put_var
25 use phyetat0_m, only: rlat, rlon, itau_phy
26
27 ! Local:
28
29 INTEGER idim2, idim3, dimid_nbsrf, dimid_nsoilmx
30 integer varid, varid_rlon, varid_rlat
31
32 !------------------------------------------------------------
33
34 PRINT *, 'Call sequence information: phyredem0'
35 CALL nf95_create("restartphy.nc", nf90_clobber, ncid_restartphy)
36
37 call nf95_put_att(ncid_restartphy, nf90_global, 'title', &
38 'start file for the physics code')
39 call nf95_put_att(ncid_restartphy, nf90_global, "itau_phy", &
40 itau_phy + nday * lmt_pas)
41
42 call nf95_def_dim(ncid_restartphy, 'points_physiques', klon, idim2)
43 call nf95_def_dim(ncid_restartphy, 'klev', klev, idim3)
44 call nf95_def_dim(ncid_restartphy, 'nbsrf', nbsrf, dimid_nbsrf)
45 call nf95_def_dim(ncid_restartphy, 'nsoilmx', nsoilmx, dimid_nsoilmx)
46
47 call nf95_def_var(ncid_restartphy, 'longitude', nf90_float, idim2, &
48 varid_rlon)
49 call nf95_def_var(ncid_restartphy, 'latitude', nf90_float, idim2, &
50 varid_rlat)
51
52 call nf95_def_var(ncid_restartphy, 'masque', nf90_float, idim2, varid)
53 call nf95_put_att(ncid_restartphy, varid, 'title', 'masque terre mer')
54
55 ! Fractions de chaque sous-surface
56
57 call nf95_def_var(ncid_restartphy, 'FTER', nf90_float, idim2, varid)
58 call nf95_put_att(ncid_restartphy, varid, 'title', 'fraction de continent')
59
60 call nf95_def_var(ncid_restartphy, 'FLIC', nf90_float, idim2, varid)
61 call nf95_put_att(ncid_restartphy, varid, 'title', &
62 'fraction glace de terre')
63
64 call nf95_def_var(ncid_restartphy, 'FOCE', nf90_float, idim2, varid)
65 call nf95_put_att(ncid_restartphy, varid, 'title', 'fraction ocean')
66
67 call nf95_def_var(ncid_restartphy, 'FSIC', nf90_float, idim2, varid)
68 call nf95_put_att(ncid_restartphy, varid, 'title', 'fraction glace mer')
69
70 call nf95_def_var(ncid_restartphy, 'TS', nf90_float, &
71 (/idim2, dimid_nbsrf/), varid)
72 call nf95_put_att(ncid_restartphy, varid, 'title', 'surface temperature')
73
74 call nf95_def_var(ncid_restartphy, 'Tsoil', nf90_float, &
75 (/idim2, dimid_nsoilmx, dimid_nbsrf/), varid)
76 call nf95_put_att(ncid_restartphy, varid, 'title', 'soil temperature')
77
78 call nf95_def_var(ncid_restartphy, 'QS', nf90_float, &
79 (/idim2, dimid_nbsrf/), varid)
80 call nf95_put_att(ncid_restartphy, varid, 'title', 'Humidite de surface')
81
82 call nf95_def_var(ncid_restartphy, 'QSOL', nf90_float, idim2, varid)
83 call nf95_put_att(ncid_restartphy, varid, 'title', 'Eau dans le sol (mm)')
84
85 call nf95_def_var(ncid_restartphy, 'ALBE', nf90_float, &
86 (/idim2, dimid_nbsrf/), varid)
87 call nf95_put_att(ncid_restartphy, varid, 'title', 'albedo de surface')
88
89 call nf95_def_var(ncid_restartphy, 'SNOW', nf90_float, &
90 (/idim2, dimid_nbsrf/), varid)
91 call nf95_put_att(ncid_restartphy, varid, 'title', 'Neige de surface')
92
93 call nf95_def_var(ncid_restartphy, 'RADS', nf90_float, idim2, varid)
94 call nf95_put_att(ncid_restartphy, varid, 'title', &
95 'Rayonnement net a la surface')
96
97 call nf95_def_var(ncid_restartphy, 'solsw', nf90_float, idim2, varid)
98 call nf95_put_att(ncid_restartphy, varid, 'title', &
99 'Rayonnement solaire a la surface')
100
101 call nf95_def_var(ncid_restartphy, 'sollw', nf90_float, idim2, varid)
102 call nf95_put_att(ncid_restartphy, varid, 'title', &
103 'Rayonnement IR a la surface')
104
105 call nf95_def_var(ncid_restartphy, 'fder', nf90_float, idim2, varid)
106 call nf95_put_att(ncid_restartphy, varid, 'title', 'Derive de flux')
107
108 call nf95_def_var(ncid_restartphy, 'rain_f', nf90_float, idim2, varid)
109 call nf95_put_att(ncid_restartphy, varid, 'title', 'precipitation liquide')
110
111 call nf95_def_var(ncid_restartphy, 'snow_f', nf90_float, idim2, varid)
112 call nf95_put_att(ncid_restartphy, varid, 'title', 'precipitation solide')
113
114 call nf95_def_var(ncid_restartphy, 'RUG', nf90_float, &
115 (/idim2, dimid_nbsrf/), varid)
116 call nf95_put_att(ncid_restartphy, varid, 'title', 'rugosite de surface')
117
118 call nf95_def_var(ncid_restartphy, 'AGESNO', nf90_float, &
119 (/idim2, dimid_nbsrf/), varid)
120 call nf95_put_att(ncid_restartphy, varid, 'title', &
121 'Age de la neige surface')
122
123 call nf95_def_var(ncid_restartphy, 'ZMEA', nf90_float, idim2, varid)
124 call nf95_def_var(ncid_restartphy, 'ZSTD', nf90_float, idim2, varid)
125 call nf95_def_var(ncid_restartphy, 'ZSIG', nf90_float, idim2, varid)
126 call nf95_def_var(ncid_restartphy, 'ZGAM', nf90_float, idim2, varid)
127 call nf95_def_var(ncid_restartphy, 'ZTHE', nf90_float, idim2, varid)
128 call nf95_def_var(ncid_restartphy, 'ZPIC', nf90_float, idim2, varid)
129 call nf95_def_var(ncid_restartphy, 'ZVAL', nf90_float, idim2, varid)
130 call nf95_def_var(ncid_restartphy, 'TANCIEN', nf90_float, &
131 (/idim2, idim3/), varid)
132 call nf95_def_var(ncid_restartphy, 'QANCIEN', nf90_float, &
133 (/idim2, idim3/), varid)
134
135 call nf95_def_var(ncid_restartphy, 'RUGMER', nf90_float, idim2, varid)
136 call nf95_put_att(ncid_restartphy, varid, 'title', &
137 'Longueur de rugosite sur mer')
138
139 call nf95_def_var(ncid_restartphy, 'CLWCON', nf90_float, idim2, varid)
140 call nf95_put_att(ncid_restartphy, varid, 'title', 'Eau liquide convective')
141
142 call nf95_def_var(ncid_restartphy, 'RNEBCON', nf90_float, idim2, varid)
143 call nf95_put_att(ncid_restartphy, varid, 'title', 'Nebulosite convective')
144
145 call nf95_def_var(ncid_restartphy, 'RATQS', nf90_float, idim2, varid)
146 call nf95_put_att(ncid_restartphy, varid, 'title', 'Ratqs')
147
148 call nf95_def_var(ncid_restartphy, 'RUNOFFLIC0', nf90_float, idim2, varid)
149 call nf95_put_att(ncid_restartphy, varid, 'title', 'Runofflic0')
150
151 call nf95_def_var(ncid_restartphy, 'sig1', nf90_float, (/idim2, idim3/), &
152 varid)
153 call nf95_put_att(ncid_restartphy, varid, 'long_name', &
154 'section adiabatic updraft')
155
156 call nf95_def_var(ncid_restartphy, 'w01', nf90_float, (/idim2, idim3/), &
157 varid)
158 call nf95_put_att(ncid_restartphy, varid, 'long_name', &
159 'vertical velocity within adiabatic updraft')
160
161 call nf95_def_var(ncid_restartphy, 'trs', nf90_float, idim2, varid)
162 call nf95_put_att(ncid_restartphy, varid, 'long_name', &
163 'radon concentation in soil')
164
165 call nf95_enddef(ncid_restartphy)
166
167 call nf95_put_var(ncid_restartphy, varid_rlon, rlon)
168 call nf95_put_var(ncid_restartphy, varid_rlat, rlat)
169
170 END SUBROUTINE phyredem0
171
172 end module phyredem0_m

  ViewVC Help
Powered by ViewVC 1.1.21