Changeset 1340 for palm/trunk/SOURCE/plant_canopy_model.f90
- Timestamp:
- Mar 25, 2014 7:45:13 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model.f90
r1321 r1340 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants defined as wp-kind 23 23 ! 24 24 ! Former revisions: … … 100 100 v(k,j+1,i) + & 101 101 v(k,j+1,i-1) ) & 102 / 4.0 )**2+ &102 / 4.0_wp )**2 + & 103 103 ( ( w(k-1,j,i-1) + & 104 104 w(k-1,j,i) + & 105 105 w(k,j,i-1) + & 106 106 w(k,j,i) ) & 107 / 4.0 )**2 )* &107 / 4.0_wp )**2 ) * & 108 108 u(k,j,i) 109 109 ENDDO … … 123 123 u(k,j,i) + & 124 124 u(k,j,i+1) ) & 125 / 4.0 )**2+ &125 / 4.0_wp )**2 + & 126 126 v(k,j,i)**2 + & 127 127 ( ( w(k-1,j-1,i) + & … … 129 129 w(k,j-1,i) + & 130 130 w(k,j,i) ) & 131 / 4.0 )**2 )* &131 / 4.0_wp )**2 ) * & 132 132 v(k,j,i) 133 133 ENDDO … … 147 147 u(k+1,j,i) + & 148 148 u(k+1,j,i+1) ) & 149 / 4.0 )**2+ &149 / 4.0_wp )**2 + & 150 150 ( ( v(k,j,i) + & 151 151 v(k,j+1,i) + & 152 152 v(k+1,j,i) + & 153 153 v(k+1,j+1,i) ) & 154 / 4.0 )**2+ &154 / 4.0_wp )**2 + & 155 155 w(k,j,i)**2 ) * & 156 156 w(k,j,i) … … 183 183 SQRT( ( ( u(k,j,i) + & 184 184 u(k,j,i+1) ) & 185 / 2.0 )**2+ &185 / 2.0_wp )**2 + & 186 186 ( ( v(k,j,i) + & 187 187 v(k,j+1,i) ) & 188 / 2.0 )**2+ &188 / 2.0_wp )**2 + & 189 189 ( ( w(k-1,j,i) + & 190 190 w(k,j,i) ) & 191 / 2.0 )**2 )* &191 / 2.0_wp )**2 ) * & 192 192 ( q(k,j,i) - sls(k,j,i) ) 193 193 ENDDO … … 202 202 DO k = nzb_s_inner(j,i)+1, pch_index 203 203 tend(k,j,i) = tend(k,j,i) - & 204 2.0 * cdc(k,j,i) * lad_s(k,j,i) * &204 2.0_wp * cdc(k,j,i) * lad_s(k,j,i) * & 205 205 SQRT( ( ( u(k,j,i) + & 206 206 u(k,j,i+1) ) & 207 / 2.0 )**2+ &207 / 2.0_wp )**2 + & 208 208 ( ( v(k,j,i) + & 209 209 v(k,j+1,i) ) & 210 / 2.0 )**2+ &210 / 2.0_wp )**2 + & 211 211 ( ( w(k,j,i) + & 212 212 w(k+1,j,i) ) & 213 / 2.0 )**2 )* &213 / 2.0_wp )**2 ) * & 214 214 e(k,j,i) 215 215 ENDDO … … 267 267 v(k,j+1,i) + & 268 268 v(k,j+1,i-1) ) & 269 / 4.0 )**2+ &269 / 4.0_wp )**2 + & 270 270 ( ( w(k-1,j,i-1) + & 271 271 w(k-1,j,i) + & 272 272 w(k,j,i-1) + & 273 273 w(k,j,i) ) & 274 / 4.0 )**2 ) *&274 / 4.0_wp )**2 ) * & 275 275 u(k,j,i) 276 276 ENDDO … … 286 286 u(k,j,i) + & 287 287 u(k,j,i+1) ) & 288 / 4.0 )**2+ &288 / 4.0_wp )**2 + & 289 289 v(k,j,i)**2 + & 290 290 ( ( w(k-1,j-1,i) + & … … 292 292 w(k,j-1,i) + & 293 293 w(k,j,i) ) & 294 / 4.0 )**2 ) *&294 / 4.0_wp )**2 ) * & 295 295 v(k,j,i) 296 296 ENDDO … … 306 306 u(k+1,j,i) + & 307 307 u(k+1,j,i+1) ) & 308 / 4.0 )**2+ &308 / 4.0_wp )**2 + & 309 309 ( ( v(k,j,i) + & 310 310 v(k,j+1,i) + & 311 311 v(k+1,j,i) + & 312 312 v(k+1,j+1,i) ) & 313 / 4.0 )**2+ &313 / 4.0_wp )**2 + & 314 314 w(k,j,i)**2 ) * & 315 315 w(k,j,i) … … 336 336 SQRT( ( ( u(k,j,i) + & 337 337 u(k,j,i+1) ) & 338 / 2.0 )**2+ &338 / 2.0_wp )**2 + & 339 339 ( ( v(k,j,i) + & 340 340 v(k,j+1,i) ) & 341 / 2.0 )**2+ &341 / 2.0_wp )**2 + & 342 342 ( ( w(k-1,j,i) + & 343 343 w(k,j,i) ) & 344 / 2.0 )**2 )* &344 / 2.0_wp )**2 ) * & 345 345 ( q(k,j,i) - sls(k,j,i) ) 346 346 ENDDO … … 351 351 DO k = nzb_s_inner(j,i)+1, pch_index 352 352 tend(k,j,i) = tend(k,j,i) - & 353 2.0 * cdc(k,j,i) * lad_s(k,j,i) * &353 2.0_wp * cdc(k,j,i) * lad_s(k,j,i) * & 354 354 SQRT( ( ( u(k,j,i) + & 355 355 u(k,j,i+1) ) & 356 / 2.0 )**2+ &356 / 2.0_wp )**2 + & 357 357 ( ( v(k,j,i) + & 358 358 v(k,j+1,i) ) & 359 / 2.0 )**2+ &359 / 2.0_wp )**2 + & 360 360 ( ( w(k,j,i) + & 361 361 w(k+1,j,i) ) & 362 / 2.0 )**2 )* &362 / 2.0_wp )**2 ) * & 363 363 e(k,j,i) 364 364 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.