Changeset 989

Timestamp:

11/08/19 17:23:11 (5 years ago)

Author:

jisesh

Message:

devel: Added nudging outside the boundary of the domain in addidion to inside nudging

File:

• codes/icosagcm/devel/src/dissip/nudging_mod.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/dissip/nudging_mod.f90

@@ -12 +12 @@
   PRIVATE
   ! nudging will be active outside a disc of radius 'radius' centered at 'center_lon', 'center lat'.
-  REAL(rstd) :: center_lon, center_lat, nudging_radius, nudging_relaxation_time
-  !$OMP THREADPRIVATE(center_lon, center_lat, nudging_radius, nudging_relaxation_time)
+  REAL(rstd) :: center_lon, center_lat, nudging_radius, nudging_relaxation_time, &
+       nudging_relaxation_time_in,nudging_relaxation_time_out
+  !$OMP THREADPRIVATE(center_lon, center_lat, nudging_radius, nudging_relaxation_time,nudging_relaxation_time_in,nudging_relaxation_time_out)
   TYPE(t_field),POINTER :: f_relax_coef_e(:), f_target_ue(:), &
        f_relax_coef_i(:), f_target_theta_rhodz(:),f_target_ps(:)
@@ -49 +50 @@
     !nudging_time=0
 
-    nudging_relaxation_time = 0.
-    CALL getin('nudging_relaxation_time', nudging_relaxation_time)
-    nudging_relaxation_time = nudging_relaxation_time/scale_factor
+    nudging_relaxation_time_in = 0.
+    nudging_relaxation_time_out = 0.
+    CALL getin('nudging_relaxation_time_in', nudging_relaxation_time_in)
+    CALL getin('nudging_relaxation_time_out', nudging_relaxation_time_out)
+    nudging_relaxation_time_in = nudging_relaxation_time_in/scale_factor
+    nudging_relaxation_time_out = nudging_relaxation_time_out/scale_factor
+    IF(nudging_relaxation_time_in .EQ. 0. .AND. nudging_relaxation_time_out .EQ. 0.) THEN
+        PRINT*," nudging_relaxation_time_in and nudging_relaxation_time_out values can not be 0. at the same time"
+    ENDIF
 
     CALL getin("guided_nudging_field",guided_nudging_field)
@@ -114 +121 @@
        IF (c>0.99) c=1
        IF (c<-0.99) c=-1
-       relax_coef = .5*(1.+c) ! rescale to [0,1] range ; ! 1 inside circle, 0 outside
+       c = .5*(1.+c) ! rescale to [0,1] range ; ! 1 inside circle, 0 outside
     ELSE
        c = 0.
     END IF
 
-    relax_coef = 1. - (1.-relax_coef) * dt / (dt + nudging_relaxation_time)
+    nudging_relaxation_time = c*nudging_relaxation_time_in + (1.-c)*nudging_relaxation_time_out 
+    relax_coef = 1. - dt/(dt + nudging_relaxation_time)
   END FUNCTION relax_coef