Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90
r5086 r5087 21 21 epmax_diag & ! epmax_cape 22 22 #ifdef ISO 23 &,xt,xt_wake,d_xt,xtrain,xtsnow &24 &,xtVprecip,xtVprecipi &25 &,xtclw,dd_xt,xtevap,xtwdtrainA &23 ,xt,xt_wake,d_xt,xtrain,xtsnow & 24 ,xtVprecip,xtVprecipi & 25 ,xtclw,dd_xt,xtevap,xtwdtrainA & 26 26 #ifdef DIAGISO 27 & , qlp,xtlp,qvp,xtvp & ! juste diagnostique28 &, fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip &29 &, fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip &30 &, f_detrainement,q_detrainement,xt_detrainement &27 , qlp,xtlp,qvp,xtvp & ! juste diagnostique 28 , fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip & 29 , fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip & 30 , f_detrainement,q_detrainement,xt_detrainement & 31 31 #endif 32 32 #endif 33 &) ! **************************************************************33 ) ! ************************************************************** 34 34 ! * 35 35 ! CONCVL * … … 482 482 do i = 1, klon 483 483 call iso_verif_egalite_choix(xt(iso_eau,i,k),q(i,k), & 484 & 'concvl 174',errmax,errmaxrel)484 'concvl 174',errmax,errmaxrel) 485 485 enddo !do i = 1, klon 486 486 enddo !do k = 1, klev … … 520 520 epmax_diag & ! epmax_cape 521 521 #ifdef ISO 522 &, xt,d_xt &523 &, xtrain,xtVprecip &524 &, xtevap,xtclw,xtwdtrainA &522 , xt,d_xt & 523 , xtrain,xtVprecip & 524 , xtevap,xtclw,xtwdtrainA & 525 525 #ifdef DIAGISO 526 & , qlp,xtlp,qvp,xtvp & ! juste diagnostique527 &, fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip &528 &, fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip &529 &, f_detrainement,q_detrainement,xt_detrainement &526 , qlp,xtlp,qvp,xtvp & ! juste diagnostique 527 , fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip & 528 , fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip & 529 , f_detrainement,q_detrainement,xt_detrainement & 530 530 #endif 531 531 #endif 532 &)532 ) 533 533 ! print *, 'cv_driver ->' !jyg 534 534 … … 543 543 call iso_verif_egalite(xt_wake(iso_eau,i,k),q_wake(i,k),'concvl 204') 544 544 call iso_verif_egalite(d_xt(iso_eau,i,k),d_q(i,k), & 545 &'concvl 452')545 'concvl 452') 546 546 endif !if (iso_eau.gt.0) then 547 547 #ifdef DIAGISO … … 557 557 call iso_verif_traceur_vect(xt,klon,klev,'concvl 218') 558 558 call iso_verif_trac_masse_vect(d_xt,klon,klev, & 559 &'concvl 464',errmax,errmaxrel)559 'concvl 464',errmax,errmaxrel) 560 560 #endif 561 561 #endif … … 604 604 epmax_diag & ! epmax_cape 605 605 #ifdef ISO 606 &,xt,xt_wake,d_xt, xtrain &607 &,xtvprecip,xtvprecipi &608 &,xtclw,dd_xt,xtevap,xtwdtrainA &606 ,xt,xt_wake,d_xt, xtrain & 607 ,xtvprecip,xtvprecipi & 608 ,xtclw,dd_xt,xtevap,xtwdtrainA & 609 609 #ifdef DIAGISO 610 &, qlp,xtlp,qvp,xtvp &611 &, fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip &612 &, fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip &613 &, f_detrainement,q_detrainement,xt_detrainement &610 , qlp,xtlp,qvp,xtvp & 611 , fq_detrainement,fq_ddft,fq_fluxmasse,fq_evapprecip & 612 , fxt_detrainement,fxt_ddft,fxt_fluxmasse,fxt_evapprecip & 613 , f_detrainement,q_detrainement,xt_detrainement & 614 614 #endif 615 615 #endif 616 & )616 ) 617 617 END IF 618 618 ! ------------------------------------------------------------------ … … 658 658 if (q(i,k).gt.ridicule) then 659 659 call iso_verif_aberrant_encadre((xt(iso_HDO,i,k) & 660 &+d_xt(iso_HDO,i,k))/(q(i,k)+d_q(i,k)),'concvl 250')660 +d_xt(iso_HDO,i,k))/(q(i,k)+d_q(i,k)),'concvl 250') 661 661 endif !if (q_seri(i,k).gt.ridicule) then 662 662 endif !if (iso_HDO.gt.0) then 663 663 if (iso_eau.gt.0) then 664 664 call iso_verif_egalite_choix(d_xt(iso_eau,i,k), & 665 &d_q(i,k),'concvl 530',errmax*dtime,errmaxrel)665 d_q(i,k),'concvl 530',errmax*dtime,errmaxrel) 666 666 endif !if (iso_HDO.gt.0) then 667 667 #ifdef ISOTRAC … … 682 682 #ifdef ISOTRAC 683 683 call iso_verif_traceur_jbid_vect(d_xt, & 684 & klon,klev)684 klon,klev) 685 685 #endif 686 686 endif !if ((iso_eau.gt.0).and.(bidouille_anti_divergence)) then
Note: See TracChangeset
for help on using the changeset viewer.