- Timestamp:
- Jul 19, 2024, 7:54:50 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90
r3491 r5086 222 222 retrievedTau(i) = R_UNDEF 223 223 end if 224 end do224 END DO 225 225 where((retrievedSize(1:nSubCols) < 0.).and.(retrievedSize(1:nSubCols) /= R_UNDEF)) & 226 226 retrievedSize(1:nSubCols) = 1.0e-06_wp*re_fill … … 455 455 end if 456 456 if(totalTau >= tauLimit) exit 457 end do457 END DO 458 458 459 459 if (totalTau > 0._wp) then … … 489 489 end if 490 490 if(totalTau >= tauLimit) exit 491 end do491 END DO 492 492 493 493 if (totalTau > 0._wp) then … … 715 715 do i = 1, size(cloudIndicies) 716 716 call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i)) 717 end do717 END DO 718 718 719 719 call adding_doubling(count(tau(1:nLevels) > 0),Refl(:), Trans(:), Refl_tot, Trans_tot) … … 897 897 Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i)) 898 898 Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1._wp - Refl_cumulative(i-1) * Refl(i)) 899 end do899 END DO 900 900 901 901 Refl_tot = Refl_cumulative(size(Refl))
Note: See TracChangeset
for help on using the changeset viewer.