- Timestamp:
- Sep 11, 2024, 4:27:07 PM (9 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spartacus_sw.F90
r5159 r5185 493 493 end if 494 494 495 if (config%do_3d_effects . and. &496 & allocated(cloud%inv_cloud_effective_size) . and. &497 & .not. (nreg == 2 . and. cloud%fraction(jcol,jlev) &495 if (config%do_3d_effects .AND. & 496 & allocated(cloud%inv_cloud_effective_size) .AND. & 497 & .not. (nreg == 2 .AND. cloud%fraction(jcol,jlev) & 498 498 & > 1.0-config%cloud_fraction_threshold)) then 499 499 if (cloud%inv_cloud_effective_size(jcol,jlev) > 0.0_jprb) then … … 662 662 ! 3D effects for any further g-points 663 663 if (ng3D == ng & 664 & . and. od_region(jg,1) > config%max_gas_od_3D) then664 & .AND. od_region(jg,1) > config%max_gas_od_3D) then 665 665 ng3D = jg-1 666 666 end if … … 935 935 if ((config%i_3d_sw_entrapment == IEntrapmentExplicitNonFractal & 936 936 & .or. config%i_3d_sw_entrapment == IEntrapmentExplicit) & 937 & . and. jlev >= i_cloud_top) then937 & .AND. jlev >= i_cloud_top) then 938 938 #else 939 939 if (config%i_3d_sw_entrapment == IEntrapmentExplicitNonFractal & … … 969 969 ! Account for cloud overlap when converting albedo and source 970 970 ! below a layer interface to the equivalent values just above 971 if (is_clear_sky_layer(jlev) . and. is_clear_sky_layer(jlev-1)) then971 if (is_clear_sky_layer(jlev) .AND. is_clear_sky_layer(jlev-1)) then 972 972 ! If both layers are cloud free, this is trivial... 973 973 total_albedo(:,:,:,jlev) = 0.0_jprb … … 1217 1217 & / max(config%cloud_fraction_threshold, region_fracs(jreg,jlev,jcol)) 1218 1218 DO jreg4 = 1,nreg ! VIA first lower region (jreg2 is second lower region) 1219 if (.not. (jreg4 == jreg . and. jreg4 /= jreg2)) then1219 if (.not. (jreg4 == jreg .AND. jreg4 /= jreg2)) then 1220 1220 albedo_part(:,jreg3,jreg) = albedo_part(:,jreg3,jreg) + entrapment(:,jreg3,jreg) & 1221 1221 & * v_matrix(jreg4,jreg,jlev,jcol) * total_albedo_below(:,jreg2,jreg4) … … 1305 1305 & / max(config%cloud_fraction_threshold, region_fracs(jreg,jlev,jcol)) 1306 1306 DO jreg4 = 1,nreg 1307 if (.not. (jreg4 == jreg . and. jreg4 /= jreg2)) then1307 if (.not. (jreg4 == jreg .AND. jreg4 /= jreg2)) then 1308 1308 albedo_part(:,jreg3,jreg) = albedo_part(:,jreg3,jreg) + entrapment(:,jreg3,jreg) & 1309 1309 & * v_matrix(jreg4,jreg,jlev,jcol) * total_albedo_below_direct(:,jreg2,jreg4) … … 1329 1329 if ((config%i_3d_sw_entrapment == IEntrapmentExplicitNonFractal & 1330 1330 & .or. config%i_3d_sw_entrapment == IEntrapmentExplicit) & 1331 & . and. .not. (is_clear_sky_layer(jlev) .and. is_clear_sky_layer(jlev-1))) then1331 & .AND. .not. (is_clear_sky_layer(jlev) .AND. is_clear_sky_layer(jlev-1))) then 1332 1332 ! Horizontal migration distances are averaged when 1333 1333 ! applying overlap rules, so equation is … … 1525 1525 ! Account for overlap rules in translating fluxes just above 1526 1526 ! a layer interface to the values just below 1527 if (is_clear_sky_layer(jlev) . and. is_clear_sky_layer(jlev+1)) then1527 if (is_clear_sky_layer(jlev) .AND. is_clear_sky_layer(jlev+1)) then 1528 1528 ! Regions in current layer map directly on to regions in 1529 1529 ! layer below
Note: See TracChangeset
for help on using the changeset viewer.