Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/ener_conserv.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ener_conserv.F90
r5144 r5158 179 179 fluxt(:, klev + 1) = 0. 180 180 181 dok = klev, 1, -1181 DO k = klev, 1, -1 182 182 fluxu(:, k) = fluxu(:, k + 1) + masse(:, k) * d_u(:, k) 183 183 fluxv(:, k) = fluxv(:, k + 1) + masse(:, k) * d_v(:, k) … … 189 189 dddt(:, 1) = (exner(:, 1) - 1.) * fluxt(:, 1) 190 190 191 dok = 2, klev191 DO k = 2, klev 192 192 dddu(:, k) = (zu(:, k) - zu(:, k - 1)) * fluxu(:, k) 193 193 dddv(:, k) = (zv(:, k) - zv(:, k - 1)) * fluxv(:, k) … … 198 198 dddt(:, klev + 1) = 0. 199 199 200 dok = 1, klev200 DO k = 1, klev 201 201 d_t_ech(:, k) = -(rcpd * (dddt(:, k) + dddt(:, k + 1))) / (2. * rcpd * masse(:, k)) 202 202 d_t_ec(:, k) = -(dddu(:, k) + dddu(:, k + 1) + dddv(:, k) + dddv(:, k + 1)) / (2. * rcpd * masse(:, k)) + d_t_ech(:, k)
Note: See TracChangeset
for help on using the changeset viewer.