Changeset 102 in lmdz_wrf

Timestamp:

Jul 25, 2014, 4:02:35 PM (10 years ago)

Author:

lfita

Message:

fm NaN from thermcell_flux2

Location:

branches/LMDZ_WRFmeas/WRFV3/lmdz

Files:

• thermcell_flux2.F90 (modified) (7 diffs)
• thermcell_main.F90 (modified) (5 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90

@@ -55 +55 @@
       character (len=80) :: abort_message
 
+! Lluis
+      INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfname
+      REAL                                               :: largest
+
+      llp = 734
+      lfname = 'physiq'
+      largest = 10.e5
+
       fomass_max=0.5
       alphamax=0.7
@@ -179 +188 @@
 ! les flux avant et apres modif
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      lfname='thermcell_flux2 before first loop'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
 
       do l=1,klev
@@ -193 +209 @@
          enddo
 
+      lfname='thermcell_flux2 after first loop'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
 
 !-------------------------------------------------------------------------
@@ -234 +257 @@
       endif
 
+      lfname='thermcell_flux2 after fraca croissant'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
+
       if (prt_level.ge.10) &
      &   write(lunout1,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
@@ -254 +285 @@
 !     enddo
       endif
+
+      lfname='thermcell_flux2 after massa croissant'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
 
       if (prt_level.ge.10) &
@@ -415 +454 @@
 ! Fin de la grande boucle sur les niveaux verticaux
       enddo
+      lfname='thermcell_flux2 after big loop'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'rhobarz'
+      CALL check_var3D(lfname, lvarname, rhobarz, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'zw2'
+      CALL check_var3D(lfname, lvarname, zw2, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
 
 !      if (prt_level.ge.10) &
@@ -461 +511 @@
          enddo
       enddo
+
+      lfname='thermcell_flux2 after quantile'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      lvarname = 'entr'
+      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr'
+      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
+
       if (labort_gcm) then
                          ig=igout

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -548 +548 @@
 !     print*,"ZMAX ",lalim,lmin,linter,lmix,lmax,zmax,zmax0,zmix,wmax
 
-      lfname='after thermcell_height'
+      lfname='thermcell_main after thermcell_height'
       lvarname = 'zw2'
       CALL check_var3D(lfname, lvarname, zw2, ngrid, nlay+1, largest, .FALSE.)
@@ -618 +618 @@
 !deduction des flux
 !-------------------------------------------------------------------------------
-      lfname='before flux'
+      lfname='thermcell_main before flux'
       lvarname = 'fm'
       CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
@@ -635 +635 @@
      &       detr,zqla,lev_out,lunout1,igout)
 !IM 060508    &       detr,zqla,zmax,lev_out,lunout,igout)
-      lfname='after flux'
+      lfname='thermcell_main after flux'
       lvarname = 'fm'
       CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
@@ -666 +666 @@
 !------------------------------------------------------------------
 
-      lfname='before transport_vertical'
+      lfname='thermcell_main before transport_vertical'
       lvarname = 'zdthladj'
       CALL check_var3D(lfname, lvarname, zdthladj, ngrid, nlay, largest, .FALSE.)
@@ -674 +674 @@
       call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse,  &
      &                   po,pdoadj,zoa,lev_out)
-      lfname='after transport_vertical'
+      lfname='thermcell_main after transport_vertical'
       lvarname = 'zdthladj'
       CALL check_var3D(lfname, lvarname, zdthladj, ngrid, nlay, largest, .FALSE.)