Ignore:
Timestamp:
Nov 29, 2018 7:16:36 PM (3 years ago)
Author:
suehring
Message:

Merge branch salsa with trunk

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/chem_emissions_mod.f90

    r3570 r3582  
    2222! Current revisions:
    2323! ------------------
    24 !
     24! - Removed salsa dependency.
     25! - Enabled PARAMETRIZED mode for default surfaces when LSM is not applied but
     26!   salsa is (M. Kurppa)
    2527!
    2628! Former revisions:
     
    890892   USE indices,                                                               &
    891893       ONLY: nnx,nny,nnz
    892    USE salsa_mod,                                                             &
    893        ONLY:  salsa
    894894   USE surface_mod,                                                           &
    895895       ONLY:  surf_lsm_h,surf_def_h,surf_usm_h
     
    14621462
    14631463                ENDDO
    1464              ELSEIF ( salsa ) THEN
    1465                 DO  m = 1, surf_def_h(0)%ns
    1466                    i = surf_def_h(0)%i(m)
    1467                    j = surf_def_h(0)%j(m)
    1468                    k = surf_def_h(0)%k(m)
    1469 
    1470 
    1471                    IF ( street_type_f%var(j,i) >= main_street_id  .AND.        &
    1472                         street_type_f%var(j,i) < max_street_id )               &
    1473                    THEN
    1474 
    1475                       !> Cycle over already matched species
    1476                       DO  ispec=1,nspec_out
    1477 
    1478                          !> PMs are already in mass units:micrograms: have to be converted to kilograms
    1479                          IF ( TRIM(spc_names(match_spec_model(ispec)))=="PM1"       &
    1480                               .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" &
    1481                               .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10")&
    1482                          THEN
    1483                                  
    1484                             surf_def_h(0)%cssws(match_spec_model(ispec),m) =   &
    1485                                   emiss_factor_main(match_spec_input(ispec)) * &
    1486                                   emis_distribution(1,j,i,ispec) * rho_air(k) /&
    1487                                   time_factor(1)
    1488                          ELSE
    1489 
    1490                          !> Other Species: inputs are micromoles: have to be converted             
    1491                             surf_def_h(0)%cssws(match_spec_model(ispec),m) =   &
    1492                                emiss_factor_main(match_spec_input(ispec)) *    &
    1493                                emis_distribution(1,j,i,ispec) *                &
    1494                                conv_to_ratio(k,j,i) *  rho_air(k) / time_factor(1)
    1495                          ENDIF
    1496                       ENDDO
    1497 
    1498                    ELSEIF ( street_type_f%var(j,i) >= side_street_id  .AND.    &
    1499                             street_type_f%var(j,i) < main_street_id )          &
    1500                    THEN
    1501 
    1502                    !> Cycle over already matched species
    1503                       DO  ispec=1,nspec_out
    1504 
    1505                          !> PMs are already in mass units: micrograms
    1506                          IF ( TRIM(spc_names(match_spec_model(ispec)))=="PM1"   &
    1507                               .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" &
    1508                               .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10")&
    1509                          THEN
    1510 
    1511                             surf_def_h(0)%cssws(match_spec_model(ispec),m) =   &
    1512                                emiss_factor_side(match_spec_input(ispec)) *    &
    1513                                emis_distribution(1,j,i,ispec) * rho_air(k) /   &
    1514                                time_factor(1) 
    1515                          ELSE
    1516                
    1517                             surf_def_h(0)%cssws(match_spec_model(ispec),m) =   &
    1518                                emiss_factor_side(match_spec_input(ispec)) *    &
    1519                                emis_distribution(1,j,i,ispec) *                &
    1520                                conv_to_ratio(k,j,i) * rho_air(k) / time_factor(1)
    1521                          ENDIF
    1522 
    1523                       ENDDO
    1524 
    1525                    ELSE
    1526 
    1527                    !> If no street type is defined, then assign null emissions to all the species
    1528                       surf_def_h(0)%cssws(:,m) = 0.0_wp
    1529 
    1530                    ENDIF
    1531 
    1532                 ENDDO
    1533 
    15341464             ENDIF
    15351465
Note: See TracChangeset for help on using the changeset viewer.