Changeset 2541 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jul 1, 2021, 8:36:13 AM (4 years ago)

Author:

emillour

Message:

Mars GCM:
Add possibility of additional tests (NaNs?, but also of unrealistic values) of
fields at the begining of physics (i.e. coming from the dynamics) and at the
end of the physics integration. These are respectively triggered by seting
flags "check_physics_inputs" and "check_physics_outputs" to .true.
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

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

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -74 +74 @@
       use compute_dtau_mod, only: compute_dtau
       use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
+      use check_fields_mod, only: check_physics_fields
 #ifdef MESOSCALE
       use comsoil_h, only: mlayer,layer
@@ -96 +97 @@
 #endif
       USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured
+      use ioipsl_getin_p_mod, only: getin_p
 
       IMPLICIT NONE
@@ -507 +509 @@
       integer iloop
 
-!      LOGICAL startphy_file
+      ! flags to trigger extra sanity checks
+      logical,save :: check_physics_inputs=.false.
+      logical,save :: check_physics_outputs=.false.
 
 c=======================================================================
@@ -519 +523 @@
 
       IF (firstcall) THEN
+
+         call getin_p("check_physics_inputs",check_physics_inputs)
+         call getin_p("check_physics_outputs",check_physics_outputs)
 
 c        variables set to 0
@@ -735 +742 @@
       ENDIF        !  (end of "if 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)
+      endif
+
 c ---------------------------------------------------
 c 1.2   Initializations done at every physical timestep:
@@ -3888 +3900 @@
 #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,zplev)
+      endif
 
       icount=icount+1