Ignore:
Timestamp:
05/28/19 22:44:58 (5 years ago)
Author:
jisesh
Message:

devel: moved DYSL into compute_caldyn_slow_NH.F90 and compute_caldyn_Coriolis.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r858 r876  
    127127END_BLOCK 
    128128 
    129 KERNEL(coriolis) 
    130 ! 
    131   DO iq=1,nqdyn 
    132     FORALL_CELLS_EXT() 
    133       ON_EDGES 
    134         Ftheta(EDGE) = .5*(theta(CELL1,iq)+theta(CELL2,iq))*hflux(EDGE) 
    135       END_BLOCK 
    136     END_BLOCK 
    137     FORALL_CELLS() 
    138       ON_PRIMAL 
    139         divF=0. 
    140         FORALL_EDGES 
    141           divF = divF + Ftheta(EDGE)*SIGN 
    142         END_BLOCK 
    143         dtheta_rhodz(CELL,iq) = -divF / AI 
    144       END_BLOCK 
    145     END_BLOCK 
    146   END DO ! iq 
    147 ! 
    148   FORALL_CELLS() 
    149     ON_PRIMAL 
    150       divF=0. 
    151       FORALL_EDGES 
    152         divF = divF + hflux(EDGE)*SIGN 
    153       END_BLOCK 
    154       convm(CELL) = -divF / AI 
    155     END_BLOCK 
    156   END_BLOCK 
    157 ! 
    158   FORALL_CELLS()  
    159     ON_EDGES 
    160       du_trisk=0. 
    161       FORALL_TRISK 
    162         du_trisk = du_trisk + WEE*hflux(EDGE_TRISK)*(qu(EDGE)+qu(EDGE_TRISK)) 
    163       END_BLOCK 
    164       du(EDGE) = du(EDGE) + .5*du_trisk 
    165     END_BLOCK 
    166   END_BLOCK 
    167  
    168 END_BLOCK 
Note: See TracChangeset for help on using the changeset viewer.