- Timestamp:
- Nov 5, 2018, 3:24:59 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
/LMDZ6/trunk removed
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cosp/lidar_simulator.F90
r3356 r3411 121 121 !--------------------------------------------------------------------------------- 122 122 123 USE MOD_COSP_CONSTANTS, only : ok_debug_cosp124 123 IMPLICIT NONE 125 124 REAL :: SRsat … … 202 201 REAL pnorm_perp_liq(npoints,nlev) ! perpendicular lidar backscattered signal power for liq 203 202 204 REAL :: seuil205 206 203 ! Output variable 207 204 REAL pnorm_perp_tot (npoints,nlev) ! perpendicular lidar backscattered signal power … … 238 235 Gamma = 0.009 239 236 240 if (ok_debug_cosp) then241 seuil=1.e-18242 else243 seuil=0.0244 endif245 237 !------------------------------------------------------------ 246 238 !---- 1. Preliminary definitions and calculations : … … 502 494 503 495 DO k= nlev, 1, -1 504 505 496 pnorm_perp_ice(:,k) = Alpha * pnorm_ice(:,k) ! Ice particles 497 pnorm_perp_liq(:,k) = 1000*Beta * pnorm_liq(:,k)**2 + Gamma * pnorm_liq(:,k) ! Liquid particles 506 498 ENDDO 507 499 … … 531 523 tautot_lay_liq(:) = tautot_liq(:,k)-tautot_liq(:,k+1) 532 524 WHERE (tautot_lay_liq(:).GT.0.) 533 beta_perp_liq(:,k) = pnorm_perp_liq(:,k)/ max(seuil,EXP(-2.0*tautot_liq(:,k+1))) &534 & * (2.*tautot_lay_liq(:))/ (1.-exp(-2.0*tautot_lay_liq(:)))525 beta_perp_liq(:,k) = pnorm_perp_liq(:,k)/ EXP(-2.0*tautot_liq(:,k+1)) * (2.*tautot_lay_liq(:)) & 526 & / (1.-exp(-2.0*tautot_lay_liq(:))) 535 527 536 528 ELSEWHERE
Note: See TracChangeset
for help on using the changeset viewer.