Ignore:
Timestamp:
05/03/19 13:40:15 (5 years ago)
Author:
dubos
Message:

devel : Cp(T) thermodynamics (TBC)

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r685 r837  
    151151      END_BLOCK 
    152152    END_BLOCK 
     153  CASE(thermo_variable_Cp) 
     154    ! thermodynamics with variable Cp 
     155    !           Cp(T) = Cp0 * (T/T0)^nu 
     156    ! =>            h = Cp(T).T/(nu+1) 
     157 
     158    FORALL_CELLS() 
     159      ON_PRIMAL 
     160        berni(CELL) = .5*(geopot(CELL)+geopot(UP(CELL))) 
     161        cp_ik = cpp*(pk(CELL)/Treff)**nu 
     162        berni(CELL) = berni(CELL) + pk(CELL)*(cp_ik/(nu+1.)-theta(CELL,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s) 
     163      END_BLOCK 
     164    END_BLOCK 
    153165  CASE DEFAULT 
    154166    PRINT *, 'Unsupported value of caldyn_thermo : ',caldyn_thermo  ! FIXME 
Note: See TracChangeset for help on using the changeset viewer.