Changeset 3864 for palm/trunk/SOURCE/advec_ws.f90

Timestamp:

Apr 5, 2019 9:01:56 AM (6 years ago)

Author:

monakurppa

Message:

major changes in salsa: data input, format and performance

• Time-dependent emissions enabled: lod=1 for yearly PM emissions that are normalised depending on the time, and lod=2 for preprocessed emissions (similar to the chemistry module).
• Additionally, 'uniform' emissions allowed. This emission is set constant on all horisontal upward facing surfaces and it is created based on parameters surface_aerosol_flux, aerosol_flux_dpg/sigmag/mass_fracs_a/mass_fracs_b.
• All emissions are now implemented as surface fluxes! No 3D sources anymore.
• Update the emission information by calling salsa_emission_update if skip_time_do_salsa >= time_since_reference_point and next_aero_emission_update <= time_since_reference_point
• Aerosol background concentrations read from PIDS_DYNAMIC. The vertical grid must match the one applied in the model.
• Gas emissions and background concentrations can be also read in in salsa_mod if the chemistry module is not applied.
• In deposition, information on the land use type can be now imported from the land use model
• Use SI units in PARIN, i.e. n_lognorm given in #/m3 and dpg in metres.
• Apply 100 character line limit
• Change all variable names from capital to lowercase letter
• Change real exponents to integer if possible. If not, precalculate the value of exponent
• Rename in1a to start_subrange_1a, fn2a to end_subrange_1a etc.
• Rename nbins --> nbins_aerosol, ncc_tot --> ncomponents_mass and ngast --> ngases_salsa
• Rename ibc to index_bc, idu to index_du etc.
• Renamed loop indices b, c and sg to ib, ic and ig
• run_salsa subroutine removed
• Corrected a bud in salsa_driver: falsely applied ino instead of inh
• Call salsa_tendency within salsa_prognostic_equations which is called in module_interface_mod instead of prognostic_equations_mod
• Removed tailing white spaces and unused variables
• Change error message to start by PA instead of SA

File:

• palm/trunk/SOURCE/advec_ws.f90 (modified) (11 diffs)

palm/trunk/SOURCE/advec_ws.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove tailing white spaces
+! 
+! 3696 2019-01-24 16:37:35Z suehring
 ! Bugfix in degradation height
 ! 
@@ -349 +352 @@
 
        USE salsa_util_mod,                                                     &
-           ONLY:  sums_salsa_ws_l   
-       
+           ONLY:  sums_salsa_ws_l
+
        USE statistics,                                                         &
            ONLY:  sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wsncs_ws_l,& 
@@ -357 +360 @@
                   sums_wssas_ws_l,  sums_wsss_ws_l, sums_wsus_ws_l,            &
                   sums_wsvs_ws_l
-  
+
 
 !
@@ -418 +421 @@
              sums_wssas_ws_l = 0.0_wp
           ENDIF        
-          
+
           IF ( salsa )  THEN
              ALLOCATE( sums_salsa_ws_l(nzb:nzt+1,0:threads_per_task-1) )
-             sums_salsa_ws_l = 0.0_wp   
+             sums_salsa_ws_l = 0.0_wp
           ENDIF
 
@@ -1143 +1146 @@
 
        USE salsa_util_mod,                                                     &
-           ONLY:  sums_salsa_ws_l            
-           
+           ONLY:  sums_salsa_ws_l
+
        USE statistics,                                                         &
            ONLY:  sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wsncs_ws_l,&
@@ -1221 +1224 @@
 
        USE salsa_util_mod,                                                     &
-           ONLY:  sums_salsa_ws_l   
+           ONLY:  sums_salsa_ws_l
 
        USE statistics,                                                         &
@@ -1824 +1827 @@
                     ) * weight_substep(intermediate_timestep_count)
              ENDDO
-             
+
           CASE ( 's' )
-          
+
              DO  k = nzb, nzt
                 sums_wsss_ws_l(k,tn)  = sums_wsss_ws_l(k,tn) +                 &
@@ -1835 +1838 @@
                     ) * weight_substep(intermediate_timestep_count)
              ENDDO
-                  
+
          CASE ( 'aerosol_mass', 'aerosol_number', 'salsa_gas' )
-         
+
              DO  k = nzb, nzt
                 sums_salsa_ws_l(k,tn)  = sums_salsa_ws_l(k,tn) +               &
@@ -1845 +1848 @@
                                 *   ABS( w(k,j,i) - hom(k,1,3,0)            )  &
                     ) * weight_substep(intermediate_timestep_count)
-             ENDDO             
+             ENDDO
 
 !          CASE ( 'kc' )
@@ -1852 +1855 @@
 
          END SELECT
-         
+
     END SUBROUTINE advec_s_ws_ij
 
@@ -3595 +3598 @@
            ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzb_max,    &
                   nzt, advc_flags_1
-           
+
        USE kinds
-       
+
        USE salsa_util_mod,                                                     &
-           ONLY:  sums_salsa_ws_l      
-       
+           ONLY:  sums_salsa_ws_l
+
        USE statistics,                                                         &
            ONLY:  hom, sums_wspts_ws_l, sums_wsqs_ws_l, sums_wssas_ws_l,       &
                   sums_wsqcs_ws_l, sums_wsqrs_ws_l, sums_wsncs_ws_l,           &
                   sums_wsnrs_ws_l, sums_wsss_ws_l, weight_substep  
-                  
+