Changeset 969 for codes

Timestamp:

09/02/19 16:19:52 (5 years ago)

Author:

jisesh

Message:

devel : added select case construct for nudging options

File:

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

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

@@ -186 +186 @@
     
     IF (abs(MOD(tt,REAL(nudging_time))-dt) < 1.0D-2) THEN
-
       CALL allocate_field(f_ps_guided, field_t, type_real, name='nudging_ps')
       CALL xios_read_field("ps_guided_read",f_ps_guided)
-      !f_ps = f_ps_guided  !================================================>FIX
       CALL transfert_request(f_ps_guided,req_i0)
       CALL transfert_request(f_ps_guided,req_i1)
@@ -200 +198 @@
       CALL allocate_field(f_pmid_target,field_t,type_real,llm,name='target_pressure') 
       CALL xios_read_field("T_guided_read",f_T_guided)
-      CALL pression_mid(f_ps, f_pmid_target) ! call ps nudging before this
+      CALL pression_mid(f_ps, f_pmid_target) 
       CALL vertical_remap_extdata(f_T_guided,f_pmid_target,f_T_guided_interp)
       CALL transfert_request(f_T_guided,req_i0)
@@ -220 +218 @@
       CALL transfert_request(f_ulat_guided_interp,req_i1)
       !CALL xios_write_field("ulat_guided_out",f_ulat_guided)
-      CALL xios_write_field("ulat_guided_out",f_T_guided_interp)
+      CALL xios_write_field("ulat_guided_out",f_T_guided_interp) ! FIX
       CALL xios_write_field("ulon_guided_out",f_ulon_guided)
       CALL xios_write_field("T_guided_out",f_T_guided)
-      CALL xios_write_field("PS_guided_out",f_ps_guided)
+      CALL xios_write_field("PS_guided_out",f_ps_guided) 
       CALL ulonlat2un(f_ulon_guided_interp, f_ulat_guided,f_target_ue) 
       CALL deallocate_field(f_T_guided)
@@ -232 +230 @@
       CALL deallocate_field(f_ulat_guided_interp)
       CALL deallocate_field(f_pmid_target)
-      !CALL deallocate_field(f_ps_guided)
-      print*,"-----------------------nudging end------------------"
-
     ENDIF
 
@@ -242 +237 @@
        CALL swap_dimensions(ind)
        CALL swap_geometry(ind) 
-       coef_e = f_relax_coef_e(ind)
-       target_ue = f_target_ue(ind)
-       ue = f_u(ind)
-       CALL compute_guided_u(coef_e, target_ue, ue)
-       coef_i = f_relax_coef_i(ind)
-       target_theta_rhodz = f_target_theta_rhodz(ind)
-       theta_rhodz = f_theta_rhodz(ind)
-       !CALL compute_guided_center(coef_i, target_theta_rhodz, theta_rhodz)
-       target_ps = f_ps_guided(ind)
-       ps2 = f_ps(ind)
-       CALL compute_guided_center2(coef_i, target_ps, ps2)
+       
+       SELECT CASE(TRIM(guided_nudging_field))
+         CASE ('wind')
+           coef_e = f_relax_coef_e(ind)
+           target_ue = f_target_ue(ind)
+           ue = f_u(ind)
+           CALL compute_guided_u(coef_e, target_ue, ue)
+         CASE ('temperature')
+           coef_i = f_relax_coef_i(ind)
+           target_theta_rhodz = f_target_theta_rhodz(ind)
+           theta_rhodz = f_theta_rhodz(ind)
+           CALL compute_guided_center(coef_i, target_theta_rhodz, theta_rhodz)
+         CASE ('ps')
+           coef_i = f_relax_coef_i(ind)
+           target_ps = f_ps_guided(ind)
+           ps2 = f_ps(ind)
+           CALL compute_guided_center2(coef_i, target_ps, ps2)
+         CASE ('wind_temperature_ps')
+           coef_e = f_relax_coef_e(ind)
+           target_ue = f_target_ue(ind)
+           ue = f_u(ind)
+           CALL compute_guided_u(coef_e, target_ue, ue)
+           coef_i = f_relax_coef_i(ind)
+           target_theta_rhodz = f_target_theta_rhodz(ind)
+           theta_rhodz = f_theta_rhodz(ind)
+           CALL compute_guided_center(coef_i, target_theta_rhodz, theta_rhodz)
+           target_ps = f_ps_guided(ind)
+           ps2 = f_ps(ind)
+           CALL compute_guided_center2(coef_i, target_ps, ps2)
+         END SELECT
     END DO