Changeset 4033


Ignore:
Timestamp:
Nov 30, 2021, 2:00:13 PM (3 years ago)
Author:
crisi
Message:

nettoyage de sorties ecran

Location:
LMDZ6/trunk/libf/phylmdiso
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmdiso/calwake.F90

    r4004 r4033  
    140140  REAL                                               :: rdcp
    141141
    142 #ifdef ISOVERIF       
    143         write(*,*) 'calwake 143 tmp: wake_deltaq(419,1)=',wake_deltaq(419,1)
    144         write(*,*) 'wake_deltaxt(iso_eau,419,1)=',wake_deltaxt(iso_eau,419,1)
    145 #endif
    146142  IF (prt_level >= 10) THEN
    147143    print *, '-> calwake, wake_s, wgen input ', wake_s(1), wgen(1)

  • LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90

    r4004 r4033  
    21232123#endif 
    21242124#ifdef ISOVERIF
    2125         write(*,*) 'cv3_routine 1259: avant condiso'
     2125        !write(*,*) 'cv3_routine 1259: avant condiso'
    21262126        do i=1,ncum
    21272127           if (iso_HDO.gt.0) then           
    … …  
    60076007     &           fxt(iso_HDO,il,i)/fr(il,i), &
    60086008     &           'cv3_yield 3662').eq.1) then
    6009                 write(*,*) 'il,i,icb(il),inb(il)=',il,i,icb(il),inb(il)
    6010                 write(*,*) 'fr(il,i),delt=',fr(il,i),delt
     6009!                write(*,*) 'il,i,icb(il),inb(il)=',il,i,icb(il),inb(il)
     6010!                write(*,*) 'fr(il,i),delt=',fr(il,i),delt
    60116011#ifdef DIAGISO
    60126012                if (fq_ddft(il,i).ne.0.0) then

  • LMDZ6/trunk/libf/phylmdiso/cv3p_mixing.F90

    r4004 r4033  
    12901290#ifdef ISO
    12911291#ifdef ISOVERIF
    1292        write(*,*) 'cv3p_mixing 2540: ', &
    1293         'verif finale en sortant de cv3p_mixing'
    1294        write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1)
     1292!       write(*,*) 'cv3p_mixing 2540: ', &
     1293!        'verif finale en sortant de cv3p_mixing'
     1294!       write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1)
    12951295       do im = 1, nd
    12961296       do jm = 1, nd

  • LMDZ6/trunk/libf/phylmdiso/cva_driver.F90

    r4004 r4033  
    5555#ifdef ISO
    5656  USE infotrac_phy, ONLY: ntraciso,niso,niso,index_trac,ntraceurs_zone
    57   USE isotopes_mod, ONLY: iso_eau,iso_HDO,ridicule,bidouille_anti_divergence
     57  USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,ridicule,bidouille_anti_divergence
    5858#ifdef ISOVERIF
    5959    use isotopes_verif_mod
    … …  
    963963#ifdef ISO
    964964#ifdef ISOVERIF
    965        write(*,*) 'cva_driver 621: après cv3_undilute1'
     965       write(*,*) 'cva_driver 621: apres cv3_undilute1'
    966966       do k=1,nd
    967967        do i = 1, len
    … …  
    11211121#ifdef ISO
    11221122#ifdef ISOVERIF
    1123        write(*,*) 'cva_driver 720: après cv3_compress'
    1124        write(*,*) 'len, nloc, ncum,nd=',len, nloc, ncum,nd
     1123       write(*,*) 'cva_driver 720: apres cv3_compress'
     1124!       write(*,*) 'len, nloc, ncum,nd=',len, nloc, ncum,nd
    11251125       do k=1,nd
    11261126        do i = 1, ncum
    … …  
    11491149         call iso_verif_positif(qnk(i),'cva_driver 966b') 
    11501150       enddo !do i = 1, ncum
    1151        write(*,*) 'cva_driver 1142: après cv3_compress OK'
     1151!       write(*,*) 'cva_driver 1142: apres cv3_compress OK'
    11521152#endif
    11531153#endif
    … …  
    13571357#ifdef ISO
    13581358#ifdef ISOVERIF
    1359        write(*,*) 'cva_driver 837: après cv3_mixing'
    1360        write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1)
     1359       write(*,*) 'cva_driver 837: apres cv3_mixing'
     1360!       write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1)
    13611361       do k=1,nd
    13621362       do j = 1, nd
    … …  
    15971597            call iso_verif_aberrant( &
    15981598     &          (xt(iso_HDO,i,k)+delt*fxt(iso_HDO,i,k)) &
    1599      &          /(q(i,k)+delt*fq(i,k)),'cva_driver 855')
     1599     &          /(q(i,k)+delt*fq(i,k)),'cva_driver 855a')
     1600                if (iso_O18.gt.0) then
     1601            call iso_verif_O18_aberrant( &
     1602     &          (xt(iso_HDO,i,k)+delt*fxt(iso_HDO,i,k)) &
     1603     &          /(q(i,k)+delt*fq(i,k)), &
     1604     &          (xt(iso_O18,i,k)+delt*fxt(iso_O18,i,k)) &
     1605     &          /(q(i,k)+delt*fq(i,k)),'cva_driver 855b')
     1606                endif
    16001607          endif
    16011608         endif

  • LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90

    r3940 r4033  
    345345    snow_evap = 0.
    346346 
    347 #ifdef ISOVERIF
    348         write(*,*) 'klon,snow_evap(413)=',klon,snow_evap(413)
    349 #endif
    350347
    351348    IF (.NOT. ok_lic_cond) THEN
    … …  
    358355         snow = MAX(0.0, snow)                   !---just in case
    359356      END WHERE
    360 #ifdef ISOVERIF
    361         write(*,*) 'fonte_neige 342: snow_evap(413)=',snow_evap(413)
    362 #endif
    363357
    364358    ELSE
    … …  
    367361      snow = snow - snow_evap * dtime         !---snow that remains or deposits on the ground
    368362      snow = MAX(0.0, snow)                   !---just in case
    369 #ifdef ISOVERIF
    370         write(*,*) 'fonte_neige 351: snow_evap(413)=',snow_evap(413)
    371         write(*,*) 'evap(413)=',evap(413)
    372         write(*,*) 'snow(413),dtime=',snow(413),dtime
    373 #endif
    374363
    375364   ENDIF
    … …  
    380369        snow_evap_diag(:)=snow_evap(:)
    381370        coeff_rel_diag=coeff_rel
    382 #ifdef ISOVERIF
    383         write(*,*) 'fonte neige 350: snow_evap_diag(1)=',snow_evap_diag(1)
    384         write(*,*) 'klon,snow_evap_diag(413)=',klon,snow_evap_diag(413)
    385         write(*,*) 'snow_evap(413)=',snow_evap(413)
    386 #endif
    387371#endif
    388372
    … …  
    681665      IF (nisurf == is_lic) THEN
    682666!         coeff_rel = dtime/(tau_calv * rday)
    683 #ifdef ISOVERIF
    684         j=61
    685         write(*,*) 'fonte_neige 636:'
    686         write(*,*) 'run_off_lic_0(j)=',run_off_lic_0(j)
    687         write(*,*) 'xtrun_off_lic_0(:,j)=',xtrun_off_lic_0(:,j)
    688 #endif
     667
    689668         DO i = 1, knon
    690669          j = knindex(i)

  • LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90

    r3927 r4033  
    63666366              else !if (qp(il,i).gt.0) then
    63676367                  ! si qp est négatif, on met les isos dedans à 0
    6368                 write(*,*) 'appel_stewart temporaire 230: qp=', &
    6369      &        qp(il,i)
    63706368                do ixt=1,ntraciso
    63716369                 xtp_avantevap(ixt,il)=0.0
    … …  
    95989596              else !if (qp(il,i).gt.0) then
    95999597                  ! si qp est négatif, on met les isos dedans à 0
    9600                 write(*,*) 'appel_stewart_np temporaire 230: qp=', &
    9601      &                         qp(il,i)
    96029598                do ixt=1,ntraciso
    96039599                 xtp_avantevap(ixt,il)=0.0
    … …  
    99819977#ifdef ISOVERIF
    99829978      ! vérif de la compression
    9983       write(*,*) 'appel_stewart_np tmp 506: ', &
    9984      &          'après compress_condensation_facftmr'
     9979!      write(*,*) 'appel_stewart_np tmp 506: ', &
     9980!     &          'après compress_condensation_facftmr'
    99859981!      write(*,*) 'sigd_cas(1:3)=',sigd_cas(1:3)
    9986       if (ncas_condensation_facftmr.ge.4) then
    9987           write(*,*) 'cas_condensation_facftmr(4)=', &
    9988      &          cas_condensation_facftmr(4)
    9989       endif
     9982!      if (ncas_condensation_facftmr.ge.4) then
     9983!          write(*,*) 'cas_condensation_facftmr(4)=', &
     9984!     &          cas_condensation_facftmr(4)
     9985!      endif
    99909986      do il=1,ncas_condensation_facftmr
    99919987        call iso_verif_egalite_choix((Pqisup_cas(il)), &
    … …  
    1015610152#ifdef ISOVERIF
    1015710153      ! vérif de la compression
    10158       write(*,*) 'appel_stewart_np tmp 616: ', &
    10159      &          'apres compress condensation_nofacftmr'
    10160       write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3)
     10154!      write(*,*) 'appel_stewart_np tmp 616: ', &
     10155!     &          'apres compress condensation_nofacftmr'
     10156!      write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3)
    1016110157      do il=1,ncas_condensation_nofacftmr
    1016210158        call iso_verif_egalite_choix((Pqisup_cas(il)), &
    … …  
    1024410240       enddo !do izone=1,ntraceurs_zone
    1024510241#ifdef ISOVERIF
    10246        write(*,*) 'appel_stewart_np tmp 690: ', &
    10247      &          'fin du cas condensation_nofacftmr'
     10242!       write(*,*) 'appel_stewart_np tmp 690: ', &
     10243!     &          'fin du cas condensation_nofacftmr'
    1024810244            do il=1,ncas_condensation_nofacftmr
    1024910245               call iso_verif_traceur(xtp &
    … …  
    1072710723
    1072810724#ifdef ISOVERIF
    10729       write(*,*) 'appel_stewart_np tmp 898 après compress glace'
     10725!      write(*,*) 'appel_stewart_np tmp 898 apres compress glace'
    1073010726!      write(*,*) 'qp_avantevap_cas(1),qp_avantevap(cas(1))=',
    1073110727!     &   qp_avantevap_cas(1),qp_avantevap(cas_evap_glace(1))   
    10732       write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
     10728      !write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
    1073310729      ! vérif de la compression
    1073410730      do il=1,ncas_evap_glace
    … …  
    1131111307     &       fac_ftmr_cas(1))
    1131211308        else !if (no_pce.eq.1) then
    11313 #ifdef ISOVERIF                   
    11314             write(*,*) 'appel_stewart_np 1957 tmp'
    11315 #endif           
     11309           
    1131611310      call stewart_explicite_vectall(ncas, &
    1131711311     &       qp_avantevap_cas(1),xtp_avantevap_cas(1,1), &
    … …  
    1185711851     &       T_cas(1))
    1185811852      else !if (frac_sublim.eq.1) then
    11859 #ifdef ISOVERIF
    11860             write(*,*) 'appel_stewart_explicite_np 2269'
    11861             write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
    11862             write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1)
    11863             write(*,*) 'Eqi_cas(1)=',Eqi_cas(1)
    11864             write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1)
    11865             write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1)
    11866 #endif         
     11853!#ifdef ISOVERIF
     11854!            write(*,*) 'appel_stewart_explicite_np 2269'
     11855!            write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1)
     11856!            write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1)
     11857!            write(*,*) 'Eqi_cas(1)=',Eqi_cas(1)
     11858!            write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1)
     11859!            write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1)
     11860!#endif         
    1186711861            call stewart_sublim_nofrac_vectall( &
    1186811862     &        ncas,qp_avantevap_cas(1), &
    … …  
    1405014044
    1405114045#ifdef ISOVERIF
    14052         write(*,*) 'calcul_iso_surf_sic 175: entree'
    14053 #endif 
    14054 #ifdef ISOVERIF
    1405514046        do i=1,knon
    1405614047         do ixt=1,ntraciso
    … …  
    1471814709
    1471914710#ifdef ISOVERIF   
    14720       write(*,*) 'calcul_iso_surf_ter 494'
     14711!      write(*,*) 'calcul_iso_surf_ter 494'
    1472114712      do i=1,knon         
    1472214713        if (iso_eau.gt.0) then 
    … …  
    1476614757#endif     
    1476714758#ifdef ISOVERIF
    14768         write(*,*) 'calcul_iso_surf_ter 910'
     14759!        write(*,*) 'calcul_iso_surf_ter 910'
    1476914760        do i=1,knon
    1477014761          if (iso_eau.gt.0) then

  • LMDZ6/trunk/libf/phylmdiso/isotopes_verif_mod.F90

    r3927 r4033  
    528528            write(*,*) 'o17excess=',o17excess(R17,R18)
    529529            write(*,*) 'deltaO17=',(R17/tnat(iso_o17)-1.0)*1000.0
    530             write(*,*) 'deltaO18=',(R18/tnat(iso_o18)-1.0)*1000.0
     530            write(*,*) 'deltaO18=',(R18/tnat(iso_O18)-1.0)*1000.0
    531531            ! attention, vérifier que la ligne suivante est bien activée
    532532            iso_verif_aberrant_o17_nostop=1
    … …  
    998998
    999999       
    1000         subroutine iso_verif_o18_aberrant(Rd,Ro,err_msg)
     1000        subroutine iso_verif_O18_aberrant(Rd,Ro,err_msg)
    10011001        implicit none
    10021002
    … …  
    10091009
    10101010        ! local
    1011         !integer iso_verif_o18_aberrant_nostop
    1012 
    1013         if (iso_verif_o18_aberrant_nostop(Rd,Ro,err_msg).eq.1) then
     1011        !integer iso_verif_O18_aberrant_nostop
     1012
     1013        if (iso_verif_O18_aberrant_nostop(Rd,Ro,err_msg).eq.1) then
    10141014            stop
    10151015        endif
    10161016
    1017         end subroutine iso_verif_o18_aberrant
    1018        
    1019         function iso_verif_o18_aberrant_nostop(Rd,Ro,err_msg)
     1017        end subroutine iso_verif_O18_aberrant
     1018       
     1019        function iso_verif_O18_aberrant_nostop(Rd,Ro,err_msg)
    10201020        USE isotopes_mod, ONLY: tnat, iso_HDO, iso_O18
    10211021        implicit none
    … …  
    10301030
    10311031        ! outputs
    1032         integer iso_verif_o18_aberrant_nostop
     1032        integer iso_verif_O18_aberrant_nostop
    10331033
    10341034        !locals
    … …  
    10391039        dexcess=deltaD-8*deltao
    10401040
    1041         iso_verif_o18_aberrant_nostop=0
     1041        iso_verif_O18_aberrant_nostop=0
    10421042        if ((deltaD.lt.deltaDmin).or.(deltao.lt.deltaDmin/2.0).or. &
    10431043     &        (deltaD.gt.deltalim).or.(deltao.gt.deltalim/8.0).or. &
    10441044     &        ((deltaD.gt.-500.0).and.((dexcess.lt.dexcess_min) &
    10451045     &        .or.(dexcess.gt.dexcess_max)))) then
    1046             write(*,*) 'erreur detectee par iso_verif_o18_aberrant:'
     1046            write(*,*) 'erreur detectee par iso_verif_O18_aberrant:'
    10471047            write(*,*) err_msg
    10481048            write(*,*) 'delta180=',deltao
    … …  
    10501050            write(*,*) 'Dexcess=',dexcess
    10511051!            stop
    1052             iso_verif_o18_aberrant_nostop=1
     1052            iso_verif_O18_aberrant_nostop=1
    10531053          endif
    10541054
    … …  
    10601060
    10611061        return
    1062         end function iso_verif_o18_aberrant_nostop
     1062        end function iso_verif_O18_aberrant_nostop
    10631063
    10641064
    … …  
    11461146           
    11471147           o17excess=1e6*(log(R17/tnat(iso_o17)) &
    1148      &           -0.528*log(R18/tnat(iso_o18)))
     1148     &           -0.528*log(R18/tnat(iso_O18)))
    11491149!           write(*,*) 'o17excess=',o17excess
    11501150        else
    … …  
    14981498        end subroutine iso_verif_aberrant_vect2Dch     
    14991499
    1500         subroutine iso_verif_o18_aberrant_enc_vect2D( &
     1500        subroutine iso_verif_O18_aberrant_enc_vect2D( &
    15011501     &           xt,q,err_msg,ni,n,m)
    15021502
    … …  
    15501550        endif !if (iso_HDO.gt.0) then
    15511551
    1552         end subroutine iso_verif_o18_aberrant_enc_vect2D   
     1552        end subroutine iso_verif_O18_aberrant_enc_vect2D   
    15531553
    15541554

  • LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90

    r3975 r4033  
    467467#ifdef ISO
    468468! isotopes: tout est externalisé
    469 #ifdef ISOVERIF
    470         write(*,*) 'ocean_forced_mod 377: call calcul_iso_surf_sic_vectall'
    471         write(*,*) 'klon,knon=',klon,knon
    472 #endif
     469!#ifdef ISOVERIF
     470!        write(*,*) 'ocean_forced_mod 377: call calcul_iso_surf_sic_vectall'
     471!        write(*,*) 'klon,knon=',klon,knon
     472!#endif
    473473         call calcul_iso_surf_sic_vectall(klon,knon, &
    474474             &   evap,snow_evap_diag,Tsurf_new,Roce,snow, &
    … …  
    480480     &  )   
    481481#ifdef ISOVERIF
    482         write(*,*) 'ocean_forced_mod 391: sortie calcul_iso_surf_sic_vectall'
     482        !write(*,*) 'ocean_forced_mod 391: sortie calcul_iso_surf_sic_vectall'
    483483          if (iso_eau.gt.0) then
    484484           do i=1,knon 

  • LMDZ6/trunk/libf/phylmdiso/phys_output_write_mod.F90

    r3940 r4033  
    25302530
    25312531#ifdef ISO
    2532     write(*,*) 'phys_output_write_mod 1745'
    2533         write(*,*) 'o_precip=',o_precip
    25342532    do ixt=1,ntraciso
    25352533        write(*,*) 'ixt'
    2536         write(*,*) 'o_xtprecip(ixt)=',ixt,o_xtprecip(ixt)
    25372534        IF (vars_defined) THEN
    25382535             zx_tmp_fi2d(:) = xtrain_fall(ixt,:) + xtsnow_fall(ixt,:)
    … …  
    25402537        CALL histwrite_phy(o_xtprecip(ixt), zx_tmp_fi2d)
    25412538
    2542     write(*,*) 'phys_output_write_mod 2409'
    25432539        IF (vars_defined) THEN
    25442540             zx_tmp_fi2d(:) = xtrain_lsc(ixt,:) + xtsnow_lsc(ixt,:)
    … …  
    25602556        CALL histwrite_phy(o_xtcond(ixt), zx_tmp_fi3d)
    25612557
    2562         write(*,*) 'phys_output_write_mod 2429'
    25632558        CALL histwrite_phy(o_dxtdyn(ixt), d_xt_dyn(ixt,:,:))
    25642559        CALL histwrite_phy(o_dxtldyn(ixt), d_xtl_dyn(ixt,:,:))
    2565         write(*,*) 'phys_output_write_mod 2432'
    25662560
    25672561        IF (vars_defined) THEN
    … …  
    25692563        ENDIF
    25702564        CALL histwrite_phy(o_dxtcon(ixt), zx_tmp_fi3d)
    2571         write(*,*) 'phys_output_write_mod 2438'
    25722565
    25732566        IF (vars_defined) THEN

  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r4009 r4033  
    24212421       do itr=1,nqtottr
    24222422         iq=itr_indice(itr)
    2423 #ifdef ISOVERIF
    2424          write(*,*) 'physiq 1973: itr,iq=',itr,iq
    2425          write(*,*) 'qx(1,1,iq)=',qx(1,1,iq)
    2426 #endif
     2423!#ifdef ISOVERIF
     2424!         write(*,*) 'physiq 1973: itr,iq=',itr,iq
     2425!         write(*,*) 'qx(1,1,iq)=',qx(1,1,iq)
     2426!#endif
    24272427         DO  k = 1, klev
    24282428             DO  i = 1, klon
    … …  
    35233523     &           'physiq 1456, avant concvl')
    35243524            endif   
    3525 #endif      
     3525#endif     
    35263526      enddo !do k=1,nlev
    3527    enddo  !do i=1,klon
    3528         if (iso_eau.gt.0) then
    3529         i=1
    3530         k=1
    3531         write(*,*) 'physic 2376: xt_seri(iso_eau,i,k),q_seri(i,k)=',xt_seri(iso_eau,i,k),q_seri(i,k)
    3532         write(*,*) 'xt_seri(:,i,k)=',xt_seri(:,i,k)
    3533         write(*,*) 'physic 2376: xt_x(iso_eau,i,k),q_x(i,k)=',xt_x(iso_eau,i,k),q_x(i,k)
    3534         write(*,*) 'xt_x(:,i,k)=',xt_x(:,i,k)
    3535         endif
    3536 #endif
     3527   enddo  !do i=1,klon   
     3528#endif
    35373529!ISOVERIF         
    35383530          if ((bidouille_anti_divergence).and. &
    … …  
    36103602#ifdef ISOVERIF
    36113603!          write(*,*) 'q_detrainement(1,:)=',q_detrainement(1,:)
    3612           call iso_verif_noNaN_vect2D(d_xt_con, &
     3604      call iso_verif_noNaN_vect2D(d_xt_con, &
    36133605     &           'physiq 3203a apres conv',ntraciso,klon,klev)
    3614           call iso_verif_noNaN_vect2D(xt_seri, &
    3615      &           'physiq 3203b apres conv',ntraciso,klon,klev)         
     3606      call iso_verif_noNaN_vect2D(xt_seri, &
     3607     &           'physiq 3203b apres conv',ntraciso,klon,klev)
     3608      if (iso_HDO.gt.0) then
     3609      call iso_verif_aberrant_enc_vect2D( &
     3610     &           xt_seri,q_seri, &
     3611     &           'physiq 3619a',ntraciso,klon,klev)
     3612      endif
     3613      if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
     3614      call iso_verif_O18_aberrant_enc_vect2D( &
     3615     &           xt_seri,q_seri, &
     3616     &           'physiq 3619b',ntraciso,klon,klev)
     3617      endif         
    36163618#endif
    36173619#ifdef ISOVERIF
    … …  
    39063908#ifdef ISO
    39073909#ifdef ISOVERIF
    3908      write(*,*) 'physiq 3425'
     3910     write(*,*) 'physiq 3425: apres convection'
    39093911      if (iso_HDO.gt.0) then
    39103912      call iso_verif_aberrant_enc_vect2D( &
    … …  
    39153917      call iso_verif_O18_aberrant_enc_vect2D( &
    39163918     &           xt_seri,q_seri, &
    3917      &           'physiq 3691a',ntraciso,klon,klev)
     3919     &           'physiq 3691b',ntraciso,klon,klev)
    39183920      endif
    39193921#ifdef ISOTRAC

  • LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90

    r3975 r4033  
    128128#ifdef ISO
    129129#ifdef ISOVERIF
    130         write(*,*) 'surf_land_bucket 152'
     130        !write(*,*) 'surf_land_bucket 152'
    131131        do i=1,knon
    132132          if (iso_eau.gt.0) then
    … …  
    146146         enddo !do ixt=1,niso
    147147        enddo !do i=1,knon
    148         write(*,*) 'surf_land_bucket 152'
     148        !write(*,*) 'surf_land_bucket 152'
    149149#endif
    150150#endif
    … …  
    211211   ! verif
    212212#ifdef ISOVERIF
    213     write(*,*) 'surf_land_bucket 211'
     213    !write(*,*) 'surf_land_bucket 211'
    214214    do i=1,knon
    215215      if (iso_eau.gt.0) then
    … …  
    236236   enddo !do i=1,knon
    237237#ifdef ISOVERIF
    238         write(*,*) 'surf_land_bucket 235'
     238       ! write(*,*) 'surf_land_bucket 235'
    239239       do i=1,knon 
    240240        if (iso_eau.gt.0) then
    … …  
    243243        endif
    244244      enddo !do i=1,knon
    245         write(*,*) 'snow_prec(1)=',snow_prec(1)
    246         write(*,*) 'xtsnow(:,1)=',xtsnow(:,1)
    247245#endif   
    248246#endif   
    … …  
    261259#ifdef ISO
    262260#ifdef ISOVERIF
    263         write(*,*) 'surf_land_bucket 258'
    264         write(*,*) 'snow_prec(1)=',snow_prec(1)
    265         write(*,*) 'xtsnow(:,1)=',xtsnow(:,1)
    266261        do i=1,knon
    267262         do ixt=1,niso
    … …  
    271266#endif
    272267#ifdef ISOVERIF
    273         write(*,*) 'surf_land_bucket 235'
     268        !write(*,*) 'surf_land_bucket 235'
    274269        do i=1,knon
    275270          if (iso_eau.gt.0) then
    … …  
    297292          endif !if ((iso_O17.gt.0).and.(iso_O18.gt.0)) then
    298293        enddo  !do i=1,knon
    299         write(*,*) 'surf_land_mod 291'
    300         write(*,*) 'snow_evap_diag(1)=',snow_evap_diag(1)
     294        !write(*,*) 'surf_land_mod 291'
     295        !write(*,*) 'snow_evap_diag(1)=',snow_evap_diag(1)
    301296#endif         
    302297        call calcul_iso_surf_ter_vectall(klon,knon, &

  • LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90

    r3927 r4033  
    166166#ifdef ISO
    167167#ifdef ISOVERIF
    168         write(*,*) 'surf_land_mod 162'
     168!        write(*,*) 'surf_land_mod 162'
    169169        do i=1,knon
    170170          if (iso_eau.gt.0) then
    … …  
    179179#endif
    180180#ifdef ISOVERIF
    181        write(*,*) 'surf_land 169: ok_veget=',ok_veget
     181!       write(*,*) 'surf_land 169: ok_veget=',ok_veget
    182182        do i=1,knon
    183183         do ixt=1,ntraciso
    … …  
    262262#ifdef ISO
    263263#ifdef ISOVERIF
    264        write(*,*) 'surf_land 247'
     264!       write(*,*) 'surf_land 247'
    265265        call iso_verif_egalite_vect1D( &
    266266     &           xtsnow,snow,'surf_land_mod 207',niso,klon)
    … …  
    306306#ifdef ISO
    307307#ifdef ISOVERIF
    308      write(*,*) 'surf_land 237: sortie'   
     308!     write(*,*) 'surf_land 237: sortie'   
    309309     do i=1,knon
    310310        if (iso_eau.gt.0) then

  • LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90

    r3975 r4033  
    375375   ! verif
    376376#ifdef ISOVERIF
    377      write(*,*) 'surf_land_ice 1499'   
     377     !write(*,*) 'surf_land_ice 1499'   
    378378     do i=1,knon
    379379       if (iso_eau.gt.0) then
Note: See TracChangeset for help on using the changeset viewer.