Ignore:
Timestamp:
Oct 30, 2018 7:05:21 PM (3 years ago)
Author:
suehring
Message:

Branch salsa @3446 re-integrated into trunk

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/advec_ws.f90

    r3302 r3467  
    2525! -----------------
    2626! $Id$
     27! Implementation of a new aerosol module salsa.
     28!
     29! 3302 2018-10-03 02:39:40Z raasch
    2730! Stokes drift velocity added in scalar advection
    2831!
     
    4750! 2731 2018-01-09 17:44:02Z suehring
    4851! Enable loop vectorization by splitting the k-loop
    49 ! 
     52!
    5053! 2718 2018-01-02 08:49:38Z maronga
    5154! Corrected "Former revisions" section
     
    315318       USE pegrid
    316319
     320       USE salsa_util_mod,                                                     &
     321           ONLY:  sums_salsa_ws_l   
     322       
    317323       USE statistics,                                                         &
    318324           ONLY:  sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wsncs_ws_l,&
     
    381387             ALLOCATE( sums_wssas_ws_l(nzb:nzt+1,0:threads_per_task-1) )
    382388             sums_wssas_ws_l = 0.0_wp
    383           ENDIF
     389          ENDIF       
     390!         
     391!--       Mona: Now always allocated
     392          ALLOCATE( sums_salsa_ws_l(nzb:nzt+1,0:threads_per_task-1) )
     393          sums_salsa_ws_l = 0.0_wp         
    384394
    385395       ENDIF
     
    10961106                  ws_scheme_sca
    10971107
    1098        USE kinds
     1108       USE kinds 
    10991109
    11001110       USE bulk_cloud_model_mod,                                               &
    11011111           ONLY: bulk_cloud_model, microphysics_morrison, microphysics_seifert
    11021112
     1113       USE salsa_util_mod,                                                     &
     1114           ONLY:  sums_salsa_ws_l           
     1115           
    11031116       USE statistics,                                                         &
    11041117           ONLY:  sums_us2_ws_l, sums_vs2_ws_l, sums_ws2_ws_l, sums_wsncs_ws_l,&
    11051118                  sums_wsnrs_ws_l, sums_wspts_ws_l, sums_wsqcs_ws_l,           &
    11061119                  sums_wsqrs_ws_l, sums_wsqs_ws_l, sums_wsss_ws_l,             &
    1107                   sums_wssas_ws_l, sums_wsus_ws_l, sums_wsvs_ws_l     
     1120                  sums_wssas_ws_l, sums_wsus_ws_l, sums_wsvs_ws_l
    11081121                   
    11091122
     
    11341147          ENDIF
    11351148          IF ( ocean_mode )  sums_wssas_ws_l = 0.0_wp
    1136 
     1149       
     1150          sums_salsa_ws_l = 0.0_wp
     1151         
     1152         
    11371153       ENDIF
    11381154
     
    11671183       USE pegrid
    11681184
     1185       USE salsa_util_mod,                                                     &
     1186           ONLY:  sums_salsa_ws_l   
     1187
    11691188       USE statistics,                                                         &
    11701189           ONLY:  hom, sums_wsncs_ws_l, sums_wsnrs_ws_l, sums_wspts_ws_l,      &
    11711190                  sums_wsqcs_ws_l,  sums_wsqrs_ws_l, sums_wsqs_ws_l,           &
    1172                   sums_wssas_ws_l, sums_wsss_ws_l, weight_substep   
     1191                  sums_wssas_ws_l, sums_wsss_ws_l, weight_substep
    11731192                 
    11741193
     
    17041723                    ) * weight_substep(intermediate_timestep_count)
    17051724             ENDDO
     1725                 
     1726         CASE ( 'aerosol_mass', 'aerosol_number', 'salsa_gas' )
     1727         
     1728             DO  k = nzb, nzt
     1729                sums_salsa_ws_l(k,tn)  = sums_salsa_ws_l(k,tn) +               &
     1730                    ( flux_t(k) / ( w(k,j,i) + SIGN( 1.0E-20_wp, w(k,j,i) ) )  &
     1731                                * ( w(k,j,i) - hom(k,1,3,0)                 )  &
     1732                    + diss_t(k) / ( ABS(w(k,j,i)) + 1.0E-20_wp              )  &
     1733                                *   ABS( w(k,j,i) - hom(k,1,3,0)            )  &
     1734                    ) * weight_substep(intermediate_timestep_count)
     1735             ENDDO             
    17061736
    17071737!          CASE ( 'kc' )
     
    32113241           
    32123242       USE kinds
     3243       
     3244       USE salsa_util_mod,                                                     &
     3245           ONLY:  sums_salsa_ws_l     
    32133246       
    32143247       USE statistics,                                                         &
     
    37313764                                *   ABS(w(k,j,i) - hom(k,1,3,0)             )  &
    37323765                            ) * weight_substep(intermediate_timestep_count)
    3733                     ENDDO   
     3766                    ENDDO 
     3767     
     3768                 CASE ( 'aerosol_mass', 'aerosol_number', 'salsa_gas' )
     3769                     DO  k = nzb, nzt
     3770                        sums_salsa_ws_l(k,tn)  = sums_salsa_ws_l(k,tn)         &
     3771                          + ( flux_t(k)                                        &
     3772                                / ( w(k,j,i) + SIGN( 1.0E-20_wp, w(k,j,i) ) )  &
     3773                                * ( w(k,j,i) - hom(k,1,3,0)                 )  &
     3774                            + diss_t(k)                                        &
     3775                                / ( ABS(w(k,j,i)) + 1.0E-20_wp              )  &
     3776                                *   ABS(w(k,j,i) - hom(k,1,3,0)             )  &
     3777                            ) * weight_substep(intermediate_timestep_count)
     3778                     ENDDO                   
    37343779                                   
    37353780
Note: See TracChangeset for help on using the changeset viewer.