Changeset 628 for codes/icosagcm/devel/src/dynamics/caldyn_gcm.F90

Timestamp:

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

Author:

dubos

Message:

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

File:

• codes/icosagcm/devel/src/dynamics/caldyn_gcm.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)