setfluxlook_mod Module Reference

Functions/Subroutines
type(component_descriptor_type) function, public	setfluxlook_get_descriptor ()
	Descriptor of this component for registration. More...
subroutine	initialisation_callback (current_state)
subroutine	timestep_callback (current_state)
subroutine	set_look (current_state)
subroutine	change_look (current_state)
subroutine	set_flux (current_state)
subroutine	read_configuration (current_state)
subroutine	read_dimensions (ncid, time_dim)
	Reads the dimensions from the NetCDF file. More...
subroutine	read_variables (filename, ncid, time_dim, time, surface_temperatures, surface_humidities, surface_latent_heat_flux, surface_sensible_heat_flux)
	Reads the variables from the NetCDF KiD model file. More...
subroutine	read_single_variable (ncid, key, data1d, data3d)
	Reads a single variable out of a NetCDF file. More...
subroutine	check_status (status)
	Will check a NetCDF status and write to log_log error any decoded statuses. More...

Variables
character(len= *), parameter	time_key = "time"
	NetCDF data time key. More...
character(len= *), parameter	surface_temperatures_key = "surface_temperature"
	NetCDF data surface_temperatures. More...
character(len= *), parameter	surface_humidities_key = "surface_humidity"
	NetCDF data surface_humidities. More...
character(len= *), parameter	surface_shf_key = "surface_sensible_heat_flux"
	NetCDF data surface_sensible_heat_flux. More...
character(len= *), parameter	surface_lhf_key = "surface_latent_heat_flux"
	NetCDF data surface_latent_heat_flux. More...
integer, parameter	lookup_entries = 80
	Number of entries for MO lookup tables. More...
integer, parameter	max_file_len =200
	Maximum length of surface condition input filename. More...
integer, parameter	max_surface_inputs =750
	Specifies the maximum number of surface inputs through configuration file Inputs through netcdf files are not limitted by this. More...
character(max_file_len)	input_file
character(len=string_length)	units_surface_temp ='unset'
real(kind=default_precision)	max_change_buoyancy_flux
real(kind=default_precision), dimension(:), allocatable	surface_boundary_input_times
real(kind=default_precision), dimension(:), allocatable	surface_sensible_heat_flux
real(kind=default_precision), dimension(:), allocatable	surface_latent_heat_flux
real(kind=default_precision), dimension(:), allocatable	surface_temperatures
real(kind=default_precision), dimension(:), allocatable	surface_humidities
integer	iqv

Function/Subroutine Documentation

change_look()

subroutine setfluxlook_mod::change_look ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 186 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
 
    real(kind=default_precision):: crelax,&     ! Relaxation parameter to move to new USTR
         dfb,&        ! Iterative difference in buoyancy flux
         rnit,&    ! real number of iterations required
         velnew,&
         dvelustr(current_state%lookup_table_entries), ob, x1, x0
    integer n, ik, &   ! Loop counters
         nit     ! Number of iterations
    if (current_state%fbuoynew .le. 0.0_default_precision) then
      current_state%fbuoy=current_state%fbuoynew
      return
    end if
    if (max_change_buoyancy_flux .gt. 0.0_default_precision) then
      rnit = abs(current_state%fbuoynew-current_state%fbuoy)/max_change_buoyancy_flux
      nit = 1+int(rnit)
    else
      nit=1
    end if              
    if (nit .gt. 10 .or. (current_state%fbuoynew .gt. 0.0_default_precision .and. &
         current_state%fbuoy .le. 0.0_default_precision)) then
      current_state%fbuoy=current_state%fbuoynew
      call set_look(current_state)
      return                                                                   
    end if
    crelax=1./sqrt(real(nit))    ! maybe better with crelax=1.
    dfb=(current_state%fbuoynew-current_state%fbuoy)/real(nit)
    do n=1, nit
      current_state%fbuoy=current_state%fbuoy+dfb                                                          
      do ik=2, current_state%lookup_table_entries-1
        dvelustr(ik)=log(current_state%lookup_table_velocity(ik+1)/current_state%lookup_table_velocity(ik-1))&
             /log(current_state%lookup_table_ustr(ik+1)/current_state%lookup_table_ustr(ik-1))
      end do
      dvelustr(1)=log(current_state%lookup_table_velocity(2)/current_state%lookup_table_velocity(1))/&
           log(current_state%lookup_table_ustr(2)/current_state%lookup_table_ustr(1))
      dvelustr(current_state%lookup_table_entries)=log(current_state%lookup_table_velocity(current_state%lookup_table_entries)/&
           current_state%lookup_table_velocity(current_state%lookup_table_entries-1))/&
           log(current_state%lookup_table_ustr(current_state%lookup_table_entries)/&
           current_state%lookup_table_ustr(current_state%lookup_table_entries-1))
      do ik=1, current_state%lookup_table_entries        
        ! compute new mean vel. based on old ustar and new fbuoy        
        ob=-current_state%lookup_table_ustr(ik)**3/(von_karman_constant*current_state%fbuoy)
        x1=sqrt(sqrt(1.-gammam*(current_state%global_grid%configuration%vertical%zn(2)+z0)/ob))
        x0=sqrt(sqrt(1.-gammam*z0/ob))
        velnew=(current_state%lookup_table_ustr(ik)/von_karman_constant)*(current_state%global_grid%configuration%vertical%zlogm-&
             (2.0_default_precision*log((x1+1.0_default_precision)/(x0+1.0_default_precision)) + &
             log((x1*x1+1.0_default_precision)/(x0*x0+1.0_default_precision)) + 2.0_default_precision*atan(x0) &
             -2.0_default_precision*atan(x1)))        
        ! relax to new ustar        
        current_state%lookup_table_ustr(ik)=current_state%lookup_table_ustr(ik)/&
             (velnew/current_state%lookup_table_velocity(ik))**(crelax/dvelustr(ik))
      end do
    end do
    current_state%cneut=current_state%lookup_table_ustr(current_state%lookup_table_entries)/&
         current_state%lookup_table_velocity(current_state%lookup_table_entries)
    current_state%cfc=current_state%lookup_table_ustr(1)*current_state%lookup_table_velocity(1)**convective_limit  ! _Businger-Dyer
    current_state%fbuoy=current_state%fbuoynew

Here is the call graph for this function:

Here is the caller graph for this function:

check_status()

subroutine setfluxlook_mod::check_status ( integer, intent(in)  status )

private

Will check a NetCDF status and write to log_log error any decoded statuses.

Parameters

status

The NetCDF status flag

Definition at line 494 of file setfluxlook.F90.

    integer, intent(in) :: status
 
    if (status /= nf90_noerr) then
      call log_log(log_error, "NetCDF returned error code of "//trim(nf90_strerror(status)))
    end if

Here is the caller graph for this function:

initialisation_callback()

subroutine setfluxlook_mod::initialisation_callback ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 62 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
    
    current_state%lookup_table_entries=lookup_entries
    call read_configuration(current_state)
 
    if (.not. current_state%passive_q .and. current_state%number_q_fields > 0)then
      if (.not. allocated(current_state%cq))then
        allocate(current_state%cq(current_state%number_q_fields))
        current_state%cq=0.0_default_precision
      end if
      iqv = get_q_index(standard_q_names%VAPOUR, 'setfluxlook')
      current_state%cq(iqv) = ratio_mol_wts-1.0
    end if
 
    if (current_state%use_surface_boundary_conditions)then
    if (current_state%type_of_surface_boundary_conditions == prescribed_surface_fluxes) then
       if (current_state%use_time_varying_surface_values) then
          call set_flux(current_state)
       else
          current_state%surface_temperature_flux=surface_sensible_heat_flux(1) &
               /(current_state%global_grid%configuration%vertical%rho(1)*cp)
          current_state%surface_vapour_flux = surface_latent_heat_flux(1) &
               /(current_state%global_grid%configuration%vertical%rho(1)*rlvap)
       end if
       
       current_state%fbuoy=0.                                                                 
       if(.not. current_state%passive_th) current_state%fbuoy=&
            current_state%global_grid%configuration%vertical%buoy_co(1)*current_state%surface_temperature_flux
       if(.not. current_state%passive_q .and. current_state%number_q_fields > 0)then                                                      
          current_state%fbuoy=current_state%fbuoy+current_state%cq(iqv)*current_state%surface_vapour_flux*g
       end if
       call set_look(current_state)  ! _set M-O lookup table
       current_state%theta_surf=0.0_default_precision
       current_state%surface_vapour_mixing_ratio=0.0_default_precision
       !
    else if (current_state%type_of_surface_boundary_conditions == prescribed_surface_values) then
       !
       ! Prescribed surface temperatures
       if (current_state%use_time_varying_surface_values) then
          call set_flux(current_state)
       else
       ! If surface_values are constant then surface_temperatures prescribed in 
       ! config and read in read_configuration but if humidity not set then
       ! surface vapour (surface_vapour_mixing_ratio) set to saturated value (see read_config)
          if (current_state%saturated_surface)then
             current_state%surface_vapour_mixing_ratio = qsaturation(surface_temperatures(1),current_state%surface_pressure*0.01)
          else 
             current_state%surface_vapour_mixing_ratio =  &
                  options_get_real(current_state%options_database, "surface_vapour_mixing_ratio")
          endif
          
       ! The code below copied from set_flux as these values need to be 
       ! set for both time varying and constant surface values
       ! Set theta_v
          current_state%theta_surf = surface_temperatures(1)*&
               (current_state%surface_reference_pressure/current_state%surface_pressure)**r_over_cp
          current_state%theta_virtual_surf = current_state%theta_surf
          if (.not. current_state%passive_q .and. current_state%number_q_fields > 0) then
             current_state%theta_virtual_surf = current_state%theta_surf +  &
                  current_state%global_grid%configuration%vertical%thref(2)*  &
                  current_state%cq(iqv)*current_state%surface_vapour_mixing_ratio
          end if          
 
          ! Finally set up new values of THVSURF dependent constants
          current_state%cmbc=betam*current_state%global_grid%configuration%vertical%zn(2)*g*&
               von_karman_constant/current_state%theta_virtual_surf
          current_state%rcmbc=1.0_default_precision/current_state%cmbc
          current_state%ellmocon=current_state%theta_virtual_surf/(g*von_karman_constant)
       end if
    end if
    end if
 

Here is the call graph for this function:

Here is the caller graph for this function:

read_configuration()

subroutine setfluxlook_mod::read_configuration ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 315 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
 
 
 
    integer :: ncid, time_dim
    integer :: number_input_humidities
 
 
    current_state%use_surface_boundary_conditions= & 
       options_get_logical(current_state%options_database, "use_surface_boundary_conditions")
 
    if (current_state%use_surface_boundary_conditions)then
    current_state%type_of_surface_boundary_conditions=options_get_integer(current_state%options_database, &
       "type_of_surface_boundary_conditions")
    current_state%use_time_varying_surface_values=options_get_logical(current_state%options_database, &
       "use_time_varying_surface_values")
  
    current_state%saturated_surface = .true. ! We will change this if we find some humidity data
 
    input_file=options_get_string(current_state%options_database, "surface_conditions_file")
    ! Read in the input_file
    if (trim(input_file)=='' .or. trim(input_file)=='None')then
      if (current_state%use_time_varying_surface_values)then 
        allocate(surface_boundary_input_times(max_surface_inputs))
        surface_boundary_input_times=0.0
        call options_get_real_array(current_state%options_database, "surface_boundary_input_times", surface_boundary_input_times)
      end if
      if (current_state%type_of_surface_boundary_conditions == prescribed_surface_fluxes)then
        allocate(surface_sensible_heat_flux(max_surface_inputs),   & 
           surface_latent_heat_flux(max_surface_inputs)            &
           )
        surface_sensible_heat_flux=0.0
        surface_latent_heat_flux=0.0
        call options_get_real_array(current_state%options_database, "surface_sensible_heat_flux", surface_sensible_heat_flux)
        call options_get_real_array(current_state%options_database, "surface_latent_heat_flux", surface_latent_heat_flux)
        number_input_humidities=0
      else if (current_state%type_of_surface_boundary_conditions == prescribed_surface_values) then 
        allocate(surface_temperatures(max_surface_inputs),         &
           surface_humidities(max_surface_inputs)                  &
           )
        surface_temperatures=0.0
        surface_humidities=0.0
        call options_get_real_array(current_state%options_database, "surface_temperatures", surface_temperatures)
        units_surface_temp=options_get_string(current_state%options_database, "units_surface_temp")
        
        call options_get_real_array(current_state%options_database, "surface_humidities", surface_humidities)
        number_input_humidities=options_get_array_size(current_state%options_database, "surface_humidities")
      end if
    else
      call check_status(nf90_open(path = trim(input_file), mode = nf90_nowrite, ncid = ncid))
      if (log_get_logging_level() .ge. log_debug) then
        call log_master_log(log_debug, "Reading in surface boundary conditions from:"//trim(input_file) )
      end if
      call read_dimensions(ncid, time_dim)
      call read_variables(trim(input_file), ncid, time_dim, &
         surface_boundary_input_times, surface_temperatures, surface_humidities, &
         surface_latent_heat_flux, surface_sensible_heat_flux)
      if (.not. allocated(surface_humidities))then
        number_input_humidities=0
      else
        number_input_humidities=size(surface_humidities)
      end if
      call check_status(nf90_close(ncid))
    end if
    if (number_input_humidities>0)current_state%saturated_surface=.false.
 
    max_change_buoyancy_flux=options_get_real(current_state%options_database, "max_change_buoyancy_flux")
  end if
 
    ! Allocate lookup_tables
    allocate(current_state%lookup_table_velocity(current_state%lookup_table_entries), &
         current_state%lookup_table_ustr(current_state%lookup_table_entries))
 

Here is the call graph for this function:

Here is the caller graph for this function:

read_dimensions()

subroutine setfluxlook_mod::read_dimensions ( integer, intent(in)  ncid, integer, intent(out)  time_dim  )

private

Reads the dimensions from the NetCDF file.

Parameters

ncid	The NetCDF file id
time_dim	Number of elements in the time dimension

Definition at line 396 of file setfluxlook.F90.

    integer, intent(in) :: ncid
    integer, intent(out) ::  time_dim
    integer ::  time_dimid
 
    call check_status(nf90_inq_dimid(ncid, time_key, time_dimid))
 
    call check_status(nf90_inquire_dimension(ncid, time_dimid, len=time_dim))
 

Here is the call graph for this function:

Here is the caller graph for this function:

read_single_variable()

subroutine setfluxlook_mod::read_single_variable ( integer, intent(in)  ncid, character(len=*), intent(in)  key, real(kind=default_precision), dimension(:), intent(inout), optional  data1d, real(kind=default_precision), dimension(:,:,:), intent(inout), optional  data3d  )

private

Reads a single variable out of a NetCDF file.

Parameters

ncid	The NetCDF file id
key	The variable key (name) to access
data1d	Optional one dimensional data to read into
data3d	Optional three dimensional data to read into

Definition at line 467 of file setfluxlook.F90.

    integer, intent(in) :: ncid
    character(len=*), intent(in) :: key
    real(kind=default_precision), dimension(:), intent(inout), optional :: data1d
    real(kind=default_precision), dimension(:,:,:), intent(inout), optional :: data3d
 
    integer :: variable_id
    real(kind=default_precision), dimension(:,:,:), allocatable :: sdata
 
    call check_status(nf90_inq_varid(ncid, key, variable_id))
 
    if (.not. present(data1d) .and. .not. present(data3d)) return
 
    if (present(data1d)) then
      call check_status(nf90_get_var(ncid, variable_id, data1d))
    else
      ! 3D will reshape the data to take account of the column-row major C-F transposition
      allocate(sdata(size(data3d,1),size(data3d,3), size(data3d,2)))
      call check_status(nf90_get_var(ncid, variable_id, sdata))
      data3d(:,:,:)=reshape(sdata(:,:,:),(/size(data3d,1),size(data3d,2),size(data3d,3)/))
      deallocate(sdata)
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

read_variables()

subroutine setfluxlook_mod::read_variables ( character(*), intent(in)  filename, integer, intent(in)  ncid, integer, intent(in)  time_dim, real(kind=default_precision), dimension(:), intent(inout), allocatable  time, real(kind=default_precision), dimension(:), intent(inout), allocatable  surface_temperatures, real(kind=default_precision), dimension(:), intent(inout), allocatable  surface_humidities, real(kind=default_precision), dimension(:), intent(inout), allocatable  surface_latent_heat_flux, real(kind=default_precision), dimension(:), intent(inout), allocatable  surface_sensible_heat_flux  )

private

Reads the variables from the NetCDF KiD model file.

Parameters

ncid	The id of the NetCDF file
time_dim	The number of elements in the time dimension
time	The time data field that is to be read
surface_temperatures	surface temperatures (K)

Definition at line 412 of file setfluxlook.F90.

    character(*), intent(in) :: filename
    integer, intent(in) :: ncid, time_dim
    real(kind=default_precision), dimension(:), allocatable, intent(inout) :: time
    real(kind=default_precision), dimension(:), allocatable, intent(inout) :: surface_temperatures
    real(kind=default_precision), dimension(:), allocatable, intent(inout) :: surface_humidities
    real(kind=default_precision), dimension(:), allocatable, intent(inout) :: surface_latent_heat_flux
    real(kind=default_precision), dimension(:), allocatable, intent(inout) :: surface_sensible_heat_flux
 
    integer :: status, variable_id
 
    ! Do some checking on the variable contents so that we can deal with different 
    ! variable names or missing variables
    
    ! time...
    status=nf90_inq_varid(ncid, time_key, variable_id)
    if (status==nf90_noerr)then
      allocate(time(time_dim))
      call read_single_variable(ncid, time_key, data1d=time)
    else
      call log_log(log_error, "No recognized time variable found in"//trim(filename))
    end if
    
    status=nf90_inq_varid(ncid, surface_temperatures_key, variable_id)
    if (status==nf90_noerr)then
      allocate(surface_temperatures(time_dim))
      call read_single_variable(ncid, surface_temperatures_key, data1d=surface_temperatures)
    end if
    
    status=nf90_inq_varid(ncid, surface_humidities_key, variable_id)
    if (status==nf90_noerr)then
      allocate(surface_humidities(time_dim))
      call read_single_variable(ncid, surface_humidities_key, data1d=surface_humidities)
    end if
    
    status=nf90_inq_varid(ncid, surface_lhf_key, variable_id)
    if (status==nf90_noerr)then
      allocate(surface_latent_heat_flux(time_dim))
      call read_single_variable(ncid, surface_lhf_key, data1d=surface_latent_heat_flux)
    end if
    
    status=nf90_inq_varid(ncid, surface_shf_key, variable_id)
    if (status==nf90_noerr)then
      allocate(surface_sensible_heat_flux(time_dim))
      call read_single_variable(ncid, surface_shf_key, data1d=surface_sensible_heat_flux)
    end if
 

Here is the call graph for this function:

Here is the caller graph for this function:

set_flux()

subroutine setfluxlook_mod::set_flux ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 246 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
 
    real(kind=default_precision) :: surface_temp   ! Surface temperature
  
    if (current_state%type_of_surface_boundary_conditions == prescribed_surface_fluxes) then  ! Prescribed surface fluxes
      
      ! Linear interpolation of input data...
      call interpolate_point_linear_1d(surface_boundary_input_times, & 
         surface_sensible_heat_flux/(current_state%global_grid%configuration%vertical%rho(1)*cp), &
         current_state%time, current_state%surface_temperature_flux, &
         extrapolate='constant') 
      
      call interpolate_point_linear_1d(surface_boundary_input_times, & 
         surface_latent_heat_flux/(current_state%global_grid%configuration%vertical%rho(1)*rlvap), &
         current_state%time, current_state%surface_vapour_flux, &
         extrapolate='constant')
 
      ! Update buoyancy flux...
      current_state%fbuoynew=0.0_default_precision
      if (.not. current_state%passive_th) current_state%fbuoynew=&
         current_state%global_grid%configuration%vertical%buoy_co(1)*current_state%surface_temperature_flux
      if (.not. current_state%passive_q) then                                                      
        current_state%fbuoynew=current_state%fbuoynew+current_state%cq(iqv)*current_state%surface_vapour_flux*g
      end if
 
    else if (current_state%type_of_surface_boundary_conditions == prescribed_surface_values) then   ! Prescribed surface temperatures
 
      ! Linear interpolation of input data...
      call interpolate_point_linear_1d(surface_boundary_input_times, & 
         surface_temperatures, &
         current_state%time, surface_temp, &
         extrapolate='constant') 
      
      select case(trim(units_surface_temp))
      case(degc) ! degrees C
         surface_temp = surface_temp + 273.15_default_precision
      case default ! kelvin
      end select
 
      if (current_state%saturated_surface)then
        current_state%surface_vapour_mixing_ratio = qsaturation(surface_temp,current_state%surface_pressure*0.01)
      else
        call interpolate_point_linear_1d(surface_boundary_input_times, & 
           surface_humidities, &
           current_state%time, current_state%surface_vapour_mixing_ratio, &
           extrapolate='constant') 
      end if
 
      ! Set theta_v
      current_state%theta_surf = surface_temp*&
         (current_state%surface_reference_pressure/current_state%surface_pressure)**r_over_cp
      current_state%theta_virtual_surf = current_state%theta_surf
      if (.not. current_state%passive_q) then
          current_state%theta_virtual_surf = current_state%theta_surf +  &
             current_state%global_grid%configuration%vertical%thref(2)*  &
             current_state%cq(iqv)*current_state%surface_vapour_mixing_ratio
      end if
 
      ! Finally set up new values of THVSURF dependent constants
      current_state%cmbc=betam*current_state%global_grid%configuration%vertical%zn(2)*g*&
         von_karman_constant/current_state%theta_virtual_surf
      current_state%rcmbc=1.0_default_precision/current_state%cmbc
      current_state%ellmocon=current_state%theta_virtual_surf/(g*von_karman_constant)
 
    end if
 

Here is the caller graph for this function:

set_look()

subroutine setfluxlook_mod::set_look ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 149 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
 
    integer I, ik, iters(current_state%lookup_table_entries)    ! Loop counters
    real(kind=default_precision) :: smth, &       ! Relaxation parameter for unstable case
         ob, x1, x0
 
    if (current_state%fbuoy .le. 0.0_default_precision) return
    current_state%velmax=100.0_default_precision
    current_state%velmin=0.1_default_precision
    current_state%aloginv=1.0_default_precision/log(current_state%velmax/current_state%velmin)
    smth=0.1_default_precision  ! _relaxation parameter for unstable case
    do ik=1, current_state%lookup_table_entries
      current_state%lookup_table_velocity(ik)=current_state%velmin*(current_state%velmax/current_state%velmin)**&
           (real(ik-1)/real(current_state%lookup_table_entries-1))            
      current_state%lookup_table_ustr(ik)=current_state%lookup_table_velocity(ik)*&
           current_state%global_grid%configuration%vertical%vk_on_zlogm
      if (current_state%fbuoy .gt. 0.0_default_precision) then         ! _unstable                                  
        iters(ik)=0
        do i=1, 30      ! @check how many iterations needed!!
          iters(ik)=iters(ik)+1
          ob=-current_state%lookup_table_ustr(ik)**3/(von_karman_constant*current_state%fbuoy)
          x1=sqrt(sqrt(1.-gammam*(current_state%global_grid%configuration%vertical%zn(2)+z0)/ob))
          x0=sqrt(sqrt(1.-gammam*z0/ob))
          current_state%lookup_table_ustr(ik)=(1.-smth)*current_state%lookup_table_ustr(ik) + smth*&
               (von_karman_constant*current_state%lookup_table_velocity(ik)) / &
               (current_state%global_grid%configuration%vertical%zlogm-(2.0_default_precision*&
               log((x1+1.0_default_precision)/(x0+1.0_default_precision)) + log((x1*x1+1.0_default_precision)/&
               (x0*x0+1.0_default_precision)) + 2.0_default_precision*atan(x0) -2.0_default_precision*atan(x1)))
        end do
      end if
    end do
    current_state%cneut=current_state%lookup_table_ustr(current_state%lookup_table_entries)/&
         current_state%lookup_table_velocity(current_state%lookup_table_entries)
    current_state%cfc=current_state%lookup_table_ustr(1)*current_state%lookup_table_velocity(1)**convective_limit  ! _Businger-Dyer

Here is the caller graph for this function:

setfluxlook_get_descriptor()

type(component_descriptor_type) function, public setfluxlook_mod::setfluxlook_get_descriptor

Descriptor of this component for registration.

Returns

The diverr component descriptor

Definition at line 55 of file setfluxlook.F90.

    setfluxlook_get_descriptor%name="setfluxlook"
    setfluxlook_get_descriptor%version=0.1
    setfluxlook_get_descriptor%initialisation=>initialisation_callback
    setfluxlook_get_descriptor%timestep=>timestep_callback

Here is the call graph for this function:

timestep_callback()

subroutine setfluxlook_mod::timestep_callback ( type(model_state_type), intent(inout), target  current_state )

private

Definition at line 137 of file setfluxlook.F90.

    type(model_state_type), intent(inout), target :: current_state
 
    if (current_state%use_surface_boundary_conditions)then
      if (current_state%use_time_varying_surface_values)then
        call set_flux(current_state)
        call change_look(current_state)
      end if
    end if

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

input_file

character(max_file_len) setfluxlook_mod::input_file

private

Definition at line 36 of file setfluxlook.F90.

  character(MAX_FILE_LEN) :: input_file

iqv

integer setfluxlook_mod::iqv

private

Definition at line 47 of file setfluxlook.F90.

  integer :: iqv  ! index for vapour

lookup_entries

integer, parameter setfluxlook_mod::lookup_entries = 80

private

Number of entries for MO lookup tables.

Definition at line 31 of file setfluxlook.F90.

  integer, parameter :: LOOKUP_ENTRIES = 80

max_change_buoyancy_flux

real(kind=default_precision) setfluxlook_mod::max_change_buoyancy_flux

private

Definition at line 40 of file setfluxlook.F90.

  real(kind=default_precision) :: max_change_buoyancy_flux

max_file_len

integer, parameter setfluxlook_mod::max_file_len =200

private

Maximum length of surface condition input filename.

Definition at line 32 of file setfluxlook.F90.

  integer, parameter :: MAX_FILE_LEN=200

max_surface_inputs

integer, parameter setfluxlook_mod::max_surface_inputs =750

private

Specifies the maximum number of surface inputs through configuration file Inputs through netcdf files are not limitted by this.

Definition at line 33 of file setfluxlook.F90.

  integer, parameter :: MAX_SURFACE_INPUTS=750

surface_boundary_input_times

real(kind=default_precision), dimension(:), allocatable setfluxlook_mod::surface_boundary_input_times

private

Definition at line 41 of file setfluxlook.F90.

  real(kind=default_precision), allocatable :: surface_boundary_input_times(:)

surface_humidities

real(kind=default_precision), dimension(:), allocatable setfluxlook_mod::surface_humidities

private

Definition at line 45 of file setfluxlook.F90.

  real(kind=default_precision), allocatable :: surface_humidities(:)

surface_humidities_key

character(len=*), parameter setfluxlook_mod::surface_humidities_key = "surface_humidity"

private

NetCDF data surface_humidities.

Definition at line 24 of file setfluxlook.F90.

surface_latent_heat_flux

real(kind=default_precision), dimension(:), allocatable setfluxlook_mod::surface_latent_heat_flux

private

Definition at line 43 of file setfluxlook.F90.

  real(kind=default_precision), allocatable :: surface_latent_heat_flux(:)

surface_lhf_key

character(len=*), parameter setfluxlook_mod::surface_lhf_key = "surface_latent_heat_flux"

private

NetCDF data surface_latent_heat_flux.

Definition at line 24 of file setfluxlook.F90.

surface_sensible_heat_flux

real(kind=default_precision), dimension(:), allocatable setfluxlook_mod::surface_sensible_heat_flux

private

Definition at line 42 of file setfluxlook.F90.

  real(kind=default_precision), allocatable :: surface_sensible_heat_flux(:)

surface_shf_key

character(len=*), parameter setfluxlook_mod::surface_shf_key = "surface_sensible_heat_flux"

private

NetCDF data surface_sensible_heat_flux.

Definition at line 24 of file setfluxlook.F90.

surface_temperatures

real(kind=default_precision), dimension(:), allocatable setfluxlook_mod::surface_temperatures

private

Definition at line 44 of file setfluxlook.F90.

  real(kind=default_precision), allocatable :: surface_temperatures(:)

surface_temperatures_key

character(len=*), parameter setfluxlook_mod::surface_temperatures_key = "surface_temperature"

private

NetCDF data surface_temperatures.

Definition at line 24 of file setfluxlook.F90.

time_key

character(len=*), parameter setfluxlook_mod::time_key = "time"

private

NetCDF data time key.

Definition at line 24 of file setfluxlook.F90.

  character(len=*), parameter ::                              &
     TIME_KEY =                  "time",                      & !<  NetCDF data time key
     surface_temperatures_key =  "surface_temperature",       & 
     surface_humidities_key   =  "surface_humidity",          & 
     surface_shf_key          =  "surface_sensible_heat_flux",& 
     surface_lhf_key          =  "surface_latent_heat_flux"     

units_surface_temp

character(len=string_length) setfluxlook_mod::units_surface_temp ='unset'

private

Definition at line 38 of file setfluxlook.F90.

  character(len=STRING_LENGTH) :: units_surface_temp='unset'  ! units of theta variable forcing