Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_old.F90
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_old.F90
r5082 r5087 709 709 zdqs(i) = FOEDE(zt(i),zdelta,zcvm5,zqs(i),zcor) 710 710 zdqsdT_raw(i) = zdqs(i)* & 711 &RCPD*(1.0+RVTMP2*zq(i)) / (RLVTT*(1.-zdelta) + RLSTT*zdelta)711 RCPD*(1.0+RVTMP2*zq(i)) / (RLVTT*(1.-zdelta) + RLSTT*zdelta) 712 712 ENDDO 713 713 ELSE … … 966 966 zfice(i) = zfice(i)**exposant_glace_old 967 967 dzfice(i)= exposant_glace_old * zfice(i)**(exposant_glace_old-1) & 968 &/ (t_glace_min_old - RTT)968 / (t_glace_min_old - RTT) 969 969 endif 970 970 971 971 if (iflag_t_glace>=1.and.zfice(i)>0.) then 972 972 dzfice(i)= exposant_glace * zfice(i)**(exposant_glace-1) & 973 &/ (t_glace_min - t_glace_max)973 / (t_glace_min - t_glace_max) 974 974 endif 975 975 … … 987 987 if (fl_cor_ebil > 0) then 988 988 num = -Tbef(i)+zt(i)+rneb(i,k)*((1-zfice(i))*RLVTT & 989 &+zfice(i)*RLSTT)/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))*qlbef(i)989 +zfice(i)*RLSTT)/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))*qlbef(i) 990 990 denom = 1.+rneb(i,k)*((1-zfice(i))*RLVTT+zfice(i)*RLSTT)/cste*zdqs(i) & 991 991 -(RLSTT-RLVTT)/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))*rneb(i,k) & 992 &*qlbef(i)*dzfice(i)992 *qlbef(i)*dzfice(i) 993 993 else 994 994 num = -Tbef(i)+zt(i)+rneb(i,k)*((1-zfice(i))*RLVTT & 995 &+zfice(i)*RLSTT)/RCPD/(1.0+RVTMP2*zq(i))*qlbef(i)995 +zfice(i)*RLSTT)/RCPD/(1.0+RVTMP2*zq(i))*qlbef(i) 996 996 denom = 1.+rneb(i,k)*((1-zfice(i))*RLVTT+zfice(i)*RLSTT)/cste*zdqs(i) & 997 997 -(RLSTT-RLVTT)/RCPD/(1.0+RVTMP2*zq(i))*rneb(i,k)*qlbef(i)*dzfice(i) … … 1151 1151 if (fl_cor_ebil > 0) then 1152 1152 zt(i) = zt(i) + (1.-zfice(i))*zcond(i) & 1153 &* RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)+zcond(i))) &1153 * RLVTT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)+zcond(i))) & 1154 1154 +zfice(i)*zcond(i) * RLSTT/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)+zcond(i))) 1155 1155 else
Note: See TracChangeset
for help on using the changeset viewer.