Changeset 240 in lmdz_wrf for trunk/WRFV3

Timestamp:

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

Author:

lfita

Message:

Adding map factor on tendencies

Location:

trunk/WRFV3

Files:

• dyn_em/solve_em.F (modified) (2 diffs)
• lmdz/diagphy_mod.F90 (modified) (7 diffs)
• phys/module_lmdz_phys.F (modified) (3 diffs)

trunk/WRFV3/dyn_em/solve_em.F

@@ -915 +915 @@
         &        WRF_PERMASS=grid%mu_1, WRF_BASEMASS=grid%mub,                       &
         &        WRF_MUT=grid%mut, WRF_MUU=grid%muu, WRF_MUV=grid%muv,               &
+        &        WRF_MAPFAC=grid%msft,                                      &
 !!        &        WRF_UTEND=grid%ru_tend, WRF_VTEND=grid%rv_tend,                     &
 !!        &        WRF_TTEND=t_tend,                                                   &
@@ -3702 +3703 @@
    ENDIF
 
+!!L. Fita, January 2105, LMD. We have mp_physics == 0, thus, WRF is not doing
+!  this...
    IF (config_flags%mp_physics /= 0)  then
+!   IF (config_flags%mp_physics /= 0)  then
 
      !$OMP PARALLEL DO   &

trunk/WRFV3/lmdz/diagphy_mod.F90

@@ -164 +164 @@
         PRINT *,'    fs_bound = stops-stopl - (ssols+ssoll)+ssens+sfront'
         PRINT *,'    airetot= ',airetot
-        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e6) THEN
+        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e16) THEN
           varname = 'airephy'
-          CALL check_var(fname, varname, airephy, klon, largest*10.e12, .FALSE.)
+          CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         END IF
         PRINT *,'    stops= ',stops
@@ -223 +223 @@
         PRINT *,'    evap_tot= ',evap_tot
         PRINT *,'    airetot= ',airetot
-        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e5) THEN
+        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e15) THEN
           varname = 'airephy'
-          CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+          CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         END IF
         IF (evap_tot .NE. evap_tot .OR. ABS(evap_tot) > largest) THEN
@@ -673 +673 @@
           PRINT *,'      h_vcol_tot = h_dair_tot+h_qw_tot+h_ql_tot+h_qs_tot'
           PRINT *,'    airetot= ',airetot
-          IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e5) THEN
+          IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e15) THEN
             varname = 'airephy'
-            CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+            CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
           END IF
           PRINT *,'    h_dair_tot= ',h_dair_tot
@@ -718 +718 @@
         PRINT *,'    qw_tot = f(q,ql,qs,airephy,t,paprs)'
         varname = 'airephy'
-        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         varname = 'q'
         CALL check_var3D(fname, varname, q, klon, klev, largest, .FALSE.)
@@ -738 +738 @@
         PRINT *,'    ql_tot = f(ql,airephy,t,paprs)'
         varname = 'airephy'
-        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         varname = 'ql'
         CALL check_var3D(fname, varname, ql, klon, klev, largest, .FALSE.)
@@ -754 +754 @@
         PRINT *,'    qs_tot = f(qs,airephy,t,paprs)'
         varname = 'airephy'
-        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         varname = 'qs'
         CALL check_var3D(fname, varname, qs, klon, klev, largest, .FALSE.)
@@ -770 +770 @@
         PRINT *,'    ec_tot = f(u,v,airephy,t,paprs)'
         varname = 'airephy'
-        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        CALL check_var(fname, varname, airephy, klon, largest*10.e15, .FALSE.)
         varname = 'u'
         CALL check_var3D(fname, varname, u, klon, klev, largest, .FALSE.)

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