Changeset 282 for codes/icosagcm/branches

Timestamp:

10/17/14 18:26:34 (10 years ago)

Author:

milmd

Message:

Sponge layer in dynamics.

Location:

codes/icosagcm/branches/SATURN_DYNAMICO

Files:

• ICOSAGCM/src/sponge.f90 (added)
• ICOSAGCM/src/timeloop_gcm.f90 (modified) (4 diffs)
• TEST/run_icosa.def (modified) (1 diff)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/timeloop_gcm.f90

@@ -39 +39 @@
   USE output_field_mod
   USE write_field
+  USE sponge_mod
   IMPLICIT NONE
 
@@ -119 +120 @@
 
     CALL init_dissip
+    CALL init_sponge
     CALL init_caldyn
     CALL init_guided
@@ -141 +143 @@
   USE icosa
   USE dissip_gcm_mod
+  USE sponge_mod
   USE disvert_mod
   USE caldyn_mod
@@ -269 +272 @@
 !       CALL wait_message(req_mass)  
        CALL euler_scheme(.FALSE.)  ! update only u, theta
+       IF (iflag_sponge > 0) THEN
+         CALL sponge(f_u,f_du,f_theta_rhodz,f_dtheta_rhodz)
+         CALL euler_scheme(.FALSE.)  ! update only u, theta
+       ENDIF
     END IF
 

codes/icosagcm/branches/SATURN_DYNAMICO/TEST/run_icosa.def

@@ -142 +142 @@
 rayleigh_friction_type=none
 
+#Sponge layer
+### iflag_sponge=0 for no sponge
+### iflag_sponge=1 for sponge over 4 topmost layers
+### iflag_sponge=2 for sponge from top to ~1% of top layer pressure
+### tau_sponge --> damping frequency at last layer
+### e-5 medium / e-4 strong yet reasonable / e-3 very strong
+### mode_sponge=1 for u,v --> 0
+### mode_sponge=2 for u,v --> zonal mean (NOT IMPLEMENTED)
+### mode_sponge=3 for u,v,h --> zonal mean (NOT IMPLEMENTED)
+iflag_sponge = 1
+tau_sponge = 1.e-4
+mode_sponge = 1
+
 # -------------------------------- Physics -------------------------------------