Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Timestamp:
- Jul 24, 2024, 1:17:08 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmdiso
- Files:
-
- 9 edited
-
cv30_routines.F90 (modified) (11 diffs)
-
cv3_routines.F90 (modified) (1 diff)
-
cv_routines.F90 (modified) (1 diff)
-
isotopes_routines_mod.F90 (modified) (63 diffs)
-
isotopes_verif_mod.F90 (modified) (78 diffs)
-
isotrac_routines_mod.F90 (modified) (41 diffs)
-
phyaqua_mod.F90 (modified) (1 diff)
-
phyredem.F90 (modified) (1 diff)
-
physiq_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv30_routines.F90
r5112 r5113 355 355 icbmax = 2 356 356 DO i = 1, len 357 ! !icbmax=max(icbmax,icb(i))357 ! icbmax=max(icbmax,icb(i)) 358 358 IF (iflag(i)<7) icbmax = max(icbmax, icb(i)) ! sb Jun7th02 359 359 END DO … … 489 489 490 490 ! Re-compute icbsmax (icbsmax2): !convect3 491 ! !convect3491 !convect3 492 492 icbsmax2 = 2 !convect3 493 493 DO i = 1, len !convect3 … … 1691 1691 END DO 1692 1692 1693 ! !if(inb.lt.(nl-1))then1694 ! !do 85 i=inb+1,nl-11695 ! !sig(i)=beta*sig(i)+2.*alpha*buoy(inb)*1696 ! !1 abs(buoy(inb))1697 ! !sig(i)=amax1(sig(i),0.0)1698 ! !w0(i)=beta*w0(i)1699 ! !85 continue1700 ! !end if1701 1702 ! !do 87 i=1,icb1703 ! !sig(i)=beta*sig(i)-2.*alpha*buoy(icb)*buoy(icb)1704 ! !sig(i)=amax1(sig(i),0.0)1705 ! !w0(i)=beta*w0(i)1706 ! !87 continue1693 ! if(inb.lt.(nl-1))then 1694 ! do 85 i=inb+1,nl-1 1695 ! sig(i)=beta*sig(i)+2.*alpha*buoy(inb)* 1696 ! 1 abs(buoy(inb)) 1697 ! sig(i)=amax1(sig(i),0.0) 1698 ! w0(i)=beta*w0(i) 1699 ! 85 continue 1700 ! end if 1701 1702 ! do 87 i=1,icb 1703 ! sig(i)=beta*sig(i)-2.*alpha*buoy(icb)*buoy(icb) 1704 ! sig(i)=amax1(sig(i),0.0) 1705 ! w0(i)=beta*w0(i) 1706 ! 87 continue 1707 1707 1708 1708 ! ------------------------------------------------------------- … … 1787 1787 1788 1788 1789 ! !cape=0.01790 ! !do 98 i=icb+1,inb1791 ! !deltap = min(pbase,ph(i-1))-min(pbase,ph(i))1792 ! !cape=cape+rrd*buoy(i-1)*deltap/p(i-1)1793 ! !dcape=rrd*buoy(i-1)*deltap/p(i-1)1794 ! !dlnp=deltap/p(i-1)1795 ! !cape=amax1(0.0,cape)1796 ! !sigold=sig(i)1797 1798 ! !dtmin=100.01799 ! !do 97 j=icb,i-11800 ! !dtmin=amin1(dtmin,buoy(j))1801 ! !97 continue1802 1803 ! !sig(i)=beta*sig(i)+alpha*dtmin*abs(dtmin)1804 ! !sig(i)=amax1(sig(i),0.0)1805 ! !sig(i)=amin1(sig(i),0.01)1806 ! !fac=amin1(((dtcrit-dtmin)/dtcrit),1.0)1807 ! !w=(1.-beta)*fac*sqrt(cape)+beta*w0(i)1808 ! !amu=0.5*(sig(i)+sigold)*w1809 ! !m(i)=amu*0.007*p(i)*(ph(i)-ph(i+1))/tv(i)1810 ! !w0(i)=w1811 ! !98 continue1812 ! !w0(icb)=0.5*w0(icb+1)1813 ! !m(icb)=0.5*m(icb+1)*(ph(icb)-ph(icb+1))/(ph(icb+1)-ph(icb+2))1814 ! !sig(icb)=sig(icb+1)1815 ! !sig(icb-1)=sig(icb)1789 ! cape=0.0 1790 ! do 98 i=icb+1,inb 1791 ! deltap = min(pbase,ph(i-1))-min(pbase,ph(i)) 1792 ! cape=cape+rrd*buoy(i-1)*deltap/p(i-1) 1793 ! dcape=rrd*buoy(i-1)*deltap/p(i-1) 1794 ! dlnp=deltap/p(i-1) 1795 ! cape=amax1(0.0,cape) 1796 ! sigold=sig(i) 1797 1798 ! dtmin=100.0 1799 ! do 97 j=icb,i-1 1800 ! dtmin=amin1(dtmin,buoy(j)) 1801 ! 97 continue 1802 1803 ! sig(i)=beta*sig(i)+alpha*dtmin*abs(dtmin) 1804 ! sig(i)=amax1(sig(i),0.0) 1805 ! sig(i)=amin1(sig(i),0.01) 1806 ! fac=amin1(((dtcrit-dtmin)/dtcrit),1.0) 1807 ! w=(1.-beta)*fac*sqrt(cape)+beta*w0(i) 1808 ! amu=0.5*(sig(i)+sigold)*w 1809 ! m(i)=amu*0.007*p(i)*(ph(i)-ph(i+1))/tv(i) 1810 ! w0(i)=w 1811 ! 98 continue 1812 ! w0(icb)=0.5*w0(icb+1) 1813 ! m(icb)=0.5*m(icb+1)*(ph(icb)-ph(icb+1))/(ph(icb+1)-ph(icb+2)) 1814 ! sig(icb)=sig(icb+1) 1815 ! sig(icb-1)=sig(icb) 1816 1816 1817 1817 … … 2798 2798 ! enddo 2799 2799 ! enddo 2800 ! !RomP >>>2800 ! RomP >>> 2801 2801 DO i = 1, nd 2802 2802 DO il = 1, ncum … … 2805 2805 END DO 2806 2806 END DO 2807 ! !RomP <<<2807 ! RomP <<< 2808 2808 2809 2809 ! *** check whether ep(inb)=0, if so, skip precipitating *** … … 3683 3683 ! *** NE PAS UTILISER POUR L'INSTANT *** 3684 3684 3685 ! !do il=1,ncum3686 ! !wd(il)=betad*abs(mp(il,icb(il)))*0.01*rrd*t(il,icb(il))3687 ! !: /(sigd*p(il,icb(il)))3688 ! !enddo3685 ! do il=1,ncum 3686 ! wd(il)=betad*abs(mp(il,icb(il)))*0.01*rrd*t(il,icb(il)) 3687 ! : /(sigd*p(il,icb(il))) 3688 ! enddo 3689 3689 3690 3690 … … 5922 5922 ! *** diagnose the in-cloud mixing ratio *** ! cld 5923 5923 ! *** of condensed water *** ! cld 5924 ! !cld5924 ! cld 5925 5925 5926 5926 DO i = 1, nd ! cld … … 5984 5984 END SUBROUTINE cv30_yield 5985 5985 5986 ! !RomP >>>5986 !RomP >>> 5987 5987 SUBROUTINE cv30_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi, phi2, & 5988 5988 d1a, dam, ep, vprecip, elij, clw, epmlmmm, eplamm, icb, inb) … … 6042 6042 DO i = 1, ncum 6043 6043 IF (k>=icb(i) .AND. k<=inb(i) .AND. j<=inb(i)) THEN 6044 ! !jyg epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/elij(i,k,j)6044 !jyg epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/elij(i,k,j) 6045 6045 epm(i, j, k) = 1. - (1.-ep(i,j))*clw(i, j)/max(elij(i,k,j), 1.E-16) 6046 ! !6046 ! 6047 6047 epm(i, j, k) = max(epm(i,j,k), 0.0) 6048 6048 END IF … … 6353 6353 ,cape,ep,hp,icb,inb,clw,nk,t,h,lv & 6354 6354 ,epmax_diag) 6355 implicit none6355 IMPLICIT NONE 6356 6356 6357 6357 ! On fait varier epmax en fn de la cape -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90
r5112 r5113 7600 7600 , pbase, p, ph, tv, buoy, sig, w0,iflag & 7601 7601 , epmax_diag) 7602 implicit none7602 IMPLICIT NONE 7603 7603 7604 7604 ! On fait varier epmax en fn de la cape -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv_routines.F90
r5112 r5113 1631 1631 ! *** diagnose the in-cloud mixing ratio *** ! cld 1632 1632 ! *** of condensed water *** ! cld 1633 ! !cld1633 ! cld 1634 1634 DO ij = 1, ncum ! cld 1635 1635 DO i = 1, nd ! cld -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/isotopes_routines_mod.F90
r5112 r5113 29 29 & uncompress_ilp_zone 30 30 #endif 31 implicit none 31 IMPLICIT NONE 32 32 33 33 ! inputs … … 866 866 USE isotopes_verif_mod 867 867 #endif 868 implicit none 868 IMPLICIT NONE 869 869 870 870 ! calcul de Rsol et Revap lors de l'évaporation de l'eau du sol … … 1315 1315 SUBROUTINE calcul_kcin(Vsurf,KCIN) 1316 1316 USE isotopes_mod, ONLY: tv0cin,tkcin0,tkcin1,tkcin2 1317 implicit none 1317 IMPLICIT NONE 1318 1318 ! calcul de kcin en fonction de Vsurf 1319 1319 … … 1348 1348 USE isotopes_verif_mod 1349 1349 #endif 1350 implicit none1350 IMPLICIT NONE 1351 1351 1352 1352 !C------------------------------------------------------------------------- … … 1475 1475 USE isotopes_verif_mod 1476 1476 #endif 1477 implicit none1477 IMPLICIT NONE 1478 1478 1479 1479 !C------------------------------------------------------------------------- … … 1538 1538 USE isotopes_verif_mod 1539 1539 #endif 1540 implicit none1540 IMPLICIT NONE 1541 1541 1542 1542 !C------------------------------------------------------------------------- … … 1646 1646 USE isotopes_verif_mod 1647 1647 #endif 1648 implicit none1648 IMPLICIT NONE 1649 1649 1650 1650 !C------------------------------------------------------------------------- … … 1817 1817 USE isotopes_verif_mod 1818 1818 #endif 1819 implicit none1819 IMPLICIT NONE 1820 1820 1821 1821 !C------------------------------------------------------------------------- … … 1895 1895 USE isotopes_verif_mod 1896 1896 #endif 1897 implicit none1897 IMPLICIT NONE 1898 1898 1899 1899 !C------------------------------------------------------------------------- … … 2034 2034 USE isotopes_verif_mod 2035 2035 #endif 2036 implicit none 2036 IMPLICIT NONE 2037 2037 2038 2038 ! calcul de Rsol … … 2246 2246 USE isotrac_mod, ONLY: index_iso 2247 2247 #endif 2248 implicit none2248 IMPLICIT NONE 2249 2249 ! inputs 2250 2250 integer i … … 2454 2454 USE isotopes_verif_mod 2455 2455 #endif 2456 implicit none 2456 IMPLICIT NONE 2457 2457 2458 2458 ! version véctorisée: sur les isotopes et sur les points de … … 4411 4411 USE isotopes_verif_mod 4412 4412 #endif 4413 implicit none 4413 IMPLICIT NONE 4414 4414 4415 4415 ! idem que stewart_loop, mais pour la rrévap de la glace. … … 4568 4568 USE isotopes_verif_mod 4569 4569 #endif 4570 implicit none 4570 IMPLICIT NONE 4571 4571 4572 4572 ! rrévap de la glace. … … 4809 4809 USE isotopes_verif_mod 4810 4810 #endif 4811 implicit none4811 IMPLICIT NONE 4812 4812 4813 4813 ! compression des variables en tableaux spécifiques pour le … … 4967 4967 USE isotopes_verif_mod 4968 4968 #endif 4969 implicit none4969 IMPLICIT NONE 4970 4970 4971 4971 ! des compressions des cas de calcul de J dans stewart_explicit … … 5012 5012 USE isotopes_verif_mod 5013 5013 #endif 5014 implicit none5014 IMPLICIT NONE 5015 5015 5016 5016 ! decompression des outputs communs à tous les cas dans … … 5076 5076 USE isotopes_verif_mod 5077 5077 #endif 5078 implicit none5078 IMPLICIT NONE 5079 5079 5080 5080 ! compression dans le cas condensation_facftmr … … 5151 5151 USE isotopes_verif_mod 5152 5152 #endif 5153 implicit none5153 IMPLICIT NONE 5154 5154 5155 5155 ! compression dans le cas condensation_facftmr … … 5224 5224 USE isotopes_verif_mod 5225 5225 #endif 5226 implicit none5226 IMPLICIT NONE 5227 5227 5228 5228 ! compression dans le cas condensation_facftmr … … 5293 5293 USE isotopes_verif_mod 5294 5294 #endif 5295 implicit none5295 IMPLICIT NONE 5296 5296 5297 5297 ! compression dans le cas condensation_facftmr … … 5465 5465 #endif 5466 5466 #endif 5467 implicit none5467 IMPLICIT NONE 5468 5468 5469 5469 ! compression dans le cas condensation_facftmr … … 5594 5594 USE isotopes_mod, ONLY: iso_eau, iso_HDO 5595 5595 5596 implicit none5596 IMPLICIT NONE 5597 5597 5598 5598 ! inputs … … 5628 5628 5629 5629 USE isotopes_mod, ONLY: iso_eau, iso_HDO 5630 implicit none5630 IMPLICIT NONE 5631 5631 5632 5632 integer ncas … … 5665 5665 5666 5666 USE isotopes_mod, ONLY: iso_eau, iso_HDO 5667 implicit none5667 IMPLICIT NONE 5668 5668 5669 5669 integer ncas … … 5719 5719 USE isotopes_verif_mod 5720 5720 #endif 5721 implicit none5721 IMPLICIT NONE 5722 5722 5723 5723 ! inputs … … 5794 5794 USE isotopes_verif_mod 5795 5795 #endif 5796 implicit none5796 IMPLICIT NONE 5797 5797 5798 5798 ! version vectorisée en ncas. La vesrion _bak23fev2008 était … … 5817 5817 5818 5818 ! **locals 5819 ! !nombre max d'itération dans integrale rieman5819 ! nombre max d'itération dans integrale rieman 5820 5820 integer j 5821 5821 integer il,ixt,k … … 6002 6002 #endif 6003 6003 #endif 6004 implicit none6004 IMPLICIT NONE 6005 6005 6006 6006 !*inputs et outputs … … 6503 6503 #endif 6504 6504 endif !if (T(il,i).ge.ztglace_kelvin) then 6505 endif ! !if (Eqi_prime.lt.0.0) then6505 endif !if (Eqi_prime.lt.0.0) then 6506 6506 endif !if (i.le.inb(il) .and. lwork(il)) then 6507 6507 enddo !do il=1,ncum … … 7567 7567 USE isotopes_verif_mod 7568 7568 #endif 7569 implicit none7569 IMPLICIT NONE 7570 7570 7571 7571 ! inputs … … 7934 7934 USE isotopes_verif_mod 7935 7935 #endif 7936 implicit none7936 IMPLICIT NONE 7937 7937 7938 7938 ! inputs … … 8138 8138 USE isotopes_verif_mod 8139 8139 #endif 8140 implicit none8140 IMPLICIT NONE 8141 8141 8142 8142 ! inputs … … 8276 8276 USE isotrac_mod, only: ridicule_trac 8277 8277 #endif 8278 implicit none8278 IMPLICIT NONE 8279 8279 8280 8280 ! inputs … … 8810 8810 USE isotopes_verif_mod 8811 8811 #endif 8812 implicit none8812 IMPLICIT NONE 8813 8813 8814 8814 ! inputs … … 9211 9211 #endif 9212 9212 #endif 9213 implicit none9213 IMPLICIT NONE 9214 9214 9215 9215 !*inputs et outputs … … 9786 9786 !#endif 9787 9787 endif !if (T(il,i).ge.ztglace_kelvin) then 9788 endif ! !if (Eqi_prime.lt.0.0) then9788 endif !if (Eqi_prime.lt.0.0) then 9789 9789 endif !if (i.le.inb(il) .and. lwork(il)) then 9790 9790 enddo !do il=1,ncum … … 10902 10902 USE isotopes_verif_mod 10903 10903 #endif 10904 implicit none10904 IMPLICIT NONE 10905 10905 10906 10906 ! inputs … … 11054 11054 #endif 11055 11055 #endif 11056 implicit none11056 IMPLICIT NONE 11057 11057 11058 11058 ! inputs … … 11648 11648 USE isotopes_verif_mod 11649 11649 #endif 11650 implicit none11650 IMPLICIT NONE 11651 11651 11652 11652 ! inputs … … 12076 12076 USE isotopes_verif_mod 12077 12077 #endif 12078 implicit none12078 IMPLICIT NONE 12079 12079 12080 12080 ! version vectorisée de condiso_liq_ice … … 12313 12313 #endif 12314 12314 #endif 12315 implicit none12315 IMPLICIT NONE 12316 12316 12317 12317 ! version vectorisée de condiso_liq_ice … … 12768 12768 USE isotopes_verif_mod 12769 12769 #endif 12770 implicit none12770 IMPLICIT NONE 12771 12771 12772 12772 ! on s'interresse à l'isotope ixt. … … 12940 12940 12941 12941 USE isotopes_mod, ONLY: pxtmelt,pxtice 12942 implicit none12942 IMPLICIT NONE 12943 12943 12944 12944 ! inputs … … 12971 12971 USE isotrac_mod, ONLY: izone_cont,index_zone,index_iso 12972 12972 #endif 12973 implicit none12973 IMPLICIT NONE 12974 12974 12975 12975 ! gestion de la neige: on precipte dessus, sublime, effondre, … … 13541 13541 & bassin_map 13542 13542 #endif 13543 implicit none13543 IMPLICIT NONE 13544 13544 13545 13545 include "YOETHF.h" … … 13864 13864 & bassin_map 13865 13865 #endif 13866 implicit none13866 IMPLICIT NONE 13867 13867 13868 13868 ! inputs … … 14116 14116 & bassin_map 14117 14117 #endif 14118 implicit none14118 IMPLICIT NONE 14119 14119 14120 14120 ! inputs … … 14456 14456 #endif 14457 14457 14458 implicit none14458 IMPLICIT NONE 14459 14459 14460 14460 include "YOETHF.h" … … 15867 15867 USE isotopes_verif_mod 15868 15868 #endif 15869 implicit none15869 IMPLICIT NONE 15870 15870 15871 15871 ! equivalent de phyetat0 pour les isotopes … … 16044 16044 USE isotrac_mod, ONLY: index_iso,index_zone,izone_init 16045 16045 #endif 16046 implicit none16046 IMPLICIT NONE 16047 16047 16048 16048 include "dimsoil.h" … … 16633 16633 USE isotopes_verif_mod 16634 16634 #endif 16635 implicit none16635 IMPLICIT NONE 16636 16636 16637 16637 ! input … … 17952 17952 USE isotopes_verif_mod 17953 17953 #endif 17954 implicit none17954 IMPLICIT NONE 17955 17955 17956 17956 ! Arguments … … 18080 18080 use phys_cal_mod ! pour le calendrier 18081 18081 use isotopes_mod, only: nessai, day_nucl, month_nucl, year_nucl 18082 implicit none18082 IMPLICIT NONE 18083 18083 18084 18084 ! Arguments … … 18152 18152 USE isotopes_verif_mod 18153 18153 #endif 18154 implicit none18154 IMPLICIT NONE 18155 18155 18156 18156 ! Arguments … … 18230 18230 USE isotopes_verif_mod 18231 18231 #endif 18232 implicit none18232 IMPLICIT NONE 18233 18233 18234 18234 include "dimensions.h" ! pour avoir iim et jjm … … 18372 18372 USE isotopes_verif_mod 18373 18373 #endif 18374 implicit none18374 IMPLICIT NONE 18375 18375 18376 18376 ! Arguments … … 18516 18516 USE isotopes_verif_mod 18517 18517 #endif 18518 implicit none18518 IMPLICIT NONE 18519 18519 18520 18520 ! version vectorisée de condiso_liq_ice … … 18646 18646 USE isotopes_verif_mod 18647 18647 #endif 18648 implicit none18648 IMPLICIT NONE 18649 18649 18650 18650 ! version vectorisée de condiso_liq_ice … … 18780 18780 USE isotopes_verif_mod 18781 18781 #endif 18782 implicit none18782 IMPLICIT NONE 18783 18783 18784 18784 ! inputs … … 18840 18840 USE isotopes_verif_mod 18841 18841 #endif 18842 implicit none 18842 IMPLICIT NONE 18843 18843 18844 18844 … … 18892 18892 18893 18893 use infotrac_phy, ONLY: nqtot,nqo,ivap,iliq,isol,iqIsoPha,ntraciso=>ntiso 18894 implicit none 18894 IMPLICIT NONE 18895 18895 18896 18896 ! inputs … … 18932 18932 18933 18933 use infotrac_phy, ONLY: nqtot,nqo,ivap,iliq,isol,iqIsoPha,ntraciso=>ntiso 18934 implicit none 18934 IMPLICIT NONE 18935 18935 18936 18936 ! inputs -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/isotopes_verif_mod.F90
r5112 r5113 9 9 !#endif 10 10 USE infotrac_phy, ONLY: ntraciso=>ntiso, niso, itZonIso, nzone 11 implicit none 11 IMPLICIT NONE 12 12 save 13 13 … … 95 95 use lmdz_ioipsl_getin_p, ONLY: getin_p 96 96 use isotopes_mod, ONLY: iso_O17, iso_O18, iso_HDO 97 implicit none97 IMPLICIT NONE 98 98 99 99 write(*,*) 'iso_verif_init 99: entree' … … 137 137 138 138 SUBROUTINE iso_verif_egalite(a,b,err_msg) 139 implicit none139 IMPLICIT NONE 140 140 ! compare a et b. Si pas egal à errmax près, on affiche message 141 141 ! d'erreur et stoppe … … 166 166 167 167 function iso_verif_egalite_nostop(a,b,err_msg) 168 implicit none168 IMPLICIT NONE 169 169 ! compare a et b. Si pas egal à errmax près, on affiche message 170 170 ! d'erreur et stoppe … … 197 197 SUBROUTINE iso_verif_aberrant(R,err_msg) 198 198 use isotopes_mod, ONLY: ridicule, iso_HDO 199 implicit none199 IMPLICIT NONE 200 200 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 201 201 ! d'erreur et stoppe … … 227 227 SUBROUTINE iso_verif_aberrant_encadre(R,err_msg) 228 228 use isotopes_mod, ONLY: ridicule, iso_HDO 229 implicit none229 IMPLICIT NONE 230 230 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 231 231 ! d'erreur et stoppe … … 262 262 SUBROUTINE iso_verif_aberrant_choix(xt,q,qmin,deltaDmax,err_msg) 263 263 use isotopes_mod, ONLY: iso_HDO 264 implicit none264 IMPLICIT NONE 265 265 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 266 266 ! d'erreur et stoppe … … 296 296 function iso_verif_aberrant_nostop(R,err_msg) 297 297 use isotopes_mod, ONLY: ridicule,iso_HDO 298 implicit none298 IMPLICIT NONE 299 299 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 300 300 ! d'erreur et stoppe … … 327 327 function iso_verif_aberrant_enc_nostop(R,err_msg) 328 328 use isotopes_mod, ONLY: ridicule,iso_HDO 329 implicit none329 IMPLICIT NONE 330 330 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 331 331 ! d'erreur et stoppe … … 362 362 363 363 use isotopes_mod, ONLY: iso_HDO 364 implicit none364 IMPLICIT NONE 365 365 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 366 366 ! d'erreur et stoppe … … 423 423 qmin,deltaDmax,err_msg) 424 424 use isotopes_mod, ONLY: iso_HDO 425 implicit none425 IMPLICIT NONE 426 426 ! si le rapprot iso R est plus grand que deltaDlim, on affiche message 427 427 ! d'erreur et stoppe … … 467 467 468 468 SUBROUTINE iso_verif_aberrant_o17(R17,R18,err_msg) 469 implicit none469 IMPLICIT NONE 470 470 ! si l'O17-excess est aberrant, on affiche un message 471 471 ! d'erreur et stoppe … … 500 500 function iso_verif_aberrant_o17_nostop(R17,R18,err_msg) 501 501 USE isotopes_mod, ONLY: tnat,iso_O17,iso_O18 502 implicit none502 IMPLICIT NONE 503 503 ! si l'O17-excess est aberrant, on affiche un message 504 504 ! d'erreur et renvoit 1 … … 541 541 542 542 SUBROUTINE iso_verif_noNaN(x,err_msg) 543 implicit none543 IMPLICIT NONE 544 544 ! si x est NaN, on affiche message 545 545 ! d'erreur et stoppe … … 566 566 567 567 function iso_verif_noNaN_nostop(x,err_msg) 568 implicit none568 IMPLICIT NONE 569 569 ! si x est NaN, on affiche message 570 570 ! d'erreur et return 1 si erreur … … 596 596 597 597 SUBROUTINE iso_verif_noNaN_vect2D(x,err_msg,ni,n,m) 598 implicit none598 IMPLICIT NONE 599 599 ! si x est NaN, on affiche message 600 600 ! d'erreur et return 1 si erreur … … 635 635 636 636 SUBROUTINE iso_verif_noNaN_vect1D(x,err_msg,ni,n) 637 implicit none637 IMPLICIT NONE 638 638 ! si x est NaN, on affiche message 639 639 ! d'erreur et return 1 si erreur … … 674 674 !************************ 675 675 SUBROUTINE iso_verif_egalite_choix(a,b,err_msg,erabs,errel) 676 implicit none676 IMPLICIT NONE 677 677 ! compare a et b. Si pas egal, on affiche message 678 678 ! d'erreur et stoppe … … 703 703 function iso_verif_egalite_choix_nostop & 704 704 (a,b,err_msg,erabs,errel) 705 implicit none705 IMPLICIT NONE 706 706 ! compare a et b. Si pas egal, on affiche message 707 707 ! d'erreur et stoppe … … 758 758 SUBROUTINE iso_verif_positif(x,err_msg) 759 759 use isotopes_mod, ONLY: ridicule 760 implicit none760 IMPLICIT NONE 761 761 ! si x<0, on plante. 762 762 ! si très limite, on le met à 0. … … 786 786 SUBROUTINE iso_verif_positif_vect(x,n,err_msg) 787 787 use isotopes_mod, ONLY: ridicule 788 implicit none788 IMPLICIT NONE 789 789 ! si x<0, on plante. 790 790 … … 816 816 817 817 SUBROUTINE iso_verif_positif_choix_vect(x,n,ridic,err_msg) 818 implicit none818 IMPLICIT NONE 819 819 ! si x<0, on plante. 820 820 … … 848 848 !****************** 849 849 SUBROUTINE iso_verif_positif_strict(x,err_msg) 850 implicit none850 IMPLICIT NONE 851 851 ! si x<0, on plante. 852 852 ! si très limite, on le met à 0. … … 869 869 870 870 function iso_verif_positif_strict_nostop(x,err_msg) 871 implicit none871 IMPLICIT NONE 872 872 ! si x<0, on plante. 873 873 ! si très limite, on le met à 0. … … 896 896 897 897 SUBROUTINE iso_verif_positif_choix(x,ridic,err_msg) 898 implicit none898 IMPLICIT NONE 899 899 ! si x<0, on plante. 900 900 ! si très limite, on le met à 0. … … 925 925 function iso_verif_positif_nostop(x,err_msg) 926 926 use isotopes_mod, ONLY: ridicule 927 implicit none927 IMPLICIT NONE 928 928 ! si x<0, on plante. 929 929 ! si très limite, on le met à 0. … … 955 955 !****************** 956 956 function iso_verif_positif_choix_nostop(x,ridic,err_msg) 957 implicit none957 IMPLICIT NONE 958 958 ! si x<0, on plante. 959 959 ! si très limite, on le met à 0. … … 992 992 993 993 SUBROUTINE iso_verif_O18_aberrant(Rd,Ro,err_msg) 994 implicit none994 IMPLICIT NONE 995 995 996 996 ! vérifie que: … … 1012 1012 function iso_verif_O18_aberrant_nostop(Rd,Ro,err_msg) 1013 1013 USE isotopes_mod, ONLY: tnat, iso_HDO, iso_O18 1014 implicit none1014 IMPLICIT NONE 1015 1015 1016 1016 ! vérifie que: … … 1060 1060 function deltaD(R) 1061 1061 USE isotopes_mod, ONLY: tnat,iso_HDO 1062 implicit none1062 IMPLICIT NONE 1063 1063 real R,deltaD 1064 1064 … … 1076 1076 function deltaO(R) 1077 1077 USE isotopes_mod, ONLY: tnat,iso_O18 1078 implicit none1078 IMPLICIT NONE 1079 1079 real R,deltaO 1080 1080 … … 1091 1091 function dexcess(RD,RO) 1092 1092 USE isotopes_mod, ONLY: tnat,iso_O18,iso_HDO 1093 implicit none1093 IMPLICIT NONE 1094 1094 real RD,RO,deltaD,deltaO,dexcess 1095 1095 … … 1108 1108 function delta_all(R,ixt) 1109 1109 USE isotopes_mod, ONLY: tnat 1110 implicit none1110 IMPLICIT NONE 1111 1111 real R,delta_all 1112 1112 integer ixt … … 1119 1119 function delta_to_R(delta,ixt) 1120 1120 USE isotopes_mod, ONLY: tnat 1121 implicit none1121 IMPLICIT NONE 1122 1122 real delta,delta_to_R 1123 1123 integer ixt … … 1130 1130 function o17excess(R17,R18) 1131 1131 USE isotopes_mod, ONLY: tnat,iso_O18,iso_O17 1132 implicit none1132 IMPLICIT NONE 1133 1133 real R17,R18,o17excess 1134 1134 … … 1151 1151 1152 1152 USE isotopes_mod, ONLY: iso_eau 1153 implicit none1153 IMPLICIT NONE 1154 1154 1155 1155 ! inputs … … 1202 1202 1203 1203 USE isotopes_mod, ONLY: iso_eau 1204 implicit none1204 IMPLICIT NONE 1205 1205 1206 1206 ! inputs … … 1242 1242 a,b,err_msg,n,m,errmax,errmaxrel) 1243 1243 1244 implicit none1244 IMPLICIT NONE 1245 1245 1246 1246 ! inputs … … 1285 1285 xt,q,err_msg,ni,n,m) 1286 1286 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO 1287 implicit none1287 IMPLICIT NONE 1288 1288 1289 1289 ! inputs … … 1333 1333 1334 1334 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO 1335 implicit none1335 IMPLICIT NONE 1336 1336 1337 1337 ! inputs … … 1386 1386 1387 1387 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO 1388 implicit none1388 IMPLICIT NONE 1389 1389 1390 1390 ! inputs … … 1436 1436 1437 1437 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO 1438 implicit none1438 IMPLICIT NONE 1439 1439 1440 1440 … … 1486 1486 1487 1487 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO,iso_O18 1488 implicit none1488 IMPLICIT NONE 1489 1489 1490 1490 ! inputs … … 1538 1538 SUBROUTINE select_dim23_from4D(n1,n2,n3,n4, & 1539 1539 var,vec,i1,i4) 1540 implicit none1540 IMPLICIT NONE 1541 1541 1542 1542 ! inputs … … 1561 1561 SUBROUTINE select_dim4_from4D(ntime,nlev,nlat,nlon, & 1562 1562 var,vec,itime,ilev,ilat) 1563 implicit none1563 IMPLICIT NONE 1564 1564 1565 1565 ! inputs … … 1581 1581 SUBROUTINE select_dim5_from5D(n1,n2,n3,n4,n5, & 1582 1582 var,vec,i1,i2,i3,i4) 1583 implicit none1583 IMPLICIT NONE 1584 1584 1585 1585 ! inputs … … 1601 1601 SUBROUTINE select_dim3_from3D(ntime,nlat,nlon, & 1602 1602 var,vec,itime,ilat) 1603 implicit none1603 IMPLICIT NONE 1604 1604 1605 1605 ! inputs … … 1621 1621 SUBROUTINE select_dim23_from3D(n1,n2,n3, & 1622 1622 var,vec,i1) 1623 implicit none1623 IMPLICIT NONE 1624 1624 1625 1625 ! inputs … … 1642 1642 SUBROUTINE putinto_dim23_from4D(n1,n2,n3,n4, & 1643 1643 var,vec,i1,i4) 1644 implicit none1644 IMPLICIT NONE 1645 1645 1646 1646 ! inputs … … 1664 1664 SUBROUTINE putinto_dim12_from4D(n1,n2,n3,n4, & 1665 1665 var,vec,i3,i4) 1666 implicit none1666 IMPLICIT NONE 1667 1667 1668 1668 ! inputs … … 1685 1685 SUBROUTINE putinto_dim23_from3D(n1,n2,n3, & 1686 1686 var,vec,i1) 1687 implicit none1687 IMPLICIT NONE 1688 1688 1689 1689 ! inputs … … 1707 1707 1708 1708 SUBROUTINE iso_verif_noNaN_par2D(x,err_msg,ni,n,m,ib,ie) 1709 implicit none1709 IMPLICIT NONE 1710 1710 ! si x est NaN, on affiche message 1711 1711 ! d'erreur et return 1 si erreur … … 1750 1750 1751 1751 use isotopes_mod, ONLY: ridicule,tnat,iso_HDO 1752 implicit none1752 IMPLICIT NONE 1753 1753 1754 1754 ! inputs … … 1800 1800 1801 1801 USE isotopes_mod, ONLY: iso_eau 1802 implicit none1802 IMPLICIT NONE 1803 1803 1804 1804 ! inputs … … 1845 1845 errmax,errmaxrel,ridicule_trac,deltalimtrac) 1846 1846 use isotopes_mod, ONLY: iso_HDO 1847 implicit none1847 IMPLICIT NONE 1848 1848 ! vérifier des choses sur les traceurs 1849 1849 ! * toutes les zones donne t l'istope total … … 1896 1896 errmax,errmaxrel,ridicule_trac,deltalimtrac) 1897 1897 USE isotopes_mod, ONLY: iso_HDO 1898 implicit none1898 IMPLICIT NONE 1899 1899 ! vérifier des choses sur les traceurs 1900 1900 ! * toutes les zones donne t l'istope total … … 1941 1941 function iso_verif_tracpos_choix_nostop(x,err_msg,seuil) 1942 1942 use isotopes_mod, only: isoName 1943 implicit none1943 IMPLICIT NONE 1944 1944 1945 1945 ! inputs … … 1970 1970 function iso_verif_traceur_noNaN_nostop(x,err_msg) 1971 1971 use isotopes_mod, only: isoName 1972 implicit none1972 IMPLICIT NONE 1973 1973 1974 1974 ! on vérifie juste que pas NaN … … 2003 2003 use isotopes_mod, ONLY: ridicule,isoName 2004 2004 ! on vérifie juste bilan de masse 2005 implicit none2005 IMPLICIT NONE 2006 2006 2007 2007 ! inputs … … 2060 2060 use isotrac_mod, only: strtrac 2061 2061 ! on vérifie juste deltaD 2062 implicit none2062 IMPLICIT NONE 2063 2063 2064 2064 ! inputs … … 2140 2140 SUBROUTINE iso_verif_traceur(x,err_msg) 2141 2141 use isotrac_mod, only: ridicule_trac 2142 implicit none2142 IMPLICIT NONE 2143 2143 ! vérifier des choses sur les traceurs 2144 2144 ! * toutes les zones donne t l'istope total … … 2168 2168 use isotrac_mod, only: ridicule_trac 2169 2169 2170 implicit none2170 IMPLICIT NONE 2171 2171 ! vérifier des choses sur les traceurs 2172 2172 ! * toutes les zones donne t l'istope total … … 2200 2200 use isotrac_mod, only: ridicule_trac 2201 2201 2202 implicit none2202 IMPLICIT NONE 2203 2203 ! vérifier des choses sur les traceurs 2204 2204 ! * toutes les zones donne t l'istope total … … 2232 2232 i1,i2,err_msg) 2233 2233 use isotrac_mod, only: ridicule_trac 2234 implicit none2234 IMPLICIT NONE 2235 2235 ! vérifier des choses sur les traceurs 2236 2236 ! * toutes les zones donne t l'istope total … … 2262 2262 SUBROUTINE iso_verif_traceur_vect(x,n,m,err_msg) 2263 2263 USE isotopes_mod, ONLY: iso_HDO 2264 implicit none2264 IMPLICIT NONE 2265 2265 ! vérifier des choses sur les traceurs 2266 2266 ! * toutes les zones donne t l'istope total … … 2297 2297 SUBROUTINE iso_verif_tracnps_vect(x,n,m,err_msg) 2298 2298 USE isotopes_mod, ONLY: iso_HDO 2299 implicit none2299 IMPLICIT NONE 2300 2300 ! vérifier des choses sur les traceurs 2301 2301 ! * toutes les zones donne t l'istope total … … 2329 2329 2330 2330 SUBROUTINE iso_verif_traceur_noNaN_vect(x,n,m,err_msg) 2331 implicit none2331 IMPLICIT NONE 2332 2332 2333 2333 ! inputs … … 2373 2373 errmax,errmaxrel) 2374 2374 use isotopes_mod, only: isoName 2375 implicit none2375 IMPLICIT NONE 2376 2376 2377 2377 ! inputs … … 2413 2413 use isotopes_mod, only: iso_HDO,iso_eau 2414 2414 use isotrac_mod, only: strtrac 2415 implicit none2415 IMPLICIT NONE 2416 2416 2417 2417 ! inputs … … 2453 2453 2454 2454 SUBROUTINE iso_verif_tracpos_vect(x,n,m,err_msg,seuil) 2455 implicit none2455 IMPLICIT NONE 2456 2456 2457 2457 ! inputs … … 2496 2496 use isotrac_mod, only: ridicule_trac 2497 2497 2498 implicit none2498 IMPLICIT NONE 2499 2499 ! vérifier des choses sur les traceurs 2500 2500 ! * toutes les zones donne t l'istope total … … 2520 2520 2521 2521 SUBROUTINE iso_verif_tracpos_choix(x,err_msg,seuil) 2522 implicit none2522 IMPLICIT NONE 2523 2523 ! vérifier des choses sur les traceurs 2524 2524 ! * toutes les zones donne t l'istope total … … 2545 2545 SUBROUTINE iso_verif_traceur_choix(x,err_msg, & 2546 2546 errmax,errmaxrel,ridicule_trac_loc,deltalimtrac) 2547 implicit none2547 IMPLICIT NONE 2548 2548 ! vérifier des choses sur les traceurs 2549 2549 ! * toutes les zones donne t l'istope total … … 2569 2569 use isotrac_mod, only: ridicule_trac 2570 2570 !use isotopes_verif, only: errmax,errmaxrel,deltalimtrac 2571 implicit none2571 IMPLICIT NONE 2572 2572 ! vérifier des choses sur les traceurs 2573 2573 ! * toutes les zones donne t l'istope total … … 2595 2595 2596 2596 SUBROUTINE iso_verif_traceur_justmass(x,err_msg) 2597 implicit none2597 IMPLICIT NONE 2598 2598 ! on vérifie que noNaN et masse 2599 2599 … … 2623 2623 2624 2624 function iso_verif_traceur_jm_nostop(x,err_msg) 2625 implicit none2625 IMPLICIT NONE 2626 2626 ! on vérifie que noNaN et masse 2627 2627 … … 2657 2657 USE isotopes_mod, ONLY: tnat,iso_eau, ridicule,iso_HDO 2658 2658 use isotrac_mod, only: option_traceurs,nzone_temp 2659 implicit none2659 IMPLICIT NONE 2660 2660 2661 2661 ! inputs … … 2716 2716 USE isotopes_mod, ONLY: tnat,iso_eau,iso_HDO,ridicule 2717 2717 use isotrac_mod, only: option_traceurs,nzone_temp 2718 implicit none2718 IMPLICIT NONE 2719 2719 2720 2720 ! inputs -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/isotrac_routines_mod.F90
r5112 r5113 20 20 USE isotopes_mod, ONLY: ridicule,iso_eau 21 21 22 implicit none22 IMPLICIT NONE 23 23 24 24 ! decompression des outputs communs à tous les cas dans … … 70 70 #endif 71 71 72 implicit none72 IMPLICIT NONE 73 73 74 74 ! decompression des outputs communs à tous les cas dans … … 235 235 #endif 236 236 237 implicit none237 IMPLICIT NONE 238 238 239 239 ! compression dans le cas condensation_facftmr … … 323 323 #endif 324 324 325 implicit none325 IMPLICIT NONE 326 326 327 327 ! compression dans le cas condensation_facftmr … … 409 409 USE isotopes_verif_mod 410 410 #endif 411 implicit none411 IMPLICIT NONE 412 412 413 413 ! compression dans le cas condensation_facftmr … … 476 476 USE isotopes_verif_mod 477 477 #endif 478 implicit none478 IMPLICIT NONE 479 479 480 480 ! compression dans le cas condensation_facftmr … … 644 644 USE isotopes_verif_mod 645 645 #endif 646 implicit none646 IMPLICIT NONE 647 647 648 648 ! compression dans le cas condensation_facftmr … … 804 804 #endif 805 805 806 implicit none806 IMPLICIT NONE 807 807 808 808 ! inputs … … 891 891 USE isotopes_verif_mod 892 892 #endif 893 implicit none893 IMPLICIT NONE 894 894 895 895 ! inputs … … 978 978 #endif 979 979 980 implicit none980 IMPLICIT NONE 981 981 982 982 ! inputs … … 1042 1042 #endif 1043 1043 USE isotrac_mod, only: izone_revap 1044 implicit none1044 IMPLICIT NONE 1045 1045 1046 1046 ! ajoute xtrevap_tag (evaps des différents traceurs d'isotopes) … … 1114 1114 & bassin_MerArabie,bassin_BengalGolf,bassin_SouthIndian, & 1115 1115 & bassin_tropics,bassin_midlats,bassin_HighLats 1116 implicit none1116 IMPLICIT NONE 1117 1117 ! répond true si lat,lon se trouve dans le bassin numéroté bassin 1118 1118 … … 1320 1320 use isotopes_verif_mod 1321 1321 #endif 1322 implicit none1322 IMPLICIT NONE 1323 1323 1324 1324 ! inputs … … 1378 1378 use isotopes_verif_mod 1379 1379 #endif 1380 implicit none1380 IMPLICIT NONE 1381 1381 real presnivs(klev) 1382 1382 … … 1403 1403 #endif 1404 1404 1405 implicit none1405 IMPLICIT NONE 1406 1406 1407 1407 ! inputs … … 1424 1424 1425 1425 function is_in_rectangle(x,y,x1,y1,x2,y2) 1426 implicit none1426 IMPLICIT NONE 1427 1427 ! inputs 1428 1428 real x,y … … 1451 1451 1452 1452 function is_in_triangle(x,y,x1,y1,x2,y2,x3,y3) 1453 implicit none1453 IMPLICIT NONE 1454 1454 ! inputs 1455 1455 real x,y … … 1494 1494 USE isotopes_verif_mod 1495 1495 #endif 1496 implicit none1496 IMPLICIT NONE 1497 1497 1498 1498 ! inout … … 1578 1578 USE isotopes_verif_mod 1579 1579 #endif 1580 implicit none1580 IMPLICIT NONE 1581 1581 1582 1582 ! recolorise selon la température minimum, mais les tags de … … 1633 1633 USE isotopes_verif_mod 1634 1634 #endif 1635 implicit none1635 IMPLICIT NONE 1636 1636 1637 1637 ! recolorise selon la température minimum, mais les tags de … … 1697 1697 USE isotopes_verif_mod 1698 1698 #endif 1699 implicit none1699 IMPLICIT NONE 1700 1700 1701 1701 ! SUBROUTINE écrite à la base pour tagguer 3 boites AMMA. … … 1750 1750 USE isotopes_verif_mod 1751 1751 #endif 1752 implicit none1752 IMPLICIT NONE 1753 1753 1754 1754 ! SUBROUTINE écrite pour l'option de taggage 20 … … 1798 1798 USE isotopes_verif_mod 1799 1799 #endif 1800 implicit none1800 IMPLICIT NONE 1801 1801 1802 1802 ! SUBROUTINE écrite pour l'option de taggage 20 … … 1867 1867 #endif 1868 1868 USE isotrac_mod, only: izone_aej,izone_mousson,izone_harmattan 1869 implicit none1869 IMPLICIT NONE 1870 1870 1871 1871 real lat(klon),lon(klon) … … 1920 1920 USE isotopes_verif_mod 1921 1921 #endif 1922 implicit none1922 IMPLICIT NONE 1923 1923 1924 1924 real lat(klon),lon(klon) … … 1954 1954 function index_zone_lat(lat) 1955 1955 use isotrac_mod, only: lattag_min,dlattag,nzone_lat 1956 implicit none1956 IMPLICIT NONE 1957 1957 1958 1958 ! inputs … … 1974 1974 function index_zone_pres(pres) 1975 1975 use isotrac_mod, only: nzone_pres,zone_pres 1976 implicit none1976 IMPLICIT NONE 1977 1977 1978 1978 ! inputs … … 1992 1992 1993 1993 function find_index(pres,nzone_pres,zone_pres) 1994 implicit none1994 IMPLICIT NONE 1995 1995 1996 1996 ! inputs … … 2029 2029 function index_zone_latpres(lat,pres) 2030 2030 use isotrac_mod, only: nzone_lat 2031 implicit none2031 IMPLICIT NONE 2032 2032 2033 2033 ! inputs … … 2058 2058 USE isotopes_verif_mod 2059 2059 #endif 2060 implicit none2060 IMPLICIT NONE 2061 2061 2062 2062 ! on recolorise la vapeur résiduelle selon la température de condensation … … 2263 2263 USE isotopes_verif_mod 2264 2264 #endif 2265 implicit none2265 IMPLICIT NONE 2266 2266 2267 2267 ! inputs … … 2288 2288 USE dimphy, only: klon,klev 2289 2289 use isotrac_mod, only: option_traceurs,boite_map 2290 implicit none2290 IMPLICIT NONE 2291 2291 2292 2292 ! inputs … … 2319 2319 !use isotrac_mod, only: ntraceurs_zone=>nzone 2320 2320 USE infotrac_phy, ONLY: ntraceurs_zone=>nzone 2321 implicit none2321 IMPLICIT NONE 2322 2322 2323 2323 ! version vectrisée de iso_verif_traceur_jbidouille … … 2382 2382 SUBROUTINE iso_verif_traceur_jbidouille(x) 2383 2383 USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule 2384 implicit none2384 IMPLICIT NONE 2385 2385 2386 2386 ! on réajuste aussi les valeurs des traceurs pour la … … 2424 2424 ! use isotopes_verif_mod, only: iso_verif_traceur_pbidouille 2425 2425 !#endif 2426 implicit none2426 IMPLICIT NONE 2427 2427 2428 2428 ! on réajuste les valeurs des traceurs pour qu'il n'y ai pas de … … 2497 2497 USE isotopes_verif_mod 2498 2498 #endif 2499 implicit none2499 IMPLICIT NONE 2500 2500 2501 2501 ! inputs … … 2577 2577 use isotopes_verif_mod 2578 2578 #endif 2579 implicit none2579 IMPLICIT NONE 2580 2580 2581 2581 ! même but que iso_verif_traceur_jbid_pos, mais n'induit … … 2644 2644 SUBROUTINE iso_verif_traceur_jbid_vect1D(x,n) 2645 2645 USE isotopes_mod, ONLY: bidouille_anti_divergence,iso_eau,ridicule 2646 implicit none2646 IMPLICIT NONE 2647 2647 2648 2648 ! version vectrisée de iso_verif_traceur_jbidouille … … 2687 2687 SUBROUTINE iso_verif_traceur_pbidouille(x,err_msg) 2688 2688 use isotopes_verif_mod 2689 implicit none2689 IMPLICIT NONE 2690 2690 ! vérifier des choses sur les traceurs 2691 2691 ! * toutes les zones donne t l'istope total … … 2714 2714 use isotrac_mod, only: ridicule_trac 2715 2715 use isotopes_verif_mod 2716 implicit none2716 IMPLICIT NONE 2717 2717 ! vérifier des choses sur les traceurs 2718 2718 ! * toutes les zones donne t l'istope total … … 2775 2775 use isotopes_mod, ONLY: iso_HDO,bidouille_anti_divergence 2776 2776 use isotopes_verif_mod 2777 implicit none2777 IMPLICIT NONE 2778 2778 2779 2779 ! inputs -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyaqua_mod.F90
r5112 r5113 630 630 ierr = nf90_put_att(nid, nf90_global, 'title', & 631 631 'Fichier conditions aux limites') 632 ! !ierr = nf90_def_dim (nid, "points_physiques", klon, ndim)632 ! ierr = nf90_def_dim (nid, "points_physiques", klon, ndim) 633 633 ierr = nf90_def_dim(nid, 'points_physiques', klon_glo, ndim) 634 634 ierr = nf90_def_dim(nid, 'time', nf90_unlimited, ntim) -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyredem.F90
r5112 r5113 509 509 !USE phyredem, ONLY: put_field_srf1 510 510 511 implicit none511 IMPLICIT NONE 512 512 513 513 ! equivalent isotopique de phyredem -
LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90
r5112 r5113 24 24 USE add_phys_tend_mod, ONLY: add_pbl_tend, add_phys_tend, diag_phys_tend, prt_enerbil, & 25 25 fl_ebil, fl_cor_ebil 26 USE assert_m, only: assert26 USE lmdz_assert, only: assert 27 27 USE change_srf_frac_mod 28 28 USE conf_phys_m, only: conf_phys … … 679 679 680 680 ! REAL,allocatable,save :: run_off_lic_0(:) 681 ! !$OMP THREADPRIVATE(run_off_lic_0)681 !$OMP THREADPRIVATE(run_off_lic_0) 682 682 !ym SAVE run_off_lic_0 683 683 !KE43
Note: See TracChangeset
for help on using the changeset viewer.
