Changeset 3570 for palm/trunk


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

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

Location:
palm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/palmrun

    r3555 r3570  
    2727# -----------------
    2828# $Id$
     29# Rename job-id by run_id
     30#
     31# 3555 2018-11-22 13:50:37Z raasch
    2932# sed-command to remove repeated space-characters replaced by tr-command, to
    3033# guarantee mac compatible
     
    30153018
    30163019             # SUBMIT THE JOB
    3017           printf " *** submit the job (output of submit command, e.g. the job-id, may follow) \n"
     3020          printf " *** submit the job (output of submit command, e.g. the run_id, may follow) \n"
    30183021          if [[ $do_trace = true ]]
    30193022          then
  • 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
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r3543 r3570  
    2121!
    2222! Current revisions:
    23 ! -----------------
     23! ------------------
    2424!
    2525!
     
    2727! -----------------
    2828! $Id$
     29! resler:
     30! Break lines at 132 characters
     31!
     32! 3543 2018-11-20 17:06:15Z suehring
    2933! Remove tabs
    3034!
     
    323327    !
    324328    ! Set alpha for f_light (4.57 is conversion factor from 1./(mumol m-2 s-1) naar W m-2
    325     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: alpha   =(/0.009,0.009, 0.009,0.006,0.006, -999., -999.,0.009,-999.,-999.,0.009,0.006,-999.,0.009,0.008/)*4.57
     329    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: alpha   = &
     330      (/0.009,0.009, 0.009,0.006,0.006, -999., -999.,0.009,-999.,-999.,0.009,0.006,-999.,0.009,0.008/)*4.57
    326331    !
    327332    ! Set temperatures per land use for F_temp
    328     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmin =   (/12.0, 12.0,  12.0,  0.0,  0.0, -999., -999., 12.0, -999., -999., 12.0,  0.0, -999., 12.0,  8.0/)
    329     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Topt =   (/26.0, 26.0,  26.0, 18.0, 20.0, -999., -999., 26.0, -999., -999., 26.0, 20.0, -999., 26.0, 24.0 /)
    330     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmax =   (/40.0, 40.0,  40.0, 36.0, 35.0, -999., -999., 40.0, -999., -999., 40.0, 35.0, -999., 40.0, 39.0 /)
     333    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmin = &
     334      (/12.0, 12.0,  12.0,  0.0,  0.0, -999., -999., 12.0, -999., -999., 12.0,  0.0, -999., 12.0,  8.0/)
     335    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Topt = &
     336      (/26.0, 26.0,  26.0, 18.0, 20.0, -999., -999., 26.0, -999., -999., 26.0, 20.0, -999., 26.0, 24.0 /)
     337    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: Tmax = &
     338      (/40.0, 40.0,  40.0, 36.0, 35.0, -999., -999., 40.0, -999., -999., 40.0, 35.0, -999., 40.0, 39.0 /)
    331339    !
    332340    ! Set F_min:
    333     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: F_min   =(/0.01, 0.01,  0.01, 0.1,  0.1,   -999., -999.,0.01, -999.,-999.,0.01,0.1,-999.,0.01, 0.04/)
     341    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: F_min = &
     342      (/0.01, 0.01,  0.01, 0.1,  0.1,   -999., -999.,0.01, -999.,-999.,0.01,0.1,-999.,0.01, 0.04/)
    334343   
    335344    ! Set maximal conductance (m/s)
    336345    ! (R T/P) = 1/41000 mmol/m3 is given for 20 deg C to go from  mmol O3/m2/s to m/s
    337346    ! Could be refined to a function of T and P. in Jones
    338     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: g_max   =(/270., 300.,  300., 140., 150.,  -999., -999.,270., -999.,-999.,270., 150.,-999.,300., 422./)/41000
     347    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: g_max = &
     348      (/270., 300.,  300., 140., 150.,  -999., -999.,270., -999.,-999.,270., 150.,-999.,300., 422./)/41000
    339349    !
    340350    ! Set max, min for vapour pressure deficit vpd;
    341     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_max =(/1.3,  0.9,   0.9,  0.5,  1.0,   -999., -999.,1.3,  -999.,-999.,1.3,1.0,  -999.,0.9, 2.8/)
    342     REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_min =(/3.0,  2.8,   2.8,  3.0,  3.25,  -999., -999.,3.0,  -999.,-999.,3.0,3.25, -999.,2.8, 4.5/)
     351    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_max = &
     352      (/1.3,  0.9,   0.9,  0.5,  1.0,   -999., -999.,1.3,  -999.,-999.,1.3,1.0,  -999.,0.9, 2.8/)
     353    REAL(wp), DIMENSION(nlu_dep), PARAMETER :: vpd_min = &
     354      (/3.0,  2.8,   2.8,  3.0,  3.25,  -999., -999.,3.0,  -999.,-999.,3.0,3.25, -999.,2.8, 4.5/)
    343355    !
    344356    ! 
     
    49584970
    49594971
    4960       REAL(kind=wp), PARAMETER     :: rsoil_wet(ncmp)    = (/2000.,      10.,    2000.,    -999.,      10.,    -999., -999., -999., -999., -999./)
    4961       REAL(kind=wp), PARAMETER     :: rsoil_frozen(ncmp) = (/2000.,     500.,    2000.,    -999.,    1000.,    -999., -999., -999., -999., -999./)
     4972      REAL(kind=wp), PARAMETER     :: rsoil_wet(ncmp)    = &
     4973          (/2000.,      10.,    2000.,    -999.,      10.,    -999., -999., -999., -999., -999./)
     4974      REAL(kind=wp), PARAMETER     :: rsoil_frozen(ncmp) = &
     4975          (/2000.,     500.,    2000.,    -999.,    1000.,    -999., -999., -999., -999., -999./)
    49624976
    49634977
     
    50275041      ! b = empirical constant for computation of rinc (in canopy resistance) (= 14 m-1 or -999 IF not applicable)
    50285042      ! h = vegetation height (m)                             grass arabl  crops conIF decid   water  urba   othr  desr   ice  sav  trf   wai   med semi
    5029       REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: b = (/ -999,  14,    14,   14,   14,    -999,  -999, -999, -999, -999, -999, 14, -999, 14, 14 /)
    5030       REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: h = (/ -999,  1,     1,    20,   20,    -999,  -999, -999, -999, -999, -999, 20,  -999, 1 ,  1 /)
     5043      REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: b = &
     5044        (/ -999,  14,    14,   14,   14,    -999,  -999, -999, -999, -999, -999, 14, -999, 14, 14 /)
     5045      REAL(kind=wp), DIMENSION(nlu_dep), PARAMETER :: h = &
     5046        (/ -999,  1,     1,    20,   20,    -999,  -999, -999, -999, -999, -999, 20,  -999, 1 ,  1 /)
    50315047
    50325048      ! Compute Rinc only for arable land, perm. crops, forest; otherwise Rinc = 0:
Note: See TracChangeset for help on using the changeset viewer.