Changeset 1826 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Apr 7, 2016 12:01:39 PM (8 years ago)

Author:

maronga

Message:

further modularization of radiation model and plant canopy model

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Cleanup after modularization
 ! 
 ! Former revisions:
@@ -532 +532 @@
     PUBLIC lsm_check_data_output, lsm_check_data_output_pr,                    &
            lsm_check_parameters, lsm_energy_balance, lsm_header, lsm_init,     &
-           lsm_init_arrays, lsm_parin, lsm_soil_model 
+           lsm_init_arrays, lsm_parin, lsm_soil_model, lsm_swap_timelevel 
 !
 !-- Public parameters, constants and initial values
-    PUBLIC alpha_vangenuchten, c_surface, canopy_resistance_coefficient,       &
-           conserve_water_content, field_capacity,                             &
-           f_shortwave_incoming, hydraulic_conductivity,                       &
-           lambda_surface_stable, lambda_surface_unstable,                     &
-           land_surface, leaf_area_index,                                      &
-           lsm_swap_timelevel, l_vangenuchten,                                 &
-           min_canopy_resistance, min_soil_resistance, n_vangenuchten,         &
-           pave_heat_capacity, pave_depth, pave_heat_conductivity,             &
-           residual_moisture, rho_cp, rho_lv, root_fraction,                   &
-           saturation_moisture, skip_time_do_lsm, soil_moisture,               &
-           soil_temperature, soil_type, soil_type_name, vegetation_coverage,   &
-           veg_type, veg_type_name, wilting_point, z0_eb, z0h_eb, z0q_eb
+    PUBLIC land_surface, skip_time_do_lsm
 
 !
@@ -644 +633 @@
              unit = 'K'   
              
-          CASE ( 'c_liq*', 'c_soil*', 'c_veg*', 'ghf_eb*',      &
-                 'm_liq_eb*', 'qsws_eb*',      &
-                 'qsws_liq_eb*', 'qsws_soil_eb*', 'qsws_veg_eb*', &
+          CASE ( 'lai*', 'c_liq*', 'c_soil*', 'c_veg*', 'ghf_eb*', 'm_liq_eb*',&
+                 'qsws_eb*', 'qsws_liq_eb*', 'qsws_soil_eb*', 'qsws_veg_eb*',  &
                  'r_a*', 'r_s*', 'shf_eb*' )
              IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
@@ -654 +642 @@
                 CALL message( 'check_parameters', 'PA0111', 1, 2, 0, 6, 0 )
              ENDIF
-
+             IF ( TRIM( var ) == 'lai*'  .AND.  .NOT.  land_surface )  THEN
+                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                                 'res land_surface = .TRUE.'
+                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
+             ENDIF
              IF ( TRIM( var ) == 'c_liq*'  .AND.  .NOT.  land_surface )  THEN
                 message_string = 'output of "' // TRIM( var ) // '" requi' //  &
@@ -716 +708 @@
              ENDIF
 
+             IF ( TRIM( var ) == 'lai*'   )  unit = 'none' 
              IF ( TRIM( var ) == 'c_liq*' )  unit = 'none'
              IF ( TRIM( var ) == 'c_soil*')  unit = 'none'
@@ -1564 +1557 @@
 1   FORMAT (//' Land surface model information:'/                              &
               ' ------------------------------'/)
-2   FORMAT (' --> Soil bottom is closed (water content is conserved, default)')
-3   FORMAT (' --> Soil bottom is open (water content is not conserved)')         
-4   FORMAT (' --> Land surface type  : ',A,/                                   &
-            ' --> Soil porosity type : ',A)
+2   FORMAT ('    --> Soil bottom is closed (water content is conserved',       &
+            ', default)')
+3   FORMAT ('    --> Soil bottom is open (water content is not conserved)')         
+4   FORMAT ('    --> Land surface type  : ',A,/                                &
+            '    --> Soil porosity type : ',A)
 5   FORMAT (/'    Initial soil temperature and moisture profile:'//            &
             '       Height:        ',A,'  m'/                                  &
@@ -1573 +1567 @@
             '       Moisture:      ',A,'  m**3/m**3'/                          &
             '       Root fraction: ',A,'  '/                                   &
-            '       Gridpoint:     ',A)
-
-
+            '       Grid point:    ',A)
 
     END SUBROUTINE lsm_header
@@ -1627 +1619 @@
 
        IF ( humidity )  THEN
-          qsws_eb = rho_l * l_v * qsws
+          qsws_eb = rho_lv * qsws
        ELSE
           qsws_eb = 0.0_wp