Changeset 1691 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Oct 26, 2015 4:17:44 PM (9 years ago)

Author:

maronga

Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (10 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -14 +14 @@
 ! PALM. If not, see <http://www.gnu.org/licenses/>.
 !
-! Copyright 1997-2014 Leibniz Universitaet Hannover
+! Copyright 1997-2015 Leibniz Universitaet Hannover
 !--------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! ------------------
-! 
+! Call to init_surface_layer added. rif is replaced by ol and zeta.
 ! 
 ! Former revisions:
@@ -266 +266 @@
                sums_l_l, sums_up_fraction_l, sums_wsts_bc_l, ts_value,         &
                var_d, weight_pres, weight_substep 
-    
+ 
+    USE surface_layer_fluxes_mod,                                              &
+        ONLY:  init_surface_layer_fluxes
+   
     USE transpose_indices 
 
@@ -313 +316 @@
     ALLOCATE( ptdf_x(nxlg:nxrg), ptdf_y(nysg:nyng) )
 
-    ALLOCATE( rif(nysg:nyng,nxlg:nxrg), shf(nysg:nyng,nxlg:nxrg),     &
+    ALLOCATE( ol(nysg:nyng,nxlg:nxrg), shf(nysg:nyng,nxlg:nxrg),      &
               ts(nysg:nyng,nxlg:nxrg), tswst(nysg:nyng,nxlg:nxrg),    &
               us(nysg:nyng,nxlg:nxrg), usws(nysg:nyng,nxlg:nxrg),     &
@@ -714 +717 @@
              hom(:,1,25,:) = SPREAD( l1d, 2, statistic_regions+1 )
 
-             IF ( prandtl_layer )  THEN
-                rif  = rif1d(nzb+1)
+             IF ( constant_flux_layer )  THEN
+                ol   = ( zu(nzb+1) - zw(nzb) ) / rif1d(nzb+1)
                 ts   = 0.0_wp  ! could actually be computed more accurately in the
                                ! 1D model. Update when opportunity arises.
@@ -723 +726 @@
              ELSE
                 ts   = 0.0_wp  ! must be set, because used in
-                rif  = 0.0_wp  ! flowste
+                ol   = ( zu(nzb+1) - zw(nzb) ) / zeta_min  ! flowste
                 us   = 0.0_wp
                 usws = 0.0_wp
@@ -731 +734 @@
           ELSE
              e    = 0.0_wp  ! must be set, because used in
-             rif  = 0.0_wp  ! flowste
+             ol   = ( zu(nzb+1) - zw(nzb) ) / zeta_min  ! flowste
              ts   = 0.0_wp
              us   = 0.0_wp
@@ -886 +889 @@
              e    = 0.0_wp
           ENDIF
-          rif   = 0.0_wp
+          ol    = ( zu(nzb+1) - zw(nzb) ) / zeta_min
           ts    = 0.0_wp
           us    = 0.0_wp
@@ -1095 +1098 @@
 !
 !--    Initialize Prandtl layer quantities
-       IF ( prandtl_layer )  THEN
+       IF ( constant_flux_layer )  THEN
 
           z0 = roughness_length
@@ -1103 +1106 @@
 !
 !--          Surface temperature is prescribed. Here the heat flux cannot be
-!--          simply estimated, because therefore rif, u* and theta* would have
+!--          simply estimated, because therefore ol, u* and theta* would have
 !--          to be computed by iteration. This is why the heat flux is assumed
 !--          to be zero before the first time step. It approaches its correct
@@ -1624 +1627 @@
        CALL location_message( 'initializing land surface model', .FALSE. )
        CALL init_lsm
+       CALL location_message( 'finished', .TRUE. )
+    ENDIF
+
+!
+!-- Initialize surface layer (done after LSM as roughness length are required
+!-- for initialization
+    IF ( constant_flux_layer )  THEN
+       CALL location_message( 'initializing surface layer', .FALSE. )
+       CALL init_surface_layer_fluxes
        CALL location_message( 'finished', .TRUE. )
     ENDIF