Ignore:
Timestamp:
12/30/17 02:00:38 (7 years ago)
Author:
dubos
Message:

devel/unstructured : updated kernels

File:
1 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/devel/src/kernels_unst/caldyn_vert.k90

    r624 r658  
    44      !$OMP DO SCHEDULE(STATIC) 
    55      DO ij = 1, primal_num 
     6         !DIR$ SIMD 
    67         DO l = 2, llm 
    78            dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) + 0.5*(theta(l,ij,iq)+theta(l-1,ij,iq))*wflux(l,ij) 
     
    1112      !$OMP DO SCHEDULE(STATIC) 
    1213      DO ij = 1, primal_num 
     14         !DIR$ SIMD 
    1315         DO l = 1, llm-1 
    1416            dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) - 0.5*(theta(l,ij,iq)+theta(l+1,ij,iq))*wflux(l+1,ij) 
     
    2325         ij_left = left(edge) 
    2426         ij_right = right(edge) 
     27         !DIR$ SIMD 
    2528         DO l = 2, llm 
    2629            wwuu(l,edge) = .25*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)+u(l-1,edge)) ! Fu 
     
    3437         ij_left = left(edge) 
    3538         ij_right = right(edge) 
     39         !DIR$ SIMD 
    3640         DO l = 1, llm 
    3741            dFu_deta = wwuu(l+1,edge)-wwuu(l,edge) ! d/deta (F*u) 
     
    4650         ij_left = left(edge) 
    4751         ij_right = right(edge) 
     52         !DIR$ SIMD 
    4853         DO l = 2, llm 
    4954            wwuu(l,edge) = .5*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)-u(l-1,edge)) 
     
    5762         ij_left = left(edge) 
    5863         ij_right = right(edge) 
     64         !DIR$ SIMD 
    5965         DO l = 1, llm 
    6066            du(l,edge) = du(l,edge) - (wwuu(l,edge)+wwuu(l+1,edge)) / (rhodz(l,ij_left)+rhodz(l,ij_right)) 
Note: See TracChangeset for help on using the changeset viewer.