Module land_model_mod

OVERVIEW

This module contains calls to the sub-models to calculate quantities on the fast and slow time scales and update the boundary conditions and tiling structure where necessary.

OTHER MODULES USED

  time_manager_mod
   mpp_domains_mod
           fms_mod
    vegetation_mod
          soil_mod
    land_types_mod
     constants_mod
  diag_manager_mod
 field_manager_mod
tracer_manager_mod

PUBLIC INTERFACE

PUBLIC DATA

PUBLIC ROUTINES

1. land_model_init

   subroutine init_land_with_xgrid & (atmos2land, bnd, time_init, time, dt_fast, dt_slow, atmos_domain) type(atmos_land_boundary_type), intent(inout) :: atmos2land ! land boundary data type(land_data_type), intent(inout) :: bnd ! land boundary data type(time_type), intent(in) :: time_init ! initial time of simulation (?) type(time_type), intent(in) :: time ! current time type(time_type), intent(in) :: dt_fast ! fast time step type(time_type), intent(in) :: dt_slow ! slow time step type(domain2d), intent(in), optional :: atmos_domain ! domain of computations

   subroutine init_land_with_area & ( atmos2land, bnd, gblon, gblat, gcellarea, gfrac, time, dt_fast, dt_slow, domain, & glon, glat ) type(atmos_land_boundary_type), intent(inout) :: atmos2land ! land boundary data type(land_data_type), intent(inout) :: bnd ! state to update real, intent(in) :: gblon(:)! lon boundaries of the grid cells real, intent(in) :: gblat(:)! lat boundaries of the grid cells real, intent(in) :: gcellarea(:,:) ! full area of the grid cells real, intent(in) :: gfrac(:,:) ! fraction of land in the grid cell type(time_type), intent(in) :: time ! current time type(time_type), intent(in) :: dt_fast ! fast time step type(time_type), intent(in) :: dt_slow ! slow time step type(domain2d), intent(in), optional :: domain ! domain of computations real, intent(in), optional :: glon(:), glat(:) ! centers ! of the grid cells

   
   

   DESCRIPTION

   Initializes the land state in three ways. 1) The grid description file may be used as an input. The land grid boundaries and area of land are read from a grid description file. 2) A logical land mask may be used to calculate the area of the grid cells where the mask is true and calls init_land_with_area. Naturally, in this case there are no partly covered land cells -- every cell is either land or ocean. 3) The land state may be initialized from the land area for each of the grid points.

   
   
   

2. land_model_end

   subroutine land_model_end ( atmos2land, bnd )
     type(atmos_land_boundary_type), intent(inout) :: atmos2land
     type(land_data_type), intent(inout) :: bnd        ! state to update

   call land_model_end ( bnd )

   DESCRIPTION

   Destructs the land model and its sub-models. Deallocates the arrays and the boundary exchange data

   
   
   

3. update_land_model_fast

   subroutine update_land_model_fast ( atmos2land, bnd )
   
     type(atmos_land_boundary_type), intent(inout)    :: atmos2land ! quantities exchanged between
                                                                 ! the atmosphere and the land
     type(land_data_type),  intent(inout) :: bnd ! state to update

   call update_land_model_fast ( atmos2land, bnd )

   DESCRIPTION

   Updates the state of the land model on the fast time scale. To implement implicit calculation scheme, first the land sub-models are called in top-bottom general order to evaluate the derivatives; then they are called in bottom-top order to finish implicit calculations. On the upward pass they have a chance to update boundary values returned to the atmosphere. After calculations, update the boundary conditions visible to the atmosphere.

   
   
   

4. update_land_model_slow

   subroutine update_land_model_slow ( atmos2land, bnd )
     type(atmos_land_boundary_type), intent(inout) :: atmos2land
     type(land_data_type), intent(inout) :: bnd        ! state to update

   call update_land_model_slow ( bnd )

   DESCRIPTION

   Updates land, and boundary conditions if necessary, on the slow time scale.

   
   
   

5. update_land_bnd_fast

   subroutine update_land_bnd_fast ( bnd )
   
     type(land_data_type), intent(inout) :: bnd        ! state to update

   call update_land_bnd_fast ( bnd )

   DESCRIPTION

   Updates land boundary data (the ones that atmosphere sees) on the fast time scale. This routine does not update the tiling structure because it is assumed that the tiling does not change on fast time scale.

   
   
   

6. update_land_bnd_slow

   subroutine update_land_bnd_slow ( bnd )
   
     type(land_data_type), intent(inout) :: bnd        ! state to update

   call update_land_bnd_slow ( bnd )

   DESCRIPTION

   Updates the land boundary data for the atmosphere on the slow time scale. This subroutine is responsible for changing of the tiling structure, if necessary, as well as for changing albedo, drag coefficients and such.

   
   
   

   NOTE

   If the tiling structure has been modified, then probably the distribution of other boundary values, such as temperature and surface humidity, should be modified too, not yet clear how.

   
   
   

7. init_land_diag

   subroutine init_land_diag(glonb, glatb, domain, id_lon, id_lat, glon, glat)
   
     real, intent(in)           :: glonb(:), glatb(:) ! longitude and latitude boundaries of grid
                                                      ! cells, specified for the global grid
                                                      ! (not only domain)
     type(domain2d), intent(in) :: domain ! the domain of operations
     integer, intent(out)       :: id_lon, id_lat     ! IDs of respective diag. manager axes
     real, intent(in),optional  :: glon(:), glat(:)   ! coordinates of grid cell centers

   call init_land_diag (glonb, glatb, domain, id_lon, id_lat)

   DESCRIPTION

   Initialize the horizontal axes for the land grid so that all submodules use them, instead of creating their own.

   
   
   

8. err

   subroutine err(ierr, file, line)
   
     integer,          intent(in) :: ierr ! error code
     character(len=*), intent(in) :: file ! name of the file
     integer,          intent(in) :: line ! number of line in the file

   call err (ierr, file, line)

   DESCRIPTION

   Prints out NetCDF library error message, including file name and line number.

   
   
   

PUBLIC TYPES

DIAGNOSTIC FIELDS

Diagnostic fields for module name identifier: 

field name   
----------   
lonb   
lon   
latb   
lat   

field description (units)
-------------------------
Longitude cell boundaries (degrees_E)
Diagnostic manager longitude axis (degrees_E)
Latitude cell boundaries (degrees_N)
Diagnostic manager latitude axis (degrees_N)

DATA SETS

ERROR MESSAGES