Ignore:
Timestamp:
05/29/19 20:33:00 (5 years ago)
Author:
jisesh
Message:

devel: moved DYSL into compute_caldyn_solver.F90,compute_theta.F90 and compute_NH_geopot.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r876 r878  
    9696END_BLOCK 
    9797 
    98 KERNEL(theta) 
    99   IF(caldyn_eta==eta_mass) THEN ! Compute mass 
    100     ! FIXME : here mass_col is computed from rhodz 
    101     ! so that the DOFs are the same whatever caldyn_eta 
    102     ! in DYNAMICO mass_col is prognosed rather than rhodz 
    103     SEQUENCE_C1 
    104       PROLOGUE(0) 
    105         mass_col(HIDX(CELL))=0. 
    106       END_BLOCK 
    107       BODY('1,llm') 
    108           mass_col(HIDX(CELL)) = mass_col(HIDX(CELL)) + rhodz(CELL) 
    109       END_BLOCK 
    110     END_BLOCK 
    111     FORALL_CELLS_EXT() 
    112       ON_PRIMAL 
    113         ! FIXME : formula below (used in DYNAMICO) is for dak, dbk based on pressure rather than mass 
    114         !        m = mass_dak(CELL)+(mass_col(HIDX(CELL))*g+ptop)*mass_dbk(CELL) 
    115         !        rhodz(CELL) = m/g 
    116         rhodz(CELL) = mass_dak(CELL) + mass_col(HIDX(CELL))*mass_dbk(CELL) 
    117       END_BLOCK 
    118     END_BLOCK 
    119   END IF 
    120   DO iq=1,nqdyn 
    121     FORALL_CELLS_EXT() 
    122       ON_PRIMAL 
    123         theta(CELL,iq) = theta_rhodz(CELL,iq)/rhodz(CELL) 
    124       END_BLOCK 
    125     END_BLOCK 
    126   END DO 
    127 END_BLOCK 
    12898 
Note: See TracChangeset for help on using the changeset viewer.