- 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/radiation/radiation_mcica_lw.F90
r4946 r5158 146 146 147 147 ! Loop through columns 148 dojcol = istartcol,iendcol148 DO jcol = istartcol,iendcol 149 149 150 150 ! Clear-sky calculation … … 203 203 is_clear_sky_layer = .true. 204 204 i_cloud_top = nlev+1 205 dojlev = 1,nlev205 DO jlev = 1,nlev 206 206 ! Compute combined gas+aerosol+cloud optical properties 207 207 if (cloud%fraction(jcol,jlev) >= config%cloud_fraction_threshold) then … … 212 212 end if 213 213 214 dojg = 1,ng214 DO jg = 1,ng 215 215 od_cloud_new(jg) = od_scaling(jg,jlev) & 216 216 & * od_cloud(config%i_band_from_reordered_g_lw(jg),jlev,jcol) … … 228 228 ! case that od_total > 0.0 and ssa_total > 0.0 but 229 229 ! od_total*ssa_total == 0 due to underflow 230 dojg = 1,ng230 DO jg = 1,ng 231 231 if (od_total(jg) > 0.0_jprb) then 232 232 scat_od_total(jg) = ssa(jg,jlev,jcol)*od(jg,jlev,jcol) & … … 247 247 else 248 248 249 dojg = 1,ng249 DO jg = 1,ng 250 250 if (od_total(jg) > 0.0_jprb) then 251 251 scat_od = ssa_cloud(config%i_band_from_reordered_g_lw(jg),jlev,jcol) &
Note: See TracChangeset
for help on using the changeset viewer.