- Timestamp:
- Jul 20, 2024, 12:00:23 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/climb_hq_mod.F90
r5082 r5087 574 574 if (iso_eau.gt.0) then 575 575 call iso_verif_egalite(xt_old(iso_eau,i,k), & 576 &q_old(i,k),'climb_hq_mod 421')576 q_old(i,k),'climb_hq_mod 421') 577 577 endif 578 578 #ifdef ISOTRAC … … 692 692 if (iso_eau.gt.0) then 693 693 call iso_verif_egalite_vect2D( & 694 &xt_new,q_new, &695 &'climb_hq_mod 504',ntraciso,klon,klev)694 xt_new,q_new, & 695 'climb_hq_mod 504',ntraciso,klon,klev) 696 696 endif !if (iso_eau.gt.0) then 697 697 if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then … … 700 700 if (q_new(i,k).gt.ridicule) then 701 701 if (iso_verif_o18_aberrant_nostop( & 702 &xt_new(iso_HDO,i,k)/q_new(i,k), &703 &xt_new(iso_O18,i,k)/q_new(i,k), &704 &'climb_hq_mod 690').eq.1) then702 xt_new(iso_HDO,i,k)/q_new(i,k), & 703 xt_new(iso_O18,i,k)/q_new(i,k), & 704 'climb_hq_mod 690').eq.1) then 705 705 write(*,*) 'i,k,q_new(i,k)=',i,k,q_new(i,k) 706 706 stop … … 773 773 if (iso_eau.gt.0) then 774 774 call iso_verif_egalite(d_xt(iso_eau,i,k), & 775 & d_q(i,k),'climb_hq_mod 503')775 d_q(i,k),'climb_hq_mod 503') 776 776 call iso_verif_egalite(xt_new(iso_eau,i,k), & 777 &q_new(i,k),'climb_hq_mod 503b')777 q_new(i,k),'climb_hq_mod 503b') 778 778 endif 779 779 #ifdef ISOTRAC
Note: See TracChangeset
for help on using the changeset viewer.