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

Annotation of /trunk/phylmd/phyredem0.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 304 - (hide annotations)
Thu Sep 6 15:51:09 2018 UTC (5 years, 8 months ago) by guez
File size: 7133 byte(s)
Variable fevap of physiq is not used. Remove it from physiq and from
the restart file. Remove the corresponding argument evap of
pbl_surface.

Use directly yqsurf instead of qairsol in pbl_surface.

1 guez 157 module phyredem0_m
2    
3     IMPLICIT NONE
4    
5     INTEGER ncid_restartphy
6    
7     contains
8    
9 guez 202 SUBROUTINE phyredem0
10 guez 157
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 guez 202 use conf_gcm_m, only: nday, lmt_pas
19 guez 157 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 guez 191 use phyetat0_m, only: rlat, rlon, itau_phy
26 guez 157
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 IF 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