Ignore:
Timestamp:
Jun 4, 2007 8:07:41 AM (17 years ago)
Author:
raasch
Message:

more preliminary uncomplete changes for ocean version

File:
1 edited

Legend:

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

    r95 r96  
    14851485             hom(:,2,22,:) = SPREAD( zw, 2, statistic_regions+1 )
    14861486
     1487          CASE ( 'sa', '#sa' )
     1488             IF ( .NOT. ocean )  THEN
     1489                IF ( myid == 0 )  THEN
     1490                   PRINT*, '+++ check_parameters:  data_output_pr = ', &
     1491                           data_output_pr(i),                          &
     1492                           '    is not implemented for ocean = FALSE'
     1493                ENDIF
     1494                CALL local_stop
     1495             ELSE
     1496                dopr_index(i) = 23
     1497                dopr_unit(i)  = 'psu'
     1498                hom(:,2,23,:) = SPREAD( zu, 2, statistic_regions+1 )
     1499                IF ( data_output_pr(i)(1:1) == '#' )  THEN
     1500                   dopr_initial_index(i) = 26
     1501                   hom(:,2,26,:)         = SPREAD( zu, 2, statistic_regions+1 )
     1502                   hom(nzb,2,26,:)       = 0.0    ! weil zu(nzb) negativ ist
     1503                   data_output_pr(i)     = data_output_pr(i)(2:)
     1504                ENDIF
     1505             ENDIF
     1506
    14871507          CASE ( 'u*2' )
    14881508             dopr_index(i) = 30
     
    18461866             dopr_unit(i)  = 'K m/s'
    18471867             hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 )
     1868
     1869          CASE ( 'rho' )
     1870             dopr_index(i) = 64
     1871             dopr_unit(i)  = 'kg/m3'
     1872             hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
     1873
     1874          CASE ( 'w"sa"' )
     1875             IF ( .NOT. ocean ) THEN
     1876                IF ( myid == 0 )  THEN
     1877                   PRINT*, '+++ check_parameters:  data_output_pr = ', &
     1878                           data_output_pr(i),                          &
     1879                           '    is not implemented for ocean = FALSE'
     1880                ENDIF
     1881                CALL local_stop
     1882             ELSE
     1883                dopr_index(i) = 65
     1884                dopr_unit(i)  = 'psu m/s'
     1885                hom(:,2,65,:) = SPREAD( zw, 2, statistic_regions+1 )
     1886             ENDIF
     1887
     1888          CASE ( 'w*sa*' )
     1889             IF ( .NOT. ocean ) THEN
     1890                IF ( myid == 0 )  THEN
     1891                   PRINT*, '+++ check_parameters:  data_output_pr = ', &
     1892                           data_output_pr(i),                          &
     1893                           '    is not implemented for ocean = FALSE'
     1894                ENDIF
     1895                CALL local_stop
     1896             ELSE
     1897                dopr_index(i) = 66
     1898                dopr_unit(i)  = 'psu m/s'
     1899                hom(:,2,66,:) = SPREAD( zw, 2, statistic_regions+1 )
     1900             ENDIF
     1901
     1902          CASE ( 'wsa' )
     1903             IF ( .NOT. ocean ) THEN
     1904                IF ( myid == 0 )  THEN
     1905                   PRINT*, '+++ check_parameters:  data_output_pr = ', &
     1906                           data_output_pr(i),                          &
     1907                           '    is not implemented for ocean = FALSE'
     1908                ENDIF
     1909                CALL local_stop
     1910             ELSE
     1911                dopr_index(i) = 67
     1912                dopr_unit(i)  = 'psu m/s'
     1913                hom(:,2,67,:) = SPREAD( zw, 2, statistic_regions+1 )
     1914             ENDIF
    18481915
    18491916
     
    20682135             unit = 'kg/kg'
    20692136
     2137          CASE ( 'rho' )
     2138             IF ( .NOT. ocean )  THEN
     2139                IF ( myid == 0 )  THEN
     2140                   PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
     2141                                '" requires ocean = .TRUE.'
     2142                ENDIF
     2143                CALL local_stop
     2144             ENDIF
     2145             unit = 'kg/m3'
     2146
    20702147          CASE ( 's' )
    20712148             IF ( .NOT. passive_scalar )  THEN
     
    20772154             ENDIF
    20782155             unit = 'conc'
     2156
     2157          CASE ( 'sa' )
     2158             IF ( .NOT. ocean )  THEN
     2159                IF ( myid == 0 )  THEN
     2160                   PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
     2161                                '" requires ocean = .TRUE.'
     2162                ENDIF
     2163                CALL local_stop
     2164             ENDIF
     2165             unit = 'psu'
    20792166
    20802167          CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*' )
Note: See TracChangeset for help on using the changeset viewer.