Changeset 5185 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spectral_definition.F90
- Timestamp:
- Sep 11, 2024, 4:27:07 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spectral_definition.F90
r5159 r5185 202 202 find_wavenumber = 1 203 203 DO while (wavenumber > this%wavenumber2(find_wavenumber) & 204 & . and. find_wavenumber < this%nwav)204 & .AND. find_wavenumber < this%nwav) 205 205 find_wavenumber = find_wavenumber + 1 206 206 end do … … 290 290 ! will be applicable 291 291 if (wavenumber(jwav) >= this%wavenumber1_band(jband) & 292 & . and. wavenumber(jwav) <= this%wavenumber2_band(jband)) then292 & .AND. wavenumber(jwav) <= this%wavenumber2_band(jband)) then 293 293 if (jwav > 1) then 294 294 wavenum1 = max(this%wavenumber1_band(jband), & … … 432 432 & / (this%wavenumber2(isd1)-this%wavenumber1(isd1)) 433 433 else 434 if (isd2 >= 1 . and. isd2 <= this%nwav) then434 if (isd2 >= 1 .AND. isd2 <= this%nwav) then 435 435 ! Right part of triangle 436 436 weight(isd2) = weight(isd2) + 0.5_jprb * (wavenum2-this%wavenumber1(isd2))**2 & … … 696 696 wavenumber2_bound = 0.01_jprb / wavelength_bound(jint-1) 697 697 where (wavenumber_mid > wavenumber1_bound & 698 & . and. wavenumber_mid <= wavenumber2_bound)698 & .AND. wavenumber_mid <= wavenumber2_bound) 699 699 i_input = i_intervals(jint) 700 700 end where
Note: See TracChangeset
for help on using the changeset viewer.