Ignore:
Timestamp:
Jun 14, 2015, 9:13:32 PM (9 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes -r2237:2291 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/dyn3dmem/vlspltgen_loc.F

    r1910 r2298  
    2727      USE Write_Field_loc
    2828      USE VAMPIR
    29       USE infotrac, ONLY : nqtot
     29      ! CRisi: on rajoute variables utiles d'infotrac 
     30      USE infotrac, ONLY : nqtot,nqperes,nqdesc,nqfils,iqfils,
     31     &    ok_iso_verif
    3032      USE vlspltgen_mod
    3133      IMPLICIT NONE
    … …  
    6466      REAL ptarg,pdelarg,foeew,zdelta
    6567      REAL tempe(ijb_u:ije_u)
    66       INTEGER ijb,ije,iq
     68      INTEGER ijb,ije,iq,iq2,ifils
    6769      LOGICAL, SAVE :: firstcall=.TRUE.
    6870!$OMP THREADPRIVATE(firstcall)
    … …  
    150152      ije=ij_end
    151153
     154      DO iq=1,nqtot
    152155c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
    153156      DO l=1,llm
    154157         DO ij=ijb,ije
    155             mw(ij,l)=w(ij,l) * zzw
     158            mw(ij,l,iq)=w(ij,l) * zzw
    156159         ENDDO
    157160      ENDDO
    158161c$OMP END DO NOWAIT
    159 
     162      ENDDO
     163
     164      DO iq=1,nqtot 
    160165c$OMP MASTER
    161166      DO ij=ijb,ije
    162          mw(ij,llm+1)=0.
     167         mw(ij,llm+1,iq)=0.
    163168      ENDDO
    164169c$OMP END MASTER
     170      ENDDO
    165171
    166172c      CALL SCOPY(ijp1llm,q,1,zq,1)
    … …  
    170176       ije=ij_end
    171177
    172       DO iq=1,nqtot
     178      DO iq=1,nqtot       
    173179c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    174180        DO l=1,llm
    … …  
    179185      ENDDO
    180186
    181 #ifdef DEBUG_IO   
     187#ifdef DEBUG_IO     
    182188       CALL WriteField_u('mu',mu)
    183189       CALL WriteField_v('mv',mv)
    … …  
    186192#endif
    187193
     194      ! verif temporaire
     195      ijb=ij_begin
     196      ije=ij_end 
     197      if (ok_iso_verif) then
     198        call check_isotopes(zq,ijb,ije,'vlspltgen_loc 191')
     199      endif !if (ok_iso_verif) then   
     200
    188201c$OMP BARRIER           
    189       DO iq=1,nqtot
    190 
     202!      DO iq=1,nqtot
     203      DO iq=1,nqperes ! CRisi: on ne boucle que sur les pères= ceux qui sont transportés directement par l'air
     204       !write(*,*) 'vlspltgen 192: iq,iadv=',iq,iadv(iq)
     205#ifdef DEBUG_IO   
     206       CALL WriteField_u('zq',zq(:,:,iq))
     207       CALL WriteField_u('zm',zm(:,:,iq))
     208#endif
     209        if(iadv(iq) == 0) then
     210       
     211          cycle
     212       
     213        else if (iadv(iq)==10) then
     214
     215#ifdef _ADV_HALO       
     216! CRisi: on ajoute les nombres de fils et tableaux des fils
     217! On suppose qu'on ne peut advecter les fils que par le schéma 10. 
     218          call vlx_loc(zq,pente_max,zm,mu,
     219     &               ij_begin,ij_begin+2*iip1-1,iq)
     220          call vlx_loc(zq,pente_max,zm,mu,
     221     &               ij_end-2*iip1+1,ij_end,iq)
     222#else
     223          call vlx_loc(zq,pente_max,zm,mu,
     224     &               ij_begin,ij_end,iq)
     225#endif
     226
     227c$OMP MASTER
     228          call VTb(VTHallo)
     229c$OMP END MASTER
     230          call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest1)
     231          call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest1)
     232! CRisi
     233          do ifils=1,nqdesc(iq)
     234            iq2=iqfils(ifils,iq)
     235            call Register_Hallo_u(zq(:,:,iq2),llm,2,2,2,2,MyRequest1)
     236            call Register_Hallo_u(zm(:,:,iq2),llm,1,1,1,1,MyRequest1)
     237          enddo
     238
     239c$OMP MASTER
     240          call VTe(VTHallo)
     241c$OMP END MASTER
     242        else if (iadv(iq)==14) then
     243
     244#ifdef _ADV_HALO           
     245          call vlxqs_loc(zq,pente_max,zm,mu,
     246     &                   qsat,ij_begin,ij_begin+2*iip1-1,iq)
     247          call vlxqs_loc(zq,pente_max,zm,mu,
     248     &                   qsat,ij_end-2*iip1+1,ij_end,iq)
     249#else
     250          call vlxqs_loc(zq,pente_max,zm,mu,
     251     &                   qsat,ij_begin,ij_end,iq)
     252#endif
     253
     254c$OMP MASTER
     255          call VTb(VTHallo)
     256c$OMP END MASTER
     257
     258          call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest1)
     259          call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest1)
     260          do ifils=1,nqdesc(iq)
     261            iq2=iqfils(ifils,iq)
     262            call Register_Hallo_u(zq(:,:,iq2),llm,2,2,2,2,MyRequest1)
     263            call Register_Hallo_u(zm(:,:,iq2),llm,1,1,1,1,MyRequest1)
     264          enddo
     265
     266c$OMP MASTER
     267          call VTe(VTHallo)
     268c$OMP END MASTER
     269        else
     270       
     271          stop 'vlspltgen_p : schema non parallelise'
     272     
     273        endif
     274     
     275      enddo !DO iq=1,nqperes
     276     
     277     
     278c$OMP BARRIER     
     279c$OMP MASTER     
     280      call VTb(VTHallo)
     281c$OMP END MASTER
     282
     283      call SendRequest(MyRequest1)
     284
     285c$OMP MASTER
     286      call VTe(VTHallo)
     287c$OMP END MASTER       
     288c$OMP BARRIER
     289
     290      ! verif temporaire
     291      ijb=ij_begin-2*iip1
     292      ije=ij_end+2*iip1 
     293      if (pole_nord) ijb=ij_begin
     294      if (pole_sud)  ije=ij_end 
     295      if (ok_iso_verif) then
     296           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 280')
     297      endif !if (ok_iso_verif) then
     298
     299      do iq=1,nqperes
     300        !write(*,*) 'vlspltgen 279: iq=',iq
     301
     302        if(iadv(iq) == 0) then
     303       
     304          cycle
     305       
     306        else if (iadv(iq)==10) then
     307
     308#ifdef _ADV_HALLO
     309          call vlx_loc(zq,pente_max,zm,mu,
     310     &                 ij_begin+2*iip1,ij_end-2*iip1,iq)
     311#endif       
     312        else if (iadv(iq)==14) then
     313#ifdef _ADV_HALLO
     314          call vlxqs_loc(zq,pente_max,zm,mu,
     315     &                    qsat,ij_begin+2*iip1,ij_end-2*iip1,iq)
     316#endif   
     317        else
     318       
     319          stop 'vlspltgen_p : schema non parallelise'
     320     
     321        endif
     322     
     323      enddo
     324c$OMP BARRIER     
     325c$OMP MASTER
     326      call VTb(VTHallo)
     327c$OMP END MASTER
     328
     329!      call WaitRecvRequest(MyRequest1)
     330!      call WaitSendRequest(MyRequest1)
     331c$OMP BARRIER
     332       call WaitRequest(MyRequest1)
     333
     334
     335c$OMP MASTER
     336      call VTe(VTHallo)
     337c$OMP END MASTER
     338c$OMP BARRIER
     339
     340     
     341      if (ok_iso_verif) then
     342           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 326')
     343      endif !if (ok_iso_verif) then       
     344      if (ok_iso_verif) then
     345           ijb=ij_begin-2*iip1
     346           ije=ij_end+2*iip1
     347           if (pole_nord) ijb=ij_begin
     348           if (pole_sud)  ije=ij_end
     349           call check_isotopes(zq,ijb,ije,'vlspltgen_loc 336')
     350      endif !if (ok_iso_verif) then 
     351
     352      do iq=1,nqperes
     353       !write(*,*) 'vlspltgen 321: iq=',iq
    191354#ifdef DEBUG_IO   
    192355       CALL WriteField_u('zq',zq(:,:,iq))
    193356       CALL WriteField_u('zm',zm(:,:,iq))
    194357#endif
     358
    195359        if(iadv(iq) == 0) then
    196360       
    … …  
    198362       
    199363        else if (iadv(iq)==10) then
    200 
    201 #ifdef _ADV_HALO       
    202           call vlx_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    203      &               ij_begin,ij_begin+2*iip1-1)
    204           call vlx_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    205      &               ij_end-2*iip1+1,ij_end)
    206 #else
    207           call vlx_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    208      &               ij_begin,ij_end)
    209 #endif
    210 
    211 c$OMP MASTER
    212           call VTb(VTHallo)
    213 c$OMP END MASTER
    214           call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest1)
    215           call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest1)
    216 
    217 c$OMP MASTER
    218           call VTe(VTHallo)
    219 c$OMP END MASTER
     364       
     365          call vly_loc(zq,pente_max,zm,mv,iq)
     366 
    220367        else if (iadv(iq)==14) then
    221 
    222 #ifdef _ADV_HALO           
    223           call vlxqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    224      &                   qsat,ij_begin,ij_begin+2*iip1-1)
    225           call vlxqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    226      &                   qsat,ij_end-2*iip1+1,ij_end)
    227 #else
    228 
    229           call vlxqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    230      &                   qsat,ij_begin,ij_end)
    231 #endif
    232 
    233 c$OMP MASTER
    234           call VTb(VTHallo)
    235 c$OMP END MASTER
    236 
    237           call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest1)
    238           call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest1)
    239 
    240 c$OMP MASTER
    241           call VTe(VTHallo)
    242 c$OMP END MASTER
     368     
     369          call vlyqs_loc(zq,pente_max,zm,mv,
     370     &                   qsat,iq)
     371 
    243372        else
    244373       
    … …  
    246375     
    247376        endif
    248      
    249       enddo
    250      
    251      
    252 c$OMP BARRIER     
    253 c$OMP MASTER     
    254       call VTb(VTHallo)
    255 c$OMP END MASTER
    256 
    257       call SendRequest(MyRequest1)
    258 
    259 c$OMP MASTER
    260       call VTe(VTHallo)
    261 c$OMP END MASTER       
    262 c$OMP BARRIER
    263       do iq=1,nqtot
    264 
    265         if(iadv(iq) == 0) then
    266        
    267           cycle
    268        
    269         else if (iadv(iq)==10) then
    270 
    271 #ifdef _ADV_HALLO
    272           call vlx_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    273      &                 ij_begin+2*iip1,ij_end-2*iip1)
    274 #endif       
    275         else if (iadv(iq)==14) then
    276 #ifdef _ADV_HALLO
    277           call vlxqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    278      &                    qsat,ij_begin+2*iip1,ij_end-2*iip1)
    279 #endif   
    280         else
    281        
    282           stop 'vlspltgen_p : schema non parallelise'
    283      
    284         endif
    285      
    286       enddo
    287 c$OMP BARRIER     
    288 c$OMP MASTER
    289       call VTb(VTHallo)
    290 c$OMP END MASTER
    291 
    292 !      call WaitRecvRequest(MyRequest1)
    293 !      call WaitSendRequest(MyRequest1)
    294 c$OMP BARRIER
    295        call WaitRequest(MyRequest1)
    296 
    297 
    298 c$OMP MASTER
    299       call VTe(VTHallo)
    300 c$OMP END MASTER
    301 c$OMP BARRIER
    302 
    303 
    304       do iq=1,nqtot
     377       
     378       enddo
     379
     380      if (ok_iso_verif) then
     381           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 357')
     382      endif !if (ok_iso_verif) then
     383
     384      do iq=1,nqperes
     385      !write(*,*) 'vlspltgen 349: iq=',iq
    305386#ifdef DEBUG_IO   
    306387       CALL WriteField_u('zq',zq(:,:,iq))
    307388       CALL WriteField_u('zm',zm(:,:,iq))
    308389#endif
     390        if(iadv(iq) == 0) then
     391         
     392          cycle
     393       
     394        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     395
     396c$OMP BARRIER       
     397#ifdef _ADV_HALLO
     398          call vlz_loc(zq,pente_max,zm,mw,
     399     &               ij_begin,ij_begin+2*iip1-1,iq)
     400          call vlz_loc(zq,pente_max,zm,mw,
     401     &               ij_end-2*iip1+1,ij_end,iq)
     402#else
     403          call vlz_loc(zq,pente_max,zm,mw,
     404     &               ij_begin,ij_end,iq)
     405#endif
     406c$OMP BARRIER
     407
     408c$OMP MASTER
     409          call VTb(VTHallo)
     410c$OMP END MASTER
     411
     412          call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest2)
     413          call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest2)
     414          ! CRisi
     415          do ifils=1,nqdesc(iq)
     416            iq2=iqfils(ifils,iq)
     417            call Register_Hallo_u(zq(:,:,iq2),llm,2,2,2,2,MyRequest2)
     418            call Register_Hallo_u(zm(:,:,iq2),llm,1,1,1,1,MyRequest2)
     419          enddo     
     420c$OMP MASTER
     421          call VTe(VTHallo)
     422c$OMP END MASTER       
     423c$OMP BARRIER
     424        else
     425       
     426          stop 'vlspltgen_p : schema non parallelise'
     427     
     428        endif
     429     
     430      enddo
     431c$OMP BARRIER     
     432
     433c$OMP MASTER       
     434      call VTb(VTHallo)
     435c$OMP END MASTER
     436
     437      call SendRequest(MyRequest2)
     438
     439c$OMP MASTER
     440      call VTe(VTHallo)
     441c$OMP END MASTER       
     442
     443
     444      if (ok_iso_verif) then
     445           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 429')
     446      endif !if (ok_iso_verif) then
     447
     448c$OMP BARRIER
     449      do iq=1,nqperes
     450      !write(*,*) 'vlspltgen 409: iq=',iq
     451
    309452        if(iadv(iq) == 0) then
    310        
     453         
    311454          cycle
    312455       
    313         else if (iadv(iq)==10) then
    314        
    315           call vly_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mv)
    316  
    317         else if (iadv(iq)==14) then
    318       
    319           call vlyqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mv,
    320      &                   qsat)
    321 
     456        else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
     457c$OMP BARRIER       
     458
     459#ifdef _ADV_HALLO
     460          call vlz_loc(zq,pente_max,zm,mw,
     461     &               ij_begin+2*iip1,ij_end-2*iip1,iq)
     462#endif
     463
     464c$OMP BARRIER       
    322465        else
    323466       
    … …  
    325468     
    326469        endif
    327        
    328        enddo
    329 
    330 
    331       do iq=1,nqtot
     470     
     471      enddo
     472      !write(*,*) 'vlspltgen_loc 476'
     473
     474c$OMP BARRIER
     475      !write(*,*) 'vlspltgen_loc 477'
     476c$OMP MASTER
     477      call VTb(VTHallo)
     478c$OMP END MASTER
     479
     480!      call WaitRecvRequest(MyRequest2)
     481!      call WaitSendRequest(MyRequest2)
     482c$OMP BARRIER
     483       CALL WaitRequest(MyRequest2)
     484
     485c$OMP MASTER
     486      call VTe(VTHallo)
     487c$OMP END MASTER
     488c$OMP BARRIER
     489
     490
     491      !write(*,*) 'vlspltgen_loc 494'
     492      if (ok_iso_verif) then
     493           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 461')
     494      endif !if (ok_iso_verif) then
     495
     496      do iq=1,nqperes
     497      !write(*,*) 'vlspltgen 449: iq=',iq
    332498#ifdef DEBUG_IO   
    333499       CALL WriteField_u('zq',zq(:,:,iq))
    334500       CALL WriteField_u('zm',zm(:,:,iq))
    335501#endif
     502        if(iadv(iq) == 0) then
     503       
     504          cycle
     505       
     506        else if (iadv(iq)==10) then
     507       
     508          call vly_loc(zq,pente_max,zm,mv,iq)
     509 
     510        else if (iadv(iq)==14) then
     511     
     512          call vlyqs_loc(zq,pente_max,zm,mv,
     513     &                   qsat,iq)
     514 
     515        else
     516       
     517          stop 'vlspltgen_p : schema non parallelise'
     518     
     519        endif
     520       
     521       enddo !do iq=1,nqperes
     522
     523      if (ok_iso_verif) then
     524           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 493')
     525      endif !if (ok_iso_verif) then
     526
     527      do iq=1,nqperes
     528      !write(*,*) 'vlspltgen 477: iq=',iq
     529#ifdef DEBUG_IO   
     530       CALL WriteField_u('zq',zq(:,:,iq))
     531       CALL WriteField_u('zm',zm(:,:,iq))
     532#endif
    336533        if(iadv(iq) == 0) then
    337534         
    338535          cycle
    339536       
    340         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    341 
    342 c$OMP BARRIER       
    343 #ifdef _ADV_HALLO
    344           call vlz_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mw,
    345      &               ij_begin,ij_begin+2*iip1-1)
    346           call vlz_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mw,
    347      &               ij_end-2*iip1+1,ij_end)
    348 #else
    349           call vlz_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mw,
    350      &               ij_begin,ij_end)
    351 #endif
    352 c$OMP BARRIER
    353 
    354 c$OMP MASTER
    355           call VTb(VTHallo)
    356 c$OMP END MASTER
    357 
    358           call Register_Hallo_u(zq(:,:,iq),llm,2,2,2,2,MyRequest2)
    359           call Register_Hallo_u(zm(:,:,iq),llm,1,1,1,1,MyRequest2)
    360 
    361 c$OMP MASTER
    362           call VTe(VTHallo)
    363 c$OMP END MASTER       
    364 c$OMP BARRIER
    365         else
    366        
    367           stop 'vlspltgen_p : schema non parallelise'
    368      
    369         endif
    370      
    371       enddo
    372 c$OMP BARRIER     
    373 
    374 c$OMP MASTER       
    375       call VTb(VTHallo)
    376 c$OMP END MASTER
    377 
    378       call SendRequest(MyRequest2)
    379 
    380 c$OMP MASTER
    381       call VTe(VTHallo)
    382 c$OMP END MASTER       
    383 
    384 c$OMP BARRIER
    385       do iq=1,nqtot
    386 
    387         if(iadv(iq) == 0) then
    388          
    389           cycle
    390        
    391         else if (iadv(iq)==10 .or. iadv(iq)==14 ) then
    392 c$OMP BARRIER       
    393 
    394 #ifdef _ADV_HALLO
    395           call vlz_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mw,
    396      &               ij_begin+2*iip1,ij_end-2*iip1)
    397 #endif
    398 
    399 c$OMP BARRIER       
    400         else
    401        
    402           stop 'vlspltgen_p : schema non parallelise'
    403      
    404         endif
    405      
    406       enddo
    407 
    408 c$OMP BARRIER
    409 c$OMP MASTER
    410       call VTb(VTHallo)
    411 c$OMP END MASTER
    412 
    413 !      call WaitRecvRequest(MyRequest2)
    414 !      call WaitSendRequest(MyRequest2)
    415 c$OMP BARRIER
    416        CALL WaitRequest(MyRequest2)
    417 
    418 c$OMP MASTER
    419       call VTe(VTHallo)
    420 c$OMP END MASTER
    421 c$OMP BARRIER
    422 
    423 
    424       do iq=1,nqtot
    425 #ifdef DEBUG_IO   
    426        CALL WriteField_u('zq',zq(:,:,iq))
    427        CALL WriteField_u('zm',zm(:,:,iq))
    428 #endif
    429         if(iadv(iq) == 0) then
    430        
    431           cycle
    432        
    433537        else if (iadv(iq)==10) then
    434538       
    435           call vly_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mv)
     539          call vlx_loc(zq,pente_max,zm,mu,
     540     &               ij_begin,ij_end,iq)
    436541 
    437542        else if (iadv(iq)==14) then
    438543     
    439           call vlyqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mv,
    440      &                   qsat)
     544          call vlxqs_loc(zq,pente_max,zm,mu,
     545     &                 qsat, ij_begin,ij_end,iq)
    441546 
    442547        else
    443548       
    444           stop 'vlspltgen_p : schema non parallelise'
     549          stop 'vlspltgen_p : schema non parallelise'
    445550     
    446551        endif
    447552       
    448        enddo
    449 
    450 
    451       do iq=1,nqtot
    452 #ifdef DEBUG_IO   
    453        CALL WriteField_u('zq',zq(:,:,iq))
    454        CALL WriteField_u('zm',zm(:,:,iq))
    455 #endif
    456         if(iadv(iq) == 0) then
    457          
    458           cycle
    459        
    460         else if (iadv(iq)==10) then
    461        
    462           call vlx_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    463      &               ij_begin,ij_end)
    464  
    465         else if (iadv(iq)==14) then
    466      
    467           call vlxqs_loc(zq(ijb_u,1,iq),pente_max,zm(ijb_u,1,iq),mu,
    468      &                 qsat, ij_begin,ij_end)
    469  
    470         else
    471        
    472           stop 'vlspltgen_p : schema non parallelise'
    473      
    474         endif
    475        
    476        enddo
    477 
     553       enddo !do iq=1,nqperes
     554
     555      !write(*,*) 'vlspltgen 550: apres derniere serie de call vlx'
     556      if (ok_iso_verif) then
     557           call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 521')
     558      endif !if (ok_iso_verif) then
    478559     
    479560      ijb=ij_begin
    480561      ije=ij_end
     562      !write(*,*) 'vlspltgen_loc 557'
    481563c$OMP BARRIER     
    482564
    483 
     565      !write(*,*) 'vlspltgen_loc 559' 
    484566      DO iq=1,nqtot
     567       !write(*,*) 'vlspltgen_loc 561, iq=',iq 
    485568#ifdef DEBUG_IO   
    486569       CALL WriteField_u('zq',zq(:,:,iq))
    … …  
    495578           ENDDO
    496579        ENDDO
    497 c$OMP END DO NOWAIT         
     580c$OMP END DO NOWAIT   
     581      !write(*,*) 'vlspltgen_loc 575'     
    498582
    499583c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    … …  
    504588        ENDDO
    505589c$OMP END DO NOWAIT 
    506 
    507       ENDDO
     590      !write(*,*) 'vlspltgen_loc 583' 
     591      ENDDO !DO iq=1,nqtot
    508592       
    509      
     593      if (ok_iso_verif) then
     594           call check_isotopes(q,ij_begin,ij_end,'vlspltgen_loc 557')
     595      endif !if (ok_iso_verif) then
     596
    510597c$OMP BARRIER
    511598
    … …  
    516603cc$OMP BARRIER
    517604
     605      !write(*,*) 'vlspltgen 597: sortie' 
    518606      RETURN
    519607      END
Note: See TracChangeset for help on using the changeset viewer.