Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (5 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp.F90
r5082 r5086 731 731 modisRetrievedCloudTopPressure(i,:), & 732 732 modisRetrievedTau(i,:),modisRetrievedSize(i,:)) 733 end do733 END DO 734 734 endif 735 735 endif -
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/icarus.F90
r5082 r5086 155 155 (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev),& 156 156 (cchar(acc(ilev,ibox)+1),ilev=1,nlev) 157 end do157 END DO 158 158 close(9) 159 159 … … 322 322 do ibox=1,ncol 323 323 fluxtop(1:npoints,ibox) = fluxtop(1:npoints,ibox) + emsfc_lw*bb(1:npoints)*trans_layers_above(1:npoints,ibox) 324 end do324 END DO 325 325 326 326 ! All Sky brightness temperature … … 445 445 levmatch(1:npoints,ibox)=ilev 446 446 endwhere 447 end do447 END DO 448 448 end if 449 449 where(tau(1:npoints,ibox) <= tauchk) -
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90
r3358 r5086 209 209 exit 210 210 end if 211 end do211 END DO 212 212 213 213 if (lerror) then … … 244 244 end if 245 245 ilo = ilo + 1 246 end do246 END DO 247 247 248 248 ilo = max ( 2, ilo ) … … 254 254 end if 255 255 ihi = ihi - 1 256 end do256 END DO 257 257 258 258 ihi = min ( ihi, ntab - 1 ) … … 305 305 syl = x2 306 306 307 end do307 END DO 308 308 309 309 result = sum1 + ca * ( b**3 - syl**3 ) / 3._wp & -
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/modis_simulator.F90
r3358 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)) -
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90
r5082 r5086 539 539 if (alam < cutice) then 540 540 ! Region from 0.045 microns to 167.0 microns - no temperature depend 541 do i=2,nwl542 if(alam < wl(i)) continue543 enddo544 541 x1 = log(wl(i-1)) 545 542 x2 = log(wl(i)) -
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/quickbeam_optics.F90
r5082 r5086 1009 1009 call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), & 1010 1010 dg, xs1, xs2, dph, err) 1011 end do1011 END DO 1012 1012 1013 1013 else
Note: See TracChangeset
for help on using the changeset viewer.