Changeset 1844 for LMDZ5/trunk


Ignore:
Timestamp:
Aug 22, 2013, 11:58:56 AM (11 years ago)
Author:
yann meurdesoif
Message:

Bug fix in transport for planetar version.
=> advection was not done

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dpar/advtrac_p.F90

    r1823 r1844  
    460460
    461461        CALL qminimum_p( q, 2, finmasse )
     462     endif ! of if (planet_type=="earth")
    462463
    463464        !------------------------------------------------------------------
     
    466467        !          iadvtr=0
    467468
    468         !$OMP MASTER
    469         call VTe(VTadvection)
    470         call stop_timer(timer_vanleer)
    471         call VTb(VThallo)
    472         !$OMP END MASTER
    473 
    474         do j=1,nqtot
    475            call Register_SwapFieldHallo(q(1,1,j),q(1,1,j),ip1jmp1,llm, &
    476                 jj_nb_caldyn,0,0,Request_vanleer)
    477         enddo
    478 
    479         call Register_SwapFieldHallo(flxw,flxw,ip1jmp1,llm, &
     469      !$OMP MASTER
     470      call VTe(VTadvection)
     471      call stop_timer(timer_vanleer)
     472      call VTb(VThallo)
     473     !$OMP END MASTER
     474
     475      do j=1,nqtot
     476        call Register_SwapFieldHallo(q(1,1,j),q(1,1,j),ip1jmp1,llm, &
    480477             jj_nb_caldyn,0,0,Request_vanleer)
    481 
    482         call SendRequest(Request_vanleer)
    483         !$OMP BARRIER
    484         call WaitRequest(Request_vanleer)     
    485 
    486         !$OMP BARRIER
    487         !$OMP MASTER
    488         call SetDistrib(jj_nb_caldyn)
    489         call VTe(VThallo)
    490         call resume_timer(timer_caldyn)
     478      enddo
     479
     480      call Register_SwapFieldHallo(flxw,flxw,ip1jmp1,llm, &
     481             jj_nb_caldyn,0,0,Request_vanleer)
     482
     483      call SendRequest(Request_vanleer)
     484      !$OMP BARRIER
     485      call WaitRequest(Request_vanleer)     
     486
     487 !$OMP BARRIER
     488 !$OMP MASTER
     489      call SetDistrib(jj_nb_caldyn)
     490      call VTe(VThallo)
     491      call resume_timer(timer_caldyn)
    491492 !$OMP END MASTER
    492493 !$OMP BARRIER 
    493494        iadvtr=0
    494      endif ! of if (planet_type=="earth")
    495495  ENDIF ! if iadvtr.EQ.iapp_tracvl
    496496
Note: See TracChangeset for help on using the changeset viewer.