- 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.v1.5.1/radiation_monochromatic.F90
r4489 r5158 164 164 integer :: jlev 165 165 166 dojlev = 1,nlev166 DO jlev = 1,nlev 167 167 ! The fraction of the total optical depth in the current layer 168 168 ! is proportional to the fraction of the mass of the atmosphere … … 192 192 & = StefanBoltzmann * single_level%skin_temperature(istartcol:iendcol)**4 & 193 193 & * single_level%lw_emissivity(istartcol:iendcol,1) 194 dojlev = 1,nlev+1194 DO jlev = 1,nlev+1 195 195 planck_hl(1,jlev,istartcol:iendcol) = StefanBoltzmann * thermodynamics%temperature_hl(istartcol:iendcol,jlev)**4 196 196 end do … … 200 200 & single_level%skin_temperature(istartcol:iendcol)) & 201 201 & * single_level%lw_emissivity(istartcol:iendcol,1) 202 dojlev = 1,nlev+1202 DO jlev = 1,nlev+1 203 203 planck_hl(1,jlev,istartcol:iendcol) = Pi*planck_function(config%mono_lw_wavelength, & 204 204 & thermodynamics%temperature_hl(istartcol:iendcol,jlev)) … … 261 261 ! Convert cloud mixing ratio into liquid and ice water path 262 262 ! in each layer 263 dojlev = 1, nlev264 dojcol = istartcol, iendcol263 DO jlev = 1, nlev 264 DO jcol = istartcol, iendcol 265 265 ! Factor to convert from gridbox-mean mass mixing ratio to 266 266 ! in-cloud water path involves the pressure difference in … … 288 288 289 289 if (config%iverbose >= 4) then 290 dojcol = istartcol,iendcol290 DO jcol = istartcol,iendcol 291 291 write(*,'(a,i0,a,f7.3,a,f7.3)') 'Profile ', jcol, ': shortwave optical depth = ', & 292 292 & sum(od_sw_cloud(1,:,jcol)*cloud%fraction(jcol,:)), &
Note: See TracChangeset
for help on using the changeset viewer.