Changeset 46 in lmdz_wrf for branches/LMDZ_WRFmeas

Timestamp:

Jul 23, 2014, 4:04:25 PM (10 years ago)

Author:

lfita

Message:

Adding new checks in diagphy_mod.F90: d_etp_tot, d_qt_tot
Adding checkig PRINTS for 'd_h_vcol' in 'physiq.F90'

Location:

branches/LMDZ_WRFmeas/WRFV3/lmdz

Files:

• diagphy_mod.F90 (modified) (2 diffs)
• pbl_surface_mod.F90 (modified) (2 diffs)
• physiq.F90 (modified) (11 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -207 +207 @@
       END IF
 
+      IF (d_etp_tot .NE. d_etp_tot .OR. ABS(d_etp_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong d_etp_tot= ',d_etp_tot,' !!!'
+        PRINT *,'    d_etp_tot: Heat flux equivalent to atmospheric enthalpy' //     &
+          ' change (W/m2)'
+        PRINT *,'    d_etp_tot = input value !'
+        STOP
+      END IF
+
       IF (fq_bound .NE. fq_bound .OR. ABS(fq_bound) > 100000.) THEN
         PRINT *,TRIM(errmsg)
@@ -232 +241 @@
           CALL check_var(fname, varname, snow_fall, klon, largest, .FALSE.)
         END IF
+        STOP
+      END IF
+
+      IF (d_qt_tot .NE. d_qt_tot .OR. ABS(d_qt_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong d_qt_tot= ',d_qt_tot,' !!!'
+        PRINT *,'    d_qt_tot: Mass flux equivalent to atmospheric watter mass' //   &
+          ' change (kg/m2/s)'
+        PRINT *,'    d_qt_tot = input value !'
         STOP
       END IF

branches/LMDZ_WRFmeas/WRFV3/lmdz/pbl_surface_mod.F90

@@ -1040 +1040 @@
              flux_v(i,k,nsrf) = y_flux_v(j,k)
 
-
+             IF (y_flux_t())
           ENDDO
        ENDDO
@@ -1247 +1247 @@
     zxfluxt(:,:) = 0.0 ; zxfluxq(:,:) = 0.0
     zxfluxu(:,:) = 0.0 ; zxfluxv(:,:) = 0.0
-    PRINT *,' Lluis checking in nsrf pbl_surface_mod zxfluxt fluxt_t pctsrf 1 ___________'
+    PRINT *,' Lluis checking in nsrf pbl_surface_mod zxfluxt flux_t pctsrf 1 ___________'
     DO nsrf = 1, nbsrf
        DO k = 1, klev

branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90

@@ -1962 +1962 @@
        &      , fs_bound, fq_bound )
       END IF 
+      PRINT *,'  Lluis 1 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 
 !c Diagnostiquer la tendance dynamique
@@ -2122 +2123 @@
 !C
       END IF 
-
+      PRINT *,'  Lluis 2 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c
 !c=========================================================================
@@ -2265 +2266 @@
 
       ENDIF
+      PRINT *,'  Lluis 3 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c =================================================================== c
 !c   Calcul de Qsat
@@ -2625 +2627 @@
       END IF 
 !C
+      PRINT *,'  Lluis 4 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
       IF (check) THEN
           za = qcheck(klon,klev,paprs,q_seri,ql_seri,airephy)
@@ -2737 +2740 @@
        &      , fs_bound, fq_bound )
       END IF 
-
+      PRINT *,'  Lluis 5 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c      print*,'apres callwake iflag_cldcon=', iflag_cldcon
 !c
@@ -2962 +2965 @@
        &      , fs_bound, fq_bound )
       END IF 
-
+      PRINT *,'  Lluis 6 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 
 !c-------------------------------------------------------------------------
@@ -3029 +3032 @@
        &      , fs_bound, fq_bound )
       END IF 
+      PRINT *,'  Lluis 7 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
       if (mydebug) then
         call writefield_phy('u_seri',u_seri,llm)
@@ -3261 +3265 @@
        &      , fs_bound, fq_bound )
       END IF 
+      PRINT *,'  Lluis 8 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c
 !c Calculer l'humidite relative pour diagnostique
@@ -3633 +3638 @@
        &      , fs_bound, fq_bound )
       END IF 
+      PRINT *,'  Lluis 9 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c
 !c
@@ -3807 +3813 @@
        &      , fs_bound, fq_bound )
       END IF 
+      PRINT *,'  Lluis 10 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
 !c
 !c
@@ -3952 +3959 @@
 !C
       END IF 
+      PRINT *,'  Lluis 11 d_h_vcol: ',d_h_vcol,' d_h_vcol_phy: ',d_h_vcol_phy
       PRINT *,'Lluis Reaching the SORTIES point'