Changeset 2257 for trunk/LMDZ.MARS/libf/phymars
- Timestamp:
- Mar 12, 2020, 9:00:47 AM (5 years ago)
- Location:
- trunk/LMDZ.MARS/libf/phymars
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/co2cloud.F
r2156 r2257 840 840 No=Nccnco2*tauscaling(ig) 841 841 Rn=-dlog(riceco2(ig,l)) 842 n_derf = derf( (rb_cldco2(1)+Rn) *dev2)842 n_derf = erf( (rb_cldco2(1)+Rn) *dev2) 843 843 Qext1bins2(ig,l)=0. 844 844 do i = 1, nbinco2_cld 845 845 n_aer(i) = -0.5 * No * n_derf !! this ith previously computed 846 n_derf = derf((rb_cldco2(i+1)+Rn) *dev2)846 n_derf = erf((rb_cldco2(i+1)+Rn) *dev2) 847 847 n_aer(i) = n_aer(i) + 0.5 * No * n_derf 848 848 Qext1bins2(ig,l)=Qext1bins2(ig,l)+Qext1bins(i)*n_aer(i) -
trunk/LMDZ.MARS/libf/phymars/improvedco2clouds_mod.F
r2156 r2257 552 552 c Ic_rice Mass transfer rate (kg/s) for a rice particle >0 si croissance ! 553 553 554 if (isnan(Ic_rice) .or. Ic_rice .eq. 0.) then 554 if ( (Ic_rice.ne.Ic_rice) ! will be true if it is Nan 555 & .or. Ic_rice .eq. 0.) then 555 556 Ic_rice=0. 556 557 subpdtcloudco2(ig,l)=-sum_subpdt(ig,l)
Note: See TracChangeset
for help on using the changeset viewer.