Ignore:
Timestamp:
05/05/19 14:58:31 (5 years ago)
Author:
dubos
Message:

devel : moved DYSL for compute_pvort_only into compute_pvort_only.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/devel/src/dynamics/compute_pvort_only.F90

    r840 r852  
    1717  END SUBROUTINE check_interface 
    1818   
     19#if BEGIN_DYSL 
     20 
     21KERNEL(pvort_only) 
     22  FORALL_CELLS_EXT() 
     23    ON_DUAL 
     24      etav = 0.d0 
     25      FORALL_EDGES 
     26         etav = etav + SIGN*u(EDGE) 
     27      END_BLOCK 
     28      hv=0. 
     29      FORALL_VERTICES 
     30        hv = hv + RIV2*rhodz(VERTEX) 
     31      END_BLOCK 
     32      qv(DUAL_CELL) = (etav + FV*AV )/(hv*AV) 
     33    END_BLOCK 
     34  END_BLOCK 
     35 
     36  FORALL_CELLS() 
     37    ON_EDGES 
     38      qu(EDGE)=0.5d0*(qv(VERTEX1)+qv(VERTEX2)) 
     39    END_BLOCK 
     40  END_BLOCK 
     41END_BLOCK 
     42 
     43#endif END_DYSL 
     44 
    1945  SUBROUTINE compute_pvort_only_unst(u,rhodz,qu,qv, hv_) 
    2046    USE ISO_C_BINDING, only : C_DOUBLE, C_FLOAT 
    … …  
    92118 
    93119END MODULE compute_pvort_only_mod 
     120 
Note: See TracChangeset for help on using the changeset viewer.