- Timestamp:
- Jul 19, 2024, 5:41:58 PM (6 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_blowing_snow_sublim_sedim.F90
r4916 r5082 89 89 !---------------------- 90 90 91 IF (iflag_sedim_bs .GT.0) THEN91 IF (iflag_sedim_bs > 0) THEN 92 92 ! begin of top-down loop 93 93 DO k = nlay, 1, -1 … … 103 103 104 104 ! thermalization of blowing snow precip coming from above 105 IF (k .LE.nlay-1) THEN105 IF (k<=nlay-1) THEN 106 106 107 107 DO i = 1, ngrid … … 164 164 ! Sublimation and melting 165 165 !++++++++++++++++++++++++++++++++++++++++++++++ 166 IF (iflag_sublim_bs .GT.0) THEN166 IF (iflag_sublim_bs > 0) THEN 167 167 168 168 DO k = 1, nlay … … 189 189 190 190 191 IF (ztemp(i) .GT.RTT) THEN191 IF (ztemp(i) > RTT) THEN 192 192 193 193 ! if temperature is positive, we assume that part of the blowing snow … … 226 226 227 227 ! we do not consider deposition, only sublimation 228 IF (zqv(i) .LT.qsi(i)) THEN228 IF (zqv(i) < qsi(i)) THEN 229 229 rhoair=zpres(i)/ztemp(i)/RD 230 230 Dv=0.0001*0.211*(p0/zpres(i))*((ztemp(i)/RTT)**1.94) ! water vapor diffusivity in air, SI … … 256 256 ! see Gerber et al. 2023, JGR Atmos for the choice of qbmin 257 257 258 IF (zqb(i) .LT.qbmin) THEN258 IF (zqb(i) < qbmin) THEN 259 259 zqv(i) = zqv(i)+zqb(i) 260 IF (ztemp(i) .LT.RTT) THEN260 IF (ztemp(i) < RTT) THEN 261 261 ztemp(i) = ztemp(i) - zqb(i) * RLSTT/RCPD/(1.0+RVTMP2*(zqv(i))) 262 262 ELSE
Note: See TracChangeset
for help on using the changeset viewer.