Changeset 250

Timestamp:

07/23/14 12:28:16 (10 years ago)

Author:

dubos

Message:

Ignore dissipation if tau_XX all 0

File:

• codes/icosagcm/trunk/src/dissip_gcm.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/dissip_gcm.f90

@@ -427 +427 @@
     ENDDO
        
-    itau_dissip=INT(mintau/dt)
+    IF(mintau>0) THEN
+       itau_dissip=INT(mintau/dt)
+       dtdissip=itau_dissip*dt
+    ELSE
+       IF (is_mpi_root) PRINT *,"No dissipation time set, setting itau_dissip to 1000000000"
+       itau_dissip=100000000
+    END IF
     itau_dissip=MAX(1,itau_dissip)
-    dtdissip=itau_dissip*dt
     IF (is_mpi_root) PRINT *,"mintau ",mintau,"itau_dissip",itau_dissip," dtdissip ",dtdissip