Changeset 871 for codes/icosagcm/trunk/src/physics/physics.f90

Timestamp:

05/17/19 15:02:07 (5 years ago)

Author:

ymipsl

Message:

experimental : add smooth physics tendency for external physics.

Activated with param :

phys_smooth_tendency=y

To be tested...

YM

File:

• codes/icosagcm/trunk/src/physics/physics.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/physics/physics.f90

@@ -20 +20 @@
   CHARACTER(LEN=255),SAVE :: physics_type
 !$OMP THREADPRIVATE(physics_type)
+  TYPE(t_message),SAVE :: req_theta0, req_ue0, req_q0
 
   PUBLIC :: physics, init_physics, zero_du_phys
@@ -115 +116 @@
 
     IF(is_mpi_root) PRINT *, 'phys_type = ',phys_type
+    
+    
 !$OMP END PARALLEL
   END SUBROUTINE init_physics
@@ -159 +162 @@
     TYPE(t_field),POINTER :: f_wflux(:)
     TYPE(t_field),POINTER :: f_q(:)
-
+    
+    LOGICAL,SAVE :: first=.TRUE.
+!$OMP THREADPRIVATE(first)
+    
     LOGICAL:: firstcall,lastcall
     INTEGER :: ind
     TYPE(t_physics_inout) :: args
 
-    IF(MOD(it,itau_physics)==0 .AND. phys_type/=phys_none) THEN
-
-       ! as a result of the the two calls to add_du_phys,
-       ! du_phys increases by u(after physics) - u (before physics)
-       CALL add_du_phys(-1., f_ue)
-
-       SELECT CASE(phys_type)
-       CASE(phys_HS94)
-          CALL held_suarez(f_ps,f_theta_rhodz,f_ue) 
-       CASE (phys_lmdz_generic)
-         CALL physics_lmdz_generic(it ,f_phis, f_ps, f_theta_rhodz, f_ue, f_wflux, f_q)
-       CASE (phys_external)
-         CALL physics_external(it ,f_phis, f_ps, f_theta_rhodz, f_ue, f_wflux, f_q)
-       CASE(phys_LB2012)
-          CALL phys_venus(f_ps,f_theta_rhodz,f_ue) 
-       CASE DEFAULT
-          CALL physics_column(it, f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
-       END SELECT
-
-       CALL transfert_request(f_theta_rhodz,req_i0)
-       CALL transfert_request(f_ue,req_e0_vect)
-       CALL transfert_request(f_q,req_i0)
-
-       CALL add_du_phys(1., f_ue)
-    END IF
-
-    IF (mod(it,itau_out)==0 ) THEN
-       CALL write_physics_tendencies
-       CALL zero_du_phys
-       SELECT CASE(phys_type)
-       CASE (phys_DCMIP)
-          CALL write_physics_dcmip
-       CASE (phys_DCMIP2016)
-          CALL write_physics_dcmip2016
-       END SELECT
-    END IF
+    IF (first) THEN
+      CALL init_message(f_theta_rhodz, req_i0, req_theta0)
+      CALL init_message(f_ue, req_e0_vect, req_ue0)
+      CALL init_message(f_q, req_i0, req_q0)
+      first=.FALSE.
+    ENDIF
+
+
+    IF (phys_external) THEN
+    
+      CALL physics_external(it ,f_phis, f_ps, f_theta_rhodz, f_ue, f_wflux, f_q)
+
+    ELSE
+    
+      IF(MOD(it,itau_physics)==0 .AND. phys_type/=phys_none) THEN
+
+         ! as a result of the the two calls to add_du_phys,
+         ! du_phys increases by u(after physics) - u (before physics)
+         CALL add_du_phys(-1., f_ue)
+
+         SELECT CASE(phys_type)
+         CASE(phys_HS94)
+            CALL held_suarez(f_ps,f_theta_rhodz,f_ue) 
+         CASE (phys_lmdz_generic)
+           CALL physics_lmdz_generic(it ,f_phis, f_ps, f_theta_rhodz, f_ue, f_wflux, f_q)
+         CASE(phys_LB2012)
+            CALL phys_venus(f_ps,f_theta_rhodz,f_ue) 
+         CASE DEFAULT
+            CALL physics_column(it, f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
+         END SELECT
+
+         CALL send_message(f_theta_rhodz, req_theta0)
+         CALL send_message(f_ue, req_ue0)
+         CALL send_message(f_q, req_q0)
+         CALL wait_message(req_theta0)
+         CALL wait_message(req_ue0)
+         CALL wait_message(req_q0)
+       
+         CALL add_du_phys(1., f_ue)
+      END IF
+
+      IF (mod(it,itau_out)==0 ) THEN
+         CALL write_physics_tendencies
+         CALL zero_du_phys
+         SELECT CASE(phys_type)
+         CASE (phys_DCMIP)
+            CALL write_physics_dcmip
+         CASE (phys_DCMIP2016)
+            CALL write_physics_dcmip2016
+         END SELECT
+      END IF
+    ENDIF
     
   END SUBROUTINE physics