Changeset 79 in lmdz_wrf


Ignore:
Timestamp:
Jul 25, 2014, 10:06:36 AM (10 years ago)
Author:
lfita
Message:

Adding checks due to NaNs? in Alp_bl!

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

    r78 r79  
    621621      lvarname = 'zgla'
    622622      CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
    623       lvarname = 'fm'
    624       CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
    625623
    626624      CALL thermcell_flux2(ngrid,nlay,ptimestep,masse, &
    … …  
    1000998      lvarname = 'pdtadj'
    1001999      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
    1002 
     1000      lvarname = 'alp_bl_det'
     1001      CALL check_var(lfname, lvarname, alp_bl_det,klon, largest, .FALSE.)
     1002      lvarname = 'rhobarz0'
     1003      CALL check_var(lfname, lvarname, rhobarz0,klon, largest, .FALSE.)
     1004      lvarname = 'fraca0'
     1005      CALL check_var(lfname, lvarname, fraca0,klon, largest, .FALSE.)
     1006      lvarname = 'w_conv'
     1007      CALL check_var(lfname, lvarname, w_conv,klon, largest, .FALSE.)
     1008      lvarname = 'w0'
     1009      CALL check_var(lfname, lvarname, w0,klon, largest, .FALSE.)
     1010      lvarname = 'alp_bl_fluct_m'
     1011      CALL check_var(lfname, lvarname, alp_bl_fluct_m,klon, largest, .FALSE.)
     1012      lvarname = 'alp_bl_fluct_tke'
     1013      CALL check_var(lfname, lvarname, alp_bl_fluct_tke,klon, largest, .FALSE.)
     1014      lvarname = 'therm_tke_max0'
     1015      CALL check_var(lfname, lvarname, therm_tke_max0,klon, largest, .FALSE.)
     1016      lvarname = 'env_tke_max0'
     1017      CALL check_var(lfname, lvarname, env_tke_max0,klon, largest, .FALSE.)
     1018      lvarname = 'pbl_tke_max0'
     1019      CALL check_var(lfname, lvarname, pbl_tke_max0,klon, largest, .FALSE.)
     1020      lvarname = 'alp_bl_conv'
     1021      CALL check_var(lfname, lvarname, alp_bl_conv,klon, largest, .FALSE.)
     1022      lvarname = 'alp_bl_stat'
     1023      CALL check_var(lfname, lvarname, alp_bl_stat,klon, largest, .FALSE.)
    10031024!------------Closure------------------
    10041025
    … …  
    10271048
    10281049  ENDIF ! (iflag_clos_bl.ge.1)
     1050      lfname='after closure'
     1051      lvarname = 'pt'
     1052      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
     1053      lvarname = 'pdtadj'
     1054      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
     1055      lvarname = 'alp_bl_det'
     1056      CALL check_var(lfname, lvarname, alp_bl_det,klon, largest, .FALSE.)
     1057      lvarname = 'rhobarz0'
     1058      CALL check_var(lfname, lvarname, rhobarz0,klon, largest, .FALSE.)
     1059      lvarname = 'fraca0'
     1060      CALL check_var(lfname, lvarname, fraca0,klon, largest, .FALSE.)
     1061      lvarname = 'w_conv'
     1062      CALL check_var(lfname, lvarname, w_conv,klon, largest, .FALSE.)
     1063      lvarname = 'w0'
     1064      CALL check_var(lfname, lvarname, w0,klon, largest, .FALSE.)
     1065      lvarname = 'alp_bl_fluct_m'
     1066      CALL check_var(lfname, lvarname, alp_bl_fluct_m,klon, largest, .FALSE.)
     1067      lvarname = 'alp_bl_fluct_tke'
     1068      CALL check_var(lfname, lvarname, alp_bl_fluct_tke,klon, largest, .FALSE.)
     1069      lvarname = 'therm_tke_max0'
     1070      CALL check_var(lfname, lvarname, therm_tke_max0,klon, largest, .FALSE.)
     1071      lvarname = 'env_tke_max0'
     1072      CALL check_var(lfname, lvarname, env_tke_max0,klon, largest, .FALSE.)
     1073      lvarname = 'pbl_tke_max0'
     1074      CALL check_var(lfname, lvarname, pbl_tke_max0,klon, largest, .FALSE.)
     1075      lvarname = 'alp_bl_conv'
     1076      CALL check_var(lfname, lvarname, alp_bl_conv,klon, largest, .FALSE.)
     1077      lvarname = 'alp_bl_stat'
     1078      CALL check_var(lfname, lvarname, alp_bl_stat,klon, largest, .FALSE.)
    10291079
    10301080!!! fin nrlmd le 10/04/2012
    … …  
    10631113         enddo
    10641114      enddo
     1115      lfname='after calculation of Al[p/e]_bl'
     1116      lvarname = 'fm'
     1117      CALL check_var3D(lfname, lvarname, fm, klon, klev, largest, .FALSE.)
     1118      lvarname = 'rhobarz'
     1119      CALL check_var3D(lfname, lvarname, rhobarz, klon, klev, largest, .FALSE.)
     1120      lvarname = 'Alp_bl'
     1121      CALL check_var(lfname, lvarname, Alp_bl,klon, largest, .FALSE.)
     1122      lvarname = 'Ale_bl'
     1123      CALL check_var(lfname, lvarname, Ale_bl,klon, largest, .FALSE.)
    10651124
    10661125! assez bizarre car, si on est dans la couche d'alim et que alim_star et
Note: See TracChangeset for help on using the changeset viewer.