Changeset 1128

Timestamp:

Dec 19, 2013, 10:12:49 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE. corrected an interp bug on surface roughness. added a security check in module_lmd_driver.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars (modified) (2 diffs)
• WRFV2/phys/module_lmd_driver.F (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars

@@ -125 +125 @@
         priority = 1
         dest_type = continuous
-        ## au-dessus de 4.... probleme de valeurs negatives...
-        smooth_option = 1-2-1; smooth_passes=4
+        smooth_option = 1-2-1; smooth_passes=10
         interp_option =         32ppd:sixteen_pt+four_pt+average_4pt+average_16pt+search
         rel_path=               32ppd:rough_GCM/
@@ -133 +132 @@
         interp_option =         64ppd:sixteen_pt+four_pt+average_4pt+average_16pt+search
         rel_path=               64ppd:rough_GCM/
-        interp_option =         64ppd_noHRti:sixteen_pt+four_pt+average_4pt+average_16pt+search
+        interp_option =         64ppd_noHRti:nearest_neighbor
         rel_path=               64ppd_noHRti:rough_GCM/
         interp_option =         custom:sixteen_pt+four_pt+average_4pt+average_16pt+search

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

@@ -981 +981 @@
    z0_val = 0.01
 ENDIF
+!!!! ADDITIONAL SECURITY. INTERP+SMOOTH IN GEOGRID MIGHT YIELD NEGATIVE Z0 !!!
+IF (z0_val < 0.) THEN
+   PRINT *, 'WELL, z0 is NEGATIVE, this is impossible. better stop here.'
+   PRINT *, 'advice --> correct interpolation / smoothing of z0 in WPS'
+   PRINT *, '           -- or check the constant value set in namelist.input'
+   STOP
+ENDIF
 #endif