- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_mcica_sw.F90
r4489 r5158 142 142 143 143 ! Loop through columns 144 dojcol = istartcol,iendcol144 DO jcol = istartcol,iendcol 145 145 ! Only perform calculation if sun above the horizon 146 146 if (single_level%cos_sza(jcol) > 0.0_jprb) then … … 152 152 ! Delta-Eddington scaling has already been performed to the 153 153 ! aerosol part of od, ssa and g 154 dojlev = 1,nlev154 DO jlev = 1,nlev 155 155 call calc_two_stream_gammas_sw(ng, & 156 156 & cos_sza, ssa(:,jlev,jcol), g(:,jlev,jcol), & … … 165 165 else 166 166 ! Apply delta-Eddington scaling to the aerosol-gas mixture 167 dojlev = 1,nlev167 DO jlev = 1,nlev 168 168 od_total = od(:,jlev,jcol) 169 169 ssa_total = ssa(:,jlev,jcol) … … 219 219 if (total_cloud_cover >= config%cloud_fraction_threshold) then 220 220 ! Total-sky calculation 221 dojlev = 1,nlev221 DO jlev = 1,nlev 222 222 ! Compute combined gas+aerosol+cloud optical properties 223 223 if (cloud%fraction(jcol,jlev) >= config%cloud_fraction_threshold) then 224 dojg = 1,ng224 DO jg = 1,ng 225 225 od_cloud_new(jg) = od_scaling(jg,jlev) & 226 226 & * od_cloud(config%i_band_from_reordered_g_sw(jg),jlev,jcol)
Note: See TracChangeset
for help on using the changeset viewer.