- Timestamp:
- 10/08/19 16:34:06 (5 years ago)
- Location:
- codes/icosagcm/devel/src/dynamics
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dynamics/compute_NH_geopot.F90
r939 r980 1 1 MODULE compute_NH_geopot_mod 2 2 USE grid_param 3 USE icosa, ONLY : rstd 3 4 IMPLICIT NONE 4 5 PRIVATE 5 6 6 7 LOGICAL, SAVE :: debug_hevi_solver = .FALSE. 8 9 REAL(rstd) :: pbot=1e5, rho_bot=100. ! for NH solver 7 10 8 11 #include "../unstructured/unstructured.h90" 9 12 10 PUBLIC :: compute_NH_geopot,compute_NH_geopot_unst 13 PUBLIC :: compute_NH_geopot,compute_NH_geopot_unst, pbot, rho_bot 11 14 12 15 CONTAINS … … 190 193 191 194 SUBROUTINE compute_NH_geopot(tau, phis, m_ik, m_il, theta, W_il, Phi_il) 192 USE icosa193 195 USE disvert_mod 194 196 USE caldyn_vars_mod 195 197 USE omp_para 196 REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6197 198 REAL(rstd),INTENT(IN) :: tau ! solve Phi-tau*dPhi/dt = Phi_rhs 198 199 REAL(rstd),INTENT(IN) :: phis(iim*jjm) -
codes/icosagcm/devel/src/dynamics/compute_caldyn_solver.F90
r939 r980 1 1 MODULE compute_caldyn_solver_mod 2 2 USE grid_param 3 USE compute_NH_geopot_mod, ONLY : compute_NH_geopot, pbot, rho_bot 3 4 IMPLICIT NONE 4 5 PRIVATE … … 104 105 DECLARE_INDICES 105 106 NUM :: X_ij, rho_ij, T_ij, gamma, Cvd, vreff, Rd_preff 106 REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6 ! FIXME107 107 #define PHI_BOT(ij) Phi_bot 108 108 #include "../kernels_unst/caldyn_mil.k90" … … 122 122 USE omp_para, ONLY : ll_begin, ll_end,ll_beginp1,ll_endp1 123 123 USE disvert_mod, ONLY : ptop 124 USE compute_NH_geopot_mod125 REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6126 124 REAL(rstd),INTENT(IN) :: tau ! "solve" Phi-tau*dPhi/dt = Phi_rhs 127 125 REAL(rstd),INTENT(IN) :: phis(iim*jjm)
Note: See TracChangeset
for help on using the changeset viewer.