Changeset 1637

Timestamp:

Nov 30, 2016, 7:29:32 PM (8 years ago)

Author:

jvatant

Message:

Minor changes on dyn. diagnostics in physiq_mod.F90 - corrected zdtdyn, added zdudyn - JVO

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -344 +344 @@
       real zdtadj(ngrid,nlayer)                          ! Convadj Diagnostic.
       real zdtdyn(ngrid,nlayer)                          ! Dynamical Heating (K/s).
+      real zdudyn(ngrid,nlayer)                          ! Dynamical Zonal Wind tendency (m.s-2).
       real,allocatable,dimension(:,:),save :: ztprevious ! Previous loop Atmospheric Temperature (K)                                                         ! Useful for Dynamical Heating calculation.
-!$OMP THREADPRIVATE(ztprevious)
+      real,allocatable,dimension(:,:),save :: zuprevious ! Previous loop Zonal Wind (m.s-1)                                                                  ! Useful for Zonal Wind tendency calculation.
+!$OMP THREADPRIVATE(ztprevious,zuprevious)
 
       real reff(ngrid,nlayer)                       ! Effective dust radius (used if doubleq=T).
@@ -469 +471 @@
         ALLOCATE(q2(ngrid,nlayer+1)) 
         ALLOCATE(ztprevious(ngrid,nlayer))
+        ALLOCATE(zuprevious(ngrid,nlayer))
         ALLOCATE(cloudfrac(ngrid,nlayer))
         ALLOCATE(totcloudfrac(ngrid))
@@ -666 +669 @@
          qsurf_hist(:,:)=qsurf(:,:)
 
-!        Initialize variable for dynamical heating diagnostic.
-!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+!        Initialize variable for dynamical heating and zonal wind tendency diagnostic
+!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          ztprevious(:,:)=pt(:,:)
+         zuprevious(:,:)=pu(:,:)
 
 !        Set temperature just above condensation temperature (for Early Mars)
@@ -1634 +1638 @@
 
       ! Diagnostic.
-      zdtdyn(1:ngrid,1:nlayer)     = pt(1:ngrid,1:nlayer)-ztprevious(1:ngrid,1:nlayer)
+      zdtdyn(1:ngrid,1:nlayer)     = (pt(1:ngrid,1:nlayer)-ztprevious(1:ngrid,1:nlayer)) / ptimestep
       ztprevious(1:ngrid,1:nlayer) = zt(1:ngrid,1:nlayer)
+
+      zdudyn(1:ngrid,1:nlayer)     = (pu(1:ngrid,1:nlayer)-zuprevious(1:ngrid,1:nlayer)) / ptimestep
+      zuprevious(1:ngrid,1:nlayer) = zu(1:ngrid,1:nlayer)
 
       if(firstcall)then
          zdtdyn(1:ngrid,1:nlayer)=0.0
+         zdudyn(1:ngrid,1:nlayer)=0.0
       endif
 
       ! Dynamical heating diagnostic.
       do ig=1,ngrid
-         fluxdyn(ig)= SUM(zdtdyn(ig,:) *mass(ig,:))*cpp/ptimestep
+         fluxdyn(ig)= SUM(zdtdyn(ig,:) *mass(ig,:))*cpp
       enddo