Ignore:
Timestamp:
Jan 31, 2022, 9:20:17 PM (3 years ago)
Author:
dcugnet
Message:
  • Fix for unadvected tracers (iadv==0)
  • The key %isH2Ofamily, from the derived type "trac_type", is replaced with the more general

key %isInPhysics, which is TRUE for tracers both in "qx" and "tr_seri".

Currently, FALSE for tracers descending on H2O (isotopes and tagging tracers included). Could be set to FALSE
for interactive CO2 (type_trac=='inco') or ice supersaturated cloud content (tranfered to "rneb_seri")

Location:
LMDZ6/trunk/libf/phylmd/Dust
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/Dust/phys_output_write_spl_mod.F90

    r4056 r4071  
    16121612           itr = 0
    16131613           DO iq = 1, nqtot
    1614              IF(tracers(iq)%isH2Ofamily) CYCLE
     1614             IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    16151615             itr = itr+1
    16161616             CALL histwrite_phy(o_trac(itr), tr_seri(:,:,itr))

  • LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

    r4056 r4071  
    12621262      itr = 0
    12631263      DO iq = 1, nqtot
    1264          IF(tracers(iq)%isH2Ofamily) CYCLE
     1264         IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    12651265         itr = itr+1
    12661266         write(str2,'(i2.2)') itr
    … …  
    14191419        itr = 0
    14201420        do iq=1,nqtot
    1421           IF(tracers(iq)%isH2Ofamily) CYCLE
     1421          IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    14221422          itr = itr+1
    14231423          print *, itr, TRIM(tracers(iq)%name)

  • LMDZ6/trunk/libf/phylmd/Dust/splaeropt_5wv_rrtm.F90

    r4056 r4071  
    114114  itr = 0
    115115  DO iq = 1, nqtot
    116     IF(tracers(iq)%isH2Ofamily) CYCLE
     116    IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    117117    itr = itr+1
    118118    SELECT CASE(tracers(iq)%name)

  • LMDZ6/trunk/libf/phylmd/Dust/splaeropt_6bands_rrtm.F90

    r4056 r4071  
    167167  itr = 0
    168168  DO iq = 1, nqtot
    169     IF(tracers(iq)%isH2Ofamily) CYCLE
     169    IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    170170    itr = itr+1
    171171    SELECT CASE(tracers(iq)%name)

  • LMDZ6/trunk/libf/phylmd/Dust/splaeropt_lw_rrtm.F90

    r4056 r4071  
    5757    itr = 0
    5858    DO iq = 1, nqtot
    59       IF(tracers(iq)%isH2Ofamily) CYCLE
     59      IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    6060      itr = itr+1
    6161      SELECT CASE(tracers(iq)%name)

  • LMDZ6/trunk/libf/phylmd/Dust/splaerosol_optic_rrtm.F90

    r4056 r4071  
    5252  itr = 0
    5353  DO iq = 1, nqtot
    54     IF(tracers(iq)%isH2Ofamily) CYCLE
     54    IF(.NOT.tracers(iq)%isInPhysics) CYCLE
    5555    itr = itr+1
    5656    IF(tracers(iq)%name/='FINE') THEN
Note: See TracChangeset for help on using the changeset viewer.