LCOV - code coverage report
Current view: top level - physics/clubb/src/CLUBB_core - stats_sfc_module.F90 (source / functions) Hit Total Coverage
Test: coverage.info Lines: 0 152 0.0 %
Date: 2024-12-17 17:57:11 Functions: 0 1 0.0 %

          Line data    Source code
       1             : !-----------------------------------------------------------------------
       2             : ! $Id$
       3             : !===============================================================================
       4             : module stats_sfc_module
       5             : 
       6             :   implicit none
       7             : 
       8             :   private ! Set Default Scope
       9             : 
      10             :   public :: stats_init_sfc
      11             : 
      12             :   ! Constant parameters
      13             :   integer, parameter, public :: nvarmax_sfc = 250  ! Maximum variables allowed
      14             : 
      15             :   contains
      16             : 
      17             : !-----------------------------------------------------------------------
      18           0 :   subroutine stats_init_sfc( vars_sfc,                    & ! intent(in)
      19             :                              l_error,                     & ! intent(inout)
      20             :                              stats_metadata, stats_sfc )    ! intent(inout)
      21             : 
      22             : ! Description:
      23             : !   Initializes array indices for stats_sfc
      24             : ! References:
      25             : !   None
      26             : !-----------------------------------------------------------------------
      27             : 
      28             :     use constants_clubb, only: &
      29             :         fstderr ! Constant(s)
      30             : 
      31             :     use stats_type_utilities, only: &
      32             :         stat_assign ! Procedure
      33             : 
      34             :     use stats_type, only: &
      35             :         stats ! Type
      36             : 
      37             :     use stats_variables, only: &
      38             :         stats_metadata_type
      39             : 
      40             :     implicit none
      41             : 
      42             :     ! External
      43             :     intrinsic :: trim
      44             : 
      45             :     !--------------------- Input Variable ---------------------
      46             :     character(len= * ), dimension(nvarmax_sfc), intent(in) :: &
      47             :       vars_sfc
      48             : 
      49             :     !--------------------- InOut Variables ---------------------      
      50             :     type (stats_metadata_type), intent(inout) :: &
      51             :       stats_metadata
      52             : 
      53             :     type (stats), target, intent(inout) :: &
      54             :       stats_sfc
      55             : 
      56             :     logical, intent(inout) :: l_error
      57             : 
      58             :     !--------------------- Local Varables ---------------------
      59             :     integer :: i, k
      60             : 
      61             :     !--------------------- Begin Code ---------------------
      62             : 
      63             :     ! Default initialization for array indices for stats_sfc is zero (see module
      64             :     ! stats_variables)
      65             : 
      66             :     ! Assign pointers for statistics variables stats_sfc using stat_assign
      67             : 
      68           0 :     k = 1
      69           0 :     do i = 1, stats_sfc%num_output_fields
      70             : 
      71           0 :       select case ( trim( vars_sfc(i) ) )
      72             :       case ('soil_heat_flux')
      73           0 :         stats_metadata%isoil_heat_flux = k
      74             : 
      75             :         call stat_assign( var_index=stats_metadata%isoil_heat_flux, var_name="soil_heat_flux", &
      76             :              var_description="soil_heat_flux, soil_heat_flux", &
      77             :              var_units="W/m^2", l_silhs=.false., &
      78           0 :              grid_kind=stats_sfc )
      79           0 :         k = k + 1
      80             :       case ('ustar')
      81           0 :         stats_metadata%iustar = k
      82             : 
      83             :         call stat_assign( var_index=stats_metadata%iustar, var_name="ustar", &
      84             :              var_description="ustar, Friction velocity", var_units="m/s", l_silhs=.false., &
      85           0 :              grid_kind=stats_sfc )
      86           0 :         k = k + 1
      87             :       case ('veg_T_in_K')
      88           0 :         stats_metadata%iveg_T_in_K = k
      89             : 
      90             :         call stat_assign( var_index=stats_metadata%iveg_T_in_K, var_name="veg_T_in_K", &
      91             :              var_description="veg_T_in_K, Surface Vegetation Temperature", var_units="K", &
      92           0 :              l_silhs=.false., grid_kind=stats_sfc )
      93           0 :         k = k + 1
      94             :       case ('sfc_soil_T_in_K')
      95           0 :         stats_metadata%isfc_soil_T_in_K = k
      96             : 
      97             :         call stat_assign( var_index=stats_metadata%isfc_soil_T_in_K, var_name="sfc_soil_T_in_K", &
      98             :              var_description="sfc_soil_T_in_K, Surface soil temperature", &
      99             :              var_units="K", l_silhs=.false., &
     100           0 :              grid_kind=stats_sfc )
     101           0 :         k = k + 1
     102             :       case ('deep_soil_T_in_K')
     103           0 :         stats_metadata%ideep_soil_T_in_K = k
     104             : 
     105             :         call stat_assign( var_index=stats_metadata%ideep_soil_T_in_K, var_name="deep_soil_T_in_K", &
     106             :              var_description="deep_soil_T_in_K, Deep soil Temperature", &
     107             :              var_units="K", l_silhs=.false., &
     108           0 :              grid_kind=stats_sfc )
     109           0 :         k = k + 1
     110             : 
     111             :       case ('lh')
     112           0 :         stats_metadata%ilh = k
     113             :         call stat_assign( var_index=stats_metadata%ilh, var_name="lh", &
     114             :              var_description="lh, Surface latent heating", var_units="W/m2", l_silhs=.false., &
     115           0 :              grid_kind=stats_sfc )
     116           0 :         k = k + 1
     117             : 
     118             :       case ('sh')
     119           0 :         stats_metadata%ish = k
     120             :         call stat_assign( var_index=stats_metadata%ish, var_name="sh", &
     121             :              var_description="sh, Surface sensible heating", var_units="W/m2", &
     122           0 :              l_silhs=.false., grid_kind=stats_sfc )
     123           0 :         k = k + 1
     124             : 
     125             :       case ('cc')
     126           0 :         stats_metadata%icc = k
     127             :         call stat_assign( var_index=stats_metadata%icc, var_name="cc", var_description="cc, Cloud cover", &
     128           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     129           0 :         k = k + 1
     130             : 
     131             :       case ('lwp')
     132           0 :         stats_metadata%ilwp = k
     133             :         call stat_assign( var_index=stats_metadata%ilwp, var_name="lwp", &
     134             :              var_description="lwp, Liquid water path", var_units="kg/m2", l_silhs=.false., &
     135           0 :              grid_kind=stats_sfc )
     136           0 :         k = k + 1
     137             : 
     138             :       case ('vwp')
     139           0 :         stats_metadata%ivwp = k
     140             :         call stat_assign( var_index=stats_metadata%ivwp, var_name="vwp", &
     141             :              var_description="vwp, Vapor water path", var_units="kg/m2", l_silhs=.false., &
     142           0 :              grid_kind=stats_sfc )
     143           0 :         k = k + 1
     144             : 
     145             :       case ('iwp')
     146           0 :         stats_metadata%iiwp = k
     147             :         call stat_assign( var_index=stats_metadata%iiwp, var_name="iwp", &
     148             :              var_description="iwp, Ice water path", var_units="kg/m2", l_silhs=.false., &
     149           0 :              grid_kind=stats_sfc )
     150           0 :         k = k + 1
     151             : 
     152             :       case ('swp')
     153           0 :         stats_metadata%iswp = k
     154             :         call stat_assign( var_index=stats_metadata%iswp, var_name="swp", &
     155             :              var_description="swp, Snow water path", var_units="kg/m2", l_silhs=.false., &
     156           0 :              grid_kind=stats_sfc )
     157           0 :         k = k + 1
     158             : 
     159             :       case ('rwp')
     160           0 :         stats_metadata%irwp = k
     161             :         call stat_assign( var_index=stats_metadata%irwp, var_name="rwp", &
     162             :              var_description="rwp, Rain water path", var_units="kg/m2", l_silhs=.false., &
     163           0 :              grid_kind=stats_sfc )
     164           0 :         k = k + 1
     165             : 
     166             :       case ('z_cloud_base')
     167           0 :         stats_metadata%iz_cloud_base = k
     168             :         call stat_assign( var_index=stats_metadata%iz_cloud_base, var_name="z_cloud_base", &
     169             :              var_description="z_cloud_base, Cloud base altitude", &
     170             :              var_units="m", l_silhs=.false., &
     171           0 :              grid_kind=stats_sfc )
     172           0 :         k = k + 1
     173             : 
     174             :       case ('z_inversion')
     175           0 :         stats_metadata%iz_inversion = k
     176             :         call stat_assign( var_index=stats_metadata%iz_inversion, var_name="z_inversion", &
     177             :              var_description="z_inversion, Inversion altitude", var_units="m", l_silhs=.false., &
     178           0 :              grid_kind=stats_sfc )
     179           0 :         k = k + 1
     180             : 
     181             :       case ('precip_rate_sfc')          ! Brian
     182           0 :         stats_metadata%iprecip_rate_sfc = k
     183             :         call stat_assign( var_index=stats_metadata%iprecip_rate_sfc, var_name="precip_rate_sfc", &
     184             :              var_description="precip_rate_sfc, Surface rainfall rate", var_units="mm/day", &
     185           0 :              l_silhs=.true., grid_kind=stats_sfc )
     186           0 :         k = k + 1
     187             : 
     188             :       case ('rain_flux_sfc')         ! Brian
     189           0 :         stats_metadata%irain_flux_sfc = k
     190             : 
     191             :         call stat_assign( var_index=stats_metadata%irain_flux_sfc, var_name="rain_flux_sfc", &
     192             :              var_description="rain_flux_sfc, Surface rain flux", &
     193             :              var_units="W/m^2", l_silhs=.false., &
     194           0 :              grid_kind=stats_sfc )
     195           0 :         k = k + 1
     196             : 
     197             :       case ('rrm_sfc')       ! Brian
     198           0 :         stats_metadata%irrm_sfc = k
     199             : 
     200             :         call stat_assign( var_index=stats_metadata%irrm_sfc, var_name="rrm_sfc", &
     201             :              var_description="rrm_sfc, Surface rain water mixing ratio", var_units="kg/kg", &
     202           0 :              l_silhs=.false., grid_kind=stats_sfc )
     203           0 :         k = k + 1
     204             : 
     205             :       case ('precip_frac_tol')
     206           0 :         stats_metadata%iprecip_frac_tol = k
     207             : 
     208             :         call stat_assign( var_index=stats_metadata%iprecip_frac_tol, &
     209             :                           var_name="precip_frac_tol", &
     210             :                           var_description="Smallest allowable precipitation " &
     211             :                           // "fraction when hydrometeors are present [-]", &
     212             :                           var_units="-", &
     213           0 :                           l_silhs=.false., grid_kind=stats_sfc )
     214           0 :         k = k + 1
     215             : 
     216             :       case ( 'morr_snow_rate' )
     217           0 :         stats_metadata%imorr_snow_rate = k
     218             :         call stat_assign( var_index=stats_metadata%imorr_snow_rate, var_name="morr_snow_rate", &
     219             :              var_description="morr_snow_rate, Snow+Ice+Graupel fallout rate " &
     220             :              // "from Morrison scheme", &
     221           0 :              var_units="mm/day", l_silhs=.false., grid_kind=stats_sfc )
     222           0 :         k = k + 1
     223             : 
     224             :       case ('wpthlp_sfc')
     225           0 :         stats_metadata%iwpthlp_sfc = k
     226             : 
     227             :         call stat_assign( var_index=stats_metadata%iwpthlp_sfc, var_name="wpthlp_sfc", &
     228             :              var_description="w'thl'_sfc, wpthlp surface flux", &
     229             :              var_units="K m/s", l_silhs=.false., &
     230           0 :              grid_kind=stats_sfc )
     231           0 :         k = k + 1
     232             : 
     233             :       case ('wprtp_sfc')
     234           0 :         stats_metadata%iwprtp_sfc = k
     235             : 
     236             :         call stat_assign( var_index=stats_metadata%iwprtp_sfc, var_name="wprtp_sfc", &
     237             :              var_description="w'rt'_sfc, wprtp surface flux", var_units="(kg/kg) m/s", &
     238           0 :              l_silhs=.false., grid_kind=stats_sfc )
     239           0 :         k = k + 1
     240             : 
     241             :       case ('upwp_sfc')
     242           0 :         stats_metadata%iupwp_sfc = k
     243             : 
     244             :         call stat_assign( var_index=stats_metadata%iupwp_sfc, var_name="upwp_sfc", &
     245             :              var_description="u'w'_sfc, upwp surface flux", &
     246             :              var_units="m^2/s^2", l_silhs=.false., &
     247           0 :              grid_kind=stats_sfc )
     248           0 :         k = k + 1
     249             : 
     250             :       case ('vpwp_sfc')
     251           0 :         stats_metadata%ivpwp_sfc = k
     252             : 
     253             :         call stat_assign( var_index=stats_metadata%ivpwp_sfc, var_name="vpwp_sfc", &
     254             :              var_description="v'w'_sfc, vpwp surface flux", &
     255             :              var_units="m^2/s^2", l_silhs=.false., &
     256           0 :              grid_kind=stats_sfc )
     257           0 :         k = k + 1
     258             : 
     259             :       case ('thlm_vert_avg')
     260           0 :         stats_metadata%ithlm_vert_avg = k
     261             : 
     262             :         call stat_assign( var_index=stats_metadata%ithlm_vert_avg, var_name="thlm_vert_avg", &
     263             :              var_description="thlm_vert_avg, Vertical average (density-weighted) of thlm", &
     264             :              var_units="K", &
     265           0 :              l_silhs=.false., grid_kind=stats_sfc )
     266           0 :         k = k + 1
     267             : 
     268             :       case ('rtm_vert_avg')
     269           0 :         stats_metadata%irtm_vert_avg = k
     270             : 
     271             :         call stat_assign( var_index=stats_metadata%irtm_vert_avg, var_name="rtm_vert_avg", &
     272             :              var_description="rtm_vert_avg, Vertical average (density-weighted) of rtm", &
     273           0 :              var_units="kg/kg", l_silhs=.false., grid_kind=stats_sfc )
     274           0 :         k = k + 1
     275             : 
     276             :       case ('um_vert_avg')
     277           0 :         stats_metadata%ium_vert_avg = k
     278             : 
     279             :         call stat_assign( var_index=stats_metadata%ium_vert_avg, var_name="um_vert_avg", &
     280             :              var_description="um_vert_avg, Vertical average (density-weighted) of um", &
     281             :              var_units="m/s", &
     282           0 :              l_silhs=.false., grid_kind=stats_sfc )
     283           0 :         k = k + 1
     284             : 
     285             :       case ('vm_vert_avg')
     286           0 :         stats_metadata%ivm_vert_avg = k
     287             : 
     288             :         call stat_assign( var_index=stats_metadata%ivm_vert_avg, var_name="vm_vert_avg", &
     289             :              var_description="vm_vert_avg, Vertical average (density-weighted) of vm", &
     290             :              var_units="m/s", &
     291           0 :              l_silhs=.false., grid_kind=stats_sfc )
     292           0 :         k = k + 1
     293             : 
     294             :       case ('wp2_vert_avg')
     295           0 :         stats_metadata%iwp2_vert_avg = k
     296             : 
     297             :         call stat_assign( var_index=stats_metadata%iwp2_vert_avg, var_name="wp2_vert_avg", &
     298             :              var_description="Density-weighted vertical average of w'^2", &
     299           0 :              var_units="m^2/s^2", l_silhs=.false., grid_kind=stats_sfc )
     300           0 :         k = k + 1
     301             : 
     302             :       case ('up2_vert_avg')
     303           0 :         stats_metadata%iup2_vert_avg = k
     304             : 
     305             :         call stat_assign( var_index=stats_metadata%iup2_vert_avg, var_name="up2_vert_avg", &
     306             :              var_description="u'^2_vert_avg, Vertical average (density-weighted) of up2", &
     307           0 :              var_units="m^2/s^2", l_silhs=.false., grid_kind=stats_sfc )
     308           0 :         k = k + 1
     309             : 
     310             :       case ('vp2_vert_avg')
     311           0 :         stats_metadata%ivp2_vert_avg = k
     312             : 
     313             :         call stat_assign( var_index=stats_metadata%ivp2_vert_avg, var_name="vp2_vert_avg", &
     314             :              var_description="v'^2_vert_avg, Vertical average (density-weighted) of vp2", &
     315           0 :              var_units="m^2/s^2", l_silhs=.false., grid_kind=stats_sfc )
     316           0 :         k = k + 1
     317             : 
     318             :       case ('rtp2_vert_avg')
     319           0 :         stats_metadata%irtp2_vert_avg = k
     320             : 
     321             :         call stat_assign( var_index=stats_metadata%irtp2_vert_avg, var_name="rtp2_vert_avg", &
     322             :              var_description="rt'^2_vert_avg, Vertical average (density-weighted) of rtp2", &
     323           0 :              var_units="kg^2/kg^2", l_silhs=.false., grid_kind=stats_sfc )
     324           0 :         k = k + 1
     325             : 
     326             :       case ('thlp2_vert_avg')
     327           0 :         stats_metadata%ithlp2_vert_avg = k
     328             : 
     329             :         call stat_assign( var_index=stats_metadata%ithlp2_vert_avg, var_name="thlp2_vert_avg", &
     330             :              var_description="thl'^2_vert_avg, Vertical average (density-weighted) of thlp2", &
     331           0 :              var_units="K^2", l_silhs=.false., grid_kind=stats_sfc )
     332           0 :         k = k + 1
     333             : 
     334             :       case ('T_sfc')
     335           0 :         stats_metadata%iT_sfc = k
     336             : 
     337             :         call stat_assign( var_index=stats_metadata%iT_sfc, var_name="T_sfc", &
     338             :              var_description="T_sfc, Surface Temperature", var_units="K", l_silhs=.false., &
     339           0 :              grid_kind=stats_sfc )
     340           0 :         k = k + 1
     341             : 
     342             :       case ('wp23_matrix_condt_num')
     343           0 :         stats_metadata%iwp23_matrix_condt_num = k
     344             :         call stat_assign( var_index=stats_metadata%iwp23_matrix_condt_num, var_name="wp23_matrix_condt_num", &
     345             :              var_description="w'^23_matrix_condt_num, Estimate of the condition number for wp2/3", &
     346           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     347           0 :         k = k + 1
     348             : 
     349             :       case ('thlm_matrix_condt_num')
     350           0 :         stats_metadata%ithlm_matrix_condt_num = k
     351             :         call stat_assign( var_index=stats_metadata%ithlm_matrix_condt_num, var_name="thlm_matrix_condt_num", &
     352             :              var_description="thlm_matrix_condt_num, Estimate of the condition " &
     353             :              // "number for thlm/wpthlp", &
     354           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     355           0 :         k = k + 1
     356             : 
     357             :       case ('rtm_matrix_condt_num')
     358           0 :         stats_metadata%irtm_matrix_condt_num = k
     359             : 
     360             :         call stat_assign( var_index=stats_metadata%irtm_matrix_condt_num, var_name="rtm_matrix_condt_num", &
     361             :              var_description="rtm_matrix_condt_num, Estimate of the condition number " &
     362             :              // "for rtm/wprtp", &
     363           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     364           0 :         k = k + 1
     365             : 
     366             :       case ('thlp2_matrix_condt_num')
     367           0 :         stats_metadata%ithlp2_matrix_condt_num = k
     368             : 
     369             :         call stat_assign( var_index=stats_metadata%ithlp2_matrix_condt_num, var_name="thlp2_matrix_condt_num", &
     370             :              var_description="thl'^2_matrix_condt_num, Estimate of the condition " &
     371             :              // "number for thlp2", &
     372           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     373           0 :         k = k + 1
     374             : 
     375             :       case ('rtp2_matrix_condt_num')
     376           0 :         stats_metadata%irtp2_matrix_condt_num = k
     377             :         call stat_assign( var_index=stats_metadata%irtp2_matrix_condt_num, var_name="rtp2_matrix_condt_num", &
     378             :              var_description="rt'^2_matrix_condt_num, Estimate of the condition " &
     379             :              // "number for rtp2", &
     380           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     381           0 :         k = k + 1
     382             : 
     383             :       case ('rtpthlp_matrix_condt_num')
     384           0 :         stats_metadata%irtpthlp_matrix_condt_num = k
     385             :         call stat_assign( var_index=stats_metadata%irtpthlp_matrix_condt_num, &
     386             :              var_name="rtpthlp_matrix_condt_num", &
     387             :              var_description="rt'thl'_matrix_condt_num, Estimate of the condition " &
     388             :              // "number for rtpthlp", &
     389           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     390           0 :         k = k + 1
     391             : 
     392             :       case ('up2_vp2_matrix_condt_num')
     393           0 :         stats_metadata%iup2_vp2_matrix_condt_num = k
     394             :         call stat_assign( var_index=stats_metadata%iup2_vp2_matrix_condt_num, &
     395             :              var_name="up2_vp2_matrix_condt_num", &
     396             :              var_description="u'^2_v'^2_matrix_condt_num, Estimate of the condition " &
     397             :              // "number for up2/vp2", &
     398           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     399           0 :         k = k + 1
     400             : 
     401             :       case ('windm_matrix_condt_num')
     402           0 :         stats_metadata%iwindm_matrix_condt_num = k
     403             :         call stat_assign( var_index=stats_metadata%iwindm_matrix_condt_num, var_name="windm_matrix_condt_num", &
     404             :              var_description="windm_matrix_condt_num, Estimate of the condition " &
     405             :              // "number for the mean wind", &
     406           0 :              var_units="count", l_silhs=.false., grid_kind=stats_sfc )
     407             : 
     408           0 :         k = k + 1
     409             : 
     410             :       case ('rtm_spur_src')
     411           0 :         stats_metadata%irtm_spur_src = k
     412             : 
     413             :         call stat_assign( var_index=stats_metadata%irtm_spur_src, var_name="rtm_spur_src", &
     414             :              var_description="rtm_spur_src, rtm spurious source", var_units="kg/(m^2 s)", &
     415           0 :              l_silhs=.false., grid_kind=stats_sfc )
     416           0 :         k = k + 1
     417             : 
     418             :       case ('thlm_spur_src')
     419           0 :         stats_metadata%ithlm_spur_src = k
     420             : 
     421             :         call stat_assign( var_index=stats_metadata%ithlm_spur_src, var_name="thlm_spur_src", &
     422             :              var_description="thlm_spur_src, thlm spurious source", &
     423           0 :              var_units="(K kg) / (m^2 s)", l_silhs=.false., grid_kind=stats_sfc )
     424           0 :         k = k + 1
     425             : 
     426             :       case ('rs_sd_morr_int')
     427           0 :         stats_metadata%irsm_sd_morr_int = k
     428             : 
     429             :         call stat_assign( var_index=stats_metadata%irsm_sd_morr_int, var_name="rs_sd_morr_int", &
     430             :              var_description="rs_sd_morr_int, rsm_sd_morr vertical integral", &
     431           0 :              var_units="(kg/kg)/s", l_silhs=.true., grid_kind=stats_sfc )
     432           0 :         k = k + 1
     433             :         
     434             :       case ('tot_vartn_normlzd_rtm')
     435           0 :         stats_metadata%itot_vartn_normlzd_rtm = k
     436             : 
     437             :         call stat_assign( var_index=stats_metadata%itot_vartn_normlzd_rtm, var_name="tot_vartn_normlzd_rtm", &
     438             :              var_description="Total variation of rtm in the vertical normalized", &
     439           0 :              var_units="-", l_silhs=.false., grid_kind=stats_sfc )
     440           0 :         k = k + 1
     441             :         
     442             :       case ('tot_vartn_normlzd_thlm')
     443           0 :         stats_metadata%itot_vartn_normlzd_thlm = k
     444             : 
     445             :         call stat_assign( var_index=stats_metadata%itot_vartn_normlzd_thlm, var_name="tot_vartn_normlzd_thlm", &
     446             :              var_description="Total variation of thlm in the vertical normalized", &
     447           0 :              var_units="-", l_silhs=.false., grid_kind=stats_sfc )
     448           0 :         k = k + 1
     449             :         
     450             :       case ('tot_vartn_normlzd_wprtp')
     451           0 :         stats_metadata%itot_vartn_normlzd_wprtp = k
     452             : 
     453             :         call stat_assign( var_index=stats_metadata%itot_vartn_normlzd_wprtp, var_name="tot_vartn_normlzd_wprtp", &
     454             :              var_description="Total variation of wprtp in the vertical normalized", &
     455           0 :              var_units="-", l_silhs=.false., grid_kind=stats_sfc )
     456           0 :         k = k + 1
     457             : 
     458             :       case default
     459           0 :         write(fstderr,*) 'Error:  unrecognized variable in vars_sfc:  ',  &
     460           0 :               trim( vars_sfc(i) )
     461           0 :         l_error = .true.  ! This will stop the run.
     462             : 
     463             :       end select
     464             : 
     465             :     end do ! 1 .. stats_sfc%num_output_fields
     466             : 
     467           0 :     return
     468             : 
     469             :   end subroutine stats_init_sfc
     470             : 
     471             : 
     472             : end module stats_sfc_module

Generated by: LCOV version 1.14