Index: LMDZ6/trunk/libf/phylmdiso/calwake.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/calwake.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/calwake.F90 (revision 4033) @@ -140,8 +140,4 @@ REAL :: rdcp -#ifdef ISOVERIF - write(*,*) 'calwake 143 tmp: wake_deltaq(419,1)=',wake_deltaq(419,1) - write(*,*) 'wake_deltaxt(iso_eau,419,1)=',wake_deltaxt(iso_eau,419,1) -#endif IF (prt_level >= 10) THEN print *, '-> calwake, wake_s, wgen input ', wake_s(1), wgen(1) Index: LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90 (revision 4033) @@ -2123,5 +2123,5 @@ #endif #ifdef ISOVERIF - write(*,*) 'cv3_routine 1259: avant condiso' + !write(*,*) 'cv3_routine 1259: avant condiso' do i=1,ncum if (iso_HDO.gt.0) then @@ -6007,6 +6007,6 @@ & fxt(iso_HDO,il,i)/fr(il,i), & & 'cv3_yield 3662').eq.1) then - write(*,*) 'il,i,icb(il),inb(il)=',il,i,icb(il),inb(il) - write(*,*) 'fr(il,i),delt=',fr(il,i),delt +! write(*,*) 'il,i,icb(il),inb(il)=',il,i,icb(il),inb(il) +! write(*,*) 'fr(il,i),delt=',fr(il,i),delt #ifdef DIAGISO if (fq_ddft(il,i).ne.0.0) then Index: LMDZ6/trunk/libf/phylmdiso/cv3p_mixing.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/cv3p_mixing.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/cv3p_mixing.F90 (revision 4033) @@ -1290,7 +1290,7 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'cv3p_mixing 2540: ', & - 'verif finale en sortant de cv3p_mixing' - write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1) +! write(*,*) 'cv3p_mixing 2540: ', & +! 'verif finale en sortant de cv3p_mixing' +! write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1) do im = 1, nd do jm = 1, nd Index: LMDZ6/trunk/libf/phylmdiso/cva_driver.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/cva_driver.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/cva_driver.F90 (revision 4033) @@ -55,5 +55,5 @@ #ifdef ISO USE infotrac_phy, ONLY: ntraciso,niso,niso,index_trac,ntraceurs_zone - USE isotopes_mod, ONLY: iso_eau,iso_HDO,ridicule,bidouille_anti_divergence + USE isotopes_mod, ONLY: iso_eau,iso_HDO,iso_O18,ridicule,bidouille_anti_divergence #ifdef ISOVERIF use isotopes_verif_mod @@ -963,5 +963,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'cva_driver 621: après cv3_undilute1' + write(*,*) 'cva_driver 621: apres cv3_undilute1' do k=1,nd do i = 1, len @@ -1121,6 +1121,6 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'cva_driver 720: après cv3_compress' - write(*,*) 'len, nloc, ncum,nd=',len, nloc, ncum,nd + write(*,*) 'cva_driver 720: apres cv3_compress' +! write(*,*) 'len, nloc, ncum,nd=',len, nloc, ncum,nd do k=1,nd do i = 1, ncum @@ -1149,5 +1149,5 @@ call iso_verif_positif(qnk(i),'cva_driver 966b') enddo !do i = 1, ncum - write(*,*) 'cva_driver 1142: après cv3_compress OK' +! write(*,*) 'cva_driver 1142: apres cv3_compress OK' #endif #endif @@ -1357,6 +1357,6 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'cva_driver 837: après cv3_mixing' - write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1) + write(*,*) 'cva_driver 837: apres cv3_mixing' +! write(*,*) 'qent,xtent(1,1,1)=',qent(1,1,1),xtent(iso_eau,1,1,1) do k=1,nd do j = 1, nd @@ -1597,5 +1597,12 @@ call iso_verif_aberrant( & & (xt(iso_HDO,i,k)+delt*fxt(iso_HDO,i,k)) & - & /(q(i,k)+delt*fq(i,k)),'cva_driver 855') + & /(q(i,k)+delt*fq(i,k)),'cva_driver 855a') + if (iso_O18.gt.0) then + call iso_verif_O18_aberrant( & + & (xt(iso_HDO,i,k)+delt*fxt(iso_HDO,i,k)) & + & /(q(i,k)+delt*fq(i,k)), & + & (xt(iso_O18,i,k)+delt*fxt(iso_O18,i,k)) & + & /(q(i,k)+delt*fq(i,k)),'cva_driver 855b') + endif endif endif Index: LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/fonte_neige_mod.F90 (revision 4033) @@ -345,7 +345,4 @@ snow_evap = 0. -#ifdef ISOVERIF - write(*,*) 'klon,snow_evap(413)=',klon,snow_evap(413) -#endif IF (.NOT. ok_lic_cond) THEN @@ -358,7 +355,4 @@ snow = MAX(0.0, snow) !---just in case END WHERE -#ifdef ISOVERIF - write(*,*) 'fonte_neige 342: snow_evap(413)=',snow_evap(413) -#endif ELSE @@ -367,9 +361,4 @@ snow = snow - snow_evap * dtime !---snow that remains or deposits on the ground snow = MAX(0.0, snow) !---just in case -#ifdef ISOVERIF - write(*,*) 'fonte_neige 351: snow_evap(413)=',snow_evap(413) - write(*,*) 'evap(413)=',evap(413) - write(*,*) 'snow(413),dtime=',snow(413),dtime -#endif ENDIF @@ -380,9 +369,4 @@ snow_evap_diag(:)=snow_evap(:) coeff_rel_diag=coeff_rel -#ifdef ISOVERIF - write(*,*) 'fonte neige 350: snow_evap_diag(1)=',snow_evap_diag(1) - write(*,*) 'klon,snow_evap_diag(413)=',klon,snow_evap_diag(413) - write(*,*) 'snow_evap(413)=',snow_evap(413) -#endif #endif @@ -681,10 +665,5 @@ IF (nisurf == is_lic) THEN ! coeff_rel = dtime/(tau_calv * rday) -#ifdef ISOVERIF - j=61 - write(*,*) 'fonte_neige 636:' - write(*,*) 'run_off_lic_0(j)=',run_off_lic_0(j) - write(*,*) 'xtrun_off_lic_0(:,j)=',xtrun_off_lic_0(:,j) -#endif + DO i = 1, knon j = knindex(i) Index: LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90 (revision 4033) @@ -6366,6 +6366,4 @@ else !if (qp(il,i).gt.0) then ! si qp est négatif, on met les isos dedans à 0 - write(*,*) 'appel_stewart temporaire 230: qp=', & - & qp(il,i) do ixt=1,ntraciso xtp_avantevap(ixt,il)=0.0 @@ -9598,6 +9596,4 @@ else !if (qp(il,i).gt.0) then ! si qp est négatif, on met les isos dedans à 0 - write(*,*) 'appel_stewart_np temporaire 230: qp=', & - & qp(il,i) do ixt=1,ntraciso xtp_avantevap(ixt,il)=0.0 @@ -9981,11 +9977,11 @@ #ifdef ISOVERIF ! vérif de la compression - write(*,*) 'appel_stewart_np tmp 506: ', & - & 'après compress_condensation_facftmr' +! write(*,*) 'appel_stewart_np tmp 506: ', & +! & 'après compress_condensation_facftmr' ! write(*,*) 'sigd_cas(1:3)=',sigd_cas(1:3) - if (ncas_condensation_facftmr.ge.4) then - write(*,*) 'cas_condensation_facftmr(4)=', & - & cas_condensation_facftmr(4) - endif +! if (ncas_condensation_facftmr.ge.4) then +! write(*,*) 'cas_condensation_facftmr(4)=', & +! & cas_condensation_facftmr(4) +! endif do il=1,ncas_condensation_facftmr call iso_verif_egalite_choix((Pqisup_cas(il)), & @@ -10156,7 +10152,7 @@ #ifdef ISOVERIF ! vérif de la compression - write(*,*) 'appel_stewart_np tmp 616: ', & - & 'apres compress condensation_nofacftmr' - write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3) +! write(*,*) 'appel_stewart_np tmp 616: ', & +! & 'apres compress condensation_nofacftmr' +! write(*,*) 'iso_routines 10153: sigd_cas(1:3)=', sigd_cas(1:3) do il=1,ncas_condensation_nofacftmr call iso_verif_egalite_choix((Pqisup_cas(il)), & @@ -10244,6 +10240,6 @@ enddo !do izone=1,ntraceurs_zone #ifdef ISOVERIF - write(*,*) 'appel_stewart_np tmp 690: ', & - & 'fin du cas condensation_nofacftmr' +! write(*,*) 'appel_stewart_np tmp 690: ', & +! & 'fin du cas condensation_nofacftmr' do il=1,ncas_condensation_nofacftmr call iso_verif_traceur(xtp & @@ -10727,8 +10723,8 @@ #ifdef ISOVERIF - write(*,*) 'appel_stewart_np tmp 898 après compress glace' +! write(*,*) 'appel_stewart_np tmp 898 apres compress glace' ! write(*,*) 'qp_avantevap_cas(1),qp_avantevap(cas(1))=', ! & qp_avantevap_cas(1),qp_avantevap(cas_evap_glace(1)) - write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1) + !write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1) ! vérif de la compression do il=1,ncas_evap_glace @@ -11311,7 +11307,5 @@ & fac_ftmr_cas(1)) else !if (no_pce.eq.1) then -#ifdef ISOVERIF - write(*,*) 'appel_stewart_np 1957 tmp' -#endif + call stewart_explicite_vectall(ncas, & & qp_avantevap_cas(1),xtp_avantevap_cas(1,1), & @@ -11857,12 +11851,12 @@ & T_cas(1)) else !if (frac_sublim.eq.1) then -#ifdef ISOVERIF - write(*,*) 'appel_stewart_explicite_np 2269' - write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1) - write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1) - write(*,*) 'Eqi_cas(1)=',Eqi_cas(1) - write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1) - write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1) -#endif +!#ifdef ISOVERIF +! write(*,*) 'appel_stewart_explicite_np 2269' +! write(*,*) 'Pqiinf_stewart(1)=',Pqiinf_stewart(1) +! write(*,*) 'Pqisup_cas(1)=',Pqisup_cas(1) +! write(*,*) 'Eqi_cas(1)=',Eqi_cas(1) +! write(*,*) 'Eqi_prime_cas(1)=',Eqi_prime_cas(1) +! write(*,*) 'Eqi_stewart(1)=',Eqi_stewart(1) +!#endif call stewart_sublim_nofrac_vectall( & & ncas,qp_avantevap_cas(1), & @@ -14050,7 +14044,4 @@ #ifdef ISOVERIF - write(*,*) 'calcul_iso_surf_sic 175: entree' -#endif -#ifdef ISOVERIF do i=1,knon do ixt=1,ntraciso @@ -14718,5 +14709,5 @@ #ifdef ISOVERIF - write(*,*) 'calcul_iso_surf_ter 494' +! write(*,*) 'calcul_iso_surf_ter 494' do i=1,knon if (iso_eau.gt.0) then @@ -14766,5 +14757,5 @@ #endif #ifdef ISOVERIF - write(*,*) 'calcul_iso_surf_ter 910' +! write(*,*) 'calcul_iso_surf_ter 910' do i=1,knon if (iso_eau.gt.0) then Index: LMDZ6/trunk/libf/phylmdiso/isotopes_verif_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/isotopes_verif_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/isotopes_verif_mod.F90 (revision 4033) @@ -528,5 +528,5 @@ write(*,*) 'o17excess=',o17excess(R17,R18) write(*,*) 'deltaO17=',(R17/tnat(iso_o17)-1.0)*1000.0 - write(*,*) 'deltaO18=',(R18/tnat(iso_o18)-1.0)*1000.0 + write(*,*) 'deltaO18=',(R18/tnat(iso_O18)-1.0)*1000.0 ! attention, vérifier que la ligne suivante est bien activée iso_verif_aberrant_o17_nostop=1 @@ -998,5 +998,5 @@ - subroutine iso_verif_o18_aberrant(Rd,Ro,err_msg) + subroutine iso_verif_O18_aberrant(Rd,Ro,err_msg) implicit none @@ -1009,13 +1009,13 @@ ! local - !integer iso_verif_o18_aberrant_nostop - - if (iso_verif_o18_aberrant_nostop(Rd,Ro,err_msg).eq.1) then + !integer iso_verif_O18_aberrant_nostop + + if (iso_verif_O18_aberrant_nostop(Rd,Ro,err_msg).eq.1) then stop endif - end subroutine iso_verif_o18_aberrant - - function iso_verif_o18_aberrant_nostop(Rd,Ro,err_msg) + end subroutine iso_verif_O18_aberrant + + function iso_verif_O18_aberrant_nostop(Rd,Ro,err_msg) USE isotopes_mod, ONLY: tnat, iso_HDO, iso_O18 implicit none @@ -1030,5 +1030,5 @@ ! outputs - integer iso_verif_o18_aberrant_nostop + integer iso_verif_O18_aberrant_nostop !locals @@ -1039,10 +1039,10 @@ dexcess=deltaD-8*deltao - iso_verif_o18_aberrant_nostop=0 + iso_verif_O18_aberrant_nostop=0 if ((deltaD.lt.deltaDmin).or.(deltao.lt.deltaDmin/2.0).or. & & (deltaD.gt.deltalim).or.(deltao.gt.deltalim/8.0).or. & & ((deltaD.gt.-500.0).and.((dexcess.lt.dexcess_min) & & .or.(dexcess.gt.dexcess_max)))) then - write(*,*) 'erreur detectee par iso_verif_o18_aberrant:' + write(*,*) 'erreur detectee par iso_verif_O18_aberrant:' write(*,*) err_msg write(*,*) 'delta180=',deltao @@ -1050,5 +1050,5 @@ write(*,*) 'Dexcess=',dexcess ! stop - iso_verif_o18_aberrant_nostop=1 + iso_verif_O18_aberrant_nostop=1 endif @@ -1060,5 +1060,5 @@ return - end function iso_verif_o18_aberrant_nostop + end function iso_verif_O18_aberrant_nostop @@ -1146,5 +1146,5 @@ o17excess=1e6*(log(R17/tnat(iso_o17)) & - & -0.528*log(R18/tnat(iso_o18))) + & -0.528*log(R18/tnat(iso_O18))) ! write(*,*) 'o17excess=',o17excess else @@ -1498,5 +1498,5 @@ end subroutine iso_verif_aberrant_vect2Dch - subroutine iso_verif_o18_aberrant_enc_vect2D( & + subroutine iso_verif_O18_aberrant_enc_vect2D( & & xt,q,err_msg,ni,n,m) @@ -1550,5 +1550,5 @@ endif !if (iso_HDO.gt.0) then - end subroutine iso_verif_o18_aberrant_enc_vect2D + end subroutine iso_verif_O18_aberrant_enc_vect2D Index: LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90 (revision 4033) @@ -467,8 +467,8 @@ #ifdef ISO ! isotopes: tout est externalisé -#ifdef ISOVERIF - write(*,*) 'ocean_forced_mod 377: call calcul_iso_surf_sic_vectall' - write(*,*) 'klon,knon=',klon,knon -#endif +!#ifdef ISOVERIF +! write(*,*) 'ocean_forced_mod 377: call calcul_iso_surf_sic_vectall' +! write(*,*) 'klon,knon=',klon,knon +!#endif call calcul_iso_surf_sic_vectall(klon,knon, & & evap,snow_evap_diag,Tsurf_new,Roce,snow, & @@ -480,5 +480,5 @@ & ) #ifdef ISOVERIF - write(*,*) 'ocean_forced_mod 391: sortie calcul_iso_surf_sic_vectall' + !write(*,*) 'ocean_forced_mod 391: sortie calcul_iso_surf_sic_vectall' if (iso_eau.gt.0) then do i=1,knon Index: LMDZ6/trunk/libf/phylmdiso/phys_output_write_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/phys_output_write_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/phys_output_write_mod.F90 (revision 4033) @@ -2530,9 +2530,6 @@ #ifdef ISO - write(*,*) 'phys_output_write_mod 1745' - write(*,*) 'o_precip=',o_precip do ixt=1,ntraciso write(*,*) 'ixt' - write(*,*) 'o_xtprecip(ixt)=',ixt,o_xtprecip(ixt) IF (vars_defined) THEN zx_tmp_fi2d(:) = xtrain_fall(ixt,:) + xtsnow_fall(ixt,:) @@ -2540,5 +2537,4 @@ CALL histwrite_phy(o_xtprecip(ixt), zx_tmp_fi2d) - write(*,*) 'phys_output_write_mod 2409' IF (vars_defined) THEN zx_tmp_fi2d(:) = xtrain_lsc(ixt,:) + xtsnow_lsc(ixt,:) @@ -2560,8 +2556,6 @@ CALL histwrite_phy(o_xtcond(ixt), zx_tmp_fi3d) - write(*,*) 'phys_output_write_mod 2429' CALL histwrite_phy(o_dxtdyn(ixt), d_xt_dyn(ixt,:,:)) CALL histwrite_phy(o_dxtldyn(ixt), d_xtl_dyn(ixt,:,:)) - write(*,*) 'phys_output_write_mod 2432' IF (vars_defined) THEN @@ -2569,5 +2563,4 @@ ENDIF CALL histwrite_phy(o_dxtcon(ixt), zx_tmp_fi3d) - write(*,*) 'phys_output_write_mod 2438' IF (vars_defined) THEN Index: LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (revision 4033) @@ -2421,8 +2421,8 @@ do itr=1,nqtottr iq=itr_indice(itr) -#ifdef ISOVERIF - write(*,*) 'physiq 1973: itr,iq=',itr,iq - write(*,*) 'qx(1,1,iq)=',qx(1,1,iq) -#endif +!#ifdef ISOVERIF +! write(*,*) 'physiq 1973: itr,iq=',itr,iq +! write(*,*) 'qx(1,1,iq)=',qx(1,1,iq) +!#endif DO k = 1, klev DO i = 1, klon @@ -3523,16 +3523,8 @@ & 'physiq 1456, avant concvl') endif -#endif +#endif enddo !do k=1,nlev - enddo !do i=1,klon - if (iso_eau.gt.0) then - i=1 - k=1 - write(*,*) 'physic 2376: xt_seri(iso_eau,i,k),q_seri(i,k)=',xt_seri(iso_eau,i,k),q_seri(i,k) - write(*,*) 'xt_seri(:,i,k)=',xt_seri(:,i,k) - write(*,*) 'physic 2376: xt_x(iso_eau,i,k),q_x(i,k)=',xt_x(iso_eau,i,k),q_x(i,k) - write(*,*) 'xt_x(:,i,k)=',xt_x(:,i,k) - endif -#endif + enddo !do i=1,klon +#endif !ISOVERIF if ((bidouille_anti_divergence).and. & @@ -3610,8 +3602,18 @@ #ifdef ISOVERIF ! write(*,*) 'q_detrainement(1,:)=',q_detrainement(1,:) - call iso_verif_noNaN_vect2D(d_xt_con, & + call iso_verif_noNaN_vect2D(d_xt_con, & & 'physiq 3203a apres conv',ntraciso,klon,klev) - call iso_verif_noNaN_vect2D(xt_seri, & - & 'physiq 3203b apres conv',ntraciso,klon,klev) + call iso_verif_noNaN_vect2D(xt_seri, & + & 'physiq 3203b apres conv',ntraciso,klon,klev) + if (iso_HDO.gt.0) then + call iso_verif_aberrant_enc_vect2D( & + & xt_seri,q_seri, & + & 'physiq 3619a',ntraciso,klon,klev) + endif + if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then + call iso_verif_O18_aberrant_enc_vect2D( & + & xt_seri,q_seri, & + & 'physiq 3619b',ntraciso,klon,klev) + endif #endif #ifdef ISOVERIF @@ -3906,5 +3908,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'physiq 3425' + write(*,*) 'physiq 3425: apres convection' if (iso_HDO.gt.0) then call iso_verif_aberrant_enc_vect2D( & @@ -3915,5 +3917,5 @@ call iso_verif_O18_aberrant_enc_vect2D( & & xt_seri,q_seri, & - & 'physiq 3691a',ntraciso,klon,klev) + & 'physiq 3691b',ntraciso,klon,klev) endif #ifdef ISOTRAC Index: LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90 (revision 4033) @@ -128,5 +128,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'surf_land_bucket 152' + !write(*,*) 'surf_land_bucket 152' do i=1,knon if (iso_eau.gt.0) then @@ -146,5 +146,5 @@ enddo !do ixt=1,niso enddo !do i=1,knon - write(*,*) 'surf_land_bucket 152' + !write(*,*) 'surf_land_bucket 152' #endif #endif @@ -211,5 +211,5 @@ ! verif #ifdef ISOVERIF - write(*,*) 'surf_land_bucket 211' + !write(*,*) 'surf_land_bucket 211' do i=1,knon if (iso_eau.gt.0) then @@ -236,5 +236,5 @@ enddo !do i=1,knon #ifdef ISOVERIF - write(*,*) 'surf_land_bucket 235' + ! write(*,*) 'surf_land_bucket 235' do i=1,knon if (iso_eau.gt.0) then @@ -243,6 +243,4 @@ endif enddo !do i=1,knon - write(*,*) 'snow_prec(1)=',snow_prec(1) - write(*,*) 'xtsnow(:,1)=',xtsnow(:,1) #endif #endif @@ -261,7 +259,4 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'surf_land_bucket 258' - write(*,*) 'snow_prec(1)=',snow_prec(1) - write(*,*) 'xtsnow(:,1)=',xtsnow(:,1) do i=1,knon do ixt=1,niso @@ -271,5 +266,5 @@ #endif #ifdef ISOVERIF - write(*,*) 'surf_land_bucket 235' + !write(*,*) 'surf_land_bucket 235' do i=1,knon if (iso_eau.gt.0) then @@ -297,6 +292,6 @@ endif !if ((iso_O17.gt.0).and.(iso_O18.gt.0)) then enddo !do i=1,knon - write(*,*) 'surf_land_mod 291' - write(*,*) 'snow_evap_diag(1)=',snow_evap_diag(1) + !write(*,*) 'surf_land_mod 291' + !write(*,*) 'snow_evap_diag(1)=',snow_evap_diag(1) #endif call calcul_iso_surf_ter_vectall(klon,knon, & Index: LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90 (revision 4033) @@ -166,5 +166,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'surf_land_mod 162' +! write(*,*) 'surf_land_mod 162' do i=1,knon if (iso_eau.gt.0) then @@ -179,5 +179,5 @@ #endif #ifdef ISOVERIF - write(*,*) 'surf_land 169: ok_veget=',ok_veget +! write(*,*) 'surf_land 169: ok_veget=',ok_veget do i=1,knon do ixt=1,ntraciso @@ -262,5 +262,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'surf_land 247' +! write(*,*) 'surf_land 247' call iso_verif_egalite_vect1D( & & xtsnow,snow,'surf_land_mod 207',niso,klon) @@ -306,5 +306,5 @@ #ifdef ISO #ifdef ISOVERIF - write(*,*) 'surf_land 237: sortie' +! write(*,*) 'surf_land 237: sortie' do i=1,knon if (iso_eau.gt.0) then Index: LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90 (revision 4032) +++ LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90 (revision 4033) @@ -375,5 +375,5 @@ ! verif #ifdef ISOVERIF - write(*,*) 'surf_land_ice 1499' + !write(*,*) 'surf_land_ice 1499' do i=1,knon if (iso_eau.gt.0) then