Changeset 4491 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
- Timestamp:
- Apr 3, 2023, 10:25:44 AM (20 months ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4479 r4491 157 157 & iso_verif_aberrant_choix,iso_verif_positif, & 158 158 & iso_verif_positif_choix_vect,iso_verif_o18_aberrant_nostop, & 159 & iso_verif_init, &159 & iso_verif_init,iso_verif_aberrant_enc_choix_nostop,& 160 160 & iso_verif_positif_strict_nostop,iso_verif_O18_aberrant_enc_vect2D 161 161 #endif … … 2480 2480 #ifdef ISO 2481 2481 #ifdef ISOVERIF 2482 ! write(*,*) 'physiq 1847: qx(1,1,:)=',qx(1,1,:) 2482 write(*,*) 'physiq 1847: qx(1,1,:)=',qx(1,1,:) 2483 write(*,*) 'iqIsoPha(:,ivap)=',iqIsoPha(:,ivap) 2483 2484 write(*,*) 'physiq 1846b: ok_isotopes,ntraciso,niso=',niso>0,ntraciso,niso 2484 2485 #endif … … 2596 2597 & 'physiq 2099 ql').eq.1) then 2597 2598 write(*,*) 'i,k,ql_seri(i,k)=',i,k,ql_seri(i,k) 2598 stop2599 !stop 2599 2600 endif ! if (iso_verif_o18_aberrant_nostop 2600 2601 endif !if (q_seri(i,k).gt.errmax) then … … 2605 2606 & 'physiq 2099 qs').eq.1) then 2606 2607 write(*,*) 'i,k,qs_seri(i,k)=',i,k,qs_seri(i,k) 2607 stop2608 !stop 2608 2609 endif ! if (iso_verif_o18_aberrant_nostop 2609 2610 endif !if (q_seri(i,k).gt.errmax) then 2610 2611 enddo !k=1,klev 2611 2612 enddo !i=1,klon 2612 endif !if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then 2613 endif !if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then 2614 #ifdef ISOTRAC 2615 DO k = 1, klev 2616 DO i = 1, klon 2617 call iso_verif_traceur(xt_seri(1,i,k),'physiq 2620a') 2618 call iso_verif_traceur(xtl_seri(1,i,k),'physiq 2620b') 2619 call iso_verif_traceur(xts_seri(1,i,k),'physiq 2620c') 2620 enddo 2621 enddo 2622 #endif 2613 2623 #endif 2614 2624 ! … … 2844 2854 ENDIF 2845 2855 ENDIF 2856 2857 #ifdef ISO 2858 #ifdef ISOVERIF 2859 #ifdef ISOTRAC 2860 DO k = 1, klev 2861 DO i = 1, klon 2862 call iso_verif_traceur(xt_seri(1,i,k), & 2863 & 'physiq 2856: avant reevp') 2864 enddo 2865 enddo 2866 #endif 2867 #endif 2868 #endif 2846 2869 ! 2847 2870 ! Re-evaporer l'eau liquide nuageuse … … 3458 3481 call iso_verif_egalite(q_w(i,k),xt_w(iso_eau,i,k),'physiq 3338') 3459 3482 endif 3483 if (iso_HDO.gt.0) then 3484 if ((iso_verif_aberrant_enc_choix_nostop(xt_x(iso_hdo,i,k),q_x(i,k), & 3485 ridicule,deltalim,'physic 3462a xt_x').eq.1).or. & 3486 (iso_verif_aberrant_enc_choix_nostop(xt_w(iso_hdo,i,k),q_w(i,k), & 3487 ridicule,deltalim,'physic 3462b xt_w').eq.1)) then 3488 write(*,*) 'i,k=',i,k 3489 write(*,*) 'q_x(i,k),q_seri(i,k),wake_s(i),wake_deltaq(i,k)=', & 3490 q_x(i,k),q_seri(i,k),wake_s(i),wake_deltaq(i,k) 3491 write(*,*) 'xt_x(iso_hdo,i,k),xt_seri(iso_hdo,i,k),wake_s(i),wake_deltaxt(iso_hdo,i,k)=', & 3492 xt_x(iso_hdo,i,k),xt_seri(iso_hdo,i,k),wake_s(i),wake_deltaxt(iso_hdo,i,k) 3493 write(*,*) 'deltaD_seri,wake=',deltaD(xt_seri(iso_hdo,i,k)/q_seri(i,k)), & 3494 deltaD(wake_deltaxt(iso_hdo,i,k)/wake_deltaq(i,k)) 3495 write(*,*) 'deltaD_x,deltaD_w=',deltaD(xt_x(iso_hdo,i,k)/q_x(i,k)),deltaD(xt_w(iso_hdo,i,k)/q_w(i,k)) 3496 stop 3497 endif 3498 endif 3460 3499 #endif 3461 3500 #endif … … 3590 3629 call iso_verif_aberrant_encadre( & 3591 3630 & xt_w(iso_hdo,i,k)/q_w(i,k), & 3592 & 'physic 2657 b')3631 & 'physic 2657c') 3593 3632 endif !if (q_x(i,k).gt.ridicule) then 3594 3633 endif !if (iso_HDO.gt.0) then … … 4255 4294 endif !if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then 4256 4295 if ((iso_O17.gt.0).and.(iso_O18.gt.0)) then 4257 if ((q_seri(i,k).gt.ridicule).and.( l.lt.nlevmaxO17)) then4296 if ((q_seri(i,k).gt.ridicule).and.(k.lt.nlevmaxO17)) then 4258 4297 call iso_verif_aberrant_o17(xt_seri(iso_o17,i,k) & 4259 4298 & /q_seri(i,k),xt_seri(iso_o18,i,k) & … … 4619 4658 do i=1,klon 4620 4659 if ((iso_O17.gt.0).and.(iso_O18.gt.0)) then 4621 if ((q_seri(i,k).gt.ridicule).and.( l.lt.nlevmaxO17)) then4660 if ((q_seri(i,k).gt.ridicule).and.(k.lt.nlevmaxO17)) then 4622 4661 call iso_verif_aberrant_o17(xt_seri(iso_o17,i,k) & 4623 4662 & /q_seri(i,k),xt_seri(iso_o18,i,k) & … … 4841 4880 endif !if (q_seri(i,k).gt.errmax) then 4842 4881 if (ql_seri(i,k).gt.ridicule) then 4843 call iso_verif_aberrant (&4844 & xtl_seri(iso_HDO,i,k)/ql_seri(i,k),'physiq 2871')4882 call iso_verif_aberrant_choix(xtl_seri(iso_HDO,i,k),ql_seri(i,k), & 4883 ridicule,deltalim_snow,'physiq 2871') 4845 4884 if (iso_O18.gt.0) then 4846 4885 if (iso_verif_o18_aberrant_nostop( & … … 4854 4893 write(*,*) 'deltaO(d_ql_lsc(i,k))=',deltaO( & 4855 4894 & d_xtl_lsc(iso_O18,i,k)/d_ql_lsc(i,k)) 4856 stop4895 !stop 4857 4896 endif 4858 4897 endif ! if (iso_O18.gt.0) then … … 4864 4903 do i=1,klon 4865 4904 do k=1,nlev 4866 if ((q_seri(i,k).gt.ridicule).and.( l.lt.nlevmaxO17)) then4905 if ((q_seri(i,k).gt.ridicule).and.(k.lt.nlevmaxO17)) then 4867 4906 call iso_verif_aberrant_o17(xt_seri(iso_o17,i,k) & 4868 4907 & /q_seri(i,k),xt_seri(iso_o18,i,k) & … … 6149 6188 & xt_seri(iso_HDO,i,k)/q_seri(i,k), & 6150 6189 & xt_seri(iso_O18,i,k)/q_seri(i,k), & 6151 & 'physiq 5937 , juste apres methox, qv').eq.1) then6190 & 'physiq 5937a, juste apres methox, qv').eq.1) then 6152 6191 write(*,*) 'physic 2444: i,k,q_seri(i,k)=',i,k,q_seri(i,k) 6153 6192 write(*,*) 'xt_seri(:,i,k)=',xt_seri(:,i,k) … … 6159 6198 & xtl_seri(iso_HDO,i,k)/ql_seri(i,k), & 6160 6199 & xtl_seri(iso_O18,i,k)/ql_seri(i,k), & 6161 & 'physiq 5937 , juste apres methox, ql').eq.1) then6200 & 'physiq 5937b, juste apres methox, ql').eq.1) then 6162 6201 write(*,*) 'i,k,ql_seri(i,k)=',i,k,ql_seri(i,k) 6163 stop6202 !stop 6164 6203 endif ! if (iso_verif_o18_aberrant_nostop 6165 6204 endif !if (q_seri(i,k).gt.errmax) then … … 6168 6207 & xts_seri(iso_HDO,i,k)/qs_seri(i,k), & 6169 6208 & xts_seri(iso_O18,i,k)/qs_seri(i,k), & 6170 & 'physiq 5937 , juste apres methox, qs').eq.1) then6209 & 'physiq 5937c, juste apres methox, qs').eq.1) then 6171 6210 write(*,*) 'i,k,qs_seri(i,k)=',i,k,qs_seri(i,k) 6172 stop6211 !stop 6173 6212 endif ! if (iso_verif_o18_aberrant_nostop 6174 6213 endif !if (q_seri(i,k).gt.errmax) then
Note: See TracChangeset
for help on using the changeset viewer.
