Changeset 628 for codes/icosagcm/devel/src/dynamics

Timestamp:

12/11/17 14:39:35 (7 years ago)

Author:

dubos

Message:

devel : new "conservative" variant for vertical transport of momentum

Location:

codes/icosagcm/devel/src/dynamics

Files:

• caldyn_gcm.F90 (modified) (1 diff)
• caldyn_kernels_base.F90 (modified) (1 diff)

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

@@ -37 +37 @@
     CALL getin("dysl_caldyn_vert",dysl_caldyn_vert)
 
+    def='advective'
+    CALL getin('caldyn_vert',def)
+    SELECT CASE(TRIM(def))
+    CASE('advective')
+       caldyn_vert_variant=caldyn_vert_noncons
+    CASE('conservative')
+       IF(.NOT. dysl_caldyn_vert) THEN
+          IF (is_mpi_root) PRINT *,'caldyn_vert=cons requires dysl_caldyn_vert=.TRUE.'
+          STOP
+       END IF
+       caldyn_vert_variant=caldyn_vert_cons
+    CASE DEFAULT
+       IF (is_mpi_root) PRINT *,'Bad selector for variable caldyn_vert : <', &
+            TRIM(def),'> options are <conservative>, <advective>'
+       STOP
+    END SELECT
+    IF (is_master) PRINT *, 'caldyn_vert=',def
+    
     def='energy'
     CALL getin('caldyn_conserv',def)

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

@@ -180 +180 @@
     REAL(rstd),INTENT(INOUT) :: dtheta_rhodz(iim*jjm,llm,nqdyn)
     REAL(rstd),INTENT(OUT) :: dps(iim*jjm)
-    INTEGER, PARAMETER :: caldyn_vert_cons=1, caldyn_vert_variant=caldyn_vert_cons
 
     ! temporary variable