Changeset 281

Timestamp:

10/11/14 23:36:13 (10 years ago)

Author:

dubos

Message:

Fixed bug preventing call to physics when itau_phys>1

Location:

codes/icosagcm/trunk/src

Files:

• physics.f90 (modified) (3 diffs)
• timeloop_gcm.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/physics.f90

@@ -73 +73 @@
     TYPE(t_physics_inout) :: args
 
-    IF(MOD(it+1,itau_physics)==0) THEN
+    IF(MOD(it,itau_physics)==0) THEN
     
        SELECT CASE(phys_type)
@@ -102 +102 @@
     USE physics_interface_mod
     USE physics_dcmip_mod, ONLY : full_physics_dcmip => full_physics
+    USE mpipara
     IMPLICIT NONE
     TYPE(t_field),POINTER :: f_phis(:)
@@ -132 +133 @@
     CASE (phys_DCMIP)
        CALL full_physics_dcmip
+    CASE DEFAULT
+       IF(is_mpi_root) PRINT *,'Internal error : illegal value of phys_type', phys_type
+       STOP
     END SELECT
 

codes/icosagcm/trunk/src/timeloop_gcm.f90

@@ -290 +290 @@
     END IF
 
-
-    IF (MOD(it,itau_physics)==0) THEN
-      CALL physics(it,f_phis, f_ps, f_theta_rhodz, f_u,  f_q)
-    ENDIF
+    CALL physics(it,f_phis, f_ps, f_theta_rhodz, f_u,  f_q)
     
   ENDDO