Ignore:
Timestamp:
Jan 22, 2018 5:17:47 PM (3 years ago)
Author:
kanani
Message:

Removal of chem directive, plus minor changes

File:
1 edited

Legend:

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

    r2759 r2766  
    2626! -----------------
    2727! $Id$
     28! Removed preprocessor directive __chem
     29!
     30! 2759 2018-01-17 16:24:59Z suehring
    2831! Bugfix, consider density in vertical fluxes of passive scalar as well as
    2932! chemical species.
     
    123126               rho_air, rho_air_zw, zu, zw, waterflux_input_conversion
    124127
    125 #if defined( __chem )
    126128    USE chem_modules
    127 #endif
    128129
    129130    USE control_parameters               
     
    968969!
    969970!--    Scaling parameter (cs*) and surface flux of chemical species
    970 #if defined ( __chem )
    971971       IF ( air_chemistry )  THEN
    972972          ALLOCATE ( surfaces%css(1:nvar,1:surfaces%ns)   )   
    973973          ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    974974       ENDIF
    975 #endif
    976975!
    977976!--    Arrays for storing potential temperature and
     
    10521051!
    10531052!--    Chemical species flux
    1054 #if defined( __chem )
    10551053       IF ( air_chemistry )  THEN
    10561054          ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    10571055       ENDIF
    1058 #endif
    10591056!
    10601057!--       
     
    11671164!
    11681165!--    Scaling parameter (cs*) and surface flux of chemical species
    1169 #if defined( __chem )
    11701166       IF ( air_chemistry )  THEN
    11711167             ALLOCATE ( surfaces%css(1:nvar,1:surfaces%ns)   )   
    11721168             ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    11731169       ENDIF
    1174 #endif
    11751170!
    11761171!--    Arrays for storing potential temperature and
     
    17371732             IF ( passive_scalar )  surf%ss(num_h) = 0.0_wp
    17381733
    1739 #if defined( __chem )
    17401734             DO  lsp = 1, nvar
    17411735                IF ( air_chemistry )  surf%css(lsp,num_h) = 0.0_wp
    17421736             ENDDO
    1743 #endif
    17441737!
    17451738!--          Set initial value for surface temperature
     
    18191812                ENDIF
    18201813
    1821 #if defined( __chem )
    18221814                IF ( air_chemistry )  THEN
    18231815                   lsp_pr = 1
     
    18491841                   ENDDO
    18501842                ENDIF
    1851 #endif
    18521843
    18531844                IF ( ocean )  THEN
     
    19201911!
    19211912!--          Prescribe top chemical species' flux
    1922 #if defined( __chem )
    19231913             DO  lsp = 1, nvar
    19241914                IF ( air_chemistry  .AND.  constant_top_csflux(lsp) )  THEN
     
    19261916                ENDIF
    19271917             ENDDO
    1928 #endif
    19291918!
    19301919!--          Prescribe top salinity flux
     
    20632052             ENDIF
    20642053
    2065 #if defined( __chem )
    20662054             IF ( air_chemistry )  THEN       
    20672055                DO  lsp = 1, nvar
     
    20702058                ENDDO
    20712059             ENDIF
    2072 #endif
    20732060
    20742061!
     
    22712258                   IF ( ALLOCATED( surf_def_h(l)%ssws ) )                      &
    22722259                      surf_h(l)%ssws(mm(l))    = surf_def_h(l)%ssws(m)
    2273 #if defined( __chem )
    22742260                   IF ( ALLOCATED( surf_def_h(l)%css ) )  THEN
    22752261                      DO  lsp = 1,nvar
     
    22822268                      ENDDO
    22832269                   ENDIF
    2284 #endif
    22852270                   IF ( ALLOCATED( surf_def_h(l)%ncsws ) )                     &
    22862271                      surf_h(l)%ncsws(mm(l))   = surf_def_h(l)%ncsws(m)
     
    23282313                      IF ( ALLOCATED( surf_lsm_h%ssws ) )                      &
    23292314                         surf_h(0)%ssws(mm(0))    = surf_lsm_h%ssws(m)
    2330 #if defined( __chem )
    23312315                      IF ( ALLOCATED( surf_lsm_h%css ) )  THEN                 
    23322316                         DO  lsp = 1, nvar
     
    23392323                         ENDDO
    23402324                      ENDIF
    2341 #endif
    23422325                      IF ( ALLOCATED( surf_lsm_h%ncsws ) )                     &
    23432326                         surf_h(0)%ncsws(mm(0))   = surf_lsm_h%ncsws(m)
     
    23852368                      IF ( ALLOCATED( surf_usm_h%ssws ) )                      &
    23862369                         surf_h(0)%ssws(mm(0))    = surf_usm_h%ssws(m)
    2387 #if defined( __chem )
    23882370                      IF ( ALLOCATED( surf_usm_h%css ) )  THEN             
    23892371                         DO lsp = 1, nvar
     
    23962378                         ENDDO
    23972379                      ENDIF
    2398 #endif
    23992380                      IF ( ALLOCATED( surf_usm_h%ncsws ) )                     &
    24002381                         surf_h(0)%ncsws(mm(0))   = surf_usm_h%ncsws(m)
     
    24742455                   IF ( ALLOCATED( surf_def_v(l)%ssws ) )                      &
    24752456                      surf_v(l)%ssws(mm(l))    = surf_def_v(l)%ssws(m)
    2476 #if defined( __chem )
    24772457                   IF ( ALLOCATED( surf_def_v(l)%css ) )  THEN               
    24782458                      DO  lsp = 1, nvar
     
    24852465                      ENDDO
    24862466                   ENDIF
    2487 #endif
    24882467                   IF ( ALLOCATED( surf_def_v(l)%ncsws ) )                     &
    24892468                      surf_v(l)%ncsws(mm(l))   = surf_def_v(l)%ncsws(m)
     
    25362515                   IF ( ALLOCATED( surf_lsm_v(l)%ssws ) )                      &
    25372516                      surf_v(l)%ssws(mm(l))    = surf_lsm_v(l)%ssws(m)
    2538 #if defined( __chem )
    25392517                   IF ( ALLOCATED( surf_lsm_v(l)%css ) )  THEN             
    25402518                      DO  lsp = 1, nvar
     
    25472525                      ENDDO
    25482526                   ENDIF
    2549 #endif
    25502527                   IF ( ALLOCATED( surf_lsm_v(l)%ncsws ) )                     &
    25512528                      surf_v(l)%ncsws(mm(l))   = surf_lsm_v(l)%ncsws(m)
     
    25982575                   IF ( ALLOCATED( surf_usm_v(l)%ssws ) )                      &
    25992576                      surf_v(l)%ssws(mm(l))    = surf_usm_v(l)%ssws(m)
    2600 #if defined( __chem )
    26012577                   IF ( ALLOCATED( surf_usm_v(l)%css ) )  THEN             
    26022578                      DO  lsp = 1, nvar
     
    26092585                      ENDDO
    26102586                   ENDIF
    2611 #endif
    26122587                   IF ( ALLOCATED( surf_usm_v(l)%ncsws ) )                     &
    26132588                      surf_v(l)%ncsws(mm(l))   = surf_usm_v(l)%ncsws(m)
     
    27172692             WRITE ( 14 )  surf_h(l)%ssws
    27182693          ENDIF
    2719 #if defined ( __chem )
    27202694          WRITE ( 14 )  'surf_h(' // dum // ')%css                 '
    27212695          IF ( ALLOCATED ( surf_h(l)%css ) )  THEN
     
    27262700             WRITE ( 14 )  surf_h(l)%cssws
    27272701          ENDIF
    2728 #endif
    27292702          WRITE ( 14 )  'surf_h(' // dum // ')%qcsws               ' 
    27302703          IF ( ALLOCATED ( surf_h(l)%qcsws ) )  THEN
     
    28142787             WRITE ( 14 )  surf_v(l)%ssws
    28152788          ENDIF
    2816 #if defined( __chem )
    28172789          WRITE ( 14 )  'surf_v(' // dum // ')%css                 ' 
    28182790          IF ( ALLOCATED ( surf_v(l)%css ) )  THEN
     
    28232795             WRITE ( 14 )  surf_v(l)%cssws
    28242796          ENDIF
    2825 #endif
    28262797          WRITE ( 14 )  'surf_v(' // dum // ')%qcsws               ' 
    28272798          IF ( ALLOCATED ( surf_v(l)%qcsws ) )  THEN
     
    30593030                      IF ( ALLOCATED( surf_h(0)%ssws )  .AND.  kk == 1 )       &
    30603031                         READ ( 13 )  surf_h(0)%ssws
    3061 #if defined( __chem )
    30623032                   CASE ( 'surf_h(0)%css' )
    30633033                      IF ( ALLOCATED( surf_h(0)%css )  .AND.  kk == 1 )        &
     
    30663036                      IF ( ALLOCATED( surf_h(0)%cssws )  .AND.  kk == 1 )      &
    30673037                         READ ( 13 )  surf_h(0)%cssws
    3068 #endif
    30693038                   CASE ( 'surf_h(0)%qcsws' )         
    30703039                      IF ( ALLOCATED( surf_h(0)%qcsws )  .AND.  kk == 1 )      &
     
    31383107                      IF ( ALLOCATED( surf_h(1)%ssws )  .AND.  kk == 1 )       &
    31393108                         READ ( 13 )  surf_h(1)%ssws
    3140 #if defined( __chem )
    31413109                   CASE ( 'surf_h(1)%css' )
    31423110                      IF ( ALLOCATED( surf_h(1)%css )  .AND.  kk == 1 )        &
     
    31453113                      IF ( ALLOCATED( surf_h(1)%cssws )  .AND.  kk == 1 )      &
    31463114                         READ ( 13 )  surf_h(1)%cssws
    3147 #endif
    31483115                   CASE ( 'surf_h(1)%qcsws' )         
    31493116                      IF ( ALLOCATED( surf_h(1)%qcsws )  .AND.  kk == 1 )      &
     
    32143181                      IF ( ALLOCATED( surf_h(2)%ssws )  .AND.  kk == 1 )       &
    32153182                         READ ( 13 )  surf_h(2)%ssws
    3216 #if defined( __chem )
    32173183                   CASE ( 'surf_h(2)%css' )
    32183184                      IF ( ALLOCATED( surf_h(2)%css )  .AND.  kk == 1 )        &
     
    32213187                      IF ( ALLOCATED( surf_h(2)%cssws )  .AND.  kk == 1 )      &
    32223188                         READ ( 13 )  surf_h(2)%cssws
    3223 #endif
    32243189                   CASE ( 'surf_h(2)%qcsws' )         
    32253190                      IF ( ALLOCATED( surf_h(2)%qcsws )  .AND.  kk == 1 )      &
     
    32893254                      IF ( ALLOCATED( surf_v(0)%ssws )  .AND.  kk == 1 )       &
    32903255                         READ ( 13 )  surf_v(0)%ssws
    3291 #if defined( __chem )
    32923256                   CASE ( 'surf_v(0)%css' )
    32933257                      IF ( ALLOCATED( surf_v(0)%css )  .AND.  kk == 1 )        &
     
    32963260                      IF ( ALLOCATED( surf_v(0)%cssws )  .AND.  kk == 1 )      &
    32973261                         READ ( 13 )  surf_v(0)%cssws
    3298 #endif
    32993262                   CASE ( 'surf_v(0)%qcsws' )         
    33003263                      IF ( ALLOCATED( surf_v(0)%qcsws )  .AND.  kk == 1 )      &
     
    33713334                      IF ( ALLOCATED( surf_v(1)%ssws )  .AND.  kk == 1 )       &
    33723335                         READ ( 13 )  surf_v(1)%ssws
    3373 #if defined( __chem )
    33743336                   CASE ( 'surf_v(1)%css' )
    33753337                      IF ( ALLOCATED( surf_v(1)%css )  .AND.  kk == 1 )        &
     
    33783340                      IF ( ALLOCATED( surf_v(1)%cssws )  .AND.  kk == 1 )      &
    33793341                         READ ( 13 )  surf_v(1)%cssws
    3380 #endif
    33813342                   CASE ( 'surf_v(1)%qcsws' )         
    33823343                      IF ( ALLOCATED( surf_v(1)%qcsws )  .AND.  kk == 1 )      &
     
    34533414                      IF ( ALLOCATED( surf_v(2)%ssws )  .AND.  kk == 1 )       &
    34543415                         READ ( 13 )  surf_v(2)%ssws
    3455 #if defined( __chem )
    34563416                   CASE ( 'surf_v(2)%css' )
    34573417                      IF ( ALLOCATED( surf_v(2)%css )  .AND.  kk == 1 )        &
     
    34603420                      IF ( ALLOCATED( surf_v(2)%cssws )  .AND.  kk == 1 )      &
    34613421                         READ ( 13 )  surf_v(2)%cssws
    3462 #endif
    34633422                   CASE ( 'surf_v(2)%qcsws' )         
    34643423                      IF ( ALLOCATED( surf_v(2)%qcsws )  .AND.  kk == 1 )      &
     
    35353494                      IF ( ALLOCATED( surf_v(3)%ssws )  .AND.  kk == 1 )       &
    35363495                         READ ( 13 )  surf_v(3)%ssws
    3537 #if defined( __chem )
    35383496                   CASE ( 'surf_v(3)%css' )
    35393497                      IF ( ALLOCATED( surf_v(3)%css )  .AND.  kk == 1 )        &
     
    35423500                      IF ( ALLOCATED( surf_v(3)%cssws )  .AND.  kk == 1 )      &
    35433501                         READ ( 13 )  surf_v(3)%cssws
    3544 #endif
    35453502                   CASE ( 'surf_v(3)%qcsws' )         
    35463503                      IF ( ALLOCATED( surf_v(3)%qcsws )  .AND.  kk == 1 )      &
     
    37623719             ENDIF
    37633720
    3764 #if defined( __chem )
    37653721             IF ( INDEX( TRIM( field_chr ), '%css' ) /= 0 )  THEN
    37663722                IF ( ALLOCATED( surf_target%css )  .AND.                     &
     
    37793735                ENDIF
    37803736             ENDIF
    3781 #endif
    37823737
    37833738             IF ( INDEX( TRIM( field_chr ), '%qcs' ) /= 0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.