Changeset 2077 for LMDZ5/trunk/libf/phylmd/fisrtilp.F90
- Timestamp:
- Jul 3, 2014, 3:38:05 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/fisrtilp.F90
r2006 r2077 13 13 ! 14 14 USE dimphy 15 USE microphys_mod ! cloud microphysics (JBM 3/14)15 USE icefrac_lsc_mod ! cloud microphysics (JBM 3/14) 16 16 IMPLICIT none 17 17 !====================================================================== … … 165 165 DATA appel1er /.TRUE./ 166 166 !ym 167 167 168 ice_thermo = iflag_ice_thermo .GE. 1 168 169 zdelq=0.0 … … 745 746 ! print*,zt(i),zrfl(i),zifl(i),'temp1' 746 747 ELSE ! of IF (iflag_t_glace.EQ.0) 748 CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1), & 749 & t_glace_min,t_glace_max,exposant_glace,zfice(:)) 747 750 if (iflag_fisrtilp_qsat.lt.1) then 748 751 DO i = 1, klon 749 752 ! JBM: icefrac_lsc is now a function contained in microphys_mod 750 751 753 ! zfice(i) = icefrac_lsc(zt(i), t_glace_min, & 754 ! t_glace_max, exposant_glace) 752 755 zt(i) = zt(i) + (1.-zfice(i))*zcond(i) * RLVTT/RCPD/(1.0+RVTMP2*zq(i)) & 753 756 +zfice(i)*zcond(i) * RLSTT/RCPD/(1.0+RVTMP2*zq(i)) … … 756 759 DO i=1, klon 757 760 ! JBM: icefrac_lsc is now a function contained in microphys_mod 758 759 761 ! zfice(i) = icefrac_lsc(zt(i), t_glace_min, & 762 ! t_glace_max, exposant_glace) 760 763 !CR: ATTENTION zt different de Tbef: à corriger 761 764 zt(i) = zt(i) + (1.-zfice(i))*zcond(i) * RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zcond(i))) & … … 790 793 ENDDO 791 794 ELSE ! of IF (iflag_t_glace.EQ.0) 792 DO i = 1, klon 793 IF (rneb(i,k).GT.0.0) THEN 795 CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1), & 796 & t_glace_min,t_glace_max,exposant_glace,zfice(:)) 797 ! DO i = 1, klon 798 ! IF (rneb(i,k).GT.0.0) THEN 794 799 ! JBM: icefrac_lsc is now a function contained in microphys_mod 795 796 797 798 800 ! zfice(i) = icefrac_lsc(zt(i), t_glace_min, & 801 ! t_glace_max, exposant_glace) 802 ! ENDIF 803 ! ENDDO 799 804 ENDIF 800 805 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.