Index: LMDZ6/trunk/libf/dyn3d/dynetat0.f90 =================================================================== --- LMDZ6/trunk/libf/dyn3d/dynetat0.f90 (revision 4051) +++ LMDZ6/trunk/libf/dyn3d/dynetat0.f90 (revision 4052) @@ -137,12 +137,11 @@ ! Rayleigh distillation law. iName = tracers(iq)%iso_iName + IF(.NOT.ok_isotopes .OR. iName<=0) CYCLE iZone = tracers(iq)%iso_iZone iPhase= tracers(iq)%iso_iPhase iqParent = tracers(iq)%iqParent - IF(ok_isotopes .AND. iName>0) THEN - IF(iZone==0) q(:,:,:,iq) = q(:,:,:,iqParent)*tnat(iName) & - *(q(:,:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.) - IF(iZone==1) q(:,:,:,iq) = q(:,:,:,iqiso(iso_indnum(iq),iPhase)) - END IF + IF(iZone==0) q(:,:,:,iq) = q(:,:,:,iqParent)*tnat(iName) & + *(q(:,:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.) + IF(iZone==1) q(:,:,:,iq) = q(:,:,:,iqiso(iso_indnum(iq),iPhase)) END DO Index: LMDZ6/trunk/libf/dyn3d/iniacademic.F90 =================================================================== --- LMDZ6/trunk/libf/dyn3d/iniacademic.F90 (revision 4051) +++ LMDZ6/trunk/libf/dyn3d/iniacademic.F90 (revision 4052) @@ -278,4 +278,6 @@ do iq=1,nqtot q(:,:,iq)=0. +! IF(tracers(iq)%name == 'H2O'//phases_sep//'g') q(:,:,iq)=1.e-10 +! IF(tracers(iq)%name == 'H2O'//phases_sep//'l') q(:,:,iq)=1.e-15 IF(tracers(iq)%name == 'H2Ov') q(:,:,iq)=1.e-10 IF(tracers(iq)%name == 'H2Ol') q(:,:,iq)=1.e-15 @@ -284,13 +286,11 @@ ! distill de Rayleigh très simplifiée iName = tracers(iq)%iso_iName + if (.NOT.ok_isotopes .OR. iName <= 0) CYCLE iZone = tracers(iq)%iso_iZone iPhase= tracers(iq)%iso_iPhase iqParent = tracers(iq)%iqParent - if (ok_isotopes .AND. iName > 0) then - if (iZone == 0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName) & - *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1) - if (iZone == 1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase)) - endif !if (ok_isotopes) then - + if (iZone == 0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName) & + *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1) + if (iZone == 1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase)) enddo else Index: LMDZ6/trunk/libf/dyn3d/vlsplt.F =================================================================== --- LMDZ6/trunk/libf/dyn3d/vlsplt.F (revision 4051) +++ LMDZ6/trunk/libf/dyn3d/vlsplt.F (revision 4052) @@ -464,5 +464,5 @@ enddo enddo - do ifils=1,tracers(iq)%nqDescen + do ifils=1,tracers(iq)%nqChilds iq2=tracers(iq)%iqDescen(ifils) call vlx(Ratio,pente_max,masseq,u_mq,iq2) @@ -1015,5 +1015,5 @@ enddo - do ifils=1,tracers(iq)%nqDescen + do ifils=1,tracers(iq)%nqChilds iq2=tracers(iq)%iqDescen(ifils) call vlz(Ratio,pente_max,masseq,wq,iq2) Index: LMDZ6/trunk/libf/dyn3d/vlspltqs.F =================================================================== --- LMDZ6/trunk/libf/dyn3d/vlspltqs.F (revision 4051) +++ LMDZ6/trunk/libf/dyn3d/vlspltqs.F (revision 4052) @@ -786,5 +786,5 @@ ! CRisi: appel récursif de l'advection sur les fils. ! Il faut faire ça avant d'avoir mis à jour q et masse - !write(*,*) 'vlyqs 689: iq,nqDesc(iq)=',iq,tracers(iq)%nqDescen + !write(*,*) 'vlyqs 689: iq,nqChilds(iq)=',iq,tracers(iq)%nqChilds do ifils=1,tracers(iq)%nqDescen Index: LMDZ6/trunk/libf/dyn3d_common/infotrac.F90 =================================================================== --- LMDZ6/trunk/libf/dyn3d_common/infotrac.F90 (revision 4051) +++ LMDZ6/trunk/libf/dyn3d_common/infotrac.F90 (revision 4052) @@ -734,4 +734,7 @@ enddo !do iq=1,nqtot tracers(:)%gen0Name = ancestor(tracers) !--- Name of the first generation ancestor + tracers(:)%isAdvected = tracers(:)%iadv > 0 +! tracers(:)%isH2Ofamily = delPhase(tracers(:)%gen0Name) == 'H2O' + tracers(:)%isH2Ofamily = [(tracers(iq)%gen0Name(1:3) == 'H2O', iq=1, nqtot)] Index: LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 =================================================================== --- LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 (revision 4051) +++ LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90 (revision 4052) @@ -174,12 +174,11 @@ ! Rayleigh distillation las. iName = tracers(iq)%iso_iName + IF(.NOT.ok_isotopes .OR. iName <= 0) CYCLE iZone = tracers(iq)%iso_iZone iPhase= tracers(iq)%iso_iPhase iqParent = tracers(iq)%iqParent - IF(ok_isotopes .AND. iName>0) THEN - IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName) & - & *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.) - IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase)) - END IF + IF(iZone==0) q(:,:,iq) = q(:,:,iqParent)*tnat(iName) & + & *(q(:,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1.) + IF(iZone==1) q(:,:,iq) = q(:,:,iqiso(iso_indnum(iq),iPhase)) END DO DEALLOCATE(q_glo) Index: LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90 =================================================================== --- LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90 (revision 4051) +++ LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90 (revision 4052) @@ -284,4 +284,6 @@ do iq=1,nqtot q(ijb_u:ije_u,:,iq)=0. +! IF(tracers(iq)%name == 'H2O'//Phases_sep//'g') q(ijb_u:ije_u,:,iq)=1.e-10 +! IF(tracers(iq)%name == 'H2O'//phases_sep//'l') q(ijb_u:ije_u,:,iq)=1.e-15 if (tracers(iq)%name == 'H2Ov') q(ijb_u:ije_u,:,iq)=1.e-10 if (tracers(iq)%name == 'H2Ol') q(ijb_u:ije_u,:,iq)=1.e-15 @@ -290,13 +292,11 @@ ! distill de Rayleigh très simplifiée iName = tracers(iq)%iso_iName + if (.NOT.ok_isotopes .OR. iName <= 0) CYCLE iZone = tracers(iq)%iso_iZone iPhase= tracers(iq)%iso_iPhase iqParent = tracers(iq)%iqParent - if (ok_isotopes .AND. iName > 0) then - if (iZone == 0) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat(iName) & - *(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1) - if (iZone == 1) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqiso(iso_indnum(iq),iPhase)) - endif !if (ok_isotopes) then - + if (iZone == 0) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat(iName) & + *(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal(iName)-1) + if (iZone == 1) q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqiso(iso_indnum(iq),iPhase)) enddo else Index: LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F =================================================================== --- LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F (revision 4051) +++ LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F (revision 4052) @@ -726,5 +726,5 @@ ! CRisi: appel récursif de l'advection sur les fils. ! Il faut faire ça avant d'avoir mis à jour q et masse - !write(*,*) 'vly 689: iq,nqDesc(iq)=',iq,tracers(iq)%nqDescen + !write(*,*) 'vly 689: iq,nqChilds(iq)=',iq,tracers(iq)%nqChilds ijb=ij_begin-2*iip1 Index: LMDZ6/trunk/libf/dyn3dmem/vlspltgen_loc.F =================================================================== --- LMDZ6/trunk/libf/dyn3dmem/vlspltgen_loc.F (revision 4051) +++ LMDZ6/trunk/libf/dyn3dmem/vlspltgen_loc.F (revision 4052) @@ -346,5 +346,5 @@ SELECT CASE(tracers(iq)%iadv) CASE(0); CYCLE - CASE(10); call vly_loc(zq,pente_max,zm,mv,iq) + CASE(10); call vly_loc(zq,pente_max,zm,mv, iq) CASE(14); call vlyqs_loc(zq,pente_max,zm,mv,qsat,iq) CASE DEFAULT Index: LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F =================================================================== --- LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F (revision 4051) +++ LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F (revision 4052) @@ -264,5 +264,5 @@ c$OMP END DO NOWAIT -cym ATTENTION ICI en OpenMP reduction pas forcement nécessaire +cym ATTENTION ICI en OpenMP reduction pas forcement necessaire cym IF(n0.gt.1) THEN cym IF(n0.gt.0) THEN @@ -335,8 +335,8 @@ c$OMP END DO NOWAIT -! CRisi: appel récursif de l'advection sur les fils. -! Il faut faire ça avant d'avoir mis à jour q et masse - !write(*,*) 'vlspltqs 336: iq,ijb_x,iqDescen(iq)=', -! & iq,ijb_x,tracers(iq)%iqDescen +! CRisi: appel recursif de l'advection sur les fils. +! Il faut faire ca avant d'avoir mis a jour q et masse + !write(*,*) 'vlspltqs 336: iq,ijb_x,nqChilds(iq)=', +! & iq,ijb_x,tracers(iq)%nqChilds do ifils=1,tracers(iq)%nqDescen @@ -356,5 +356,5 @@ c$OMP END DO NOWAIT enddo - do ifils=1,tracers(iq)%nqDescen + do ifils=1,tracers(iq)%nqChilds iq2=tracers(iq)%iqDescen(ifils) !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2 @@ -727,7 +727,7 @@ c$OMP END DO NOWAIT -! CRisi: appel récursif de l'advection sur les fils. -! Il faut faire ça avant d'avoir mis à jour q et masse - !write(*,*) 'vlyqs 689: iq,iqDescen(iq)=',iq,tracers(iq)%iqDescen +! CRisi: appel recursif de l'advection sur les fils. +! Il faut faire ca avant d'avoir mis a jour q et masse + !write(*,*) 'vlyqs 689: iq,nqChilds(iq)=',iq,tracers(iq)%nqChilds ijb=ij_begin-2*iip1 @@ -748,5 +748,5 @@ DO l=1,llm ! modif des bornes: CRisi 16 nov 2020 - ! d'abord masse avec bornes corrigées + ! d'abord masse avec bornes corrigees DO ij=ijbm,ijem !MVals: veiller a ce qu'on n'ait pas de denominateur nul @@ -767,5 +767,5 @@ c$OMP END DO NOWAIT enddo - do ifils=1,tracers(iq)%nqDescen + do ifils=1,tracers(iq)%nqChilds iq2=tracers(iq)%iqDescen(ifils) !write(lunout,*) 'vly: appel recursiv vly iq2=',iq2 Index: LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 =================================================================== --- LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (revision 4051) +++ LMDZ6/trunk/libf/phylmd/infotrac_phy.F90 (revision 4052) @@ -145,5 +145,4 @@ CHARACTER(LEN=30) :: modname="init_infotrac_phy" - INTEGER :: iq nqtot=nqtot_ @@ -153,7 +152,4 @@ nqtottr=nqtottr_ ALLOCATE(tracers(nqtot)); tracers(:) = tracers_(:) - tracers(:)%isAdvected = tracers(:)%iadv > 0 -! tracers(:)%isH2Ofamily = delPhase(tracers(:)%gen0Name) == 'H2O' - tracers(:)%isH2Ofamily = [(tracers(iq)%gen0Name(1:3) == 'H2O', iq=1, nqtot)] #ifdef CPP_StratAer nbtr_bin=nbtr_bin_ @@ -209,8 +205,8 @@ ALLOCATE(index_trac(ntraceurs_zone,niso)) index_trac(:,:)=index_trac_(:,:) - - ALLOCATE(itr_indice(nqtottr)) - itr_indice(:)=itr_indice_(:) ENDIF ! of IF(ok_isotopes) + + ALLOCATE(itr_indice(nqtottr)) + itr_indice(:)=itr_indice_(:) END SUBROUTINE init_infotrac_phy Index: LMDZ6/trunk/libf/phylmdiso/infotrac_phy.F90 =================================================================== --- LMDZ6/trunk/libf/phylmdiso/infotrac_phy.F90 (revision 4051) +++ LMDZ6/trunk/libf/phylmdiso/infotrac_phy.F90 (revision 4052) @@ -205,9 +205,8 @@ ALLOCATE(index_trac(ntraceurs_zone,niso)) index_trac(:,:)=index_trac_(:,:) - ENDIF ! of IF(ok_isotopes) - ALLOCATE(itr_indice(nqtottr)) - itr_indice(:)=itr_indice_(:) + ALLOCATE(itr_indice(nqtottr)) + itr_indice(:)=itr_indice_(:) write(*,*) 'infotrac_phy 207: nqtottr=',nqtottr