Changeset 2541

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

Files:

• README (modified) (1 diff)
• deftank/callphys.def.co2clouds (modified) (2 diffs)
• deftank/callphys.def.hdo (modified) (1 diff)
• deftank/callphys.def.watercycle (modified) (1 diff)
• libf/phymars/check_fields.F90 (added)
• libf/phymars/physiq_mod.F (modified) (6 diffs)

trunk/LMDZ.MARS/README

@@ -3408 +3408 @@
 == 29/06/2021 == AB
 Make simu_MCS.F90 compatible with DYNAMICO lon-lat output files, with a check on the latitude array order when doing the interpolation
+
+== 01/07/2021 == EM
+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.

trunk/LMDZ.MARS/deftank/callphys.def.co2clouds

@@ -3 +3 @@
 #Run with or without tracer transport ?
 tracer=.true.
+
+# check input fields when entering physics? 
+check_physics_inputs=.false.
+# check fields at the end of physics?
+check_physics_outputs=.false.
 
 #Directory where external input files are:
@@ -18 +23 @@
 #Save statistics in file "stats.nc" ?
 callstats =.true.
-
-#Save EOF profiles in file "profiles" for Climate Database?
-calleofdump = .false.
 
 ## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)

trunk/LMDZ.MARS/deftank/callphys.def.hdo

@@ -16 +16 @@
 callstats =.true.
 
-#Save EOF profiles in file "profiles" for Climate Database?
-calleofdump = .false.
+# check input fields when entering physics? 
+check_physics_inputs=.false.
+# check fields at the end of physics?
+check_physics_outputs=.false.
 
 #Directory where external input files are:

trunk/LMDZ.MARS/deftank/callphys.def.watercycle

@@ -16 +16 @@
 callstats =.true.
 
-#Save EOF profiles in file "profiles" for Climate Database?
-calleofdump = .false.
+# check input fields when entering physics? 
+check_physics_inputs=.false.
+# check fields at the end of physics?
+check_physics_outputs=.false.
 
 #Directory where external input files are:

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