Changeset 980

Timestamp:

10/08/19 16:34:06 (5 years ago)

Author:

dubos

Message:

devel : reduced ground stiffness for better stability (NH)

Location:

codes/icosagcm/devel/src/dynamics

Files:

• compute_NH_geopot.F90 (modified) (2 diffs)
• compute_caldyn_solver.F90 (modified) (3 diffs)

codes/icosagcm/devel/src/dynamics/compute_NH_geopot.F90

@@ -1 +1 @@
 MODULE compute_NH_geopot_mod
   USE grid_param
+  USE icosa, ONLY : rstd
   IMPLICIT NONE
   PRIVATE
 
   LOGICAL, SAVE :: debug_hevi_solver = .FALSE.
+ 
+  REAL(rstd) :: pbot=1e5, rho_bot=100. ! for NH solver
 
 #include "../unstructured/unstructured.h90"
 
-  PUBLIC :: compute_NH_geopot,compute_NH_geopot_unst
+  PUBLIC :: compute_NH_geopot,compute_NH_geopot_unst, pbot, rho_bot
 
 CONTAINS
@@ -190 +193 @@
 
   SUBROUTINE compute_NH_geopot(tau, phis, m_ik, m_il, theta, W_il, Phi_il)
-    USE icosa
     USE disvert_mod
     USE caldyn_vars_mod
     USE omp_para
-    REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6
     REAL(rstd),INTENT(IN)    :: tau ! solve Phi-tau*dPhi/dt = Phi_rhs
     REAL(rstd),INTENT(IN)    :: phis(iim*jjm)

codes/icosagcm/devel/src/dynamics/compute_caldyn_solver.F90

@@ -1 +1 @@
 MODULE compute_caldyn_solver_mod
   USE grid_param
+  USE compute_NH_geopot_mod, ONLY : compute_NH_geopot, pbot, rho_bot
   IMPLICIT NONE
   PRIVATE
@@ -104 +105 @@
   DECLARE_INDICES
   NUM :: X_ij, rho_ij, T_ij, gamma, Cvd, vreff, Rd_preff
-  REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6  ! FIXME
 #define PHI_BOT(ij) Phi_bot
 #include "../kernels_unst/caldyn_mil.k90"
@@ -122 +122 @@
     USE omp_para, ONLY : ll_begin, ll_end,ll_beginp1,ll_endp1
     USE disvert_mod, ONLY : ptop
-    USE compute_NH_geopot_mod
-    REAL(rstd), PARAMETER :: pbot=1e5, rho_bot=1e6
     REAL(rstd),INTENT(IN) :: tau ! "solve" Phi-tau*dPhi/dt = Phi_rhs
     REAL(rstd),INTENT(IN)    :: phis(iim*jjm)