Changeset 239 for trunk/MESOSCALE/LMD_MM_MARS/makemeso

Timestamp:

Jul 21, 2011, 5:24:11 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: corrected a link bug in launch_gcm. added a test in makemeso to check if xy grid is compliant with number of processors.

File:

• trunk/MESOSCALE/LMD_MM_MARS/makemeso (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/makemeso

@@ -524 +524 @@
   #----------------
   case ${numproc} in
-    1) physx=$(expr ${lon} - 1)
-       physy=$(expr ${lat} - 1)
-       physz=$(expr ${level} - 1) ;;
-    2) physx=$(expr ${lon} - 1)
-       physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 2)
-       physz=$(expr ${level} - 1) ;;
-    4) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 2)
-       physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 2)
-       physz=$(expr ${level} - 1) ;;
-    6) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 2)
-       physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 3)
-       physz=$(expr ${level} - 1) ;;
-    8) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 2)
-       physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 4)
-       physz=$(expr ${level} - 1) ;;
-    12) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 3)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 4)
-        physz=$(expr ${level} - 1) ;;
-    16) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 4)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 4)
-        physz=$(expr ${level} - 1) ;;
-    20) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 4)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 5)
-        physz=$(expr ${level} - 1) ;;
-    24) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 4)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 6)
-        physz=$(expr ${level} - 1) ;;
-    32) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 4)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 8)
-        physz=$(expr ${level} - 1) ;;
-    64) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 8)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 8)
-        physz=$(expr ${level} - 1) ;;
-    128) physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ 8)
-        physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ 16) 
-        physz=$(expr ${level} - 1) ;;
-    *) echo not supported ; exit
+      1) divx=1 ; divy=1  ;;
+      2) divx=1 ; divy=2  ;;
+      4) divx=2 ; divy=2  ;;
+      6) divx=2 ; divy=3  ;;
+      8) divx=2 ; divy=4  ;;
+     12) divx=3 ; divy=4  ;;
+     16) divx=4 ; divy=4  ;;
+     20) divx=4 ; divy=5  ;;
+     24) divx=4 ; divy=6  ;;
+     32) divx=4 ; divy=8  ;;
+     64) divx=8 ; divy=8  ;;
+    128) divx=8 ; divy=16 ;;
+      *) echo not yet supported ; exit ;;
   esac
+  restex=$(expr ${physx}%${divx})   
+  restey=$(expr ${physy}%${divy})
+  if [[ ${restex} != 0 || ${restey} != 0 ]]
+  then
+     echo nx-1 shall be divided by ${divx}
+     echo ny-1 shall be divided by ${divy}
+     exit
+  fi
+  physx=$(expr ${lon} - 1) ; physx=$(expr ${physx} \/ ${divx})
+  physy=$(expr ${lat} - 1) ; physy=$(expr ${physy} \/ ${divy})
+  physz=$(expr ${level} - 1)
 
   # change this if you change num_soil_layers in WRF