Changeset 5185 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_spectral_definition.F90
- Timestamp:
- Sep 11, 2024, 4:27:07 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_spectral_definition.F90
r5159 r5185 171 171 find_wavenumber = 1 172 172 DO while (wavenumber > this%wavenumber2(find_wavenumber) & 173 & . and. find_wavenumber < this%nwav)173 & .AND. find_wavenumber < this%nwav) 174 174 find_wavenumber = find_wavenumber + 1 175 175 end do … … 250 250 ! will be applicable 251 251 if (wavenumber(jwav) >= this%wavenumber1_band(jband) & 252 & . and. wavenumber(jwav) <= this%wavenumber2_band(jband)) then252 & .AND. wavenumber(jwav) <= this%wavenumber2_band(jband)) then 253 253 if (jwav > 1) then 254 254 wavenum1 = max(this%wavenumber1_band(jband), & … … 388 388 & / (this%wavenumber2(isd1)-this%wavenumber1(isd1)) 389 389 else 390 if (isd2 >= 1 . and. isd2 <= this%nwav) then390 if (isd2 >= 1 .AND. isd2 <= this%nwav) then 391 391 ! Right part of triangle 392 392 weight(isd2) = weight(isd2) + 0.5_jprb * (wavenum2-this%wavenumber1(isd2))**2 & … … 647 647 wavenumber2_bound = 0.01_jprb / wavelength_bound(jint-1) 648 648 where (wavenumber_mid > wavenumber1_bound & 649 & . and. wavenumber_mid <= wavenumber2_bound)649 & .AND. wavenumber_mid <= wavenumber2_bound) 650 650 i_input = i_intervals(jint) 651 651 end where
Note: See TracChangeset
for help on using the changeset viewer.