- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_flux.F90
r4853 r5158 423 423 & config%i_band_from_reordered_g_sw, & 424 424 & this%sw_dn_surf_band, istartcol, iendcol) 425 dojcol = istartcol,iendcol425 DO jcol = istartcol,iendcol 426 426 this%sw_dn_surf_band(:,jcol) & 427 427 & = this%sw_dn_surf_band(:,jcol) & … … 429 429 end do 430 430 else 431 dojcol = istartcol,iendcol431 DO jcol = istartcol,iendcol 432 432 call indexed_sum(this%sw_dn_direct_surf_g(:,jcol), & 433 433 & config%i_band_from_reordered_g_sw, & … … 450 450 & config%i_band_from_reordered_g_sw, & 451 451 & this%sw_dn_surf_clear_band, istartcol, iendcol) 452 dojcol = istartcol,iendcol452 DO jcol = istartcol,iendcol 453 453 this%sw_dn_surf_clear_band(:,jcol) & 454 454 & = this%sw_dn_surf_clear_band(:,jcol) & … … 456 456 end do 457 457 else 458 dojcol = istartcol,iendcol458 DO jcol = istartcol,iendcol 459 459 call indexed_sum(this%sw_dn_direct_surf_clear_g(:,jcol), & 460 460 & config%i_band_from_reordered_g_sw, & … … 486 486 & this%sw_dn_diffuse_surf_canopy, istartcol, iendcol) 487 487 else 488 dojcol = istartcol,iendcol488 DO jcol = istartcol,iendcol 489 489 call indexed_sum(this%sw_dn_direct_surf_g(:,jcol), & 490 490 & config%i_albedo_from_band_sw(config%i_band_from_reordered_g_sw), & … … 502 502 this%sw_dn_diffuse_surf_canopy(:,istartcol:iendcol) = 0.0_jprb 503 503 this%sw_dn_direct_surf_canopy (:,istartcol:iendcol) = 0.0_jprb 504 dojband = 1,config%n_bands_sw505 dojalbedoband = 1,nalbedoband504 DO jband = 1,config%n_bands_sw 505 DO jalbedoband = 1,nalbedoband 506 506 if (config%sw_albedo_weights(jalbedoband,jband) /= 0.0_jprb) then 507 507 ! Initially, "diffuse" is actually "total" … … 534 534 & this%lw_dn_surf_canopy, istartcol, iendcol) 535 535 else 536 dojcol = istartcol,iendcol536 DO jcol = istartcol,iendcol 537 537 call indexed_sum(this%lw_dn_surf_g(:,jcol), & 538 538 & config%i_emiss_from_band_lw(config%i_band_from_reordered_g_lw), & … … 548 548 & lw_dn_surf_band, istartcol, iendcol) 549 549 else 550 dojcol = istartcol,iendcol550 DO jcol = istartcol,iendcol 551 551 call indexed_sum(this%lw_dn_surf_g(:,jcol), & 552 552 & config%i_band_from_reordered_g_lw, & … … 556 556 nalbedoband = size(config%lw_emiss_weights,1) 557 557 this%lw_dn_surf_canopy(:,istartcol:iendcol) = 0.0_jprb 558 dojband = 1,config%n_bands_lw559 dojalbedoband = 1,nalbedoband558 DO jband = 1,config%n_bands_lw 559 DO jalbedoband = 1,nalbedoband 560 560 if (config%lw_emiss_weights(jalbedoband,jband) /= 0.0_jprb) then 561 561 this%lw_dn_surf_canopy(jalbedoband,istartcol:iendcol) & … … 602 602 & this%sw_up_toa_band, istartcol, iendcol) 603 603 else 604 dojcol = istartcol,iendcol604 DO jcol = istartcol,iendcol 605 605 call indexed_sum(this%sw_dn_toa_g(:,jcol), & 606 606 & config%i_band_from_reordered_g_sw, & … … 618 618 & this%sw_up_toa_clear_band, istartcol, iendcol) 619 619 else 620 dojcol = istartcol,iendcol620 DO jcol = istartcol,iendcol 621 621 call indexed_sum(this%sw_up_toa_clear_g(:,jcol), & 622 622 & config%i_band_from_reordered_g_sw, & … … 634 634 & this%lw_up_toa_band, istartcol, iendcol) 635 635 else 636 dojcol = istartcol,iendcol636 DO jcol = istartcol,iendcol 637 637 call indexed_sum(this%lw_up_toa_g(:,jcol), & 638 638 & config%i_band_from_reordered_g_lw, & … … 647 647 & this%lw_up_toa_clear_band, istartcol, iendcol) 648 648 else 649 dojcol = istartcol,iendcol649 DO jcol = istartcol,iendcol 650 650 call indexed_sum(this%lw_up_toa_clear_g(:,jcol), & 651 651 & config%i_band_from_reordered_g_lw, & … … 753 753 iend = ubound(source,1) 754 754 755 dojg = istart, iend755 DO jg = istart, iend 756 756 ig = ind(jg) 757 757 dest(ig) = dest(ig) + source(jg) … … 777 777 iend = ubound(source,1) 778 778 779 dojg = istart, iend779 DO jg = istart, iend 780 780 ig = ind(jg) 781 781 dest(ig) = dest(ig) + source(jg) … … 797 797 dest = 0.0 798 798 799 dojg = lbound(source,1), ubound(source,1)799 DO jg = lbound(source,1), ubound(source,1) 800 800 ig = ind(jg) 801 dojc = ist, iend801 DO jc = ist, iend 802 802 dest(ig,jc) = dest(ig,jc) + source(jg,jc) 803 803 end do … … 820 820 nlev = size(source,2) 821 821 822 dojlev = 1,nlev823 dojg = istart, iend822 DO jlev = 1,nlev 823 DO jg = istart, iend 824 824 ig = ind(jg) 825 825 dest(ig,jlev) = dest(ig,jlev) + source(jg,jlev) … … 846 846 nlev = size(source,2) 847 847 848 dojlev = 1,nlev849 dojg = istart, iend848 DO jlev = 1,nlev 849 DO jg = istart, iend 850 850 ig = ind(jg) 851 851 dest(ig,jlev) = dest(ig,jlev) + source(jg,jlev)
Note: See TracChangeset
for help on using the changeset viewer.