Ignore:
Timestamp:
05/29/19 11:16:16 (5 years ago)
Author:
jisesh
Message:

devel: moved DYSL into compute_caldyn_slow_hydro.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/devel/Python/src/kernels_caldyn_hydro.jin

    r836 r877  
    9393END_BLOCK 
    9494 
    95 KERNEL(caldyn_slow_hydro) 
    96   FORALL_CELLS_EXT() 
    97     ON_EDGES 
    98       uu = .5*(rhodz(CELL1)+rhodz(CELL2))*u(EDGE) 
    99       hflux(EDGE) = uu*LE_DE 
    100     END_BLOCK 
    101   END_BLOCK 
    102  
    103   FORALL_CELLS() 
    104     ON_PRIMAL 
    105       ke=0.d0 
    106       FORALL_EDGES 
    107         ke = ke + LE_DE*u(EDGE)**2 
    108       END_BLOCK 
    109       BERNI(CELL)=ke*(.25/AI) 
    110     END_BLOCK 
    111   END_BLOCK 
    112   IF(zero) THEN 
    113     FORALL_CELLS() 
    114       ON_EDGES 
    115         du(EDGE) = SIGN*(berni(CELL1)-berni(CELL2)) ! minus gradient 
    116       END_BLOCK 
    117     END_BLOCK 
    118   ELSE 
    119     FORALL_CELLS() 
    120       ON_EDGES 
    121         du(EDGE) = du(EDGE) + SIGN*(berni(CELL1)-berni(CELL2)) ! minus gradient 
    122       END_BLOCK 
    123     END_BLOCK 
    124   END IF 
    125  
    126 END_BLOCK 
Note: See TracChangeset for help on using the changeset viewer.