Ignore:
Timestamp:
Mar 2, 2017, 8:46:14 PM (8 years ago)
Author:
oboucher
Message:

Adding od443aer diagnostic for CMIP6.
Removing od670aer and AI calculations as they are not used.

File:
1 edited

Legend:

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

    r2748 r2806  
    1212  USE DIMPHY
    1313  USE aero_mod
    14   USE phys_local_var_mod, ONLY: od550aer,od865aer,ec550aer,od550lt1aer
     14  USE phys_local_var_mod, ONLY: od443aer,od550aer,od865aer,ec550aer,od550lt1aer
    1515  USE YOMCST, ONLY: RD,RG
    1616
     
    366366    DO la=1,las
    367367
    368     !--only 550, 670 and 865 nm are used
    369     IF (la.NE.la550.AND.la.NE.la670.AND.la.NE.la865) CYCLE
     368    !--only 443, 550, 670 and 865 nm are used
     369    !--to save time 670 and AI are not computed for CMIP6
     370    !IF (la.NE.la443.AND.la.NE.la550.AND.la.NE.la670.AND.la.NE.la865) CYCLE
     371    IF (la.NE.la443.AND.la.NE.la550.AND.la.NE.la865) CYCLE
    370372
    371373      IF (soluble) THEN            ! For soluble aerosol
     
    433435
    434436!--AOD calculations for diagnostics
     437  od443aer(:)=SUM(tausum(:,la443,:),dim=2)
    435438  od550aer(:)=SUM(tausum(:,la550,:),dim=2)
    436   od670aer(:)=SUM(tausum(:,la670,:),dim=2)
     439  !od670aer(:)=SUM(tausum(:,la670,:),dim=2)
    437440  od865aer(:)=SUM(tausum(:,la865,:),dim=2)
    438441
     
    441444
    442445!--aerosol index
    443   ai(:)=-LOG(MAX(od670aer(:),1.e-8)/MAX(od865aer(:),1.e-8))/LOG(670./865.)
     446  ai(:)=0.0
     447  !ai(:)=-LOG(MAX(od670aer(:),1.e-8)/MAX(od865aer(:),1.e-8))/LOG(670./865.)
    444448
    445449  od550lt1aer(:)=tausum(:,la550,id_ASSO4M_phy)+tausum(:,la550,id_ASBCM_phy) +tausum(:,la550,id_AIBCM_phy)+ &
Note: See TracChangeset for help on using the changeset viewer.