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/Python/src/kernels_caldyn_hevi.jin

    r837 r852  
    178178END_BLOCK 
    179179 
    180 KERNEL(pvort_only) 
    181   FORALL_CELLS_EXT() 
    182     ON_DUAL 
    183       etav = 0.d0 
    184       FORALL_EDGES 
    185          etav = etav + SIGN*u(EDGE) 
    186       END_BLOCK 
    187       hv=0. 
    188       FORALL_VERTICES 
    189         hv = hv + RIV2*rhodz(VERTEX) 
    190       END_BLOCK 
    191       qv(DUAL_CELL) = (etav + FV*AV )/(hv*AV) 
    192     END_BLOCK 
    193   END_BLOCK 
    194  
    195   FORALL_CELLS() 
    196     ON_EDGES 
    197       qu(EDGE)=0.5d0*(qv(VERTEX1)+qv(VERTEX2)) 
    198     END_BLOCK 
    199   END_BLOCK 
    200 END_BLOCK 
    201  
    202180KERNEL(coriolis) 
    203181! 
Note: See TracChangeset for help on using the changeset viewer.