Changeset 861
- Timestamp:
- 05/09/19 01:27:19 (5 years ago)
- Location:
- codes/icosagcm/devel/src/dynamics
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dynamics/caldyn_gcm.F90
r846 r861 17 17 SUBROUTINE init_caldyn 18 18 USE compute_mod, ONLY : compute_pvort_only 19 USE compute_pvort_only_mod, ONLY : compute_pvort_only_hex 19 USE compute_pvort_only_mod, ONLY : compute_pvort_only_hex, compute_pvort_only_unst 20 20 21 21 CHARACTER(len=255) :: def … … 149 149 150 150 CALL allocate_caldyn 151 152 IF(grid_type==grid_ico) THEN 151 152 SELECT CASE(grid_type) 153 CASE(grid_ico) 153 154 ! NB : unused now, would be useful for non-traditional / deep-atmosphere 154 155 DO ind=1,ndomain … … 159 160 CALL compute_planetvel(planetvel) 160 161 END DO 161 END IF 162 163 compute_pvort_only => compute_pvort_only_hex 162 compute_pvort_only => compute_pvort_only_hex 163 CASE(grid_unst) 164 compute_pvort_only => compute_pvort_only_unst 165 END SELECT 164 166 165 167 END SUBROUTINE init_caldyn -
codes/icosagcm/devel/src/dynamics/caldyn_hevi.f90
r859 r861 150 150 qv=f_qv(ind) 151 151 hv=f_hv(ind) 152 ! PRINT *, 'caldyn_hevi : u, mass, qu qv, hv', & 153 ! SHAPE(u), SHAPE(mass), SHAPE(qu), SHAPE(qv), SHAPE(hv) 152 154 CALL compute_pvort_only(u,mass,qu,qv,hv) 153 155 IF(caldyn_kinetic==kinetic_consistent) THEN
Note: See TracChangeset
for help on using the changeset viewer.