Changeset 925

Timestamp:

06/21/19 00:11:29 (5 years ago)

Author:

dubos

Message:

devel : added safeguard for Courant criterion for horizontal dissipation

File:

• codes/icosagcm/devel/src/dissip/dissip_gcm.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/dissip/dissip_gcm.f90

@@ -376 +376 @@
        END SELECT
     END SELECT
-    
+
+    IF(is_master) PRINT *, 'vertical profile of horizontal dissipation : ', zvert(:)
+
     DO l=1,llm
        tau_graddiv(l) = tau_graddiv(l)/zvert(l)
@@ -407 +409 @@
     itau_dissip=MAX(1,itau_dissip)
     dtdissip=itau_dissip*dt
+
     IF (is_master) THEN
        PRINT *,"init_dissip: rayleigh_tau",rayleigh_tau, "mintau ",mintau
@@ -412 +415 @@
     ENDIF
     
+    IF (dtdissip>2.*mintau) THEN 
+       IF(is_master) PRINT *, 'The CFL condition for dissipation dtdissip<2*mintau is violated : dtdissip, mintau ', dtdissip, mintau
+       STOP
+    END IF
+
   END SUBROUTINE dissip_timescale