Changeset 240 in lmdz_wrf for trunk/WRFV3/phys

Timestamp:

Feb 2, 2015, 1:05:23 PM (10 years ago)

Author:

lfita

Message:

Adding map factor on tendencies

File:

• trunk/WRFV3/phys/module_lmdz_phys.F (modified) (3 diffs)

trunk/WRFV3/phys/module_lmdz_phys.F

@@ -30 +30 @@
      &        wrf_permass, wrf_basemass,                                             &
      &        wrf_mut, wrf_muu, wrf_muv,                                             &
+     &        wrf_mapfac,                                    &
      &        wrf_Utend, wrf_Vtend, wrf_Ttend,                                       &
      &        wrf_MOISTtend, wrf_psfctend,                                           &
@@ -110 +111 @@
     REAL, INTENT(IN), DIMENSION(wims:wime,wjms:wjme)     :: wrf_Lon, wrf_Lat,        &
         & wrf_basemass, wrf_permass, wrf_mapft, wrf_mapfu, wrf_mapfv, wrf_mut,       &
-        & wrf_muu, wrf_muv
+        & wrf_muu, wrf_muv, wrf_mapfac
 
     REAL, INTENT(OUT), DIMENSION(wims:wime,wjms:wjme)    :: wrf_psfctend
@@ -1289 +1290 @@
 ! All the tendencies have to be given as decoupled from the dry air mass column! (see dyn_em/module_em.F:rk_update_scalar) 
 !!
+! L. Fita, LMD January 2015. Adding map factor!
     DO iz=1,dimz
-      wrf_Ttend(:,iz,:) = wrf_Ttend(:,iz,:)*wrf_mut
-      wrf_Utend(:,iz,:) = wrf_Utend(:,iz,:)*wrf_muu
-      wrf_Vtend(:,iz,:) = wrf_Vtend(:,iz,:)*wrf_muv
+      wrf_Ttend(:,iz,:) = wrf_Ttend(:,iz,:)*wrf_mut/wrf_mapfac
+      wrf_Utend(:,iz,:) = wrf_Utend(:,iz,:)*wrf_muu/wrf_mapfac
+      wrf_Vtend(:,iz,:) = wrf_Vtend(:,iz,:)*wrf_muv/wrf_mapfac
       DO iq=1,wnum3dm
-        wrf_MOISTtend(:,iz,:,iq) = wrf_MOISTtend(:,iz,:,iq)*wrf_mut
+        wrf_MOISTtend(:,iz,:,iq) = wrf_MOISTtend(:,iz,:,iq)*wrf_mut/wrf_mapfac
       END DO
     END DO