Changeset 2287 for trunk/MESOSCALE

Timestamp:

Apr 15, 2020, 3:07:19 PM (5 years ago)

Author:

mlefevre

Message:

MESOSCALE. Corrections of local time calculations for Venus mesoscale and LES model.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F (modified) (1 diff)
• share/module_date_time.F (modified) (3 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F

@@ -48 +48 @@
     !
   IF (JULYR .ne. 9999) THEN
-    JH_cur_split = (GMT + elaps/420000.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
+    JH_cur_split = (GMT + elaps/421200.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
     JH_cur_split = MODULO(JH_cur_split,24.)   !! the two arguments of MODULO must be of the same type 
     JH_cur_split = JH_cur_split / 24.
-    JD_cur = (JULDAY - 1 + INT((420000.0*GMT+elaps)/1.008e7))
+    JD_cur = (JULDAY - 1 + INT((421200.*GMT+elaps)/10108800.))
     JD_cur = MODULO(int(JD_cur),2)
-    MY = (JULYR-2000) + (1.008e7*(JULDAY - 1)+420000.0*GMT+elaps)/2.016e7
+    MY = (JULYR-2000) + (10108800.*(JULDAY - 1)+421200.*GMT+elaps)/20217600.
     MY = INT(MY)
   ELSE
-    JH_cur_split = lct_input - lon_input / 15. + elaps/(4200.)
-    JH_cur_split = MODULO(JH_cur_split,2808.)
-    JH_cur_split = JH_cur_split / 2808.
+    JH_cur_split = lct_input - lon_input / 15. + elaps/(421200.)
+    JH_cur_split = MODULO(JH_cur_split,24.)
+    JH_cur_split = JH_cur_split / 24.
     print*,'JH_cur_split',JH_cur_split
     JD_cur = INT((3700*(lct_input - lon_input / 15.) +elaps)/88800)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/share/module_date_time.F

@@ -57 +57 @@
 !****VNEUS
       !JULDAYV=MODULO(JULDAYV,2400.)
-      !GMT=((JULDAYV/100.))+(nh+FLOAT(ni)/13.+FLOAT(ns)/117.)/420000.
+      !GMT=((JULDAYV/100.))+(nh+FLOAT(ni)/13.+FLOAT(ns)/117.)/421200.
       !GMT=MODULO(GMT,24.)
 
@@ -104 +104 @@
 !****VENUS
      !JULDAYV=MODULO(JULDAYV,2400.)
-     !GMT=((JULDAYv/100.))+(nh+FLOAT(ni)/13.+FLOAT(ns)/117.)/420000.
+     !GMT=((JULDAYv/100.))+(nh+FLOAT(ni)/13.+FLOAT(ns)/117.)/421200..
      !GMT=MODULO(GMT,24.)
 
@@ -220 +220 @@
       mday(11) = 51
       mday(12) = 56
-!     maday=24 !VENUS
+!     mday=24 !VENUS
       
       !  Break down old hdate into parts