Changeset 278 for codes/icosagcm/trunk/src/time.f90

Timestamp:

10/11/14 23:27:39 (10 years ago)

Author:

dubos

Message:

time_style cleanup

File:

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

codes/icosagcm/trunk/src/time.f90

@@ -48 +48 @@
   
   SUBROUTINE init_time
-  USE earth_const
-  USE getin_mod
-  USE mpipara
-  IMPLICIT NONE
-  REAL(rstd) :: run_length
-
-
-   time_style='dcmip'
-   CALL getin('time_style',time_style)
-
-   IF (TRIM(time_style)=='dcmip')  Then
+    USE earth_const
+    USE getin_mod
+    USE mpipara
+    IMPLICIT NONE
+    REAL(rstd) :: run_length
+    
     dt=90.
     CALL getin('dt',dt)
-
+    write_period=0
+    CALL getin('write_period',write_period)
     itaumax=100
-    CALL getin('itaumax',itaumax)
-
+    CALL getin('itaumax',itaumax)    
     run_length=dt*itaumax
     CALL getin('run_length',run_length)
     itaumax=run_length/dt
-    dt=dt/scale_factor
-
-    write_period=0
-    CALL getin('write_period',write_period)
-    write_period=write_period/scale_factor
+    
+    time_style='dcmip'
+    CALL getin('time_style',time_style)
+    SELECT CASE(TRIM(time_style))
+    CASE('none') ! do nothing
+    CASE('dcmip') ! rescale time step for small-planet experiments
+       dt=dt/scale_factor       
+       write_period=write_period/scale_factor
+       IF (is_mpi_root) PRINT *, 'Output frequency (scaled) set to ',write_period
+    CASE DEFAULT
+       IF (is_mpi_root) PRINT*,"Bad selector for variable time_style >",TRIM(time_style),"> options are <none>, <dcmip>"
+       STOP
+    END SELECT
+    
     itau_out=FLOOR(.5+write_period/dt)
-    IF (is_mpi_root) PRINT *, 'Output frequency (scaled) set to ',write_period, ' : itau_out = ',itau_out
-    ENDIF 
+    IF (is_mpi_root) PRINT *, 'Output frequency itau_out = ',itau_out
 
     itau_adv=1
     CALL getin('itau_adv',itau_adv)
-
+    
     itau_dissip=1
     CALL getin('itau_dissip',itau_dissip)
-
+    
     itau_physics=1
     CALL getin('itau_physics',itau_physics)
-
+    
     IF (is_mpi_root)  THEN
        PRINT *, 'itaumax=',itaumax