Changeset 3094 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Nov 30, 2017, 2:57:51 PM (7 years ago)
Author:
oboucher
Message:

Adding calculation of H2SO4 to particles flux diag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90

    r2752 r3094  
    44  USE aerophys
    55  USE infotrac
    6   USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, R2SO4, DENSO4, f_r_wet
     6  USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, budg_h2so4_to_part, R2SO4, DENSO4, f_r_wet
    77  USE nucleation_tstep_mod
    88  USE cond_evap_tstep_mod
     
    6464
    6565  DO ilon=1, klon
     66!
     67!--initialisation of diagnostic
     68  budg_h2so4_to_part(ilon)=0.0
     69!
    6670  DO ilev=1, klev
     71!
     72!--initialisation of diagnostic
     73  budg_3D_nucl(ilon,ilev)=0.0
     74  budg_3D_cond_evap(ilon,ilev)=0.0
     75!
    6776  ! only in the stratosphere
    6877  IF (is_strato(ilon,ilev)) THEN
    6978    ! initialize sulfur fluxes
    70     budg_3D_nucl(ilon,ilev)=0.0
    71     budg_3D_cond_evap(ilon,ilev)=0.0
    7279    H2SO4_init=tr_seri(ilon,ilev,id_H2SO4_strat)
    7380    ! adaptive timestep for nucleation and condensation
     
    141148    budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol &
    142149             & *evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG
     150    ! compute vertically integrated flux due to the net effect of nucleation and condensation/evaporation
     151    budg_h2so4_to_part(ilon)=budg_h2so4_to_part(ilon)+(H2SO4_init-tr_seri(ilon,ilev,id_H2SO4_strat)) &
     152             & *mSatom/mH2SO4mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
    143153  ENDIF
    144154  ENDDO
Note: See TracChangeset for help on using the changeset viewer.