Ignore:
Timestamp:
Dec 17, 2016, 8:48:54 PM (8 years ago)
Author:
oboucher
Message:

Complement to previous commit. Fixing an issue with mass array index.
This does not have any impact because ASSO4M and CSSO4M have same properties.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90

    r2735 r2736  
    9999  REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
    100100  !
    101   INTEGER :: id, classbc, classbc_pi
     101  INTEGER :: aerindex, classbc, classbc_pi
    102102  REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
    103103
     
    584584
    585585    !--shortname for aerosol index
    586     id=aerosol_name(m)
     586    aerindex=aerosol_name(m)
    587587
    588588    IF (soluble) THEN
     
    594594           DO i=1, klon
    595595
    596              tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
    597              tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
     596             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
     597             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
    598598
    599599             ! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols
     
    642642                             cg_MG_6bands(RH_num(i,k),inu,classbc))
    643643
    644                tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
    645                piz_ae(i,k,id,inu)    = piz_ae2b_int
    646                cg_ae(i,k,id,inu)     = cg_ae2b_int
     644               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
     645               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
     646               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
    647647
    648648               !--natural aerosols
     
    659659                             cg_MG_6bands(RH_num(i,k),inu,classbc_pi))
    660660
    661                tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
    662                piz_ae_pi(i,k,id,inu) = piz_ae2b_int
    663                cg_ae_pi(i,k,id,inu)  = cg_ae2b_int
     661               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
     662               piz_ae_pi(i,k,aerindex,inu) = piz_ae2b_int
     663               cg_ae_pi(i,k,aerindex,inu)  = cg_ae2b_int
    664664                       
    665665             ENDDO
     
    672672         DO k=1, klev
    673673           DO i=1, klon
    674              tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
    675              tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
     674             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
     675             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
    676676
    677677             DO inu=1,NSW
     
    689689                             cg_aers_6bands(RH_num(i,k),inu,spsol))
    690690
    691                tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
    692                tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
    693                piz_ae(i,k,id,inu)    = piz_ae2b_int
    694                cg_ae(i,k,id,inu)     = cg_ae2b_int
     691               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
     692               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
     693               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
     694               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
    695695                       
    696696             ENDDO
     
    700700         !--external mixture case for soluble BC
    701701         IF (spsol.EQ.1) THEN
    702            piz_ae_pi(:,:,id,:) = piz_ae(:,:,id,:)
    703            cg_ae_pi(:,:,id,:)  = cg_ae(:,:,id,:)
     702           piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:)
     703           cg_ae_pi(:,:,aerindex,:)  = cg_ae(:,:,aerindex,:)
    704704         ENDIF
    705705
     
    710710       DO k=1, klev
    711711         DO i=1, klon
    712            tmp_var=m_allaer(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
    713            tmp_var_pi=m_allaer_pi(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
     712           tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
     713           tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
    714714
    715715           DO inu=1,NSW
     
    718718             cg_ae2b_int = cg_aeri_6bands(inu,spinsol)
    719719
    720              tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
    721              tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int
    722              piz_ae(i,k,id,inu) = piz_ae2b_int
    723              cg_ae(i,k,id,inu)= cg_ae2b_int
     720             tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
     721             tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi*tau_ae2b_int
     722             piz_ae(i,k,aerindex,inu) = piz_ae2b_int
     723             cg_ae(i,k,aerindex,inu)= cg_ae2b_int
    724724           ENDDO
    725725         ENDDO
Note: See TracChangeset for help on using the changeset viewer.