Changeset 2663

Timestamp:

Apr 19, 2022, 12:06:36 PM (3 years ago)

Author:

dbardet

Message:

Adding check_fields subroutine for checking field values after and before dynamics-to-physics and physics-to-dynamics transfers

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• check_fields.F90 (added)
• physiq_mod.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -13 +13 @@
                   flxw,                    &
                   pdu,pdv,pdt,pdq,pdpsrf)
- 
+
+      use ioipsl_getin_p_mod, only: getin_p 
       use radinc_h, only : L_NSPECTI,L_NSPECTV,naerkind
       use watercommon_h, only : RLVTT, Psat_water,epsi,su_watercycle, RV, T_h2o_ice_liq
@@ -61 +62 @@
                               waterrain, global1d, szangle
       use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
+      use check_fields_mod, only: check_physics_fields
       use conc_mod, only: rnew, cpnew, ini_conc_mod
       use phys_state_var_mod
@@ -439 +441 @@
 #endif
 
+      ! flags to trigger extra sanity checks
+      logical, save ::  check_physics_inputs=.false.
+      logical, save ::  check_physics_outputs=.false.
+!$OPM THREADPRIVATE(check_physics_inputs,check_physics_outputs)     
+
       ! Misc
       character*2 :: str2
@@ -451 +458 @@
 ! --------------------------------
       if (firstcall) then
+        call getin_p("check_physics_inputs", check_physics_inputs)
+        call getin_p("check_physics_outputs", check_physics_outputs)
+
         ! Allocate saved arrays (except for 1D model, where this has already
         ! been done)
@@ -714 +724 @@
          write(*,*) "physiq: end of firstcall"
       endif ! end of 'firstcall'
+
+      if (check_physics_inputs) then
+         !check the validity of input fields coming from the dynamics
+         call check_physics_fields("begin physiq:", pt, pu, pv, pplev, pq)
+      endif
 
 ! ------------------------------------------------------
@@ -2407 +2422 @@
 #endif
       
+      if (check_physics_outputs) then
+         ! Check the validity of updated fields at the end of the physics step
+         call check_physics_fields("end of physiq:", zt, zu, zv, pplev, zq)
+      endif
+
       icount=icount+1