Changeset 176 in lmdz_wrf

Timestamp:

Oct 8, 2014, 11:23:46 AM (10 years ago)

Author:

lfita

Message:

Fixing not output of 'EVAP' in wrfout
Adding evap checks
Increasing threshold for the checking of airephy

Location:

trunk/WRFV3/lmdz

Files:

• diagphy_mod.F90 (modified) (1 diff)
• physiq.F90 (modified) (2 diffs)
• wrf_lmdz_mod.F90 (modified) (1 diff)

trunk/WRFV3/lmdz/diagphy_mod.F90

@@ -166 +166 @@
         IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e6) THEN
           varname = 'airephy'
-          CALL check_var(fname, varname, airephy, klon, largest*10.e6, .FALSE.)
+          CALL check_var(fname, varname, airephy, klon, largest*10.e12, .FALSE.)
         END IF
         PRINT *,'    stops= ',stops

trunk/WRFV3/lmdz/physiq.F90

@@ -2206 +2206 @@
 
       PRINT *,'  Lluis before pbl_surface qsol: ',qsol(llp),    &
-        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur
+        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur,' evap: ',evap(llp)
 
       CALL pbl_surface(                                                              &
@@ -2233 +2233 @@
 
       PRINT *,'  Lluis after pbl_surface qsol: ',qsol(llp),    &
-        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur
+        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur,' evap: ',evap(llp)
 
 !-----------------------------------------------------------------------------------------

trunk/WRFV3/lmdz/wrf_lmdz_mod.F90

@@ -2499 +2499 @@
         pwrf_grid%ldtsvdfg(ix,iy) = d_ts(ixy,is_lic)
         pwrf_grid%ldtsvdfi(ix,iy) = d_ts(ixy,is_sic)
-        pwrf_grid%levap(ix,iy) = fevap(ixy,1)
+! L. Fita, October 2014. with fevap there is not evaporation!
+!        pwrf_grid%levap(ix,iy) = fevap(ixy,1)
+        pwrf_grid%levap(ix,iy) = evap(ixy)
         pwrf_grid%levap_ter(ix,iy) = fevap(ixy,1)
         pwrf_grid%levap_lic(ix,iy) = fevap(ixy,2)