Changeset 4143 for LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F
- Timestamp:
- May 9, 2022, 12:35:40 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F
r4052 r4143 13 13 USE parallel_lmdz 14 14 USE infotrac, ONLY : nqtot,tracers, ! CRisi & 15 & qperemin,masseqmin,ratiomin! MVals et CRisi15 & min_qParent,min_qMass,min_ratio ! MVals et CRisi 16 16 IMPLICIT NONE 17 17 c … … 346 346 DO ij=ijb,ije 347 347 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 348 masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),m asseqmin)349 if (q(ij,l,iq).gt. qperemin) then ! modif 13 nov 2020348 masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass) 349 if (q(ij,l,iq).gt.min_qParent) then ! modif 13 nov 2020 350 350 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 351 351 else 352 Ratio(ij,l,iq2)= ratiomin352 Ratio(ij,l,iq2)=min_ratio 353 353 endif 354 354 enddo … … 370 370 DO ij=ijb+1,ije 371 371 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 372 new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),m asseqmin)372 new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),min_qMass) 373 373 q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+ 374 374 & u_mq(ij-1,l)-u_mq(ij,l)) … … 423 423 USE parallel_lmdz 424 424 USE infotrac, ONLY : nqtot,tracers, ! CRisi & 425 & qperemin,masseqmin,ratiomin! MVals et CRisi425 & min_qParent,min_qMass,min_ratio ! MVals et CRisi 426 426 USE comconst_mod, ONLY: pi 427 427 IMPLICIT NONE … … 751 751 DO ij=ijbm,ijem 752 752 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 753 masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),m asseqmin)753 masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),min_qMass) 754 754 enddo !DO ij=ijbm,ijem 755 755 … … 758 758 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 759 759 !write(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq) 760 if (q(ij,l,iq).gt. qperemin) then ! modif 13 nov 2020760 if (q(ij,l,iq).gt.min_qParent) then ! modif 13 nov 2020 761 761 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 762 762 else 763 Ratio(ij,l,iq2)= ratiomin763 Ratio(ij,l,iq2)=min_ratio 764 764 endif 765 765 enddo !DO ij=ijbm,ijem
Note: See TracChangeset
for help on using the changeset viewer.