Ignore:
Timestamp:
May 22, 2007 3:46:47 PM (14 years ago)
Author:
raasch
Message:

Preliminary update for user defined profiles

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_profiles.f90

    r4 r87  
    44! Actual revisions:
    55! -----------------
    6 !
     6! var_hom renamed pr_palm
    77!
    88! Former revisions:
     
    111111
    112112             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 /), &
    114114                                     start = (/ 1 /), count = (/ 1 /) )
    115115             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 331 )
    116116
    117117             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 /), &
    119119                                     start = (/ 1 /), count = (/ 1 /) )
    120120             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 332 )
    121121
    122122             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 /), &
    124124                                     start = (/ 1 /), count = (/ 1 /) )
    125125             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 333 )
    126126
    127127             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)    /), &
    130130                                     start = (/ 1 /), count = (/ 1 /) )
    131131             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 334 )
    132132
    133133             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 /), &
    136136                                     start = (/ 1 /), count = (/ 1 /) )
    137137             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 335 )
    138138
    139139             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) /), &
    141141                                     start = (/ 1 /), count = (/ 1 /) )
    142142             IF ( nc_stat /= NF90_NOERR )  CALL handle_netcdf_error( 336 )
     
    237237                         CASE ( 'ts2' )
    238238                            cross_normx_factor(k,j) = &
    239                              ( hom_sum(nzb+3,var_hom,normalizing_region) )**2
     239                             ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
    240240                         CASE ( 'wpt0' )
    241241                            cross_normx_factor(k,j) = &
     
    243243                         CASE ( 'wsts2' )
    244244                            cross_normx_factor(k,j) = &
    245                              hom_sum(nzb+8,var_hom,normalizing_region)  &
    246                            * ( hom_sum(nzb+3,var_hom,normalizing_region) )**2
     245                             hom_sum(nzb+8,pr_palm,normalizing_region)  &
     246                           * ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
    247247                         CASE ( 'ws2' )
    248248                            cross_normx_factor(k,j) = &
    249                              ( hom_sum(nzb+8,var_hom,normalizing_region) )**2
     249                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
    250250                         CASE ( 'ws2ts' )
    251251                            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)
    254254                         CASE ( 'ws3' )
    255255                            cross_normx_factor(k,j) = &
    256                              ( hom_sum(nzb+8,var_hom,normalizing_region) )**3
     256                             ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
    257257
    258258                      END SELECT
     
    262262                         CASE ( 'z_i' )
    263263                            cross_normy_factor(k,j) = &
    264                                     hom_sum(nzb+6,var_hom,normalizing_region)
     264                                    hom_sum(nzb+6,pr_palm,normalizing_region)
    265265
    266266                      END SELECT
     
    351351
    352352          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 /), &
    354354                                  start = (/ dopr_time_count /),               &
    355355                                  count = (/ 1 /) )
     
    357357
    358358          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 /), &
    360360                                  start = (/ dopr_time_count /),               &
    361361                                  count = (/ 1 /) )
     
    363363
    364364          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 /), &
    366366                                  start = (/ dopr_time_count /),               &
    367367                                  count = (/ 1 /) )
     
    369369
    370370          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)    /), &
    373373                                  start = (/ dopr_time_count /),               &
    374374                                  count = (/ 1 /) )
     
    376376
    377377          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 /), &
    380380                                  start = (/ dopr_time_count /),               &
    381381                                  count = (/ 1 /) )
     
    383383
    384384          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) /), &
    386386                                  start = (/ dopr_time_count /),               &
    387387                                  count = (/ 1 /) )
     
    482482                   CASE ( 'tsw2' )
    483483                      cross_normx_factor(k,j) = &
    484                             ( hom_sum(nzb+11,var_hom,normalizing_region) )**2
     484                            ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
    485485                   CASE ( 'wpt0' )
    486486                      cross_normx_factor(k,j) = &
     
    488488                   CASE ( 'wstsw2' )
    489489                      cross_normx_factor(k,j) = &
    490                               hom_sum(nzb+8,var_hom,normalizing_region)  &
    491                           * ( hom_sum(nzb+11,var_hom,normalizing_region) )**2
     490                              hom_sum(nzb+8,pr_palm,normalizing_region)  &
     491                          * ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
    492492                   CASE ( 'ws2' )
    493493                      cross_normx_factor(k,j) = &
    494                             ( hom_sum(nzb+8,var_hom,normalizing_region) )**2
     494                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
    495495                   CASE ( 'ws2tsw' )
    496496                      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)
    499499                   CASE ( 'ws3' )
    500500                      cross_normx_factor(k,j) = &
    501                             ( hom_sum(nzb+8,var_hom,normalizing_region) )**3
     501                            ( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
    502502
    503503                END SELECT
     
    506506                   CASE ( 'z_i' )
    507507                      cross_normy_factor(k,j) = &
    508                                    hom_sum(nzb+6,var_hom,normalizing_region)
     508                                   hom_sum(nzb+6,pr_palm,normalizing_region)
    509509
    510510                END SELECT
Note: See TracChangeset for help on using the changeset viewer.