Changeset 87 for palm/trunk/SOURCE/data_output_profiles.f90
- Timestamp:
- May 22, 2007 3:46:47 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_profiles.f90
r4 r87 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! var_hom renamed pr_palm 7 7 ! 8 8 ! Former revisions: … … 111 111 112 112 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2 113 (/ hom_sum(nzb+8, var_hom,normalizing_region)**2 /), &113 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), & 114 114 start = (/ 1 /), count = (/ 1 /) ) 115 115 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 331 ) 116 116 117 117 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2 118 (/ hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &118 (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 119 119 start = (/ 1 /), count = (/ 1 /) ) 120 120 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 332 ) 121 121 122 122 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3 123 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 /), &123 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), & 124 124 start = (/ 1 /), count = (/ 1 /) ) 125 125 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 333 ) 126 126 127 127 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw 128 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 * &129 hom_sum(nzb+3, var_hom,normalizing_region) /), &128 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 * & 129 hom_sum(nzb+3,pr_palm,normalizing_region) /), & 130 130 start = (/ 1 /), count = (/ 1 /) ) 131 131 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 334 ) 132 132 133 133 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2 134 (/ hom_sum(nzb+8, var_hom,normalizing_region) * &135 hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &134 (/ hom_sum(nzb+8,pr_palm,normalizing_region) * & 135 hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 136 136 start = (/ 1 /), count = (/ 1 /) ) 137 137 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 335 ) 138 138 139 139 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i 140 (/ hom_sum(nzb+6, var_hom,normalizing_region) /), &140 (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), & 141 141 start = (/ 1 /), count = (/ 1 /) ) 142 142 IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 336 ) … … 237 237 CASE ( 'ts2' ) 238 238 cross_normx_factor(k,j) = & 239 ( hom_sum(nzb+3, var_hom,normalizing_region) )**2239 ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2 240 240 CASE ( 'wpt0' ) 241 241 cross_normx_factor(k,j) = & … … 243 243 CASE ( 'wsts2' ) 244 244 cross_normx_factor(k,j) = & 245 hom_sum(nzb+8, var_hom,normalizing_region) &246 * ( hom_sum(nzb+3, var_hom,normalizing_region) )**2245 hom_sum(nzb+8,pr_palm,normalizing_region) & 246 * ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2 247 247 CASE ( 'ws2' ) 248 248 cross_normx_factor(k,j) = & 249 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2249 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 250 250 CASE ( 'ws2ts' ) 251 251 cross_normx_factor(k,j) = & 252 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2 &253 * hom_sum(nzb+3, var_hom,normalizing_region)252 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 & 253 * hom_sum(nzb+3,pr_palm,normalizing_region) 254 254 CASE ( 'ws3' ) 255 255 cross_normx_factor(k,j) = & 256 ( hom_sum(nzb+8, var_hom,normalizing_region) )**3256 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3 257 257 258 258 END SELECT … … 262 262 CASE ( 'z_i' ) 263 263 cross_normy_factor(k,j) = & 264 hom_sum(nzb+6, var_hom,normalizing_region)264 hom_sum(nzb+6,pr_palm,normalizing_region) 265 265 266 266 END SELECT … … 351 351 352 352 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2 353 (/ hom_sum(nzb+8, var_hom,normalizing_region)**2 /), &353 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), & 354 354 start = (/ dopr_time_count /), & 355 355 count = (/ 1 /) ) … … 357 357 358 358 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2 359 (/ hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &359 (/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 360 360 start = (/ dopr_time_count /), & 361 361 count = (/ 1 /) ) … … 363 363 364 364 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3 365 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 /), &365 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), & 366 366 start = (/ dopr_time_count /), & 367 367 count = (/ 1 /) ) … … 369 369 370 370 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw 371 (/ hom_sum(nzb+8, var_hom,normalizing_region)**3 * &372 hom_sum(nzb+3, var_hom,normalizing_region) /), &371 (/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 * & 372 hom_sum(nzb+3,pr_palm,normalizing_region) /), & 373 373 start = (/ dopr_time_count /), & 374 374 count = (/ 1 /) ) … … 376 376 377 377 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2 378 (/ hom_sum(nzb+8, var_hom,normalizing_region) * &379 hom_sum(nzb+3, var_hom,normalizing_region)**2 /), &378 (/ hom_sum(nzb+8,pr_palm,normalizing_region) * & 379 hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), & 380 380 start = (/ dopr_time_count /), & 381 381 count = (/ 1 /) ) … … 383 383 384 384 nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i 385 (/ hom_sum(nzb+6, var_hom,normalizing_region) /), &385 (/ hom_sum(nzb+6,pr_palm,normalizing_region) /), & 386 386 start = (/ dopr_time_count /), & 387 387 count = (/ 1 /) ) … … 482 482 CASE ( 'tsw2' ) 483 483 cross_normx_factor(k,j) = & 484 ( hom_sum(nzb+11, var_hom,normalizing_region) )**2484 ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2 485 485 CASE ( 'wpt0' ) 486 486 cross_normx_factor(k,j) = & … … 488 488 CASE ( 'wstsw2' ) 489 489 cross_normx_factor(k,j) = & 490 hom_sum(nzb+8, var_hom,normalizing_region) &491 * ( hom_sum(nzb+11, var_hom,normalizing_region) )**2490 hom_sum(nzb+8,pr_palm,normalizing_region) & 491 * ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2 492 492 CASE ( 'ws2' ) 493 493 cross_normx_factor(k,j) = & 494 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2494 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 495 495 CASE ( 'ws2tsw' ) 496 496 cross_normx_factor(k,j) = & 497 ( hom_sum(nzb+8, var_hom,normalizing_region) )**2&498 * hom_sum(nzb+11, var_hom,normalizing_region)497 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2& 498 * hom_sum(nzb+11,pr_palm,normalizing_region) 499 499 CASE ( 'ws3' ) 500 500 cross_normx_factor(k,j) = & 501 ( hom_sum(nzb+8, var_hom,normalizing_region) )**3501 ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3 502 502 503 503 END SELECT … … 506 506 CASE ( 'z_i' ) 507 507 cross_normy_factor(k,j) = & 508 hom_sum(nzb+6, var_hom,normalizing_region)508 hom_sum(nzb+6,pr_palm,normalizing_region) 509 509 510 510 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.