Changeset 1777 for trunk

Timestamp:

Aug 27, 2017, 9:41:59 AM (7 years ago)

Author:

aslmd

Message:

MESOSCALE setting the stage for maybe fixing nesting in the LMD_MM_MARS. added saving the number of previous domai computed. print details only when changing domain. minor changes to printing.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new (modified) (8 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

@@ -159 +159 @@
    REAL :: wappel_phys
    LOGICAL, SAVE :: flag_first_restart
+   INTEGER, SAVE :: previous_id = 0
 !**************************************************
 ! IMPORTANT: pour travailler avec grid nesting
@@ -188 +189 @@
 
    REAL :: tk1,tk2
+
 !==================================================================
 ! CODE
 !==================================================================
 
+print *, '** ',planet_type,'** ENTER WRF-LMD DRIVER'
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !! determine here if this is turbulence-resolving or not
 IF (JULYR .ne. 9999) THEN
@@ -214 +219 @@
      ENDIF
 ENDIF
-
-
-print *,'** ',planet_type,' ** DOMAIN',id
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
 !-------------------------!
@@ -280 +283 @@
 ENDIF
 
-IF ((elaps .eq. 0.).or.flag_first_restart) THEN
+is_first_step: IF ((elaps .eq. 0.).or.flag_first_restart) THEN
 firstcall=.true.  !! for continuity with GCM, physics are always called at the first WRF timestep
   !firstcall=.false. !! just in case you'd want to get rid of the physics 
@@ -324 +327 @@
 flag_first_restart=.false.
 #endif
-
-!! put here some general information you'd like to print just once
-    print *, 'TILES: ', i_start,i_end, j_start, j_end  ! numbers for simple runs, arrays for parallel runs
-    print *, 'DOMAIN: ', ids, ide, jds, jde
-    print *, 'MEMORY: ', ims, ime, jms, jme
-    print *, 'ADVECTED TRACERS: ', num_3d_s-1
-    print *, 'PHYSICS IS CALLED EACH...',wappel_phys
-!! put here some general information you'd like to print just once
 ELSE
 !-------------------------------------------------!
@@ -345 +340 @@
 test = 9999   
 ENDIF
-ENDIF 
+ENDIF is_first_step
+
+!------------------------------------!
+! print info about domain initially  !
+! ... and whenever domain is changed !
+!------------------------------------!
+print *,'** ',planet_type,' ** DOMAIN',id
+IF (previousid .ne. id) THEN
+    print *, '** ',planet_type,' ** ... INITIALIZE DOMAIN',id
+    print *, '** ',planet_type,' ** ... PREVIOUS DOMAIN was',previousid
+    print *, 'ITIMESTEP: ',itimestep
+    print *, 'TILES: ', i_start,i_end, j_start, j_end  ! numbers for simple runs, arrays for parallel runs
+    print *, 'DOMAIN: ', ids, ide, jds, jde
+    print *, 'MEMORY: ', ims, ime, jms, jme
+    print *, 'COMPUT: ', ips, ipe, jps, jpe
+    print *, 'SIZE OF PHYSICS GRID for this process: ',ngrid
+    print *, 'ADVECTED TRACERS: ', num_3d_s-1
+    print *, 'PHYSICS IS CALLED EACH...',wappel_phys
+    print *, '-- means: PHYSICAL TIMESTEP=',zdt_split
+ENDIF
+
 
 !!!***********!!
@@ -722 +737 @@
             RDUST,VMR_ICE,RICE)
 !!!
-print*,"update_outputs_physiq_diag"
-
+print *, '** ',planet_type,'** OUTPUT PHYSICS DONE'
 
 ENDIF call_physics
@@ -790 +804 @@
 !!*****!!
 print *,'** ',planet_type,'** END LMD PHYSICS'
+previous_id = id
 !----------------------------------------------------------------!
 ! use debug (see solve_em) if you wanna display some message ... !