Skip to content

Commit e046818

Browse files
committed
bug fix - w/ aero on, code was satisfying this GFDL condition which is meaningless for GEOS
1 parent 1b2684b commit e046818

1 file changed

Lines changed: 3 additions & 3 deletions

File tree

GEOSaana_GridComp/GSI_GridComp/crtm_interface.f90

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -603,7 +603,6 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
603603
n_aerosols_jac_wk=0
604604
Load_AerosolCoeff=.false.
605605
endif
606-
if(mype==0) print *, 'DEBUG_RT: iasi, aero ', isis, n_actual_aerosols_wk, radmod%laerosol_fwd
607606

608607
! Initialize radiative transfer
609608

@@ -888,7 +887,8 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
888887
endif
889888

890889
! Initial GFDL saturation water vapor pressure tables
891-
if (n_actual_aerosols_wk>0 .or. n_clouds_fwd_wk>0 .and. imp_physics==11) then
890+
if (imp_physics==11) then
891+
if (n_actual_aerosols_wk>0 .or. n_clouds_fwd_wk>0) then
892892

893893
if (mype==0) write(6,*)myname_,':initial and load GFDL saturation water vapor pressure tables'
894894

@@ -909,6 +909,7 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
909909
des2 (length) = des2 (length - 1)
910910
desw (length) = desw (length - 1)
911911

912+
endif
912913
endif
913914

914915
return
@@ -2150,7 +2151,6 @@ subroutine call_crtm(obstype,obstime,data_s,nchanl,nreal,ich, &
21502151

21512152
! Set aerosols for CRTM
21522153
if(n_actual_aerosols_wk>0) then
2153-
print *, 'DEBUG_RT: crtm_interface: calling aero set'
21542154
call Set_CRTM_Aerosol ( msig, n_actual_aerosols_wk, n_aerosols_fwd_wk, aerosol_names, aero_conc, auxrh, &
21552155
atmosphere(1)%aerosol )
21562156
endif

0 commit comments

Comments
 (0)