Ignore:
Timestamp:
Nov 27, 2018 5:44:21 PM (3 years ago)
Author:
kanani
Message:

Fix too long lines (chemistry_model_mod, chem_emissions_mod), correct terminal message (palmrun)

File:
1 edited

Legend:

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

    r3483 r3570  
    2727! -----------------
    2828! $Id$
     29! resler:
     30! Break lines at 132 characters
     31!
     32! 3483 2018-11-02 14:19:26Z raasch
    2933! bugfix: wrong locations of netCDF directives fixed
    3034!
     
    704708                DO  ispec = 1 , len_index
    705709
    706                    IF ( emiss_factor_main(match_spec_input(ispec)) .LT. 0 .AND. emiss_factor_side(match_spec_input(ispec)) .LT. 0 ) THEN
     710                   IF ( emiss_factor_main(match_spec_input(ispec)) .LT. 0 .AND. &
     711                        emiss_factor_side(match_spec_input(ispec)) .LT. 0 ) THEN
    707712
    708713                      message_string = 'PARAMETERIZED emissions mode selected:'            //          &
     
    11851190       
    11861191          !Assign Constant Values of time factors, diurnal time profile for traffic sector:
    1187           par_emis_time_factor( : ) = (/ 0.009, 0.004, 0.004, 0.009, 0.029, 0.039, 0.056, 0.053, 0.051, 0.051, 0.052, 0.055, &
    1188                                                  0.059, 0.061, 0.064, 0.067, 0.069, 0.069, 0.049, 0.039, 0.039, 0.029, 0.024, 0.019 /)
     1192          par_emis_time_factor( : ) = &
     1193            (/ 0.009, 0.004, 0.004, 0.009, 0.029, 0.039, 0.056, 0.053, 0.051, 0.051, 0.052, 0.055, &
     1194               0.059, 0.061, 0.064, 0.067, 0.069, 0.069, 0.049, 0.039, 0.039, 0.029, 0.024, 0.019 /)
    11891195         
    11901196          !> in this case allocate time factor each hour in a day
     
    12491255             !> NOX Compositions
    12501256                IF (TRIM(spc_names(match_spec_model(ispec)))=="NO") THEN
    1251                 !>             Kg/m2*s                   kg/m2*s                                                   
    1252                    delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%nox_comp(icat,1)*con_factor*hour_per_day
    1253                    
     1257                !>             Kg/m2*s                   kg/m2*s
     1258                   delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1259                                                 emt_att%nox_comp(icat,1)*con_factor*hour_per_day
     1260
    12541261                   emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
    12551262
    12561263                ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="NO2") THEN
    1257    
    1258                    delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%nox_comp(icat,2)*con_factor*hour_per_day
     1264
     1265                   delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1266                                                 emt_att%nox_comp(icat,2)*con_factor*hour_per_day
    12591267
    12601268                   emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
    12611269 
    12621270             !> SOX Compositions
    1263                      
     1271
    12641272                ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="SO2") THEN
    1265                    !>             Kg/m2*s                   kg/m2*s                                                                       
    1266                    delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%sox_comp(icat,1)*con_factor*hour_per_day
     1273                   !>             Kg/m2*s                   kg/m2*s
     1274                   delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1275                                                 emt_att%sox_comp(icat,1)*con_factor*hour_per_day
    12671276
    12681277                   emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
    12691278
    12701279                ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="SO4") THEN
    1271                    !>             Kg/m2*s                   kg/m2*s                                                                         
    1272                    delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%sox_comp(icat,2)*con_factor*hour_per_day
     1280                   !>             Kg/m2*s                   kg/m2*s
     1281                   delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1282                                                 emt_att%sox_comp(icat,2)*con_factor*hour_per_day
    12731283
    12741284                   emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
     
    12821292                   DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,1))
    12831293
    1284                       delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%pm_comp(icat,i_pm_comp,1)*con_factor*hour_per_day
    1285                                                                                          
    1286 
    1287                       emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
     1294                      delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1295                                                    emt_att%pm_comp(icat,i_pm_comp,1)*con_factor*hour_per_day
     1296
     1297                      emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ &
     1298                                                                 delta_emis(nys:nyn,nxl:nxr)
    12881299 
    12891300                   ENDDO
     
    12951306                   DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,2))
    12961307
    1297                       delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%pm_comp(icat,i_pm_comp,2)*con_factor*hour_per_day
    1298                                                                                          
    1299 
    1300                       emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
     1308                      delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1309                                                    emt_att%pm_comp(icat,i_pm_comp,2)*con_factor*hour_per_day
     1310
     1311                      emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ &
     1312                                                                 delta_emis(nys:nyn,nxl:nxr)
    13011313 
    13021314                   ENDDO
     
    13071319                !> Cycle over the different pm components for PM10 type
    13081320                   DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,3)) 
    1309                        
    1310                       delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)*emt_att%pm_comp(icat,i_pm_comp,3)*con_factor*hour_per_day
    1311                                                                                                  
    1312 
    1313                       emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
     1321
     1322                      delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* &
     1323                                                    emt_att%pm_comp(icat,i_pm_comp,3)*con_factor*hour_per_day
     1324
     1325                      emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ &
     1326                                                                 delta_emis(nys:nyn,nxl:nxr)
    13141327
    13151328                   ENDDO
     
    13261339                                                       emt_att%voc_comp(icat,match_spec_voc_input(ivoc))*con_factor*hour_per_day
    13271340
    1328                          emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)
     1341                         emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ &
     1342                                                                    delta_emis(nys:nyn,nxl:nxr)
    13291343
    13301344                      ENDIF                       
     
    13791393
    13801394                         !> PMs are already in mass units: kilograms
    1381                          IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25"  &
     1395                         IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" &
     1396                             .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25"  &
    13821397                             .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN
    13831398
    1384                             !              kg/(m^2*s) *kg/m^3                               
     1399                            !              kg/(m^2*s) *kg/m^3
    13851400                            surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_main(match_spec_input(ispec)) * &
    13861401                            !                                                       kg/(m^2*s)
     
    14131428
    14141429                         !> PMs are already in mass units: micrograms
    1415                          IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25"  &
     1430                         IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" &
     1431                             .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25"  &
    14161432                             .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN
    14171433
Note: See TracChangeset for help on using the changeset viewer.