- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_cloud_optics_prop.F90
r4715 r5082 175 175 reice_pi = 0. 176 176 177 IF (iflag_t_glace .EQ.0) THEN177 IF (iflag_t_glace==0) THEN 178 178 DO k = 1, klev 179 179 DO i = 1, klon … … 245 245 !--flag_aerosol=7 => MACv2SP climatology 246 246 !--in this case there is an enhancement factor 247 IF (flag_aerosol .EQ. 7) THEN247 IF (flag_aerosol == 7) THEN 248 248 249 249 !--present-day … … 300 300 zfiwp_var = 1000.*icefrac_optics(i, k)*radocond(i, k)/pclc(i, k)*rhodz(i, k) 301 301 ! Calculation of ice cloud effective radius in micron 302 IF (iflag_rei .EQ.1) THEN302 IF (iflag_rei == 1) THEN 303 303 ! when we account for precipitation in the radiation scheme, 304 304 ! It is recommended to use the rei formula from Sun and Rikkus 1999 with a revision … … 407 407 408 408 409 IF (iflag_rei .GT.0) THEN409 IF (iflag_rei > 0) THEN 410 410 411 411 ! when we account for precipitation in the radiation scheme, … … 600 600 reffclwtop(i) = 0. 601 601 cldncl(i) = 0. 602 IF (novlp .EQ.3 .OR. novlp.EQ.1) tcc(i) = 1.603 IF (novlp .EQ.2) tcc(i) = 0.602 IF (novlp==3 .OR. novlp==1) tcc(i) = 1. 603 IF (novlp==2) tcc(i) = 0. 604 604 ENDDO 605 605 … … 612 612 IF (pcltau(i,k)>thres_tau .AND. pclc(i,k)>thres_neb) THEN 613 613 614 IF (novlp .EQ.2) THEN614 IF (novlp==2) THEN 615 615 IF (first) THEN 616 616 WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM' … … 621 621 ENDIF 622 622 623 IF (novlp .EQ.3) THEN623 IF (novlp==3) THEN 624 624 IF (first) THEN 625 625 WRITE (*, *) 'Hypothese de recouvrement: RANDOM' … … 630 630 ENDIF 631 631 632 IF (novlp .EQ.1) THEN632 IF (novlp==1) THEN 633 633 IF (first) THEN 634 634 WRITE (*, *) 'Hypothese de recouvrement: MAXIMUM_ & … … 655 655 ENDDO ! loop over k 656 656 657 IF (novlp .EQ.3 .OR. novlp.EQ.1) tcc(i) = 1. - tcc(i)657 IF (novlp==3 .OR. novlp==1) tcc(i) = 1. - tcc(i) 658 658 659 659 ENDDO ! loop over i
Note: See TracChangeset
for help on using the changeset viewer.