Changeset 947 for codes/icosagcm/devel/src/dissip/nudging_mod.f90

Timestamp:

07/10/19 16:31:55 (5 years ago)

Author:

jisesh

Message:

devel: Nudging external data and vertical interpolatiion

File:

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

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

@@ -44 +44 @@
     nudging_time=0.
     CALL getin('nudging_time', nudging_time)
-
+    nudging_time = nudging_time/scale_factor
     ! we should check that radius>0
     CALL getin("guided_nudging_field",guided_nudging_field)
@@ -147 +147 @@
     USE transfert_mod
     USE time_mod
+    USE pression_mod
+    USE vertical_remap_mod
     REAL(rstd), INTENT(IN):: tt
     TYPE(t_field),POINTER :: f_ps(:)
+    TYPE(t_field),POINTER,SAVE :: f_pmid_target(:)
     TYPE(t_field),POINTER :: f_phis(:)
     TYPE(t_field),POINTER :: f_theta_rhodz(:)
     TYPE(t_field),POINTER :: f_u(:)  
     TYPE(t_field),POINTER :: f_q(:)
-    TYPE(t_field),POINTER, SAVE :: f_T_guided(:),  f_ulon_guided(:), f_ulat_guided(:),f_pressure_mid(:) 
+    TYPE(t_field),POINTER, SAVE :: f_T_guided(:),  f_ulon_guided(:), f_ulat_guided(:),f_T_guided_interp(:), &
+                                   f_ulon_guided_interp(:),f_ulat_guided_interp(:)
     REAL(rstd), POINTER :: target_ue(:,:), ue(:,:), coef_e(:)
     REAL(rstd), POINTER :: target_theta_rhodz(:,:,:), theta_rhodz(:,:,:), coef_i(:) 
     INTEGER :: ind
-
-    IF (abs(MOD(tt,REAL(nudging_time))-dt) < 1.0D-2) THEN
-
+    
+    IF (abs(MOD(tt,REAL(nudging_time))-dt) < 1.0D-2) THEN 
+       
       CALL allocate_field(f_T_guided, field_t, type_real, llm, name='nudging_T')
+      CALL allocate_field(f_T_guided_interp, field_t, type_real, llm, name='nudging_T')
+      CALL allocate_field(f_ulon_guided_interp, field_t, type_real, llm, name='nudging_T')
+      CALL allocate_field(f_ulat_guided_interp, field_t, type_real, llm, name='nudging_T')
       CALL allocate_field(f_ulon_guided, field_t, type_real, llm, name='nudging_ulon')
       CALL allocate_field(f_ulat_guided, field_t, type_real, llm, name='nudging_ulat')
-
+      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 vertical_remap_extdata(f_T_guided,f_pmid_target,f_T_guided_interp)
       CALL transfert_request(f_T_guided,req_i0)
-      CALL temperature2theta_rhodz(f_ps,f_T_guided,f_target_theta_rhodz)
+      CALL transfert_request(f_T_guided_interp,req_i0)
+      CALL transfert_request(f_T_guided,req_i1)
+      CALL transfert_request(f_T_guided_interp,req_i1)
+      CALL temperature2theta_rhodz(f_ps,f_T_guided_interp,f_target_theta_rhodz)
       CALL xios_read_field("ulon_guided_read",f_ulon_guided)
       CALL xios_read_field("ulat_guided_read",f_ulat_guided)
+      CALL vertical_remap_extdata(f_ulon_guided,f_pmid_target,f_ulon_guided_interp)
+      CALL vertical_remap_extdata(f_ulat_guided,f_pmid_target,f_ulat_guided_interp)
       CALL transfert_request(f_ulon_guided,req_i0)
+      CALL transfert_request(f_ulon_guided_interp,req_i0)
       CALL transfert_request(f_ulat_guided,req_i0)
+      CALL transfert_request(f_ulat_guided_interp,req_i0)
       CALL transfert_request(f_ulon_guided,req_i1)
+      CALL transfert_request(f_ulon_guided_interp,req_i1)
       CALL transfert_request(f_ulat_guided,req_i1)
+      CALL transfert_request(f_ulat_guided_interp,req_i1)
       CALL xios_write_field("ulat_guided_out",f_ulat_guided)
       CALL xios_write_field("ulon_guided_out",f_ulon_guided)
-      CALL ulonlat2un(f_ulon_guided, f_ulat_guided,f_target_ue)
-
+      CALL xios_write_field("T_guided_out",f_T_guided)
+      CALL ulonlat2un(f_ulon_guided_interp, f_ulat_guided,f_target_ue) 
       CALL deallocate_field(f_T_guided)
+      CALL deallocate_field(f_T_guided_interp)
       CALL deallocate_field(f_ulon_guided)
+      CALL deallocate_field(f_ulon_guided_interp)
       CALL deallocate_field(f_ulat_guided)
+      CALL deallocate_field(f_ulat_guided_interp)
+      CALL deallocate_field(f_pmid_target)
     
     ENDIF