Changeset 559 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Mar 5, 2012, 2:34:44 AM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: debugged idealized crater case to solve problems found by JBM. the problems appeared because crater case forced the model to extrapolate at initialization ; hence the crater was not actually here. solved the problem by setting the input sounding at bottom of crater. at the same time introduced coriolis terms. it is now important to use the pert_coriolis keyword at T to avoid impacting background winds. updated all setting files for TESTCRATER, which now works fine.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F (modified) (3 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F

@@ -264 +264 @@
          grid%cosa(i,j)     = 1.
          grid%e(i,j)        = 0.
-         grid%f(i,j)        = 0.         !! MARS: put coriolis here if needed
+         grid%f(i,j)        = 0.
 
       END DO
@@ -401 +401 @@
       !    grid%ht(i,j) = hm1*exp(-(( float(i-icm)/xa1)**2))   &
       !               *( (cos(pii*float(i-icm)/xal1))**2 )
-    write(6,*) '3D hill or crater '
-    write(6,*) 'height, width ',hm,xa
-    grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
+    IF (hm .gt. 0.) THEN
+      write(6,*) '3D hill. height, width: ',hm,xa
+      write(6,*) 'input sounding is out of the mountain'
+      grid%ht(i,j) = alt_input + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
+    ELSE IF (hm .lt. 0.) THEN
+      write(6,*) '3D crater. height, width: ',hm,xa
+      write(6,*) 'input sounding is at the bottom of crater'
+      grid%ht(i,j) = (alt_input - hm) + hm/(1.+(float(i-icm)/xa)**2+(float(j-jcm)/xa)**2)
+      !! AS: cannot use same formula as hill because it would force ideal.exe to extrapolate
+      !!       which is not possible given how interp_0 is written
+    ELSE
+      write(6,*) 'Nothing. Height is 0. Flat topography'
+      grid%ht(i,j) = alt_input
+    ENDIF
   ELSE
     write(6,*) 'linear slope '
     write(6,*) 'height ',hm
-    grid%ht(i,j) = alt_input + hm * float(i)  
+    IF (hm .gt. 0.) THEN  
+      grid%ht(i,j) = alt_input + hm * float(i)  
+    ELSE IF (hm .lt. 0.) THEN
+      !! see above, crater case
+      grid%ht(i,j) = (alt_input - hm) + hm * float(i)
+    ELSE
+      write(6,*) 'Nothing. Height is 0. Flat topography'
+      grid%ht(i,j) = alt_input
+    ENDIF
   ENDIF
       !!!3D crater    
@@ -437 +456 @@
     grid%mars_cice(i,j)=0.
     grid%mars_wice(i,j)=0.
+    write(6,*) 'NOTE TO SELF. slpx and slpy set to 0 which means no slope insolation.'
     grid%slpx(i,j) = 0.
     grid%slpy(i,j) = 0.
+    !!! COMMENT THE LINES BELOW IF YOU DON'T WANT CORIOLIS TERMS
+    !!! cf. doc WRF2008 page 11 for e and f expressions
+    grid%e(i,j) = 2. * EOMEG * COS(pi*lat_input/180.)
+    grid%f(i,j) = 2. * EOMEG * SIN(pi*lat_input/180.)
+    write(6,*) 'CALCULATE CORIOLIS TERM', grid%f(i,j),grid%e(i,j)
 !!MARS
   ENDDO